2012-06-20 LSmean 自分用メモ res <- glm(y~x1*x2,data=data,family=poisson) dummy <- expand.grid(x1=unique(data$x1),x2=unique(data$x2)) dummy$pres <- predict.glm(res,newdata=dummy)ares <- aggregate(dummy$pres,list(unique(data$x1)),mean) カテゴリカル変数の場合このような感じで出せます。 変数が増えても対応できます。