Refer to full name df.residual in fastLm summary (closes #115)

Author: eddelbuettel

ChangeLog

@@ -1,3 +1,8 @@
+2022-11-04  Dirk Eddelbuettel  <[email protected]>
+
+	* R/fastLm.R (summary.fastLm,print.fastLm): Refer to correct and full
+	variable name df.residual in the returned object
+
 2022-09-15  Jonah Gabry  <[email protected]>
 
 	* inst/skeleton/rcppeigen_hello_world.cpp: Correct typo

R/fastLm.R

@@ -54,7 +54,7 @@ summary.fastLm <- function(object, ...) {
     object$coefficients <- cbind(Estimate     = coef,
                                  "Std. Error" = se,
                                  "t value"    = tval,
-                                 "Pr(>|t|)"   = 2*pt(-abs(tval), df=object$df))
+                                 "Pr(>|t|)"   = 2*pt(-abs(tval), df=object$df.residual))
 
     ## cf src/stats/R/lm.R and case with no weights and an intercept
     f <- object$fitted.values
@@ -66,7 +66,7 @@ summary.fastLm <- function(object, ...) {
     object$r.squared <- mss/(mss + rss)
     df.int <- if (object$intercept) 1L else 0L
     n <- length(f)
-    rdf <- object$df
+    rdf <- object$df.residual
     object$adj.r.squared <- 1 - (1 - object$r.squared) * ((n - df.int)/rdf)
     class(object) <- "summary.fastLm"
     object
@@ -82,7 +82,7 @@ print.summary.fastLm <- function(x, ...) {
 
     printCoefmat(x$coefficients, P.values=TRUE, has.Pvalue=TRUE, ...)
     cat("\nResidual standard error: ", formatC(x$s, digits=digits), " on ",
-        formatC(x$df), " degrees of freedom\n", sep="")
+        formatC(x$df.residual), " degrees of freedom\n", sep="")
     cat("Multiple R-squared: ", formatC(x$r.squared, digits=digits),
         ",\tAdjusted R-squared: ",formatC(x$adj.r.squared, digits=digits),
         "\n", sep="")