R中更全面的summary()函数？答案

【问题标题】：More comprehensive summary() function in R?R中更全面的summary()函数？
【发布时间】：2015-01-26 04:09:27
【问题描述】：

我想知道 R 中是否有更全面的 summary() 函数，它可能包含更多模型指标，例如估计的置信区间，可能是对数似然、AIC、BIC 之类的东西。我知道调用模型对象上的其他函数来获取这些指标非常容易，但如果它们都打包在摘要函数中，那就太酷了！

【问题讨论】：

适合什么型号？ summary() 调用对象特定的方法（例如，summary.lm()、summary.glm()、summary.merMod() 等）。您也可以编写自己的函数。
如果你有一个lm 对象然后调用summary()
是的，但是您能否从 summary() 中获得参数估计和其他模型统计信息的置信区间
置信区间，使用confint()。
您的问题留下了许多未确定的问题，这就是为什么它不是默认返回。例如，您没有指定所需的置信区间。什么阿尔法水平？您想要 est 系数向量的每个元素的置信区间还是整个向量（或其子集）的置信区域？等如果您想要一组特定的返回值，那么您可以编写自己的摘要函数和打印方法。

标签： r modeling diagnostics

【解决方案1】：

不一定更全面，但我倾向于编写自己的摘要函数。然后我可以定制特定的应用程序。我似乎总是需要“n”，它不在任何可用的摘要中。

sumstats<-function(x) {
  n<-NROW(x)
  mean<-mean(x)
  max<-max(x)
  min<-min (x)
  SD<- sd(x)
  hist<-hist(x)
  summary<-list(n=n,mean=mean,max=max,min=min,SD=SD)
  return(summary)
}

【讨论】：

Cool thxs 我在模型周围做了更多，但这是我的函数：summaryRobust

【解决方案2】：

这更像是我在摘要函数中所希望的：

summaryRobust <- function(model){
    regOutput <- cbind( as.data.frame(summary(model)[4]), as.data.frame(confint(model)) )
   colnames(regOutput) <- c("Estimate", "Std. Error", "T-Value", "P-Value", "C.I. 2.5%", "C.I. 97.5%")

   rSq <- summary(model)$r.squared
   Adj_rSq <- summary(model)$adj.r.squared
   modelFit <- rbind(rSq, Adj_rSq)
   rownames(modelFit) <- c("R-squared", "Adjusted R-squared")
   colnames(modelFit) <- c("")

   metrics <- as.data.frame(cbind(AIC(model), BIC(model), logLik(model)))
   colnames(metrics) <- c("AIC", "BIC", "logLik")
   rownames(metrics) <- c("")

   par(mfrow=c(2,2))
   plot(model)
   out <- list(Model_Output = regOutput, Model_Fit = modelFit, Other_Metrics = metrics)
   return(out)
}

【讨论】：