以下是我的问题（为了可重现性而使用虚构数据）： set.seed(42) df<-data.frame("x"=rnorm(1000),"y"=rnorm(1000),"z"=rnorm(1000)) df2<-data.frame("x"=rnorm(100),"y"=rnor...

在R语言中，使用lm()函数进行留一交叉验证的命令很简单吗？ 具体来说，是否有一个简单的命令可以对下面的代码进行操作？ x <- rnorm(1000,3,2) y <- 2*x + rnorm(1000) pred_error_sq <- c(0) for(i in ...

我不太明白为什么我在计算贝叶斯信息准则时遇到了麻烦，希望有人能指点我正确的方向。我这样做是因为我正在尝试手动计算BIC（对于plm对象，它们似乎没有与之相关联的已建立例程）。我使用了维基百科页面中的公式，该公式给出了以残差平方和为基础的BIC公式，而不是对数似然函数。 y<-rnorm...

plot.lm()如何确定哪些点是异常值（即要标记哪些点）以进行残差与拟合图？在文档中，我唯一发现的内容是： “详情”： 默认情况下，子标题（即函数调用）会显示为每个图的副标题（在x轴标题下方），当图分别位于单独的页面上时，或者作为外边距中的副标题（如果有多个图在同一页上）。 “比例位置”图...

我刚开始学习R语言，希望通过使用*apply函数来改进下面的脚本（我已经阅读了关于apply的内容，但是我无法使用它）。我想在多个自变量上（这些自变量是数据框中的列）使用lm函数。 for (i in (1:3) { assign(paste0('lm.',names(data[i]))...

我有以下数据集，想要对每个国家和地区进行线性回归，然后将预测值添加到数据集中： 添加三列后的最终数据框如下： 我已经为一个国家和一个地区完成了上述操作，但现在希望能对每个国家和地区都做一遍，并通过cbind将预测值、上限值和下限值放回到数据集中： data <- data...

如果我在R中使用lm和factor命令运行固定效应模型，如何在stargazer模型中抑制因子变量系数？ 也就是说，我的模型是： m1<-lm(GDP~pop_growth + factor(city)) 我只想报告有关人口增长的截距和系数，而不是每个城市虚拟变量的系数。 编辑...

我有一个包含未来收益的向量Y和一个包含当前收益的向量X。由于最后一个当前收益也是可用系列的末尾，因此最后一个Y元素为NA。 X = { 0.1, 0.3, 0.2, 0.5 } Y = { 0.3, 0.2, 0.5, NA } Other = { 5500, 222, 523, 3677 ...

出于某种原因，当我指定GLMs（以及LMs，事实证明也是如此）时，R无法预测数据中的缺失值。以下是一个例子： y = round(runif(50)) y = c(y,rep(NA,50)) x = rnorm(100) m = glm(y~x, family=binomial(link="...

我一直在尝试弄清楚R语言中lm()函数中的subset参数是如何工作的。特别是下面的代码对我来说似乎有些可疑： data(mtcars) summary(lm(mpg ~ wt, data=mtcars)) summary(lm(mpg ~ wt, cyl, data=mtcars)...