楼主: zhaojing830712
7599 11

急!!!求教:proc mixed语句中的问题?? [推广有奖]

  • 0关注
  • 0粉丝

高中生

12%

还不是VIP/贵宾

-

威望
0
论坛币
68 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
555 点
帖子
17
精华
0
在线时间
27 小时
注册时间
2009-7-19
最后登录
2011-4-13

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
程序如下:
proc mixed  data=bone covtest ic;
class id year group;
model tbmt=group| year age  f_sex1 mc1_a mc2_a ca_a fe_a nation1 baeduc maeduc bmi height/s ;
random int year/ subject=id ;
repeated year / subject=id ;
lsmeans year*group /adjust=bon;
by sex;
run;
日志提示:

NOTE: 122 observations are not included because of missing values.
NOTE: Convergence criteria met.
NOTE: 上述消息是针对以下 BY 组的:
      sex=1
NOTE: 117 observations are not included because of missing values.
NOTE: Convergence criteria met.
NOTE: 上述消息是针对以下 BY 组的:
      sex=2
NOTE: “PROCEDURE MIXED”所用时间(总处理时间):
      实际时间          0.29 秒
      CPU 时间          0.31 秒
二维码

扫码加我 拉你入群

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

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

关键词:proc mixed Mixed mix ROC observations 求教 proc 语句 Mixed

沙发
zhaojing830712 发表于 2010-2-18 14:55:06 |只看作者 |坛友微信交流群
求高手帮忙看看问题出在哪里啊?为什么出不来结果了呢?请大家帮忙啊!

使用道具

藤椅
jingju11 发表于 2010-2-18 23:16:55 |只看作者 |坛友微信交流群
Did you check if variable year is distinct for each record within a subject? if not, a singular R matrix will be introduced.

使用道具

板凳
zhaojing830712 发表于 2010-2-19 11:30:57 |只看作者 |坛友微信交流群
检查了,没问题

使用道具

报纸
zhaojing830712 发表于 2010-2-19 11:39:42 |只看作者 |坛友微信交流群
请教还可能是什么问题啊?

使用道具

地板
醉_清风 发表于 2010-2-19 11:40:25 |只看作者 |坛友微信交流群
进来学习 等待答案
从来不需要想起 永远也不会忘记

使用道具

7
lqyrendajinji 发表于 2010-2-19 11:56:11 |只看作者 |坛友微信交流群
post your dataset

使用道具

8
zhaojing830712 发表于 2010-2-19 12:37:20 |只看作者 |坛友微信交流群
sas的文件不能上传,转成了excel,O(∩_∩)O谢谢!

bone.xls

192.5 KB

使用道具

9
jingju11 发表于 2010-2-19 22:22:29 |只看作者 |坛友微信交流群
5# zhaojing830712

I think the problem may be from the denominator degrees of freedom for the tests of fixed effects. The degrees of freedom were mostly set as 0 in this setting.
When a random statement was specified, it is the default for SAS to compute denominator degrees of freedom by the containment method, which provides correct degrees of freedom only for ‘iid’ designs. Apparently, yours is not this case.
I would suggest not use ‘DDFM = con’ as the default does but choose other settings, such as 'DDFM = KR' in the model statement.

JingJu

使用道具

10
zhaojing830712 发表于 2010-2-20 11:17:20 |只看作者 |坛友微信交流群
多谢指教,如果添加'DDFM = KR' 选项,能运行出结果,
但是日志提示:NOTE: 122 observations are not included because of missing values.
NOTE: Convergence criteria met but final hessian is not positive definite.
NOTE: Asymptotic variance matrix of covariance parameter estimates has been found to be singular and a generalized inverse was used. Covariance parameters
      with zero variance do not contribute to degrees of freedom computed by DDFM=KENWARDROGER.
NOTE: 上述消息是针对以下 BY 组的:
      sex=1
NOTE: 117 observations are not included because of missing values.
NOTE: Convergence criteria met but final hessian is not positive definite.
NOTE: Asymptotic variance matrix of covariance parameter estimates has been found to be singular and a generalized inverse was used. Covariance parameters
      with zero variance do not contribute to degrees of freedom computed by DDFM=KENWARDROGER.
NOTE: 上述消息是针对以下 BY 组的:
      sex=2

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-30 18:04