From 7e783e89b8ef28438ec9d20fec99b1125d9a25b0 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 30 Apr 2014 09:07:22 +0200 Subject: Restore t-test (1- and 2-sided) for fitted parameters --- R/mkinfit.R | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) (limited to 'R') diff --git a/R/mkinfit.R b/R/mkinfit.R index 16b4ff3..d3edd2f 100644 --- a/R/mkinfit.R +++ b/R/mkinfit.R @@ -338,14 +338,17 @@ summary.mkinfit <- function(object, data = TRUE, distimes = TRUE, alpha = 0.05, se <- sqrt(diag(covar) * resvar) lci <- param + qt(alpha/2, rdf) * se uci <- param + qt(1-alpha/2, rdf) * se - + tval <- param/se + pval1 <- 2 * pt(abs(tval), rdf, lower.tail = FALSE) + pval2 <- pt(abs(tval), rdf, lower.tail = FALSE) } names(se) <- pnames modVariance <- object$ssr / length(object$residuals) - param <- cbind(param, se, lci, uci) - dimnames(param) <- list(pnames, c("Estimate", "Std. Error", "Lower", "Upper")) + param <- cbind(param, se, lci, uci, tval, pval1, pval2) + dimnames(param) <- list(pnames, c("Estimate", "Std. Error", "Lower", "Upper", + "t value", "Pr(>|t|)", "Pr(>t)")) blci <- buci <- numeric() # Only transform boundaries of CI for one parameter at a time -- cgit v1.2.1