SAS に、連続独立変数の特定のインスタンスで、2 項従属変数 (「カリフォリダエ」という特定のハエ科が存在する (1) または存在しない (0)) 確率を出力してもらいたい (「degree_index」は、 .055 から 2.89 ですが、2.89 を超えて継続的に記録することができ、時間の経過とともに常に増加します) Proc GENMOD を使用します. コードを変更して確率を出力するにはどうすればよいですか? たとえば、カリフォリダエが 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 です。