Question

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

因变量=冷自变量=品种（类别2），性别（类别4），Dayno。（cols 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
 }

Answer 1

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

要提取p值，可以执行以下操作：

summary(x1)$coefficients[,4]

$ p.value函数等效于$ r.squared函数

1 个答案: