基于样条函数的R-svyglm预测

我想从公式中包含自然样条线的R survey glm对象进行预测。以以下内容为例：

其中mm是获得预测的模型矩阵-它不使用predict.glm，我认为这是一个很好的理由。我只是想知道是否有一个聪明的解决办法之前，我开始黑客左右。

---

谢谢。

您能在预测之前将它们分开吗。svyglm？：宾果。谢谢这很容易…我想接受你所做的编辑，但我说我需要2K的声誉才能做到这一点。

library(survey)
data(api)
dclus2 <-svydesign(id=~dnum+snum, fpc=~fpc1+fpc2, data=apiclus2)
m <- svyglm(api00 ~ ell + meals + avg.ed, design = dclus2)

m2 <- svyglm(api00 ~ ell + meals + ns(avg.ed,4), design = dclus2)

n <- data.frame(meals=2,ell=2,avg.ed=3)
predict(m1,newdata=n,type="response")
predict(m2,newdata=n,type="response")    # error

mm %*% coef(object)