R/uncommented.r

Defines functions aic.shglm aic.gaussian aic.poisson aic.inverse.gaussian aic.binomial aic.Gamma ll.speedglm

ll.speedglm <- function(family,aic.model,nvar){
  switch(family,
            binomial = -(aic.model-2 * nvar)/2,
            Gamma = -((aic.model-2 * nvar)-2)/2,
            gaussian = -((aic.model-2 * nvar)-2)/2,
            poisson = -(aic.model- 2 * nvar)/2,
            inverse.gaussian = -((aic.model-2 * nvar)-2)/2,
            quasi = NA,
            quasibinomial = NA,
            quasipoisson = NA)
}


aic.Gamma <- function(y, n, mu, wt, dev) {
  disp <- dev/n
  -2*sum(dgamma(y,1/disp,scale=mu*disp,log=TRUE)*wt)#+2
}

aic.binomial<-function(y, n, mu, wt, dev){
    m <- if (any(n > 1)) n else wt
    -2 * sum(ifelse(m > 0, (wt/m), 0) * dbinom(round(m * y),
        round(m), mu, log = TRUE))
}

aic.inverse.gaussian <- function(y, n, mu, wt, dev)
sum(wt) * (log(dev/sum(wt) * 2 * pi) + 1) + 3 * sum(log(y) *
    wt)# + 2

aic.poisson <- function(y, n, mu, wt, dev)
-2 * sum(dpois(y, mu, log = TRUE) * wt)

aic.gaussian <- function(y, n, mu, wt, dev){
    n.obs <- length(y)
    n.obs * (log(dev/n.obs * 2 * pi) + 1) - sum(log(wt)) # + 2
}


safe_pchisq <- function (q, df, ...) #from stats
{
  df[df <= 0] <- NA
  pchisq(q = q, df = df, ...)
}

safe_pf <- function (q, df1, ...) #from stats
{
  df1[df1 <= 0] <- NA
  pf(q = q, df1 = df1, ...)
}


aic.shglm <- function(family,y, n, mu, wt, dev){
  options(warn=-1)
  a<-switch(family,
            binomial = aic.binomial(y, n, mu, wt, dev),
            Gamma = aic.Gamma(y, n, mu, wt, dev),
            gaussian = aic.gaussian(y, n, mu, wt, dev),
            poisson = aic.poisson(y, n, mu, wt, dev),
            inverse.gaussian = aic.inverse.gaussian(y, n, mu, wt, dev),
            quasi = NA,
            quasibinomial = NA,
            quasipoisson = NA)
 options(warn=1)
 return(a)
}


model.frame.speedglm <- function (formula, ...) 
{
  if (is.null(formula$offset)||is.null(formula$call[["offset"]])) model.frame.default(formula) else
    if(is.null(formula$model)) cbind(model.frame.default(formula),"(offset)"=formula$offset) else
      formula$model
}

model.frame.speedlm <- model.frame.speedglm  

Try the speedglm package in your browser

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

speedglm documentation built on March 18, 2022, 6:48 p.m.