付録 M5実習5: ロジスティック回帰分析 統計ソフトRの使い方

多多多多多多多

多多 II. 多多多多多多多多多多多

M5 多多多多多多

2011 多 2 多 23 多 多多多多多多多多多多多多多多多多 多多多多多多多多多多多多 多多多多多多多多多多多多多多

多多多多多多多多多多多

) ,

, ,

,

(x₁ x₂ x₃ x₄ x₅ f

y 

多多多多多多多多多多多多多多

多多多











 ₁x_{1 2} x_{2 3} x₃

y    

,  ,

,

, x₁ x₂ x₃ y

,  ,

,

, ₁ ₂ ₃



多多多多多多多多多 多多多多多多多多

↑ 多多多多多多多多多多多多多多多多多

↑ 多多多多多多多多多多多多多多多多

多多多多多多多多多多多

) ,

, ,

,

(x₁ x₂ x₃ x₄ x₅ f

y 

多多多多多多多多多多多多多多

多多多











 ₁x_{1 2} x_{2 3} x₃

y    

,  ,

,

, x₁ x₂ x₃ y

,  ,

,

, ₁ ₂ ₃



多多多多多多多多多 多多多多多多多多

↑ 多多多多多多多多多多多多多多多多多

↑ 多多多多多多多多多多多多多多多多

↑ 多多多多多多多多多多多 多多多多多多多多_{「 」}

多多多

多多多 多多多多多多多 多多多_{「 」} 多多多多多多多

2 多多多多多多多多多多多多 多多多多多多

多多多多多 2 多多多多多多多多多多多多多多 多多多多多多多多多多多多多

多多多多多多多多多多多多

) ,

, ,

,

(x₁ x₂ x₃ x₄ x₅ f

p_y 

x1

x2

x3

x4

x5

多多多 B

0 / 1 y 

x1

x2

x3

x4

x5

多多多 A

多多多多多多 多多多多多多

多多多多多多多多多 ^{多 多多多多多}

多 多多多多多

py

多多多多多多多多多多多多

) ,

, ,

,

(x₁ x₂ x₃ x₄ x₅ f



) (  _{1 1}  _{2 2}  _{3 3}  

 f   x  x  x















 

3 3 2

2 1

1

1_{( p ) x x x}

f _y    

y y

p p 1 log

py

y

y y

p p 1

2 多多多多 1 or 0 0 ≦ 多多 ≦ 1

0 ≦ 多多多 ≦ +∞ 多∞ ≦ 多多多多 ≦ +∞

多多多多多多多

多多多多多多多多多 多多多多多

多多多多多多多多多 多多多多多

多多多

909 多多多多多多多多多多多多多多多多多多多 多多 30kg 多多多多多多「 2 多多多多 多多多」

多多多多 30kg 多多多 多多多多多多多多多多多多多_{「 」}

多多 logit 多多多 30kg 多多多多 α+β₁ 多多 +β₂ 多 多多多多多多多多多多多 1 多多多多 0

α 多 β₁ 多 β₂ 多多多多多多多多多多多多

2

1 ^2.567

205 .

0 32.626

log1 x x

p p

y

y   



多多多多多多多多多多多多多多多多多

↑

多多多多多多多多多多多

2

1 ^2.567

205 .

0 32.626

log1 x x

p p

y

y   



567 .

2 ,

205 .

0 ,

32.626 ₁  ₂ 

  



多多多 β

多多多多多多多多多多多多多

多多多多多多 K 多多多 L 多多多多多多多

K 多多多多多多多多 Tcm L 多多多多多多多多 Tcm

L 多多多多多多多多 Tcm K 多多多多多多多多 Tcm

2

1 ¹

) 1

/

log( p_yK  p_yK    T  

2

1 ⁰

) 1

/

log( p_yL  p_yL    T  

多

K 多多多 L 多多多多多多多多







 ) log( /1 ) 1

/

log( p_yK p_yK p_yL p_{yL 2}

β₂ 多多多多多多多多多多多多多多

2

1 ^2.567

205 .

0 32.626

log1 x x

p p

y

y   

 ^log1 p ^{1x1 2 x2} p

y

y     



多多多多多多多多多多多

多多多多多

yL yL yK

yK

p p p

p

 

 ^log1 log1

多多多多多多多多多













 

 

yL yL yK

yK

p p p

p

1 log 1

K 多多多 多 多多多多多多 L 多多多 多 多多多多多多 多多 多多多多

多 log ( 多多多多 )

多多多多

多多多多多多多多

多多多

β

多多多多多多多多多多多多多多多 !!

多多多多多多多多多多多

多多多多多 多多多多多多多多

多多多

β

多多多多多多多多多多多多多多多 !!

多多多多多多多多多多 30kg 多多多多多多多多多多多多多

多多多多 0.0768 多多多

多多多多多多多多多多多

多多多多多多多多多多多多多多多 !!

0.0768 )

567 .

2 exp(

)

exp(₂   

2

1 ^2.567

205 .

0 32.626

log1 x x

p p

y

y   



567 .

2 ,

205 .

0 ,

32.626 ₁  ₂ 

  



β

多多多多多多多多多多多多多

L 多多多 M 多多多多多多多

L 多多多多多多多多 Tcm

M 多多多多多多多多 T-1cm

多多多多多多多多多多多

M 多多多多多多多多 T-1cm L 多多多多多多多多 Tcm

2

1 ⁰

) 1

