Question

我希望SAS在我的连续自变量（“degree_index”）的特定实例中打印我的二元因变量发生的概率（“Calliphoridae”一个特定的苍蝇家族存在（1）或不存在（0））记录从.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是.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

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

1 个答案: