paramVertexOnly function

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

paramVertexOnly(
  InputNetwork,
  VertexStatsvec = rep(1, nvertexstats),
  maxLag,
  VertexLag = rep(1, maxLag),
  VertexLagMatrix = matrix(1, maxLag, length(VertexStatsvec)),
  dayClass = NA,
  regMethod = "bayesglm"
)

Arguments

• InputNetwork: Input network list.
• VertexStatsvec: Binary vector of size 9, indicating vertex model.
• maxLag: maximum lag.
• VertexLag: Binary vector of size maxLag, indicating Lag terms in the model.
• VertexLagMatrix: Binary matrix indicating lagged vertex statistics in the model.
• dayClass: Any network level present time attribute vector. Here used to indicate week/weekend as 0/1.
• 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
maxLag = 3
VertexLag = rep(1, maxLag)
VertexLagMatrix <- matrix(0, maxLag, nvertexstats)
VertexLagMatrix[, c(4, 7)] <- 1
VertexLagMatrix[c(2,3),7] <- 0
getWeekend <- function(z){
    weekends <- c("Saturday", "Sunday")
    if(!network::is.network(z)){
        if(is.na(z)) return(NA)
    } else {
         zDay <- get.network.attribute(z, attrname = "day")
         out <- ifelse(zDay %in% weekends, 1, 0)
         return(out)   
    }
}

## for(i in 1:31) print(getWeekend(beach[[i]]))
## generate a vector of network level exogenous variable
dayClass <- numeric(length(beach))
for(i in seq_along(dayClass)) {
    dayClass[i] <- getWeekend(beach[[i]])
}
out <- paramVertexOnly(InputNetwork = beach,
                       maxLag = 3,
                       VertexStatsvec = rep(1, nvertexstats),
                       VertexLag = rep(1, maxLag),
                       VertexLagMatrix = VertexLagMatrix,
                       dayClass = dayClass)

Author(s)

Abhirup