如何简洁地从数据框中写出具有多个变量的公式？
假设我有一个响应变量和一组包含三个协变量的数据（作为玩具示例）：y = c(1,4,6) d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2)) 我想对数据进行线性回归拟合：fit = lm(y ~ d$x1 + d$x2 ...
在R中的线性回归和按组分组
我想使用lm()函数在R中进行线性回归。我的数据是一个年度时间序列，其中一列是年份（22年），另一列是州名（50个州）。我希望为每个州都拟合一个回归模型，以便最终得到一个lm响应向量。我可以想象为每个州使用for循环，然后在循环内部进行回归，并将每个回归的结果添加到一个向量中。但这似乎不太像R...
如何从带名称的数字中提取数字（不包括名称）？
我只需要 B1(newx) 线性模型系数的值，不需要名字。我只需要0.5的值，不需要名字"newx"。newx <- c(0.5,1.5,2.5) newy <- c(2,3,4) out <- lm(newy ~ newx) out的样子：Call: lm(formul...
在使用predict.lm进行预测时，出现“'newdata'有1行，但发现32行变量”的警告。
在使用R的predict和lm函数时，我发现一个奇怪的问题。对于相同的数据，使用数据框和向量作为输入会得到不同的结果。 数据框代码：data(mtcars) fitCar<-lm(mtcars$mpg~mtcars$wt) predict(fitCar, data.fr...
提取回归系数值
我有一个用于研究药物利用的时间序列数据回归模型。目的是将样条拟合到时间序列上，并计算95% CI等参数。该模型如下： id <- ts(1:length(drug$Date)) a1 <- ts(drug$Rate) a2 <- lag(a1-1) tg <- ts....
将最大似然估计得到的系数放入stargazer表格中
Stargazer可以为lm（和其他）对象生成非常漂亮的latex表格。假设我已经通过最大似然拟合了一个模型。我想让stargazer为我的估计值生成类似于lm的表格。我该怎么做？ 虽然有点hacky，但一种方法可能是创建一个“假”的lm对象，其中包含我的估计值--只要summary（my....
Predict() - 或许我没有理解它
我今天早些时候发布了一个问题，关于使用predict函数时出现的错误。我已经纠正了这个错误，并且认为我走在了正确的道路上。 我有许多观测值（实际值），并且我想要推断或预测一些数据点。我使用lm创建了一个模型，然后尝试使用实际值作为预测输入来使用predict。 这段代码和我的先前帖子中重复，但...
使用R语言从线性模型中提取R平方值
我知道使用 summary 可以帮助我手动计算，但是我需要计算大量的R平方值。因此，我需要电脑为我提取它。以下是一个简单的例子：library(alr3) M.lm=lm(MaxSalary~Score,data=salarygov) #Here you will see the R squa...
如何调试“对具有2个或更多级别的因子才能应用对比”的错误？
这是我正在使用的所有变量：str(ad.train) $ Date : Factor w/ 427 levels "2012-03-24","2012-03-29",..: 4 7 12 14 19 21 24 29 31 ...
R：在predict()中，numeric 'envir'参数长度不为一
我正在尝试使用predict()函数通过将变量传递到模型中来预测R中的值。 但是我遇到了以下错误：Error in eval(predvars, data, env) : numeric 'envir' arg not of length one 这是我的 数据框，名称为df:df &l...