This function creates a null dataset by directly generating a Z-statistic matrix under the multivariate normal assumption used in the ARCHIE method.
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
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