Message-ID: <4B2249F8.1060501@mnhn.fr>
Date: 2009-12-11T13:32:40Z
From: Celine Teplitsky
Subject: zipoisson in MCMCglmm
Dear all,
I am trying to run a zero inflated model with MCMCglmm as the number of
days of parades in the birds is overdispersed.
I used this code:
priorA <- list(R=list(V=diag(2),n=2,fix=2),
G=list(G1=list(V=diag(2),n=2),G2=list(V=diag(2),n=2)))
NBJ_ZIP<- MCMCglmm(NBJ~trait, random=~idh(trait):animal+idh(trait):IDb,
rcov=~idh(trait):units,data=Data ,
family="zipoisson",pedigree=PedR,prior=priorA)
I got this error message:
Error in if (max(data_tmp$MCMC_y) == 1) { :
missing value where TRUE/FALSE needed
and I don't manage to understand what I should change. I'm sorry if this
was already answered but I could not find it.
Many thanks in advance
All the best
Celine