paramVertexOnly: Parameter estimation for Vertex model only for a list of...

In dnr: Simulate Dynamic Networks using Exponential Random Graph Models (ERGM) Family

Description

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

Usage

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"

Value

List of 3 elements:
VertexFit: Output from regEngine.
VertexStats: Subsetted vertex stats matrix.
VertexStatsFull: Full matrix of vertex stats.

Author(s)

Abhirup

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)

dnr documentation built on Dec. 1, 2020, 1:07 a.m.