r - 我想“引导”（读取平均值）我的回归结果

Question

问题

所以我想做的很简单。我想运行 10 次（或 100 次）回归，并以某种方式组合或平均结果。显然，通常多次运行回归不会改变结果，但对我来说确实如此（见下文）。

背景

您可以使用包sure（链接）来计算 a 的代理残差polr。resids*问题是这些代理残差需要从连续分布中采样；因此，每次调用残差的结果都会有所不同（它们不会从引导残差中收敛）。

为了规避这个问题，我想要一种“引导”整个交易，直到（与残差相反）模型收敛（这可能要求太多，但可以说我想得到一个好的平均值）。

数据

library(sure) # for residual function and sample data sets
library(MASS) # for polr function

df1 <- df1
df1$z1 <- df1$x
df1$x <- NULL
df1$x <- df2$y
df1$z2 <- df2$x
df1$y <- df3$x

mod1 <- polr(as.ordered(x) ~ z1 + z2, data=df1, method='probit')
df1$residuals <- resids(mod1, nsim=100, method="latent")
mod2 <- lm(y ~ x + residuals + z1 + z2, data=df1)
summary(mod2)

.*该软件包基于美国统计协会杂志上的这篇论文。

标签： rstatistics-bootstrap