binary - 如何打印重复测量逻辑回归的概率？

Question

我希望 SAS 在我的连续自变量的特定实例（从.055 到 2.89，但可以连续记录超过 2.89，并且随着时间的推移总是增加）使用 Proc GENMOD。 如何更改我的代码以打印概率，例如，Calliphoridae 出现在 degree_index=.1？我的例子代码是：

proc genmod data=thesis descending ;
class Body_number ;
model Calliphoridae = degree_index / dist=binomial link=logit ; 
repeated subject=Body_number/ type=cs;
estimate 'degreeindex=.1' intercept 1 degree_index 0 /exp;
estimate 'degree_index=.2' intercept 1 degree_index .1 /exp;run;

我得到对比度估计结果的输出，因为 degree_index=.1 的平均估计是 ..99；degree_index=.2 是 0.98。

Answer 1

我认为它正确地模拟了概率......我只是没有包括平方

度日指数。如果你这样做，它允许概率增加和减少。我

当我手工计算概率时意识到了这一点

(e^-1.1307x+.2119)/(1+e^-1.1307x+.2119) 验证这确实是建模

对于特定 x 值的平均估计值，当 y=1 时的概率......然后我意识到它是

拟合回归线并且不能增加和减少，因为只有

一个 x 值。 http://www.stat.sc.edu/~hansont/stat704/chapter14a.pdf