MLE.DCSBM: MLE.DCSBM

Description Usage Arguments Value Author(s) References Examples

View source: R/MLE_DCSBM.R

Description

Estimate the maximum likelihood estimators for P and delta at each time point in the list of adjacency matrices

Usage

1
MLE.DCSBM(Adjacency.list, community.array, T, k)

Arguments

Adjacency.list:

list of adjacency matrices in the observed dynamic network

community.array:

an array of length T whose tth entry is a numeric vector of length n specifying community labels at time t

T:

number of graphs in the temporal sequence

k:

number of communities (fixed accross time) Note, this is the pre-conceived idea of how many communities there are.

Value

a list containing the objects

Author(s)

James D. Wilson and Nathaniel T. Stevens

References

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
#Generate a collection of 50 networks with a change at time 25. The change is a local 
#change in connection propensity in community 1
n <- 100
P.old <- cbind(c(0.10, 0.01), c(0.02, 0.075))
P.new <- cbind(c(0.20, 0.025), c(0.02, 0.075))
P.array <- array(c(replicate(25, P.old), replicate(25, P.new)), dim = c(2, 2, 50))
community.array <- array(rep(c(rep(1, 50), rep(2, 50)), 50), dim = c(1, 100, 50))
delta.array <- array(rep(rep(0.2, 2), 50), dim = c(1, 2, 50))

dynamic.net <- dynamic.DCSBM(n = 100, T = 50, P.array = P.array,
                             community.array = community.array,
                             delta.array = delta.array, edge.list = FALSE)
image(Matrix(dynamic.net$Adjacency.list[[1]]))
image(Matrix(dynamic.net$Adjacency.list[[30]]))

#Estimate the MLEs
MLEs.example <- MLE.DCSBM(dynamic.net$Adjacency.list, community.array = community.array,
                           T = 50, k = 2)

jdwilson4/NetSurv documentation built on May 18, 2019, 11:40 p.m.