: 都是在做 cumulative logit model : model 应该是一样的 : 只是response超过两类 若加以covariates也超过两类 : 不同procedure的预设表示法可能不同 : 也就是它们可能采用不一样的 reference group 之类的 我来自我解答了 ^^a 谢谢学姊的提醒 经过了测试之後发现真的是reference group 设定的问题 以性别为例 gender=0,1 proc logistic的预设是 beta_gender0 + beta_gender1 = 0 ( 即 gender0 = 1, beta_gender1 = -1 ) 然後会列出 beta_gender0 的估计系数 proc genmod则是跟我们习惯使用的dummy设法一样 令beta_gender1 = 0 ( 即 gender0 = 1, beta_gender1 = 0 ) 然後列出 beta_gender0 的估计系数 所以如果以这个例子来说的话 这两个模型的beta_gender0刚好会是两倍的关系 ( 但是level>2的话 就不一定罗 ) ☆☆ 让 genmod 和 logistic 可以得到相同的参数估计值的方法 ☆☆ proc logistic ; class gender / PARAM = REFERENCE ; model Y = X ; run; ---------------------------------------------------- 真的笨哪~~~ 一个指令搞两天 Orz"" --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 140.109.74.44