1d:T7a0,# Parameter estimation for Vertex model only for a list of dynamic networks. ```r 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 ```r 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

paramVertexOnly function