我有比例数据,我想用GLMM分析 - 我有一些固定和随机的术语要包含在公式中,但我不明白我是如何构造语法的。 Quad
中的YEAR
和MyData
是我的交叉随机条款(调查年份(总共3个)和调查完成后的四元组(总共4个)。
library(lme4)
# Response variable (No of positives and negatives)
y<-cbind(POSPLANTS,NEGPLANTS)
# Create dataframe with variables from MyData that will become fixed terms:
data <-data.frame(PREC,ALT,INFTYPE,VARIETY,POSPLANTS,NEGPLANTS)
model<-glmer(y~PREC+ALT+INFTYPE+VARIETY+ (1|YEAR) +(1|Quad),binomial,data=data)
summary(model)
是不是?帮助将不胜感激!