全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5795 11
2012-12-19
最近在做线性混合模型,用到lme4
代码如下:

(gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),
family = binomial, data = cbpp))
## GLMM with individual-level variability (accounting for overdispersion)
cbpp$obs <- 1:nrow(cbpp)
(gm2 <- glmer(cbind(incidence, size - incidence) ~ period +
(1 | herd) + (1|obs),
family = binomial, data = cbpp))


但是怎么都找不到data = cbpp的cbpp
请问哪位大侠可以上传一份或者告诉我怎么才可以找得到?
不胜感激
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-1-10 09:33:44
同学 我不能回答你的问题,但是看到你做了混合效应,想请教下随机效应项怎么表述啊,模型如下
Y~a+b*X+c*log10(sqrt(R^2+d^2))+e*Fn+f*Fr 其中,X,R, Fn ,Fr为自变量,Y为因变量,
a,b,c,d,e,f是需要拟合的系数。
X为代表“组”的一个变量,这个模型用混合效应的话,b只能是固定效应了,其余都可以有随机效应。
想请教的是,我让c具有随机效应,用下面的表达就可以
random=c~1|X
但是让c,d,e,f都有随机效应呢,不知道怎么表达啊

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-1-11 14:10:30
我也具体不知道什么意思?
random=c+d+e+f~1
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-1-14 08:35:21
syler 发表于 2013-1-10 09:33
同学 我不能回答你的问题,但是看到你做了混合效应,想请教下随机效应项怎么表述啊,模型如下
Y~a+b* ...
你可以把程序传上来,然后大家一起修改一下
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-1-14 11:35:57
xingzhaoh 发表于 2013-1-14 08:35
你可以把程序传上来,然后大家一起修改一下
EQ<-read.table("data.txt",header=TRUE)
EQ.nlmeR002<-nlme(Y)~a+b*X+c*log10(sqrt(R^2+D^2))+e*Fn+f*Fr,data=EQ,fixed=list(a~1,b~1,c~1,D~1,e~1,f~1),random=c~1|X,start=c(-1,0.8,-2,10,-0.2,0.3),method="ML")
summary(EQ.nlmeR002)
程序很短,,X,R, Fn ,Fr为自变量,Y为因变量,
a,b,c,d,e,f是需要拟合的系数。混合效应是把组与组(不同的X之间)的区别用系数的变化来表示了,我这里只把C为可变的,理论上,a,c,D,e,f,都是可变的。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-1-15 08:49:26
试试
random=c+d+e+f~1
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群