多元回归预测区间的误差

Question

这是错误消息：

  

qt（（1-level）/ 2，df）中的错误：数学函数的非数字参数

我要做的是建立一个模型来检查SBP与年龄，性别和种族调整之间的关联。 我的代码使用uwIntroStats包：适合模型的代码有效。性别（male）编码为女性0，男性1，race编码1至4。

library(uwIntroStats)
data(mri)
model <- regress("mean", sbp~age*male+as.factor(race), data = mri)
predict(model, data.frame(age=70,male=0,race=2),interval="prediction")

任何错误发生的原因以及如何解决？谢谢！

Answer 1

您需要命名newdata参数：否则predict方法认为您正在尝试指定下一个不匹配的参数level。来自?predict.uRegress：

  

##类'uRegress'的S3方法
     预测（object，interval =“ prediction”，level = 0.95，...）

所以

predict(model, newdata=data.frame(age=70,male=0,race=2),
        interval="prediction")

有效（您实际上不需要指定interval="prediction"-这是默认值）。