[R-sig-ME] glmmADMB

Thierry Onkelinx thierry.onkelinx at inbo.be
Thu Nov 12 17:44:10 CET 2015


Dear Veronique,

It looks like your data contains missing values. Try to remove them. Your
syntax seems to be correct.

Note that Treatment + Time + Time * Treatment can be abbreviated to Time *
Treatment

Best regards,

Thierry

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2015-11-11 12:17 GMT+01:00 Veronique Dermauw <vdermauw op itg.be>:

> Dear,
>
> Recently I have been working on a dataset which seems to require a
> generalized linear mixed model with a negative binomial distribution.
> I came across the glmmADMB package which seems ideal.
> However, I am a bit stuck here (see codes and error message below).
> The dataset I am working on, involves count variables (GIS) from animals
> from different villages who were allocated a treatment (Treatment) and were
> sampled at various Time point (hence Time and Time*Treatment).
> I was wondering if you could assist me with this command?
>
> Many thanks for your response,
> Best regards,
> Veronique Dermauw.
>
>
> mod1b<-glmmadmb(GIS~Treatment+Time+Time*Treatment+(1|Animal)+(1|Village),family="nbinom",data=camb)
> summary(mod1b)
> Error in II[, ii] = II[, ii] + REmat$codes[[i]] :
>   number of items to replace is not a multiple of replacement length
> In addition: Warning messages:
> 1: In glmmadmb(GIS ~ Treatment + Time + Time * Treatment + (1 | Animal) +
> :
>   NAs removed in constructing fixed-effect model frame: you should
> probably remove them manually, e.g. with na.omit()
> 2: In II[, ii] + REmat$codes[[i]] :
>   longer object length is not a multiple of shorter object length
>
>
> Veronique Dermauw
> PhD, MSc Veterinary Medicine
> Institute of Tropical Medicine
> Nationalestraat 155
> 2000 Antwerp
> +32(0)32476447
>
>
> P Please consider the environment before printing this e-mail
>
> Disclaimer: Http://www.itg.be/disclaimer
>
> Directions to our location(s): http://g.co/maps/ua89b
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models op r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list