Reputation: 11
I want to traslate some code from R to SAS.
In this piece of R code, there is a se.fit = TRUE
theModel <- as.formula(ifelse(timeTrend, "response~seasgroup+wtime",
"response~seasgroup"))
model <- glm(theModel, family = quasipoisson(link = "log"))
pred <- predict.glm(model, data.frame(wtime = c(kk), seasgroup=factor(1)),
dispersion = model$phi, se.fit = TRUE)
eta0 <- pred$fit
seEta0 <- pred$se.fit
I think to use a GENMOD, but what is similar to: se.fit = TRUE
?
I tried this:
proc genmod data=observed_db (where=(t_inversa >= 27));
class seasgroup;
model &model / dist=poisson link=log scale=pearson;
output out=observed_storia
pred=pred
**what is similar?**
leverage=leverage;
store out=previsione;
ods output ModelFit=ModelFit;
run;
proc plm source=previsione;
score data=observed_db out=observed_db pred **what is similar?** / nooffset ilink;
run;
eta0= log(pred);
seEta0= **what is similar?**/eta0 ..?..
Upvotes: 0
Views: 65