ABCDEF
1
参考)モデル9-0　上と似た定式化｡初期値を求めるため普通のポアソン
2
#修正前#修正後
3
#Deviance Residuals: e50qspにミスがある結果Deviance Residuals:
4
# Min 1Q Median 3Q Max Min 1Q Median 3Q Max
5
#-3.5880 -0.3992 -0.1902 -0.0590 4.2178 -3.5211 -0.3987 -0.1897 -0.0590 4.2166
6
#
7
#Coefficients:Coefficients:
8
# Estimate Std. Error z value Pr(>|z|) Estimate Std. Error z value Pr(>|z|)
9
#(Intercept) -7.172030 0.660423 -10.860 < 2e-16 ***(Intercept) -7.20925 0.66315 -10.871 < 2e-16 ***
10
#mhiro 3.842157 0.813981 4.720 2.36e-06 ***mhiro 3.78676 0.82469 4.592 4.40e-06 ***
11
#mnaga 3.823731 0.814157 4.697 2.65e-06 ***mnaga 3.76685 0.82488 4.567 4.96e-06 ***
12
#fhiro -0.042810 0.031714 -1.350 0.177047 fhiro -0.04250 0.03172 -1.340 0.180290
13
#me30 0.034250 0.032232 1.063 0.287956 me30 -0.05567 0.04266 -1.305 0.191892
14
#fe30 0.237340 0.030863 7.690 1.47e-14 ***fe30 0.20624 0.03986 5.174 2.29e-07 ***
15
#me30sq 0.006345 0.020339 0.312 0.755080 me30sq -0.03060 0.02354 -1.300 0.193727
16
#fe30sq 0.043690 0.019457 2.246 0.024735 * fe30sq 0.03169 0.02216 1.430 0.152636
17
#me30qsp -0.050273 0.028880 -1.741 0.081730 . me30qsp 0.05037 0.04114 1.224 0.220820
18
#fe30qsp -0.129221 0.027760 -4.655 3.24e-06 ***fe30qsp -0.08398 0.03848 -2.183 0.029068 *
19
#me50qsp 5.743046 6.529469 0.880 0.379099 me50qsp -0.19177 0.05966 -3.214 0.001308 **
20
#fe50qsp 11.911646 6.777319 1.758 0.078820 . fe50qsp -0.07145 0.05519 -1.295 0.195408
21
#mlage70 9.191060 1.588280 5.787 7.17e-09 ***mlage70 8.96947 1.63889 5.473 4.43e-08 ***
22
#flage70 -2.516992 2.252687 -1.117 0.263854 flage70 -2.60487 2.26167 -1.152 0.249426
23
#mlage70sq 2.267354 1.178584 1.924 0.054380 . mlage70sq 2.14842 1.22350 1.756 0.079093 .
24
#flage70sq -4.797241 1.802920 -2.661 0.007795 ** flage70sq -4.86015 1.81065 -2.684 0.007270 **
25
#mlage40qsp -3.561543 1.521291 -2.341 0.019225 * mlage40qsp -3.31870 1.56647 -2.119 0.034125 *
26
#flage40qsp 5.905774 2.104683 2.806 0.005016 ** flage40qsp 5.99867 2.11297 2.839 0.004526 **
27
#mlage70qsp -7.294940 2.019405 -3.612 0.000303 ***mlage70qsp -7.33271 2.01459 -3.640 0.000273 ***
28
#flage70qsp -2.011561 1.667643 -1.206 0.227729 flage70qsp -2.12168 1.66660 -1.273 0.202997
29
#log(1 + x) 0.474301 0.075389 6.291 3.15e-10 ***log(1 + x) 0.47122 0.07538 6.251 4.08e-10 ***
30
#e30:log(1 + x) -0.141205 0.039364 -3.587 0.000334 ***e30:log(1 + x) -0.14078 0.03930 -3.582 0.000341 ***
31
#log(1 + x):lage70 -0.541013 0.277653 -1.949 0.051353 . log(1 + x):lage70 -0.53865 0.27766 -1.940 0.052385 .
32
#log(1 + x):dsex 0.268900 0.097120 2.769 0.005627 ** log(1 + x):dsex 0.27069 0.09712 2.787 0.005319 **
33
#------
34
#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
35
#(Dispersion parameter for poisson family taken to be 1)(Dispersion parameter for poisson family taken to be 1)
36
# Null deviance: 27566 on 37059 degrees of freedom Null deviance: 27566 on 37059 degrees of freedom
37
#Residual deviance: 14251 on 37036 degrees of freedomResidual deviance: 14242 on 37036 degrees of freedom
38
#AIC: 25818AIC: 25810
39
40
41
42
43
モデル9
44
#修正前#修正後
45
[[1]]
46
# beta.estimated ses t　　　 e50qspにミスがある結果 beta.estimated ses tcalc.AIC(res8)
47
#(Intercept) -7.165854978 0.66052568 -10.8487152(Intercept) -7.20114548 0.66308663 -10.860037
48
#mhiro 3.843617740 0.81395249 4.7221648mhiro 3.78935443 0.82399285 4.598771
49
#mnaga 3.825561178 0.81413257 4.6989413mnaga 3.76984999 0.82418426 4.574038
50
#fhiro -0.039656734 0.03169809 -1.2510765fhiro -0.03939173 0.03170574 -1.242416
51
#me30 0.032055022 0.03209524 0.9987468me30 -0.05684162 0.04259556 -1.334449
52
#fe30 0.242526210 0.03093974 7.8386642fe30 0.21281297 0.03998091 5.322864
53
#me30sq 0.001079688 0.02041160 0.0528958me30sq -0.03531584 0.02359363 -1.496838
54
#fe30sq 0.036261118 0.01972119 1.8386879fe30sq 0.02496813 0.02235298 1.116993
55
#me30qsp -0.045610809 0.02890742 -1.5778235me30qsp 0.05365397 0.04115247 1.303785
56
#fe30qsp -0.123707746 0.02788874 -4.4357590fe30qsp -0.08020425 0.03852640 -2.081800
57
#me50qsp 5.645137226 6.53019793 0.8644665me50qsp -0.18927995 0.05968431 -3.171352
58
#fe50qsp 11.847541323 6.77840600 1.7478359fe50qsp -0.06817465 0.05521142 -1.234792
59
#mlage70 9.153976368 1.58726141 5.7671511mlage70 8.94142250 1.63463851 5.469969
60
#flage70 -2.556711579 2.25295986 -1.1348234flage70 -2.63884525 2.26133457 -1.166942
61
#mlage70sq 2.235315865 1.17884465 1.8961921mlage70sq 2.12190266 1.22083589 1.738074
62
#flage70sq -4.853128264 1.80496636 -2.6887638flage70sq -4.91288034 1.81215840 -2.711066
63
#mlage40qsp -3.565495464 1.52069230 -2.3446528mlage40qsp -3.33159492 1.56302819 -2.131500
64
#flage40qsp 5.911118941 2.10514256 2.8079424flage40qsp 5.99872530 2.11289049 2.839108
65
#mlage70qsp -7.219684296 2.01919430 -3.5755273mlage70qsp -7.25067676 2.01427381 -3.599648
66
#flage70qsp -1.947305481 1.66837153 -1.1671894flage70qsp -2.05687939 1.66741196 -1.233576
67
#exp(e30) -0.378706946 0.09354455 -4.0484126exp(e30) -0.37602715 0.09371096 -4.012627
68
#exp(lage70) -0.687420776 0.45260803 -1.5187994exp(lage70) -0.69497460 0.45241326 -1.536150
69
#x:dsex 0.254029054 0.09712233 2.6155577x:dsex 0.25724104 0.09736845 2.641934
70
#x 0.470076049 0.06256625 7.5132529→　これが線量の係数　x 0.46774879 0.06249616 7.484440
71
#[[2]]
72
#[1] -12883.27[[2]]
73
#[1] -12879.08
74
#[[3]]
75
#[1] 25814.55 モデル9.00よりもAICは小さいのでこちらの方がよい｡｡｡｡[[3]]
76
[1] 25806.16