SAS PROC MIXED/GLIMMIX , 3-level model （memory error）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

毕业论文跑下面code，一直出错，上网查了两天（中英文网站查遍!）都没法解决啊TAT心力交瘁，求高手解答!
数据结构是这样的，3-level 重复测量数据
level_1 (participants, n=5000+) is nested within level_2 (school, n=200+), which is nested within level_3 (countries, n=12).

Study objective:
1. random effect for level_1 and level_2;
2. Fixed effect for Level_3, and test the exposure*level_3 interaction.

proc glimmix data=;
class  level_1 level_2  level_3  covariates;
model outcome = exposure  level_3 exposure*level_3 covariates;
random intercept / subject = level_2;
random intercept / subject = level_1 (level_2);
run;


SAS 一直报错：The SAS System stopped processing this step because of insufficient memory.


我的问题如下：

1. 如何解决SAS报错问题啊？

2. 我的SAS codes对吗？level_3虽也有nested, 但是我只放在了Model statement那里。

3. 我换过PROC MIXED，也是一样的memory error，然后网上说换成HPMIXED 可以，我试了，程序跑起来了！但是只能出type 3 analysis fixed effect的结果，但是没有parameter estimates啊…求问高手解决这个问题!

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：proc mixed glimmix memory Level Mixed 毕业论文 memory error 天都

奖励100论坛币！！！

youkuifeng 发表于 2017-5-9 10:13
奖励100论坛币！！！

hpmixed上，随机系数上面要加上solution
random  ....../solution  
这样应该就可以了

dogmamongo 发表于 2017-5-9 12:29
hpmixed上，随机系数上面要加上solution
random  ....../solution  
这样应该就可以了

你好！谢谢回复！

HPMIXED这样估计是可以的~可是除了这个没有别的办法吗？因为这个HPMIXED虽然跑得动，但是一个model要跑一小时啊…我几十个这样的model呢TAT

另外，我想问一下HPMIXED和GLIMMIX结果是一模一样的吗？

dogmamongo 发表于 2017-5-9 12:29
hpmixed上，随机系数上面要加上solution
random  ....../solution  
这样应该就可以了

Hi 我试着跑了你说的在/后面加一个solution，但是出来的是random effect的parameter，但是我想要的是exposure 和covariates 的parameter estimate…你知道这样咋做吗？上网查了也没查到啊TAT

model 那边 最后也加入/solution

dogmamongo 发表于 2017-5-10 13:58
model 那边 最后也加入/solution

Hi， 我加了model后的solution，但是还是很奇怪地，有一个predictor（continuous variable）没有parameter estimate…加之这个HPMIXED跑起来要一个小时TAT，这是终极方法了吗？

另外，我换了一个code：
proc mixed data=;
class  level_1 level_2  level_3  covariates;
model outcome = exposure  level_3 exposure*level_3 covariates;
random intercept / subject = level_1;
random intercept / subject = level_1*level_2;
run;

这样就能跑了，但是还是想请教一下，这样和原来的model比较有区别吗？

我真的要被这个MIXED MODEL折磨疯了…有哪位高手解救一下吗？TAT