Lag.0: Lag within subject: older less efficient version of Lag

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

Usage

Lag.0(x, id, idx, lag = 1, at = NULL, check = T)

Arguments

x
id
idx
lag
at
check

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, id, idx, lag = 1, at = NULL, check = T) 
{
    if (check) {
        comb <- paste(id, idx)
        if (any(duplicated(comb))) 
            stop("id not unique in each level of idx")
    }
    if (any(is.na(idx))) 
        stop("NAs in idx")
    if (any(round(idx) != idx)) 
        stop("Non integers in idx")
    ret <- x
    id <- as.character(id)
    names(x) <- id
    for (i in max(idx):min(idx)) {
        to.pos <- idx == i
        if (is.null(at)) 
            from <- x[idx == (i - lag)]
        else from <- x[idx == at]
        ids <- names(x[to.pos])
        ret[to.pos] <- from[ids]
    }
    ret
  }

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