**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 wou ld 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:Kotchma r.Dennis@epamail.epa.gov]

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

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

We need to have all references on file here. I am w orking 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.govPBBc 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.0716LNpbb 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.0378LNpbb 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.1884Pbbc 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