Poly: Matrix of powers

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

Poly(x, order = 2, exclude = 0)

Arguments

x
order
exclude

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, order = 2, exclude = 0) 
{
    ret <- cbind(rep(1, length(x)))
    for (i in 1:order) ret <- cbind(ret, x^i)
    nam <- deparse(substitute(x))
    powers <- paste(nam, "^", 0:order, sep = "")
    powers[1] <- "Intercept"
    powers[2] <- nam
    dimnames(ret)[[2]] <- powers
    if (!is.null(exclude)) 
        ret <- ret[, -(1:(exclude + 1))]
    ret
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.