R/calculateLagList.R

In ClusterVAR: Fitting Latent Class Vector-Autoregressive (VAR) Models

calculateLagList <- function(K,
                             HighestLag,
                             LowestLag)
{

    lagList = vector("list", length = K)

    for (lagRunner in 1:K)
    {
        lagList[[lagRunner]] = c(HighestLag:LowestLag)
    }
    lagList = expand.grid(lagList)
    lagList = t(apply(lagList, 1, sort))
    lagList = lagList[!duplicated(lagList), , drop = FALSE]

    if(K==1) lagList <- t(lagList)

    outlist <- lagList[order(rowSums(lagList), decreasing = TRUE) , , drop = FALSE]

    return(outlist)

} # eoF