$p.value 函数等效于 $r.squared 函数[重复]答案

【问题标题】：$p.value function equivalent to the $r.squared function [duplicate]$p.value 函数等效于 $r.squared 函数[重复]
【发布时间】：2019-06-05 03:24:04
【问题描述】：

我有如下代码。这运行了许多多重回归（都略有不同）。它还在一个漂亮的表格中报告了这些回归的 R2 值，我可以轻松地复制和粘贴该表格。这是在代码末尾实现的。但是，我似乎无法以类似的方式获得 p 值表，据我所知，没有与 $r.squared 函数等效的 p 值函数。对此的任何帮助将不胜感激。

因变量=冷自变量 = 品种（类别，第 2 列），性别（类别第 4 列），Dayno。（第 6-1055 列，这是每个回归之间的差异）。

r.sqr <- matrix(NA, nrow =1, ncol = length(6:ncol(rweights)))

 for(i in 6:ncol(rweights)){
   na.n <- length(which(is.na(as.numeric(unlist(rweights[,i])))))
   if (na.n == nrow(rweights)) next
 nn <- which(!is.na(as.numeric(unlist(rweights[,i]))))
   if (is.na(unlist(rweights[nn,2])) && is.na(unlist(rweights[nn,4]))) next
   fit <- lm(Cold~ Breed + Gender + as.numeric(unlist(rweights[,i])), data=rweights, , na.action=na.omit)
   r.sqr[1,(i-5)] <- summary(fit)$r.squared
 }

【问题讨论】：

尝试运行anova(fit) 并从中提取。但是不，不能保证summary.lm 中的p.value 插槽

标签： r

【解决方案1】：

样本数据：x <- lm(hp ~ wt, mtcars)

要提取您可以执行的 p 值：

summary(x1)$coefficients[,4]

【讨论】：

summary(fit) 给我一个包含 p 值的输出，但意味着每次都手动执行。当前代码在全局环境中生成一个名为 r.qr 的新对象，它是 1055 个不同回归中每个回归的 R2 值的列表/表格。
你可以考虑在你的函数中嵌入summary(x1)$coefficients[,4]，就像你对r.sqr[1,(i-5)] <- summary(fit)$r.squared所做的那样。
是的。毫无乐趣地尝试各种迭代。 x1 在您的代码中指的是什么？编辑..越来越近（我认为）
这只是我引用的对象，是 lm() 调用的结果。所以你应该改用 fit 。抱歉没有注意到我同时使用了 x 和 x1，应该是一样的。
谢谢。我使用“fit”位只是想仔细检查发生了什么问题不是因为我不小心引用了错误的东西！