如何在R中实现

proc nlmixed data=DS;
odds=exp(b0) * (1+ b1*alcohol + b2*smoke + b3*alcohol*smoke) ;
model case ~ binary( odds/(1+odds) );
run;
R中如何实现SAS这个过程？谢谢！

你是在做logistic回归么？我不太懂sas，如果你是做logistic的话，这样打：
a=read.table(file.choose())
b=glm(y~alcohol+smoke+alcohol*smoke,binomial,a)
summary(b)
odd=exp(这里面是你各项估计值)/(1+exp(这里面是你各项估计值))

justice1989 发表于 2013-4-30 23:34
你是在做logistic回归么？我不太懂sas，如果你是做logistic的话，这样打：
a=read.table(file.choose())
...

这样是不对的！

你是要做什么啊！

justice1989 发表于 2013-4-30 23:54
你是要做什么啊！

SAS中用的是非线性混合模型，R中如何实现呢？

nlme

qoiqpwqr 发表于 2013-5-2 08:34
nlme

odds=exp(b0) * (1+ b1*alcohol + b2*smoke + b3*alcohol*smoke) ;
model case ~ binary( odds/(1+odds) );
能否按照这个写得完整些呢？谢谢！