R/summary.JGee1.R In JGEE: Joint Generalized Estimating Equation Solver

Documented in summary.JGee1

```summary.JGee1 <-
function(object, correlation = TRUE, ...)
{
coef <- object\$coefficients
resid <- object\$residuals
n <- length(resid)
p <- object\$rank
if(is.null(p))
p <- sum(!is.na(coef))
if(!p) {
warning("This model has zero rank --- no summary is provided")
return(object)
}
nas <- is.na(coef)
cnames <- names(coef[!nas])
coef <- matrix(rep(coef[!nas], 5), ncol = 5)
dimnames(coef) <- list(cnames, c("Estimate",
"Naive S.E.",  "Naive z",
"Robust S.E.", "Robust z"))
rse <- sqrt(diag(object\$robust.variance))
nse <- sqrt(diag(object\$naive.variance))
coef[,2] <- nse
coef[,3] <- coef[,1]/coef[,2]
coef[,4] <- rse
coef[,5] <- coef[,1]/coef[,4]
summary <- list()
summary\$call <- object\$call
summary\$version <- object\$version
summary\$nobs <- object\$nobs
summary\$residual.summary <- quantile(as.vector(object\$residuals))
names(summary\$residual.summary) <- c("Min", "1Q", "Median", "3Q", "Max")
summary\$model<- object\$model
summary\$title <- object\$title
summary\$coefficients <- coef
summary\$working.correlation <- object\$working.correlation
summary\$scale <- object\$scale
summary\$error <- paste("Error code was", object\$error)
summary\$working.correlation1 <- object\$working.correlation1
summary\$working.correlation2 <- object\$working.correlation2
summary\$working.correlation <-  object\$working.correlation
summary\$iterations <- object\$iterations
if ( correlation ) {
##	rob.var <- object\$robust.variance
##	nai.var <- object\$naive.variance
##	summary\$robust.correlation <- rob.var /
##	outer(sqrt(diag(rob.var)),sqrt(diag(rob.var)))
##	dimnames(summary\$robust.correlation) <- list(object\$xnames,object\$xnames)
##	summary\$naive.correlation <- nai.var /
##	outer(sqrt(diag(nai.var)),sqrt(diag(nai.var)))
##	dimnames(summary\$naive.correlation) <- list(object\$xnames,object\$xnames)
}
attr(summary,"class") <- "summary.JGee1"
summary
}
```

Try the JGEE package in your browser

Any scripts or data that you put into this service are public.

JGEE documentation built on May 2, 2019, 11:51 a.m.