使用lm()函数对整个数据框进行交互回归
我知道在 R 中有一种快捷方式可以对像这样的数据框运行一个 lm() 回归： reg<-lm(y~.,data=df) 使用 df 作为自变量 x1、x2、...、x5 的解释变量，所以这与写成以下形式相同： With df having explanatory variables x...
使用`lm`的线性模型：如何获取预测值之和的预测方差
我正在对多个预测变量的线性模型进行预测值求和，就像下面的例子一样，并希望计算这个总和的组合方差、标准误差和可能的置信区间。lm.tree <- lm(Volume ~ poly(Girth,2), data = trees) 假设我有一组名为Girths的集合：newdat <- ...
使用“与平均值的差异”模型来处理分类变量上的所有系数，需要使用“对比编码”吗？
假设我们想做一个简单的“收入描述模型”。 假设我们有三个群体，北部、中部和南部（考虑美国地区）。 在比较其他相似的群体时，假设北部的平均收入为130，中部为80，南部为60。 假设群体大小相等，因此平均值为90。 在（线性回归）模型中，应该有一种方法来报告系数与总体平均数（在多元情况下，“所...
如何从线性模型（lm）中预测x值
我有这样的数据集：x <- c(0, 40, 80, 120, 160, 200) y <- c(6.52, 5.10, 4.43, 3.99, 3.75, 3.60) 我使用lm()计算了一个线性模型：model <- lm(y ~ x) 我想知道如果我有新的y值，例如yn...
在函数内使用局部协变量更新模型的update()函数
我需要在函数内更新一个回归模型。理想情况下，该函数应适用于任何类型的模型（lm、glm、multinom、clm）。更确切地说，我需要添加一个或多个在函数内定义的协变量。以下为示例：MyUpdate <- function(model){ randData <- data...
如何在不影响传递给predict()函数的情况下，最小化lm类对象的大小？
我想在一份包含2个预测变量和5000万行数据的大型数据集上运行`lm()`。这个分析是在只有10GB空间存储数据的远程服务器上运行的。我从数据中取出了一万行进行了测试，结果生成的对象大小为2GB+。 我仅需要`lm()`函数返回的"lm"类对象来生成模型的摘要统计信息(`summary(l...
在回归中仅与一个因子水平交互的R语言操作
在回归模型中，是否可以仅使用一个因子的虚拟变量与交互作用？例如，假设我有以下数据： x: numerical vector of 3 variables (1,2 and 3) y: response variable z: numerical vector 有没有可能构建像这样的模型：...
如何正确地将由`lm`拟合的线性模型通过`dput`函数保存到ASCII文件中，并在后续重新创建该模型？
我想将lm对象保存到文件中，并在另一个程序中重新加载它。 我知道可以使用saveRDS/readRDS通过写入/读取二进制文件来实现，但我希望有一个ASCII格式的文件。更一般地说，我想知道为什么我通常用于读取dput输出的习惯用法没有按照我的预期那样工作。 以下是制作简单拟合的示例以及成功...
运行lm在不同变量上的函数
我希望创建一个函数，可以在给定数据集的不同变量上运行回归模型（例如使用lm）。在该函数中，我将指定使用的数据集，因变量y和自变量x作为参数。我希望这是一个函数而不是循环，因为我想在脚本的各个位置调用代码。我的函数看起来像这样： lmfun <- function(data, y, x)...
在R中，选择因子在lm回归中的基础类别的最佳方法是什么？
假设我想使用lm和一个factor作为右侧变量运行回归。选择哪个水平作为基础类别（被排除以避免多重共线性）的最佳方法是什么？请注意，我不想排除截距，因为我有很多因素。 我还希望有一个基于公式的解决方案，而不是直接在数据框上操作的解决方案，但如果您认为您有一个非常好的解决方案，请也发表一下。 ...