Skip to contents

This function creates a null dataset by directly generating a Z-statistic matrix under the multivariate normal assumption used in the ARCHIE method.

Usage

gene_data_null3(n, p, data_h1, A_snp_cis, A_cis_trans, A_trans)

Arguments

n

Sample size

p

Number of SNPs

data_h1

A dataset generated by the gene_data function

A_snp_cis

Binary adjacency matrix indicating SNP-cis gene relationships

A_cis_trans

Binary adjacency matrix indicating cis-trans gene relationships

A_trans

Binary adjacency matrix indicating trans-trans gene relationships

Value

A matrix of Z-statistics under the null hypothesis

Examples

# Generate original data
data_h1 <- gene_data(n = 1000, p = 60, A_snp_cis, A_cis_trans, A_trans)
#> Error: object 'A_snp_cis' not found

# Generate null Z-matrix using ARCHIE setting
Z_matrix_null <- gene_data_null3(1000, 60, data_h1, A_snp_cis, A_cis_trans, A_trans)
#> Error: object 'data_h1' not found