Learn R Programming

ppsbm (version 1.0.0)

generated_Q3_n20: Example dataset

Description

Example of undirected dataset with \(n=20\) individuals in \(Q=3\) clusters and observation Time=1

Usage

generated_Q3_n20

Arguments

Format

A list of 3 components:

data

Observed data is itself a list of 3 components:

  • time.seq - Vector containing the times of the events (length M)

  • type.seq - Vector containing the types of the events (length M)

  • Time - Positive real number. [0,Time] is the total time interval of observation

z

Latent variables. A matrix with size \(Q\times n\) and entries 1 (cluster q contains node i) or 0 (else).

intens

Intensities used to simulate data. A list of \(Q(Q+1)/2\) intensity functions. Each one is given as a list of 2 components:

  • intens - a positive function. The intensity function \(\alpha^{(q,l)}\)

  • max - positive real number. An upper bound on function \(\alpha^{(q,l)}\)

Details

This random datatset was obtained using the following code


intens <- list(NULL)
intens[[1]] <- list(intens=function(x) return (rep(4,length(x))), max=4.1)
intens[[2]] <- list(intens=function(x){
   y <- rep(0,length(x))
   y[x<.25] <- 4
   y[x>.75] <- 10
   return(y)
 }, max=10.1)
intens[[3]] <- list(intens=function(x) return(8*(1/2-abs(x-1/2))), max=4.1)
intens[[4]] <- list(intens=function(x) return(100*x*exp(-8*x)), max=4.698493)
intens[[5]] <- list(intens=function(x) return(exp(3*x)*(sin(6*pi*x-pi/2)+1)/2), max=12.59369)
intens[[6]] <- list(intens=function(x) return(8.1*(exp(-6*abs(x-1/2))-.049)), max=7.8031)

generated_Q3_n20 <- generateDynppsbm(intens,Time=1,n=20,prop.groups=rep(1/3,3),directed=F)

References

MATIAS, C., REBAFKA, T. & VILLERS, F. (2018). A semiparametric extension of the stochastic block model for longitudinal networks. Biometrika. 105(3): 665-680.