Skip to main content

A Letter to the Editor Regarding the Congenital Hypothyroidism Study by Mangano and Sherman

Alfred Körblein wrote this letter to the editor of OJPed in response to the study by Mangano and Sherman.  However, he was notified by the OJPed that the journal did not accept a letter to the editor.

The letter was also sent to both authors, Mangano and Sherman.  Mangano never replied, and Sherman only sent a cursory thank you, merely acknowledging the receipt of the e-mail.


Letter to the editor of OJPed in response to:

Joseph J. Mangano, Janette D. Sherman. Elevated airborne beta levels in Pacific/West Coast US States and trends in hypothyroidism among newborn after the Fukushima nuclear meltdown. Open Journal of Pediatrics, 2013, 3, 1-9.

Alfred Körblein,

Mangano and Sherman’s idea to study congenital hypothyroidism (CH) case numbers among newborn babies is new and seems to be promising. Their approach is reasonable; they compare the CH case numbers after Fukushima with those before Fukushima in 5 Pacific/West Coast states (Hawaii, Alaska, Washington, Oregon, and California: the study region) where the fallout from the Fukushima plumes was higher than in other parts of the USA. And they also look at the 2010-2011 change of CH cases in 36 other US states (the control region).

Mangano and Sherman state:

“The 2010-2011 ratio representing the change in CH cases was 1.16 for the five Pacific/West Coast States, rising from 281 to 327 confirmed cases. The 1.16 ratio exceeded the 0.97 ratio (decline in cases from 1208 to 1167) for the 36 control states; the difference is significant at p < 0.03. Increases in ratios were observed in the exposed areas for the periods March 17-June 30 (1.28, significant at p < 0.04) and July 1-December 31 (1.10, not significant at p < 0.21).”

Unfortunately the authors (M&S) do not clearly say how they tested the significance of their result. In their method section they state:

“CH cases for births in the periods March 17 to December 31 (2010 and 2011) will be compared, for the Pacific/West Coast States and the remainder of the US Portions of this 290 day period will also be compared. Significance testing will be conducted using a t test, where n equals the number of Pacific/West Coast cases in 2010 and 2011, the observed change will be the change in the Pacific/West Coast, and the expected change will be the change for the remainder of the US.”

So their hypothesis to be tested is whether there is a significant difference in 2010-2011 ratios between the study region and the control region. To check whether the reported p-values of p<0.03 and p<0.04 are correct, I applied a Poisson regression of CH numbers from their Table 4 in the two years (2010 and 2011) and in the two regions (study and control) and used two dummy variables: “d11”, which is 1 for 2011 and zero for 2010, and “study”, which is 1 for the study region and zero for the control region. A third dummy variable is the interaction dint=d11*study. Dummy variable “dint” is used to estimate the difference in 2010-2011 ratios between the study and control region.

The result of the Poisson regression was a 20.5% increase in CH cases in 2011 for March 17 to December 31 (p=0.041) and a 35.6% increase for March 17 to June 30 (p=0.049). So both results are statistically significant, albeit not at their stated p<0.03 or p<0.04 levels.

But it turns out that the estimate of the parameter for “d11” is not significant; the p-value s are p=0.400 for March 17 to December 31 and p=0.451 for March 17 to June30. As a rule, parameters that have no meaningful influence on the goodness of fit should be omitted in the regression and, per convention, p>0.2 is considered not meaningful. There is also a criterion that helps choosing the best model, the so called Akaike criterion (AIC), see, which is a measure for the goodness of fit. This is smaller for the regression model without d11 (AIC=35.2) than for the model with d11 (AIC=36.6).

The regression without d11, yields the following result for the shorter period March 17 through June 30:

       Estimate    Std. Error  z value    Pr(>|z|)
(Intercept) 5.96229     0.03587   166.197     <2e-16 ***
study      -1.40842     0.10869   -12.958     <2e-16 ***
dint        0.25014     0.13683     1.828     0.0675 . 

Now the parameter for dint (which estimates the 2010-2011 change in CH cases in the study region) is not significant (p=0.0675). A similar result is obtained for the longer period March 17 through December 31 (p=0.0624).

