[R-sig-ME] MCMCglmm binomial model

François Dumont dumont.francois.3 at courrier.uqam.ca
Wed Mar 23 19:43:43 CET 2011


I'm trying to run a binomial model with a MCMCglmm code that I found on 
the R-sig-ME list (see below). Unfortunately, the following error 
message occurs when I send the code to R. What I am doing wrong ?

error message : Erreur dans `[.data.frame`(data, , 
match(response.names[0:nJ + nt], names(data))) :
   undefined columns selected

prior <- list(R=list(V=1, nu=0.002), G=list(G1=list(V=1, nu=1, 
alpha.mu=0, alpha.V=100),G2=list(V=1, nu=1, alpha.mu=0, alpha.V=100),
G3=list(V=1, nu=1, alpha.mu=0, alpha.V=100)))

model1 = MCMCglmm(bik~
, data=sp,nitt=503000,thin=500,burnin=3000,verbose=FALSE, prior=prior, 
pr=TRUE, family="multinomial2", DIC=TRUE)

bik is a binomial variables (an animal exhibited or not a behaviour) 
with value being 0 or 1. I need to extract the specific sample of each 
individual so that's why pr=TRUE.


François Dumont.

Candidat au doctorat en Biologie
Laboratoire de lutte biologique - Dr. Éric Lucas
Université du Québec à Montréal

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