如何使用svyolr函数获得模型预测

Question

我正在使用svyolr包中的survey函数来运行序数回归模型。我想通过检查模型对我的调查数据和假设数据集的预测来分析模型。

由于predict函数似乎没有svyolr方法，如何从拟合的svyolr模型中获得模型预测？

这里有一些示例数据和一个简单的示例模型可以使用：

library(survey)

# Load example data
data(api)

# Create a survey design object
dclus1 <- svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc)

# Add an ordinal variable to the data
dclus1 <- update(dclus1, mealcat=cut(meals,c(0,25,50,75,100)))

# Fit an ordinal regression model
m <- svyolr(formula = mealcat ~ avg.ed + stype,
            design = dclus1)

# Generate fake data to use in evaluating model predictions
fake_data <- dclus1$variables[,c("avg.ed", "stype")]
fake_data <- subset(fake_data, !is.na(avg.ed))

Answer 1

svyolr基于MASS :: polr，所以我想您可以为此破解一种预测算法。

我怀疑glm.predict :: polr.predict（）生成的置信区间是合法的

var numColumns = 9;