To:

Cc:

Bcc:

Subject: Reanalysis of Weaver et al (2008)

De: Weaver, Virginia [mailto:vweaver@jhsph.edu] 

Enviado el: Wednesday, August 18, 2010 8:00 AM

Asunto: RE: request for model replication from old publication

Reanalyis of Weaver et al. (2008): Year 1 data would be better since that was 

the original recruited group and has the least bias related to drop-out. 

Ginnie

J-test (3).rtf

________________________________________

De: Kotchmar.Dennis@epamail.epa.gov [mailto:Kotchmar.Dennis@epamail.epa.gov] 

Enviado el: Tuesday, September 28, 2010 11:01 AM

Asunto: Please email me a copy of the further analysis of Weaver et al. 2008

We need to have all references on file here. I am working thru your

text. If there are other communications I have not found yet please send

them also.

Thank you.

Dennis J. Kotchmar M.D.

National Center for Environmental Assessment

United States Environmental Protection Agency

919-541-4158

kotchmar.dennis@epa.gov
PBBc model 

Model: MODEL1

Dependent Variable: SYST

Number of Observations Read 652

Number of Observations Used 648

Number of Observations with Missing Values 4

Analysis of Variance

Source DF

Sum of

Squares

Mean

Square F Value Pr > F

Model 12 34084 2840.31829 13.19 <.0001

Error 635 136695 215.26789

Corrected Total 647 170779

Root MSE 14.67201 R-Square 0.1996

Dependent Mean 120.63742 Adj R-Sq 0.1845

Coeff Var 12.16207

Parameter Estimates

Variable DF

Parameter

Estimate

Standard

Error t Value Pr > |t|

Intercept 1 119.09183 1.50128 79.33 <.0001

agec 1 0.28396 0.08228 3.45 0.0006

agecSQR 1 0.01507 0.00568 2.65 0.0082

FEMALE 1 -0.51260 2.04009 -0.25 0.8017

bmic 1 1.20424 0.20896 5.76 <.0001

ddose1 1 -4.38060 2.16878 -2.02 0.0438

ddose2 1 -1.92612 1.82489 -1.06 0.2916

ddose3 1 1.11859 2.02409 0.55 0.5807

ddose4 1 2.95745 2.04153 1.45 0.1479

exdrnk 1 -4.93714 2.61252 -1.89 0.0592

htnmeds 1 18.12350 3.07294 5.90 <.0001

pbbc 1 0.09958 0.04043 2.46 0.0140

jobdurc 1 -0.19959 0.11061 -1.80 0.0716
LNpbb model

Model: MODEL2

Dependent Variable: SYST

Number of Observations Read 652

Number of Observations Used 648

Number of Observations with Missing Values 4

Analysis of Variance

Source DF

Sum of

Squares

Mean

Square F Value Pr > F

Model 12 35050 2920.84839 13.67 <.0001

Error 635 135729 213.74606

Corrected Total 647 170779

Root MSE 14.62006 R-Square 0.2052

Dependent Mean 120.63742 Adj R-Sq 0.1902

Coeff Var 12.11901

Parameter Estimates

Variable DF

Parameter

Estimate

Standard

Error t Value Pr > |t|

Intercept 1 107.84005 3.86040 27.93 <.0001

agec 1 0.28131 0.08190 3.43 0.0006

agecSQR 1 0.01586 0.00567 2.80 0.0053

FEMALE 1 0.08291 2.04510 0.04 0.9677

bmic 1 1.19156 0.20785 5.73 <.0001

ddose1 1 -4.25345 2.16149 -1.97 0.0495

ddose2 1 -1.86945 1.81826 -1.03 0.3043

ddose3 1 1.10117 2.01693 0.55 0.5853

ddose4 1 3.20318 2.02661 1.58 0.1145

exdrnk 1 -4.88041 2.60137 -1.88 0.0611

htnmeds 1 17.96450 3.06297 5.87 <.0001

LNpbb 1 3.37845 1.03611 3.26 0.0012

jobdurc 1 -0.23151 0.11121 -2.08 0.0378
LNpbb model with Yhat from first first model included as a covariate

Model: MODEL1

Dependent Variable: SYST

Number of Observations Read 652

Number of Observations Used 648

Number of Observations with Missing Values 4

Analysis of Variance

Source DF

Sum of

Squares

Mean

Square F Value Pr > F

Model 13 35420 2724.64702 12.76 <.0001

Error 634 135359 213.49924

Corrected Total 647 170779

Root MSE 14.61161 R-Square 0.2074

Dependent Mean 120.63742 Adj R-Sq 0.1912

Coeff Var 12.11201

Parameter Estimates

Variable DF

Parameter

Estimate

Standard

Error t Value Pr > |t|

Intercept 1 257.44155 113.67063 2.26 0.0239

agec 1 0.66817 0.30496 2.19 0.0288

agecSQR 1 0.03641 0.01660 2.19 0.0287

FEMALE 1 -0.63912 2.11618 -0.30 0.7627

bmic 1 2.78299 1.22623 2.27 0.0236

ddose1 1 -10.02199 4.88424 -2.05 0.0406

ddose2 1 -4.49176 2.69587 -1.67 0.0962

ddose3 1 2.60440 2.31655 1.12 0.2613

ddose4 1 7.59691 3.90318 1.95 0.0521

exdrnk 1 -11.33151 5.54601 -2.04 0.0414

htnmeds 1 42.24593 18.69133 2.26 0.0241

LNpbb 1 6.53731 2.61275 2.50 0.0126

jobdurc 1 -0.50824 0.23773 -2.14 0.0329

yhat1 1 -1.34328 1.02007 -1.32 0.1884
Pbbc model with Yhat from the LNpbb model included as a covariate

Model: MODEL2

Dependent Variable: SYST

Number of Observations Read 652

Number of Observations Used 648

Number of Observations with Missing Values 4

Analysis of Variance

Source DF

Sum of

Squares

Mean

Square F Value Pr > F

Model 13 35420 2724.64702 12.76 <.0001

Error 634 135359 213.49924

Corrected Total 647 170779

Root MSE 14.61161 R-Square 0.2074

Dependent Mean 120.63742 Adj R-Sq 0.1912

Coeff Var 12.11201

Parameter Estimates

Variable DF

Parameter

Estimate

Standard

Error t Value Pr > |t|

Intercept 1 -111.20302 92.05349 -1.21 0.2275

agec 1 -0.25761 0.23144 -1.11 0.2661

agecSQR 1 -0.01452 0.01311 -1.11 0.2684

FEMALE 1 -0.11098 2.03803 -0.05 0.9566

bmic 1 -1.14033 0.95988 -1.19 0.2353

ddose1 1 4.09282 4.01667 1.02 0.3086

ddose2 1 1.71294 2.32770 0.74 0.4621

ddose3 1 -1.02894 2.19088 -0.47 0.6388

ddose4 1 -2.57394 3.00348 -0.86 0.3918

exdrnk 1 4.74406 4.66265 1.02 0.3093

htnmeds 1 -16.86038 14.31291 -1.18 0.2392

pbbc 1 -0.13376 0.10158 -1.32 0.1884

jobdurc 1 0.20784 0.19659 1.06 0.2908

yhat2 1 1.93500 0.77336 2.50 0.0126