を通じてベータのMLEを取得します。 model <- glm(y~x, data=d1, family = poisson(link="log"))
summary(model)
私は次の出力を得る:私はそれのための関数を書きたい Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 1.3948
マイコード proc reg data=mydata plots=;
title "Correlation Between Estriol Levels and Birthweights";
model birthweight=estriol/clb cli clm;
run;
プロットの回帰直線と散布図から信頼バンド及び予測バンドを除去しようとするだけでなく、信頼帯及び