/

log( p_yK  p_yK    T  

2

1^{( 1) 0}

) 1

/

log( p_yL  p_yL    T   

多

L 多多多 M 多多多多多多多多







 ) log( /1 )

1 /

log( p_yL p_yL p_yM p_yM ₁ β₁ 多多多多多多 1cm 多多多多多多多多多

2

1 ^2.567

205 .

0 32.626

log1 x x

p p

y

y   

 ^log1 p ^{1x1 2 x2} p

y

y     



多多多多多多多多多多多

多多多多多 1 多多多多多多多多多多多多 β₁ 多多多多多 1 多多多多多多多多多多多 e^β1

0.815 )

205 .

0 exp(

)

exp(₁   

多多多多多多多多多多多 1 多多多多多多多多多多

多多多 30kg 多多多多多多多多多多多 0.815 多 多多多多多多多多多 10 多多多多多多

多多多 30kg 多多多多多多多多多多多多多多多多

多多多多多多多多多多多

0.129 )

05 .

2 exp(

) 10 205

. 0

exp(    

多多多多多多多多多多多多多多

多多多 kg 多 多多多多

多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多多 2 多多多多多多多多多多多多多多多多多多 多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多

多多多多多多多多 2500g 多多多多多多多

多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多 β 多 多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多 多多多多多多多多 2 多多多多多多多多多多多 多多多多多多多多多多多多多多多多

Subjects and odds ratios according to lifetime physical activity history

多多多 多多 ^Crud_{e OR} ^(95% _CI) ^Adjusted OR*

(95% CI) 多多多多

多多多 33 26 3.8 ^(2.0-_7.2) 3.7 ^(1.7-_7.7)

多多多 37 112

多多多多多多

多多多 42 58 2.1 (1.2-_3.7) 1.2 (0.5-_2.9)

多多多 28 80

*: adjusted for working status, age at first

pregnancy, smoking, alcohol, HRT, family history, lactation and BMI

2.1 1.2 Crude Adjusted

多多多多多多多多多多多多多多多 多多多多多多多多多多多多多

多多多多 多多多多多多多多多_{「 」}

多多多多多多多多多

Crude 多多多多多多多多多多多多多多 Adjusted 多多多多多

多多多多多多多多多多

多多多多多多多多多多多多

) ,

, ,

,

(x₁ x₂ x₃ x₄ x₅ f

p_y 

) ( y p

x1

x2

x3

x4

x5

多多多 B

0 / 1 y 

x1

x2

x3

x4

x5

多多多 A

多多多多多多 多多多多多多

多多多多多多多多多 ^{多 多多多多多}

多多多多多多 166cm 多多多多多多多 30kg 多多多多多多多多多

9708 .

3 567

. 2 166

205 .

0 32.6262

log1     

 _y

y

p p

0.01885 )

9708 .

3

1 _y ^exp(^{ }

y

p p

0.01850 1885

. 0 1

01885 .

0 

  py

多多

多 risk 多

多多多多 多多多

多多 166cm 多多多多多多多 30kg 多多多多多多多多 0.01885

多多多多多多多多多

) ,

, ,

,

(x₁ x₂ x₃ x₄ x₅ f

y 











 ₁x_{1 2} x_{2 3} x₃

y    

↑ 多多多多多多多多多多多 多多多多多多多多_{「 」}

多多多 多多多多多多多 多多多_{「 」} 多多多多多多多

2 多多多多多多多多多多多多 多多多多 多多多多多多多多多多多多_{「 」}

多多多多

多多多多多多多多多多多多多多多多多 risk 多

→ 多多多多多→多多多多多

多多多多多多多多多多多

)

1(^t S

)

2 (^t S

多多多多多多多多

多多多多多多多多

多多多多

)

1(^t S

)

2 (^t S

) ( 1 )

(

1 1

t S dt

t

 dS

) ( 1 )

(

2 2

t S dt

t

 dS

多多 多 多多多多多多多 多 rate 多 hazard 多多多 多多多多多多多多

多多多

多多多

多多多多多多多多

多多 多 多多多多多多多多多多多多

Constant hazard ratio

Proportional hazard

Cox’s model

2 多多多多多多多多多多多多多

多多多多多多多多多多多多多多

多多多多多多多多多多多多 β 多多多多多多多

多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多

多多 多 多多多多多多多多多多多多

多多多多多多多多多多多多多多多

多多多多多多多 vs 多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 2 多多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多多多多多

多多多多多多多 vs 多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多

多多多多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多

多多多多 2 多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多

多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多多多多多多多多多多

多多多多 2 多多多多多多多多多多多多多多多多多多多多多多多多多 多多多多多多

多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多

多多

多多多 6 多多多多多多多多多多多多 7 多多多多多多多多多 多多多多多多多多多多多多 1cm 多多多多多多多 30kg 多 多多多多多多多多多多多多多多多多多多多多多多多多多

多多多多多多多多多多多多多多多多多多多多多多多多 多多多 A4 多多多 1 多多多多多多多多多多多多