R/lmFitStats.R

lmFitStats = function(x,rsq = F,digits = 3, show.signif.stars=TRUE) {
  fitSum = summary(x)
  out = data.frame(fitSum$coefficients)
  names(out) = c("Estimate","S.E.","t-Stat","p-Value")
  out <- round(out,digits)
  if( show.signif.stars) {
    out <- cbind(out, signifStars(out[,"p-Value"]))
    names(out) = c("Estimate","S.E.","t-Stat","p-Value","")
  }
  print(out)
  if(rsq) {
    out = fitSum$r.squared
    out = data.frame(out,row.names = "R-Squared =")
    names(out) = ""
    print(round(out,digits))
  }
}
kecoli/mpo documentation built on May 20, 2019, 8:34 a.m.