基于MR分析和共定位分析的输出文件,批量进行ld check分析,作为共定位分析的补充方法。

Usage,
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_harmonised_file

mr分析输出的工具变量文件,对应的文件名"table_s1_harmonise_data.csv"。

mr_results_file

mr分析输出的结果文件(可经过p值筛选,或不经过p值筛选),对应的文件名"table_s2_mr_results.csv"或者"table_s2_mr_results_fdr_sig.csv"。 建议使用经过p值筛选的阳性结果,可减少运算量。

outcome_files

用于ld check分析的结局文件,为标准mr分析输入格式,可以是1个或者多个结局文件。

outcome_p

进行ld_check的gwas数据中的显著性SNP阈值,默认为1e-4。

outcome_number

进行ld_check的gwas数据中的显著性SNP的数量,默认为30,即p值最小的前30个SNP。如果是“all”,则选择全部显著性snp。

exposure_annotation_file

mr分析中,暴露性状的基因注释文件,用于ld check分析选择感兴趣染色体区域(region)。

up_num

指定编码基因的上游碱基范围,coloc分析常用的是1000kb,也有研究中使用100kb。参考文献,默认为100kb。

down_num

指定编码基因的下游碱基范围,coloc分析常用的是1000kb,也有研究中使用100kb。参考文献,默认为100kb。

ld_check_r2

通过ld_check的最小r2,默认为0.7。更严格一些,也可设置为0.8。

bfile

参考bfile文件,如“F:/data_ref/1kg.v3/EUR”

plink_exe

plink,默认内置最新plink工具。

out_path

指定输出文件的地址,默认是当前工作路径"./"。

out_prefix

指定输出文件的前缀,默认为"ld_r2"。

返回ld check的计算结果。