cran
本人不喜欢glimmix,因为他用的是quasi-likelihood,我比较欣赏nlmixed的true likelihood.不过好像glimmix能做几个variance components,但是nlmixed好像之可以做一个。
其实啊,作mixed model,还是ADMBRE好啊。给你们介绍一下。
graceli
我已经下载并安装了,发现是6月更新的,我前些时候下时还不行呢,看来SAS动作还是比较快
dudu
我直接问两位吧!
我的研究目的是看年龄、性别、基因型、吸烟、饮酒这些因素是否与一个疾病的发病是否相关。
数据很简单。
因变量为患病和不患病(分别赋值为1和0)。
其他为自变量,基因型分为10种(有的基因型频率较高,有的很低)。年龄为连续变量,不符合正态分布。性别为男性和女性(男性为1, 女性为0)。吸烟为吸与不吸(分别赋值为1和0);饮酒为饮与不饮(分别赋值为1和0)。
请问,这样能否为我写一个nlmixed的程序,万份感谢。
cran
这看来是一个binomial distribution了。但是你这个好像没有variance component啊,直接用genmod就好了。还有,你这个实验有replicate么?否则实验不合格。
proc nlmixed data=yourdata;
eta=b0+b1*age+b2*smoking+b3*drinking;
expeta=exp(eta);
p=expeta/(1+expeta);
MODEL y ~ BINOMIAL(1,p);
RUN;
----------
or
----------
proc genmod data=yourdata;
class smoking drinking
model y = age smoking drinking / dist= binomial link= logit;
run;
cran
还有一个gene的变量,漏了,可以补上
dudu
谢谢cran,我已经使用了您提供的程序,可以运算了,感激之意,难以言表!
cran
proc genmod data=yourdata;
class smoking drinking
model y = age smoking drinking / dist= binomial link= logit;
run;
改正自己的错误,link = logit, 不是log
dudu
请问cran,genmod是个什么模块啊?
graceli
其实还是logistic回归,只不过是换了不同的过程步实现而已
cran
准确地说是GLM
dudu
谢谢两位高人指点!
cran
ls的还没有回答我,你的数据有replicate么?
dudu
cran,昨天和graceli探讨了一下,我的数据应该没有重复测量的数据,那样的话是不是不太适合用nlmixed模块分析啊?
cran
如果没有重复,你的实验根本不是一个及格的设计,可以分析,但是结果漏洞很大。
知道么,在上世纪70-80年代,在生态学上,有18%的研究成果在实验设计上是错误的,然后有48%的结果是沿用这些错误结果的,你说影响有多大?
dudu
cran,我的数据呢不是剂量资料,只是对某个反应的阳性与阴性的问题,不是数值大小的问题。
cran
这个。。。那我就不清楚了
viviyang
我想接着dudu的问题,如果要做Panal (Pool) Data可不可以?我的情况是有Random effect的。我尝试用IML,但是我只会做线性方程的,怎样用IML做Logit方程我不会。
我使用dudu那个问题的程序,多加了一个变量后,系统提示No valid parameter points were found.无法给出估计结果,是什么原因呢,因为我换了很多次变量也不行。
还有就是在nlmixed中怎样能够输出vif, R-Square,以及做异方差的检验呢?
viviyang
或者我说一下我的问题,
我有一组联立方程
(1)Logit(y)=Ln(p/(1-p))=a1+a2X1+a3X3+u1
(2)X1=b1+b2Logit(y)+u2;
y是一个二分类变量。方程中没有标明时间,但实际是包含5年时间的PanalData。
它的非线性导致我在手册有关model程序的讲解中找不到答案。希望给我指点指点。
yueayaya
我是新手,想向高手们请教几个问题.
我用nlmixed做多水平logistic模型,具体代码如下:
proc nlmixed data=sample1;
parms b0=22 b1=-2 b2=0.01 b3=1.49 b4=0.01 b5=8.14 b6=-1.4 b7=-0.77 b8=-8.88 b9=-2.38 b10=-4.68 sigma=1;
logit=b0+u+b1*X1+b2*X2+b3*X3+b4*X4+b5*X5+b6*X6+b7*X7+b8*X8+b9*X9+b10*X10;
p=1/(1+exp(-logit));
model y~binary(p);
random u~normal(0,sigma) subject=block;
run;
其中X1-X10是影响因素,block是水平二.
请问一下,我这个过程编得对吗?
还有就是:上述编程运行之后,只得到了specification表.不是应该还有参数修正表等表格的吗? 是不是nlmixed模块是要另外下载安装才能得出所有的结果呀?
谢谢!
hustlfm
cran对NLMIXED掌握太好了,佩服啊!兄弟我目前有个东西需要运用一下这个过程。首先我对广义非线性模型的理论掌握很是不透彻,在什么条件模型和边缘模型之间晕头转向,这两者之间有何区别?此外,条件期望和方差的计算方法是否有专门的书籍讲述?如何根据具体的问题,设计合适的模型?
我把这个模型运用到设备可靠性评估方面。比如说我已经有8台设备在4年内失效的时间(以小时为单位)数据,简单一点来讲,应该属于计数数据类型,用非时齐泊松分布。cran,能否替兄弟我整个模型出来和具体的程序?有固定效应和随机效应两个东西。
此外,我自己整了个模型,cran,能否帮我实现一下?ADMBRE能做多个方差分量分析吗?
cran,兄弟我盼望着你的回复哦。如果能够在线请教,那就更好了。你的邮箱能否告知?