MCMCGlmm";对比度仅适用于具有两个或两个以上级别的因素”;在R

MCMCGlmm";对比度仅适用于具有两个或两个以上级别的因素”;在R,r,R,我正在尝试使用R中的“mcmcglmm”包安装mcmcglmm。 我有一个7个变量的数据框 df <- NNMixedModel2 $Dyad_code $Sex $Parity $Age $Breed $Relatedness_m df可能有帮助 NNMixedModel2$Dyad_code <- as.factor(NNMixedModel2$Dyad_code) NNMixedModel2$Sex <- as.factor(NNMixedModel2$Sex) NN

我正在尝试使用R中的“mcmcglmm”包安装mcmcglmm。 我有一个7个变量的数据框

df <- NNMixedModel2
$Dyad_code 
$Sex
$Parity
$Age
$Breed
$Relatedness_m
df可能有帮助
NNMixedModel2$Dyad_code <- as.factor(NNMixedModel2$Dyad_code)
NNMixedModel2$Sex <- as.factor(NNMixedModel2$Sex)
NNMixedModel2$Parity <- as.factor(NNMixedModel2$Parity)
NNMixedModel2$Breed <- as.factor(NNMixedModel2$Breed)
NNMixedModel2$Age <- as.numeric(NNMixedModel2$Age)
NNMixedModel2$Relatedness_m <- as.numeric(NNMixedModel2$Relatedness_m)
NNMixedModel2$weightlg1p <- as.numeric(NNMixedModel2$weightlg1p)
prior1.1 <- list(G = list(G1 = list(V = 1, nu = 0.002)), R = list(V = 1, nu = 0.002))

model1.1 <- MCMCglmm(weightlg1p ~ Sex+Parity+Age+Relatedness_m, random = ~Dyad_code, pedigree =NULL, rcov=~us(trait):units, family="gaussian", data = NNMixedModel2, nitt = 13000, thin = 10, burnin = 3000, prior = prior1.1, verbose = FALSE)
Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
 contrasts can be applied only to factors with 2 or more levels