Stress after Three Mile Island: Logistic Regression and One-facet Rasch Model

The nuclear accident at Three Mile Island in Spring 1979 produced more than radioactive fall-out. There was also publication of a Rasch analysis in a statistical journal. "[There is] a method for obtaining CMLE's [conditional maximum likelihood estimates] in the Rasch model with standard log-linear model programs such as GLIM and SPSS-X ... [This] method is used to analyze the data from the Three Mile Island study" (Conaway, 1989). What message does this paper have for Rasch measurement practitioners?

Stress in 267 mothers of young children living within 10 miles of the plant were surveyed over three years. 115 mothers, living within 5 miles of the reactor, were assigned to group LT5. The other 152 mothers, living 6 or more miles away, were assigned to group GT5. Each mother was interviewed four times, and her stress level rated high, medium or low. The data were originally published by Fienberg et al. (1985).

Conaway implemented a Rasch model in his logit-linear regression. From Table 1, one of his simpler tables, he concludes that "for the LT5 group, ... levels of stress did not change [significantly] over time". From other tables, he concludes that "in the GT5 group ... there are significant decreases in the levels of stress" (p.58). Even after reading both Conaway and Fienberg, however, I was unclear as to what had happened to mothers' stress levels. Fortunately the data matrix was also published. This afforded the opportunity to reanalyze the data with standard Rasch analysis software.

One approach is a fixed-effects, one-facet Rasch analysis (implemented with Facets). At each of time-point, each group of mothers can be thought of as random replications of a shared group/time-point stress. Since 2 groups are rated at 4 time-points, there are 8 time-point effects. Each stress rating is classified solely by group/time-point, ignoring information about which mother is rated. Thus there is just one facet containing 8 group/time-point elements. Each of the four LT5 time-point elements has 115 observations. Each of the four GT5 time-point elements has 152 observations. Since the survey is intended to have a uniform rating scale, only one rating scale structure is modelled. (An equivalent BIGSTEPS result is obtained by using a data matrix of 257 mothers by 8 group/time-points and anchoring all mothers at zero logits. Then the time-point and rating scale calibrations can be used as anchors for the analysis of time-serial effects.) The Facets output is shown in Table 2. A zero logit measure corresponds to medium stress.

The rating scale step calibrations placed the low stress region at -1.44 logits and below, and the high stress region at 1.44 logits and above. The mean-square variance ratio, "Mnsq", quantifies the extent of central tendency in the stress ratings. Values above 1.0 indicate unexpected heterogeneity in the stress ratings across mothers. Values less than 1.0 indicate unexpected homogeneity. We see that the ratings of the GT5 group are more homogeneous than those of the LT5 group.

Most interesting is how stress changed over time. A plot of the 8 time-point calibrations is shown in the Figure. Immediately after the accident, the two groups of mothers were equally stressed. As time passed, the more distant GT5 group relaxed more than the LT5 group. Since each plotted point has a standard error of about 0.15 logits, the only statistically remarkable points are the first and second time points of the GT5 group. They report a significant drop in stress between Winter 1979 and Spring 1980.

What is the message? Even one-facet analysis can be useful. Though Conaway's analysis is intricate, meticulous and detailed, one simple plot can capture the whole story. John Michael Linacre

Conaway MR 1989. Analysis of repeated categorical measurements with conditional likelihood methods. Journal of the American Statistical Association, 84/5, 53-62.

Fienberg SE, Bromet EJ, Follmann D, Lambert D, May SM 1985. Longitudinal analysis of categorical epidemiological data. Environmental Health Perspectives, 63, 241-248.

   Analysis of Deviance for the LT5 Group
Model                 df   Deviance       p
----------------------------------------------
S                     45   35.94
S+H2                  44   33.58          .12
S+H2+M2               43   33.02          .76
S+H2+M2+H3+M3+H4+M4   39   30.41          .63

Table 1. Conaway's analysis of variance. S represents mothers effects. Hw and Mw stand for the medium and high stress category effects at time-point w (=2,3,4) relative to time-point 1.

Group        Total  Mean  Logit   S.E.   Mnsq.
----------------------------------------------
LT5: 115 Mothers
Winter 1979  133    1.2     .49   .16     1.2
Spring 1980  121    1.1     .16   .16     1.1
Fall   1981  130    1.1     .41   .16     1.1
Fall   1982  124    1.1     .24   .16     1.2

GT5: 152 Mothers
Winter 1979  176    1.2     .49   .14     0.8
Spring 1980  141    0.9    -.23   .14     1.2
Fall   1981  154    1.0     .04   .14     0.7
Fall   1982  150    1.0    -.04   .14     0.9

Table 2. Rasch analysis of the Stress data. Output produced by the "Facets" software.

Stress after Three Mile Island: Logistic Regression and One-facet Rasch Model. Linacre JM. … Rasch Measurement Transactions, 1992, 5:4 p.188


Please help with Standard Dataset 4: Andrich Rating Scale Model



