Power.law.network: Power law network

In HeteroGGM: Gaussian Graphical Model-Based Heterogeneity Analysis

Power law network

Description

Generating three s-block power-law precision matrices

Usage

Power.law.network(p, s = 10, umin = 0.1, umax = 0.4, I2 = 0, I3 = 0)

Arguments

p	The dimensions of the precision matrix.
s	The number of sub-networks.
umin	The lower bound of non-zero elements on non-diagonal elements.
umax	The upper bound of non-zero elements on non-diagonal elements.
I2	The replacement blocks for the precision matrix of the second subgroup.
I3	The replacement blocks for the precision matrix of the third subgroup.

Value

A list including The precision matrices of three subgroups.

Examples

p <- 20               # The dimension of the precision matrix
################ The true parameters ################
# Power law network
set.seed(2)
A.list <- Power.law.network(p,s=5,I2=c(1),I3=c(2))
Theta01 <- A.list$A1
Theta02 <- A.list$A2
Theta03 <- A.list$A3
sigma01 <- solve(Theta01)
sigma02 <- solve(Theta02)
sigma03 <- solve(Theta03)
Sigma0.list <- list(sigma01,sigma02,sigma03)
Theta0.list <- list(Theta01,Theta02,Theta03)

HeteroGGM documentation built on Oct. 11, 2023, 5:14 p.m.