varLevel: Identify level of aggregation at which a variable in...

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

Usage

varLevel(x, form, ...)

Arguments

x
form
...

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, form, ...) 
{
    sel <- model.frame(form, x)
    z <- list()
    idx <- rep("", nrow(x))
    z[[1]] <- constant(x, ...)
    for (ii in 1:ncol(sel)) {
        idx <- paste(idx, as.character(sel[[ii]]), sep = ";")
        z[[ii + 1]] <- constant(x, idx, all = T, ...)
    }
    ret <- do.call("rbind", z)
    ret <- length(z) - apply(ret * 1, 2, sum)
    ret
  }

gmonette/spida documentation built on May 17, 2019, 7:25 a.m.