Rasch Publications
Rasch Measurement Transactions (free, online) Rasch Measurement research papers (free, online) Probabilistic Models for Some Intelligence and Attainment Tests, Georg Rasch Applying the Rasch Model 3rd. Ed., Bond & Fox Best Test Design, Wright & Stone
Rating Scale Analysis, Wright & Masters Introduction to Rasch Measurement, E. Smith & R. Smith Introduction to Many-Facet Rasch Measurement, Thomas Eckes Invariant Measurement: Using Rasch Models in the Social, Behavioral, and Health Sciences, George Engelhard, Jr. Statistical Analyses for Language Testers, Rita Green
Rasch Models: Foundations, Recent Developments, and Applications, Fischer & Molenaar Journal of Applied Measurement Rasch models for measurement, David Andrich Constructing Measures, Mark Wilson Rasch Analysis in the Human Sciences, Boone, Stave, Yale
in Spanish: Análisis de Rasch para todos, Agustín Tristán Mediciones, Posicionamientos y Diagnósticos Competitivos, Juan Ramón Oreja Rodríguez

To be emailed about new material on www.rasch.org
please enter your email address here:

I want to Subscribe: & click below
I want to Unsubscribe: & click below

Please set your SPAM filter to accept emails from Rasch.org

www.rasch.org welcomes your comments:

Your email address (if you want us to reply):

 

ForumRasch Measurement Forum to discuss any Rasch-related topic

Go to Top of Page
Go to index of all Rasch Measurement Transactions
AERA members: Join the Rasch Measurement SIG and receive the printed version of RMT
Some back issues of RMT are available as bound volumes
Subscribe to Journal of Applied Measurement

Go to Institute for Objective Measurement Home Page. The Rasch Measurement SIG (AERA) thanks the Institute for Objective Measurement for inviting the publication of Rasch Measurement Transactions on the Institute's website, www.rasch.org.

Coming Rasch-related Events
May 26 - June 23, 2017, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Core Topics (E. Smith, Winsteps), www.statistics.com
June 30 - July 29, 2017, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Further Topics (E. Smith, Winsteps), www.statistics.com
July 31 - Aug. 3, 2017, Mon.-Thurs. Joint IMEKO TC1-TC7-TC13 Symposium 2017: Measurement Science challenges in Natural and Social Sciences, Rio de Janeiro, Brazil, imeko-tc7-rio.org.br
Aug. 7-9, 2017, Mon-Wed. In-person workshop and research coloquium: Effect size of family and school indexes in writing competence using TERCE data (C. Pardo, A. Atorressi, Winsteps), Bariloche Argentina. Carlos Pardo, Universidad Catòlica de Colombia
Aug. 7-9, 2017, Mon-Wed. PROMS 2017: Pacific Rim Objective Measurement Symposium, Sabah, Borneo, Malaysia, proms.promsociety.org/2017/
Aug. 10, 2017, Thurs. In-person Winsteps Training Workshop (M. Linacre, Winsteps), Sydney, Australia. www.winsteps.com/sydneyws.htm
Aug. 11 - Sept. 8, 2017, Fri.-Fri. On-line workshop: Many-Facet Rasch Measurement (E. Smith, Facets), www.statistics.com
Aug. 18-21, 2017, Fri.-Mon. IACAT 2017: International Association for Computerized Adaptive Testing, Niigata, Japan, iacat.org
Sept. 15-16, 2017, Fri.-Sat. IOMC 2017: International Outcome Measurement Conference, Chicago, jampress.org/iomc2017.htm
Oct. 13 - Nov. 10, 2017, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Core Topics (E. Smith, Winsteps), www.statistics.com
Jan. 5 - Feb. 2, 2018, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Core Topics (E. Smith, Winsteps), www.statistics.com
Jan. 10-16, 2018, Wed.-Tues. In-person workshop: Advanced Course in Rasch Measurement Theory and the application of RUMM2030, Perth, Australia (D. Andrich), Announcement
Jan. 17-19, 2018, Wed.-Fri. Rasch Conference: Seventh International Conference on Probabilistic Models for Measurement, Matilda Bay Club, Perth, Australia, Website
April 13-17, 2018, Fri.-Tues. AERA, New York, NY, www.aera.net
May 25 - June 22, 2018, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Core Topics (E. Smith, Winsteps), www.statistics.com
June 29 - July 27, 2018, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Further Topics (E. Smith, Winsteps), www.statistics.com
Aug. 10 - Sept. 7, 2018, Fri.-Fri. On-line workshop: Many-Facet Rasch Measurement (E. Smith, Facets), www.statistics.com
Oct. 12 - Nov. 9, 2018, Fri.-Fri. On-line workshop: Practical Rasch Measurement - Core Topics (E. Smith, Winsteps), www.statistics.com
The HTML to add "Coming Rasch-related Events" to your webpage is:
<script type="text/javascript" src="http://www.rasch.org/events.txt"></script>

 

The URL of this page is www.rasch.org/rmt/rmt54n.htm

Website: www.rasch.org/rmt/contents.htm