Nothing
bread <- function(x, ...)
{
UseMethod("bread")
}
bread.default <- function(x, ...) {
nobs0(x) * vcov0(x, ...)
}
bread.lm <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
sx <- summary.lm(x)
return(sx$cov.unscaled * as.vector(sum(sx$df[1L:2L])))
}
bread.mlm <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
cf <- coef(x)
rval <- summary.lm(x)
rval <- kronecker(
structure(diag(ncol(cf)), .Dimnames = rep.int(list(colnames(cf)), 2L)),
structure(rval$cov.unscaled, .Dimnames = rep.int(list(rownames(cf)), 2L)) * as.vector(sum(rval$df[1L:2L])),
make.dimnames = TRUE
)
return(rval)
}
bread.glm <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
sx <- summary(x)
wres <- as.vector(residuals(x, "working")) * weights(x, "working")
dispersion <- if(substr(x$family$family, 1L, 17L) %in% c("poisson", "binomial", "Negative Binomial")) 1
else sum(wres^2)/sum(weights(x, "working"))
return(sx$cov.unscaled * as.vector(sum(sx$df[1L:2L])) * dispersion)
}
bread.nls <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
sx <- summary(x)
return(sx$cov.unscaled * as.vector(sum(sx$df[1L:2L])))
}
bread.polr <- function(x, ...)
{
vcov(x) * x$n
}
bread.clm <- function(x, ...)
{
vcov(x) * x$n
}
bread.survreg <- function(x, ...)
length(x$linear.predictors) * if(is.null(x$naive.var)) x$var else x$naive.var
bread.gam <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
sx <- summary(x)
sx$cov.unscaled * sx$n
}
bread.coxph <- function(x, ...)
{
rval <- x$n * if(is.null(x$naive.var)) x$var else x$naive.var
dimnames(rval) <- list(names(coef(x)), names(coef(x)))
return(rval)
}
bread.hurdle <- function(x, ...)
{
x$vcov * x$n
}
bread.zeroinfl <- function(x, ...)
{
x$vcov * x$n
}
bread.mlogit <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
vcov(x) * length(residuals(x))
}
bread.rlm <- function(x, ...)
{
if(!is.null(x$na.action)) class(x$na.action) <- "omit"
xmat <- model.matrix(x)
xmat <- naresid(x$na.action, xmat)
wts <- weights(x)
if(is.null(wts)) wts <- 1
res <- residuals(x)
psi_deriv <- function(z) x$psi(z, deriv = 1)
rval <- sqrt(abs(as.vector(psi_deriv(res/x$s)/x$s))) * wts * xmat
rval <- chol2inv(qr.R(qr(rval))) * nrow(xmat)
return(rval)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.