Generate null data by removing cis-to-trans effects

This function creates a null dataset by regenerating trans gene expression data without the influence of cis genes, while maintaining trans-to-trans effects.

Usage

gene_data_null2(data_h1, A_trans, gamma_trans_trans = 0.5)

Arguments

data_h1

A dataset generated by the gene_data function

A_trans

Binary adjacency matrix indicating trans-trans gene relationships

gamma_trans_trans

Effect size for trans to trans gene effects

Value

A modified dataset with new trans gene expression lacking cis-to-trans effects

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 data without cis-to-trans effects
data_null <- gene_data_null2(data_h1, A_trans, gamma_trans_trans = 0.5)
#> Error: object 'data_h1' not found