问题标签 [r-car]

我正在使用 qqplot 比较两组生成的数据，它们都适合伽马。如何将 95% 的置信范围添加到图中？

我正在尝试编写一个循环来计算方差通货膨胀因子。我知道有一些功能和包可以为我做到这一点，但我需要某种定制。

样本数据

该数据框有 30 列（预测变量）。

这是我的循环逻辑：

1) 将每个预测变量与其他预测变量进行回归并计算 R2。使用 VIF = 1/1 - R2 将 R2 转换为 VIF。这会给我 30 个 VIF 值。

2）对VIF值进行排序。如果顶部预测变量的 VIF > 10，则从cov.mat. cov.mat现在将有 29 个预测变量。

3) 重复第 1 步，即将每个预测变量与其他预测变量进行回归并再次计算 VIF（这次是 29 个 VIF）。如果最大 VIF > 10，则删除具有最高 VIF 的变量并继续执行直到最大 VIF <= 10。

但是，问题是我想保留 X4 、 X6 和 X10 ，即使它们在给定迭代中的 VIF > 10 也是如此。因此在上述过程中，如果 X4 或 X6 或 X10 在一次迭代中具有最高的 VIF（> 10），则删除具有第二高 VIF 的变量（仅当第二高的 VIF 也 > 10 并且不是 X4 或X6 或 X10)。我希望这很清楚

假设上面的循环工作正常，我有一个矩阵，第一列为变量名称，第二列为 VIF 值。

这就是我迷路的地方。

我首先需要检查具有最高 VIF 的变量是否 > 10 并且不在 X4 或 x6 和 X10 之间，然后从 dataframe 中删除该变量cov.mat。如果 VIF 最高的变量（给定 VIF > 10）是 X4 或 X6 或 X10，则转到第二行df并评估其 VIF > 10 以及是否不在 X4、X6 或 X10 之间，如果它满足条件，将其从 中删除cov.mat并重新开始迭代。

编辑

我的原始数据框有 51 列和 1458 行。当我运行上述函数时，它给了我一个错误there are aliased coefficients in the model。为什么会这样？

在 Ubuntu 的 Rstudio 中安装“car”包时，我遇到了一个错误，说没有安装依赖项“nloptr”、“lme4”等。所以我开始安装“nloptr”，但我遇到了这个错误——

系统在最后一行之后挂起，当我在一段时间后重新启动时，没有安装“nloptr”，因此我无法安装“car”包。请帮我。提前致谢。

Anova()我已经构建了一个简单的线性模型，并使用包中的函数执行了 II 型和 III 型方差分析car。在输出中，类型 IIIAnova包含截距，而类型 II 中不包含截距Anova。

Ⅱ型

第三类

这背后的原因是什么？如果我想使用 II 类测试，如何确定模型截距是否显着？

我想使用包car中的函数vif（测试多重共线性）。包“car”已成功解包，但加载library(car)后，出现错误：

错误：inDL(x, as.logical(local), as.logical(now), ...) 中“car”的包或命名空间加载失败：无法加载共享对象“C:/Users/lenovo/Documents/ R/win-library/3.4/Rcpp/libs/x64/Rcpp.dll': LoadLibrary failure: 找不到指定的程序。

我找到了类似问题的旧解决方案：升级/重新安装 R 和包。没有任何效果。它不适用于版本 3.4.4 或 3.5.0。在 3.5.0 版本中是另一个问题，没有安装 data.table 包（这是包（汽车）的一部分）

错误：loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) 中的“car”的包或命名空间加载失败：没有包称为“数据表”。

data.table的问题已在 R 3.5 中的 data.table 包中解决，但未安装。

那么，3.4.4 版中的 package(car) 是否有问题的解决方案，或者我应该使用 3.5.0 版吗？但是该功能是否可用？

我正在尝试将多个（四个）散点图组合成一个图形（但不是散点图矩阵）。scatterplot()我正在使用包的功能制作单个散点图car。我曾经能够使用layout()orpar()函数来组合这四个图。但是，现在当我尝试在 Rstudio 中执行此操作时，它只是按顺序显示四个图。我不确定这是否是因为 R 或 Rstudio 的较新版本。

以下是使用 mtcars 数据集的示例：

我在 Windows 10 上运行 R 3.4.2、RStudio 1.1.453。任何指针将不胜感激。

我的问题是如何在 2 型 Anova（混合效应模型）上运行事后处理？到目前为止，我正在使用glmer()“lme4”包中Anova()的“car”包，并尝试从“agricolae”包中运行 HSD 测试。

在搜索了一段时间后，这是我能找到的最好的，但是，这样做时我收到一条错误消息。有谁知道如何解决这个问题或我做错了什么？或者另一种方式来做到这一点？

出现的错误信息：Error in HSD.test(totaldiv.anova, "focalspecies", group = TRUE, console = TRUE, : argument "MSerror" is missing, with no default

先感谢您！

只需运行以下代码。我想显示一个坐标网格，但没有任何反应：

从文档??car::scatter3d中我意识到

在回归曲面上绘制网格线（真或假）。

因此，该grid参数不是我想要的。有没有办法获得坐标网格？对我来说，这对眼睛的指导非常有用。

在 Carles 输入后编辑：

我想保留交互式图表 - 这就是我正在寻找car::scatter3d解决方案的原因。如果你不需要这个，scatterplot3d和的组合FactoClass真的很好。以下以非交互方式工作：

我经常在库（汽车）中使用“重新编码”功能来重新编码变量中的级别。直到今天我的代码都运行良好，但现在它给我带来了错误。df 等没有任何变化，不确定发生了什么。
可能有人可以启发我！

我的数据框（示例）：

我的实际数据集中的级别：-

我的代码（直到今天都运行良好！）：-

错误：-Error: Argument 2 must be named, not unnamed

我需要一点帮助，我是 R 的新手。我正在使用 plspm 包生成 PLS-SEM (PLS-PM) 模型。该软件包不提供 VIF 函数来评估潜在变量（构造）的指标（变量）的共线性

我想创建一个可以在循环中自动执行此操作并报告 VIF 值的函数。我手动完成了整个 VIF 计算。

在下面的示例中，我有 2 个潜在变量。第一个是具有 21 个指标（X1 到 X21）的制度压力 (IP)。第二个是具有 16 个指标（X22 到 X37）的战略响应 (SR)。

谢谢你！

知识产权

SR

IP VIF

SR VIF