We now fit the random effects model with the explanatory variables we found significant in our previous analysis.
data case move age year dur ed ch1 ch2 ch3 ch4 msb mse esb ese &
osb ose mbu mrm mfm msb1 epm eoj esb1 ops osb1 msb2 esb2 osb2 osb3
read rochmig.dat
6349 observations in dataset
yvar move
transform tempage age - 30
transform trage tempage / 10
transform trage2 trage * trage
transform trage3 trage2 * trage
transform trage4 trage3 * trage
transform trage5 trage4 * trage
transform trage6 trage5 * trage
transform ldur log dur
factor ch3 fch3
factor mbu fmbu
factor mrm fmrm
factor mfm fmfm
factor ops fops
factor esb2 fesb2
factor osb3 fosb3
factor msb2 fmsb2
C fit simple logistic main effects model
lfit int ldur year trage trage2 trage3 trage4 trage5 trage6 &
fch3 fesb2 fmbu fmrm fmfm fops fmsb2 fosb3
Iteration Deviance Reduction
__________________________________________
1 8801.5829
2 2935.4172 5866.
3 2266.9758 668.4
4 2124.9723 142.0
5 2105.4389 19.53
6 2103.4563 1.983
7 2103.0820 0.3743
8 2103.0417 0.4028E-01
9 2103.0411 0.5907E-03
10 2103.0411 0.1442E-06
dis est
Parameter Estimate S. Error
___________________________________________________
int 0.98308 0.68858
ldur -0.99954 0.73925E-01
year -0.38417E-01 0.73821E-02
trage 0.44814 0.32946
trage2 -0.18073 0.61760
trage3 -0.58213 0.52585
trage4 0.35071 0.56434
trage5 0.43121 0.20699
trage6 -0.23961 0.15231
fch3 ( 1) 0. ALIASED [I]
fch3 ( 2) -1.2225 0.74642
fesb2 ( 1) 0. ALIASED [I]
fesb2 ( 2) 0.55741 0.32354
fesb2 ( 3) 1.0932 0.39493
fmbu ( 1) 0. ALIASED [I]
fmbu ( 2) 1.2346 0.64645
fmrm ( 1) 0. ALIASED [I]
fmrm ( 2) 1.5114 0.68339
fmfm ( 1) 0. ALIASED [I]
fmfm ( 2) 0.46040 0.23846
fops ( 1) 0. ALIASED [I]
fops ( 2) 1.1069 0.28439
fmsb2 ( 1) 0. ALIASED [I]
fmsb2 ( 2) -0.46453 0.17487
fosb3 ( 1) 0. ALIASED [I]
fosb3 ( 2) 0.83115 0.16524
C fit the same model with random effects
fit .
Iteration Deviance Step End-points Orthogonality
length 0 1 criterion
________________________________________________________________________
1 2135.3134 1.0000 free free 4.6620
2 2128.4000 0.2500 free free 4.4184
3 2123.0042 0.4288 free fixed 0.21340E-01
4 2122.5984 0.0078 free fixed 7.1550
5 2088.6586 1.0000 free free 4.0641
6 2079.1069 1.0000 free free 4.5894
7 2075.6823 1.0000 free free 24.604
8 2075.6458 1.0000 free free 17.341
9 2075.6458 1.0000 free free
dis est
Parameter Estimate S. Error
___________________________________________________
int 0.73017 0.89590
ldur -0.63527 0.10783
year -0.37769E-01 0.10902E-01
trage 0.17360 0.34893
trage2 -0.16613E-01 0.61773
trage3 -0.54783 0.53830
trage4 0.28880 0.56026
trage5 0.40754 0.21096
trage6 -0.22024 0.15128
fch3 ( 1) 0. ALIASED [I]
fch3 ( 2) -1.3073 0.75078
fesb2 ( 1) 0. ALIASED [I]
fesb2 ( 2) 0.31615 0.37617
fesb2 ( 3) 0.77441 0.45042
fmbu ( 1) 0. ALIASED [I]
fmbu ( 2) 1.2513 0.66612
fmrm ( 1) 0. ALIASED [I]
fmrm ( 2) 1.5259 0.71835
fmfm ( 1) 0. ALIASED [I]
fmfm ( 2) 0.45266 0.25126
fops ( 1) 0. ALIASED [I]
fops ( 2) 1.2016 0.30209
fmsb2 ( 1) 0. ALIASED [I]
fmsb2 ( 2) -0.56390 0.19405
fosb3 ( 1) 0. ALIASED [I]
fosb3 ( 2) 0.68677 0.18610
scale 0.49269 0.18099
PROBABILITY
___________
end-point 0 0.48867 0.19067 0.32760
end-point 1 0.29991E-02 0.43654E-02 0.20105E-02
stop
|
Next: Interpretation of results |
Home page | Contents | Previous |