Fatal编程技术网
  • r/
  • GLMM随机项与固定项R

GLMM随机项与固定项R

r

GLMM随机项与固定项R,r,lme4,mixed-models,R,Lme4,Mixed Models,我有我想用GLMM分析的比例数据-我有一些固定和随机的术语要纳入公式中，但我不明白我是如何构造语法的Quad和YEAR是我的交叉随机项（调查年份（共3个）和进行调查的Quad（共4个） 库（lme4） #响应变量（正、负数量） 你的代码上有几个注释（与你的问题无关）。不要附加数据。只需使用data.frame，而不是data.frame（cbind（），就好像你的任何变量都是字符，那么它将强制所有变量都是字符（包括数字）如果是随机效应，你需要解释变量是什么。可能值得一提的是什么是不正确的；-错误

我有我想用GLMM分析的比例数据-我有一些固定和随机的术语要纳入公式中，但我不明白我是如何构造语法的<
MyData
中的code>Quad和
YEAR
是我的交叉随机项（调查年份（共3个）和进行调查的Quad（共4个）

库（lme4）
#响应变量（正、负数量）

你的代码上有几个注释（与你的问题无关）。不要附加数据。只需使用
data.frame
，而不是
data.frame（cbind（
），就好像你的任何变量都是字符，那么它将强制所有变量都是字符（包括数字）如果是随机效应，你需要解释变量是什么。可能值得一提的是什么是不正确的；-错误？警告？无意义的结果？
glmer
来自
lme4
你将四元和年份建模为嵌套随机效应。它们是嵌套的还是交叉的？Also，固定效果不能同时是随机效果。谢谢你的评论，我在问题中提到了它们。随机效果是交叉的@Roland

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)