The use of CH case numbers instead of CH incidences, however, means that their inferences can only be made if the live birth numbers don’t change between 2010 and 2011, or if the 2010-2011 relative change in birth numbers is the same in both the study and control region. 
Official numbers of live births for 2011 are not available yet, but the US Centers for Disease Control and Prevention provide preliminary annual live birth data for 2011 for individual US states on their website [1]. The respective data for 2010 can be found on the statehealthfacts website [2]. Therefore I was able to determine the numbers of live births for 2010 and 2011 in the two periods (March 17-December 31 and March 17-December 31) and calculate the incidence rates. I used a logistic regression of CH incidence rates for the two time windows with only one dummy variable (dint) to estimate the excess in 2011 in the study region. This yielded a 16.6% increase in March 17-December 31 2011 p=0.0087) and a 32.5% increase in March 17-June 30 2011 (p=0.0036). Both results are significant at the p<0.01 level.

The main reason for the lower p-values compared to the regression of the CH case numbers is that in 2010 the incidences agree fairly well in the study and control region (the p-value for the dummy variable study is p=0.735) so study can be omitted. There is also good agreement between the incidences in 2010 and 2011 in the study region, so the dummy variable d11 (p=0.578) can also be omitted.

Essentially, the use of incidence rates instead of case numbers increases the statistical significance of Mangano and Sherman’s findings.

Several technical gremlins appear in the published tables. The confidence intervals in Table 1 and the ratios in the last two columns of Table 3 are incorrect, and the case numbers for the 36 control states in the last two rows do not sum up to the numbers in the first row in Table 4. Also page 3 states “Results showed that for I-131, the highest depositions, in becquerels per cubic meter....“. This should read “per square meter“. This was first spotted on the ex-skf.blogspot website [3].

However these glitches have no influence on the authors’ main finding.



Post a Comment

Popular posts from this blog

"Strict" Management of the Fukushima Health Management Survey Data by Fukushima Medical University

The Interim Report released in March 2016 by the Oversight Committee for Fukushima Health Management Survey called for the need to establish the rules of management and provision of the data so it can be widely utilized by domestic and international experts. Accordingly, the Subcommittee to Review Provision of Data for the Purpose of Academic Research met in Fukushima City, Fukushima, for the first time on May 31, 2016. Subcommittee members were selected by the Fukushima prefectural government in order to address the role of the subcommittee in establishing technical rules in provision of data and include specialists in epidemiology, information technology, law, and legal sociology. Also included are members of the Oversight Committee for Fukushima Health Management Survey, Shoichiro Tsugane and Hokuto Hoshi, and the Health Survey Support Department Head and the Department of Epidemiology Chair at Fukushima Medical University, Tetsuya Ohira, representing the Fukushima Health Managemen…

A Letter Hand-delivered to the UN and UNSCEAR, Requesting Revision of UNSCEAR Report And a New UN Mandate for UNSCEAR

On October 24, 2014, at the Fourth Committee of the UN General Assembly being held in New York City, representatives from Physicians for Social Responsibility (USA) and Human Rights Now (Tokyo, Japan) hand-delivered a letter to the chairperson of the 4th Committee and the Secretary of UNSCEAR.  The letter, co-signed by 43 civil society groups from 9 countries, including 21 Japanese groups, requested revision of the 2014 UNSCEAR report on Fukushima accident as well as a new UN mandate for UNSCEAR.

Date: 24 October 2014
To: Members of the Fourth Committee of the UN General Assembly 69th Session, Members of UNSCEAR, and Members of the UN General Assembly:

Re:Civil Society groups request revision of the recent United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) Report: “Levels and effects of radiation exposure due to the nuclear accident after the 2011 Great East-Japan Earthquake and tsunami.”

The 2011 Fukushima disaster made UN oversight of the adverse effects of …

Fukushima Thyroid Examination May 2015: 103 Thyroid Cancer Cases Confirmed, 5 in the Second Round Screening

The Nineteenth Prefectural Oversight Committee for Fukushima Health Management Survey convened in Fukushima City, Fukushima Prefecture, on May 18, 2015, releasing the latest results of thyroid examination, consisting of Initial Screening or the first round screening (originally scheduled to be conducted from October 2011 to March 2014, but actually still ongoing) and Full-Scale Screening or the second round screening (beginning April 2014). It has been 3 months since the last committee meeting on February 12, 2015, and the latest results include 3 more months worth of data confirmed as of March 31, 2015.

An official English translation of the results is now available here.

As of March 31, 2015, there are 16 more (12 from the first round and 4 from the second round) confirmed cancer cases, all papillary thyroid cancer, and 9 more (2 from the first round and 7 from the second round) newly suspicious cases. The number of confirmed cancer cases now totals 103 (98 from the first round and 5…