R包TwoSampleMR执行harmonise_data()行数不一致报错原因咨询
报错触发原因
- 核心触发原因:p值列名不匹配
代码中read_outcome_data的pval_col参数填写的'Р-value'首字母为西里尔字母Р,与结局GWAS文件表头的英文大写P字符编码不一致,无法匹配到正确的p值列。部分SNP行因p值读取为缺失值,被函数默认过滤逻辑剔除,导致读取到的out_dat对应行数为0,和暴露数据的1行做匹配合并时就会触发行数不一致的报错。
直接将该参数修改为pval_col = 'P-value'(使用英文大写P输入)即可解决该问题。 - 参数拼写错误
调用harmonise_data时填写的参数名exporsure_dat为错误拼写,正确参数名为exposure_dat,拼写错误会导致函数无法正确识别传入的暴露数据,也可能引发数据匹配异常。 - 补充核对项
你提供的代码中filename参数行末尾缺少逗号,属于语法错误,如果实际运行时没有直接报语法错误,大概率是粘贴代码时的笔误,建议核对本地代码的参数分隔符是否正确。
内容的提问来源于stack exchange,提问作者Margobel
相关产品推荐
相关产品推荐

