x*******8
发帖数: 188 | 1 一直用SAS，现在用R。 简单的longitudinal data. 两组老鼠吃高脂食物或低脂食物，
每隔一周测一次体重，看体重是否跟进食低脂和高脂有关。
怎么写这个R code，很confused， 谢谢大拿帮忙！
体重Weight ～ 饮食Diet ＋ Weeks ＋ 老鼠ID | s***c
发帖数: 1664 | 2 前阵面试碰到过这个题，我觉得是fit a linear mixed model， diet是fixed effect
， weeks是random effect （random slope and intercept）
就是每一只老鼠，都fit a linear model， random slope and intercept是normal
distributed
我可能想复杂了，weeks可以算作fixed effect。。其实一直没搞明白GEE怎么回事
library（nlme）
fit<-lme(weight～diet*weeks,data=data, random=1+weeks|id, method="ML")
Model Specification：
𝑌𝑖𝑗 = 𝛽0 + 𝛽1 ∗ diet + beta2 &
#8727; weeks + 𝛽3 ∗diet ∗ weeks +bi,0 + 𝑏&#
119894;,1∗ weeks + 𝜀𝑖𝑗
𝛽0: fixed intercept
𝛽1: fixed slope for diet
𝛽2: fixed slope for weeks
𝛽3: fixed slope for diet and weeks interaction
𝑏𝑖,0: random intercept, 𝑏𝑖,0~𝑁(0,
sigma1)
𝑏𝑖,1: random slope for weeks, 𝑏𝑖,1~𝑁
(0, 𝜎2)
𝜀𝑖𝑗: 𝜖𝑖𝑗~𝑁(0, sigma&
#120598;) | x*******8
发帖数: 188 | 3 看到你的详细清晰解答，太感动了！顺祝一定能找到个满意的工作！！
但是我的数据是60个老鼠，11个week。R 一直run，停不下来，我减少到三个week，还
是停不下来。很奇怪，不知是不是我电脑太慢，还是R的问题。我没办法，还是用sas
跑算了！proc mixed 很快，没问题。谢谢了。
effect
&
【在 s***c 的大作中提到】
: 前阵面试碰到过这个题，我觉得是fit a linear mixed model， diet是fixed effect
: ， weeks是random effect （random slope and intercept）
: 就是每一只老鼠，都fit a linear model， random slope and intercept是normal
: distributed
: 我可能想复杂了，weeks可以算作fixed effect。。其实一直没搞明白GEE怎么回事
: library（nlme）
: fit<-lme(weight～diet*weeks,data=data, random=1+weeks|id, method="ML")
: Model Specification：
: 𝑌𝑖𝑗 = 𝛽0 + 𝛽1 ∗ diet + beta2 &
: #8727; weeks + 𝛽3 ∗diet ∗ weeks +bi,0 + 𝑏&#
| d*****0
发帖数: 3 | 4 你这个 Week 是 repeated measures， 你可以用R lme4 package 的lmer 试试，一半
不会出现这个问题的。
【在 x*******8 的大作中提到】
: 看到你的详细清晰解答，太感动了！顺祝一定能找到个满意的工作！！
: 但是我的数据是60个老鼠，11个week。R 一直run，停不下来，我减少到三个week，还
: 是停不下来。很奇怪，不知是不是我电脑太慢，还是R的问题。我没办法，还是用sas
: 跑算了！proc mixed 很快，没问题。谢谢了。
:
: effect
: &
| g**a
发帖数: 2129 | 5 lmer好像不让选variance-covariance matrix structure. 话说有人知道lmer默认的vc
structure是什么吗？
[在 xiaoyu898 (xiaoyu898) 的大作中提到：]
：一直用SAS，现在用R。 简单的longitudinal data. 两组老鼠吃高脂食物或低脂食物
，每隔一周测一次体重，看体重是否跟进食低脂和高脂有关。
：怎么写这个R code，很confused， 谢谢大拿帮忙！
：........... | s***c
发帖数: 1664 | 6 lmer没有default的，得通过code的不同写法来specify
http://stats.stackexchange.com/questions/5344/how-to-choose-nlm
vc
【在 g**a 的大作中提到】
: lmer好像不让选variance-covariance matrix structure. 话说有人知道lmer默认的vc
: structure是什么吗？
: [在 xiaoyu898 (xiaoyu898) 的大作中提到：]
: ：一直用SAS，现在用R。 简单的longitudinal data. 两组老鼠吃高脂食物或低脂食物
: ，每隔一周测一次体重，看体重是否跟进食低脂和高脂有关。
: ：怎么写这个R code，很confused， 谢谢大拿帮忙！
: ：...........
|
|
