# STEPWIZE: 0
# lhs: dead
# rhs:  1 date score genet height distance density 

# formula(AIC = 3992.23): dead ~ 1 + date + score + genet + height + distance + density
# formula(AIC = 4799.13): dead ~ 1 + score + genet + height + distance + density
# formula(AIC = 4098.30): dead ~ 1 + date + genet + height + distance + density
# formula(AIC = 4034.11): dead ~ 1 + date + score + height + distance + density
# formula(AIC = 3992.51): dead ~ 1 + date + score + genet + distance + density
# formula(AIC = 3989.43): dead ~ 1 + date + score + genet + height + density
# formula(AIC = 3990.73): dead ~ 1 + date + score + genet + height + distance

# model		AIC
 distance	3989.43
  density	3990.73
   <none>	3992.23
   height	3992.51
    genet	4034.11
    score	4098.30
     date	4799.13
#-----------------------------------------------------------------------------

# STEPWIZE: 1
# lhs: dead
# rhs:  1 date score genet height distance density 

# formula(AIC = 3989.43): dead ~ 1 + date + score + genet + height + density
# formula(AIC = 4797.59): dead ~ 1 + score + genet + height + density
# formula(AIC = 4093.23): dead ~ 1 + date + genet + height + density
# formula(AIC = 4035.91): dead ~ 1 + date + score + height + density
# formula(AIC = 3993.48): dead ~ 1 + date + score + genet + density
# formula(AIC = 3988.51): dead ~ 1 + date + score + genet + height

# model		AIC
  density	3988.51
   <none>	3989.43
   height	3993.48
    genet	4035.91
    score	4093.23
     date	4797.59
#-----------------------------------------------------------------------------

# STEPWIZE: 2
# lhs: dead
# rhs:  1 date score genet height distance density 

# formula(AIC = 3988.51): dead ~ 1 + date + score + genet + height
# formula(AIC = 4795.90): dead ~ 1 + score + genet + height
# formula(AIC = 4091.14): dead ~ 1 + date + genet + height
# formula(AIC = 4041.16): dead ~ 1 + date + score + height
# formula(AIC = 3990.19): dead ~ 1 + date + score + genet 

# model		AIC
   <none>	3988.51
   height	3990.19
    genet	4041.16
    score	4091.14
     date	4795.90
#-----------------------------------------------------------------------------

# formula(AIC = 3988.51): dead ~ 1 + date + score + genet + height

Call:  glmmML(formula= as.formula(formula.str), data = data, cluster = cluster,      family = binomial(logit), start.sigma = start.sigma) 


                  coef se(coef)        z Pr(>|z|)
(Intercept) -14.050877 0.528961 -26.5632 0.00e+00
date          0.061208 0.002555  23.9598 0.00e+00
score         0.467188 0.043647  10.7038 0.00e+00
genetG        1.415387 0.227585   6.2192 5.00e-10
genetI        0.622016 0.236062   2.6350 8.41e-03
genetQ       -0.579595 0.273484  -2.1193 3.41e-02
genetc1       4.411689 0.353970  12.4634 0.00e+00
genetc2      -0.097885 0.377600  -0.2592 7.95e-01
height       -0.007128 0.002866  -2.4867 1.29e-02

Standard deviation in mixing distribution:  1.386 
Std. Error:                                 0.08303 

Residual deviance: 3969 on 7978 degrees of freedom 	AIC: 3989