[R-sig-ME] Error with glmmADMB and beta distribution
Heather Kharouba
kharouba at zoology.ubc.ca
Tue Jun 12 17:34:47 CEST 2012
Hello,
I am new to glmmADMB and would like to use the beta distribution with
glmmadmb. However, I keep running into the same error. The response
variable is AUC (area under the curve data) from species distribution
models which is continuous and ranges from 0 to 1. The fixed variables are
the number of variables used to build the model (continuous), spatial
extent of the distribution model (continuous), model type (factor) and
taxonomic group (categorical) and the random effect is study.
Here's the model:
model1<-glmmadmb(AUC~variables+log_area+model+taxa+(1|study),
family="beta", verbose=TRUE, data=auc);
A snapshot of the data:
study taxa AUC model variables log_area
Araujo et al. 2005 BIRD 0.9156878 GAM 7 16.21771
Araujo et al. 2005 BIRD 0.9288596 GAM 7 16.21771
Araujo et al. 2005 BIRD 0.9254065 GAM 7 16.21771
Araujo et al. 2005 BIRD 0.8825593 GAM 7 16.21771
Araujo et al. 2005 BIRD 0.9388894 GAM 7 16.21771
Araujo et al. 2005 BIRD 0.9061483 GAM 7 16.21771
When I run the model I get this error:
Error in glmmadmb(AUC ~ variables + log_area + model + taxa + (1 | study),
:
The function maximizer failed (couldn't find STD file)
This is the printout from verbose=TRUE:
- final statistics:
24 variables; iteration 0; function evaluation 0
Function value inf; maximum gradient component mag 0.0000e+00
Exit code = 0; converg criter 1.0000e-04
Var Value Gradient |Var Value Gradient |Var Value Gradient
1 0.00000 nan | 2 0.00000 nan | 3 0.00000
nan
4 0.00000 nan | 5 0.00000 nan | 6 0.00000
nan
7 0.00000 nan | 8 0.00000 nan | 9 0.00000
nan
10 0.00000 nan | 11 0.00000 nan | 12 0.00000
nan
13 0.00000 nan | 14 0.00000 nan | 15 0.00000
nan
16 0.00000 nan | 17 0.00000 nan | 18 0.00000
nan
19 0.00000 nan | 20 0.00000 nan | 21 0.00000
nan
22 0.00000 nan | 23 0.00000 nan | 24 0.00000
nan
- final statistics:
25 variables; iteration 0; function evaluation 0
Function value inf; maximum gradient component mag 0.0000e+00
Exit code = 0; converg criter 1.0000e-04
Var Value Gradient |Var Value Gradient |Var Value Gradient
1 0.00000 nan | 2 0.00000 nan | 3 0.00000
nan
4 0.00000 nan | 5 0.00000 nan | 6 0.00000
nan
7 0.00000 nan | 8 0.00000 nan | 9 0.00000
nan
10 0.00000 nan | 11 0.00000 nan | 12 0.00000
nan
13 0.00000 nan | 14 0.00000 nan | 15 0.00000
nan
16 0.00000 nan | 17 0.00000 nan | 18 0.00000
nan
19 0.00000 nan | 20 0.00000 nan | 21 0.00000
nan
22 0.00000 nan | 23 0.00000 nan | 24 0.00000
nan
25 0.05795 nan |
Hessian type 4
I'm using glmm ADMB version 0.7.2.2 with R version 2.14.2 on a Mac OS X
Version 10.6.8. I've tried including
admb.opts=admbControl(shess=FALSE,noinit=FALSE) and still get the same
error. I'm guessing there's something wrong with either the response
variable or with the overall data structure?
Thanks,
Heather Kharouba
More information about the R-sig-mixed-models
mailing list