ld_check_batch.Rd基于MR分析和共定位分析的输出文件,批量进行ld check分析,作为共定位分析的补充方法。
ld_check_batch(
mr_harmonised_file,
mr_results_file,
outcome_files,
outcome_p = 1e-04,
outcome_number = 30,
exposure_annotation_file,
up_num = 1e+05,
down_num = 1e+05,
ld_check_r2 = 0.7,
bfile,
plink_exe = get_plink_exe(),
out_path = "./",
out_prefix = "lc"
)mr分析输出的工具变量文件,对应的文件名"table_s1_harmonise_data.csv"。
mr分析输出的结果文件(可经过p值筛选,或不经过p值筛选),对应的文件名"table_s2_mr_results.csv"或者"table_s2_mr_results_fdr_sig.csv"。 建议使用经过p值筛选的阳性结果,可减少运算量。
用于ld check分析的结局文件,为标准mr分析输入格式,可以是1个或者多个结局文件。
进行ld_check的gwas数据中的显著性SNP阈值,默认为1e-4。
进行ld_check的gwas数据中的显著性SNP的数量,默认为30,即p值最小的前30个SNP。如果是“all”,则选择全部显著性snp。
mr分析中,暴露性状的基因注释文件,用于ld check分析选择感兴趣染色体区域(region)。
指定编码基因的上游碱基范围,coloc分析常用的是1000kb,也有研究中使用100kb。参考文献,默认为100kb。
指定编码基因的下游碱基范围,coloc分析常用的是1000kb,也有研究中使用100kb。参考文献,默认为100kb。
通过ld_check的最小r2,默认为0.7。更严格一些,也可设置为0.8。
参考bfile文件,如“F:/data_ref/1kg.v3/EUR”
plink,默认内置最新plink工具。
指定输出文件的地址,默认是当前工作路径"./"。
指定输出文件的前缀,默认为"ld_r2"。
返回ld check的计算结果。