poireg:

In musto101/wilcox_R: What the package does (short line)

Usage

poireg(x, y, xout = FALSE, outfun = outpro, plotit = FALSE, xlab = "X", ylab = "Y", varfun = var, YHAT = FALSE, STAND = TRUE, ...)

Arguments

x
y
xout
outfun
plotit
xlab
ylab
varfun
YHAT
STAND
...

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, y, xout = FALSE, outfun = outpro, plotit = FALSE, 
    xlab = "X", ylab = "Y", varfun = var, YHAT = FALSE, STAND = TRUE, 
    ...) 
{
    xy = elimna(cbind(x, y))
    x <- as.matrix(x)
    x = xy[, 1:ncol(x)]
    y = xy[, ncol(xy)]
    x <- as.matrix(x)
    if (xout) {
        flag <- outfun(x, ...)$keep
        x <- x[flag, ]
        y <- y[flag]
        x <- as.matrix(x)
    }
    temp = glm(formula = y ~ x, family = poisson)
    init = summary(temp)
    yhat = temp$coef[1]
    for (j in 1:ncol(x)) {
        j1 = j + 1
        yhat = yhat + temp$coef[j1] * x[, j]
    }
    yhat = exp(yhat)
    if (plotit) {
        x = as.matrix(x)
        if (ncol(x) > 1) 
            stop("Cannot plot with more than one predictor")
        plot(x, y, xlab = xlab, ylab = ylab)
        xord = order(x)
        lines(x[xord], yhat[xord])
        init$coef
    }
    ex = varfun(yhat)/varfun(y)
    str = sqrt(ex)
    hatv = NULL
    if (YHAT) 
        hatv = yhat
    list(results = init, Explanatory.Power = ex, Strength.Assoc = str, 
        yhat = hatv)
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.