Description Usage Arguments Value Author(s) Examples
View source: R/paramVertexOnly.R
Parameter estimation for Vertex model only for a list of dynamic networks.
1 2 3 4 5 6 7 8 9 |
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" |
List of 3 elements:
VertexFit: Output from regEngine.
VertexStats: Subsetted vertex stats matrix.
VertexStatsFull: Full matrix of vertex stats.
Abhirup
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.