paramVertexOnlyGroup function

Parameter estimation for Vertex model only for a list of dynamic networks.

paramVertexOnlyGroup(
  InputNetwork,
  VertexStatsvec = rep(1, nvertexstats),
  maxLag,
  VertexModelGroup = NA,
  VertexLag = rep(1, maxLag),
  VertexAttLag = rep(1, maxLag),
  VertexLagMatrix = matrix(1, maxLag, length(VertexStatsvec)),
  regMethod = "bayesglm"
)

Arguments

• InputNetwork: Input network list.
• VertexStatsvec: Binary vector of size 9, indicating vertex model.
• maxLag: maximum lag.
• VertexModelGroup: Group term for vertex model.
• VertexLag: Binary vector of size maxLag, indicating Lag terms in the model.
• VertexAttLag: Vertex group term lag vector.
• VertexLagMatrix: Binary matrix indicating lagged vertex statistics in the model.
• regMethod: one of "glm", "glmnet", "bayesglm"

Returns

List of 3 elements:

VertexFit: Output from regEngine.

VertexStats: Subsetted vertex stats matrix.

VertexStatsFull: Full matrix of vertex stats.

Examples

nvertexstats <- 9
InputNetwork <- beach
maxLag <- 3
VertexStatsvec <- rep(1, nvertexstats)
VertexLag <- rep(1, maxLag)
regMethod <- "bayesglm"
VertexModelGroup <- "regular"
VertexLagMatrix <- matrix(0, maxLag, nvertexstats)
VertexLagMatrix[, c(4, 7)] <- 1
VertexLagMatrix[c(2,3),7] <- 0
Vout1 <- paramVertexOnlyGroup(InputNetwork = beach,
                          maxLag = maxLag,
                          VertexStatsvec = VertexStatsvec,
                          VertexModelGroup = VertexModelGroup,
                          VertexLag = VertexLag,
                          VertexLagMatrix = VertexLagMatrix)
summary(Vout1$VertexFit$fit)

Author(s)

Abhirup