write.csv(b, file="exposure.csv") (#此步骤是最不同于在线读取的地方。我们需将写出的文件,整理为下方代码读取的格式,然后放在TwoSampleMR包所在的位置处) bmi<-system.file("exposure.csv",package="TwoSampleMR") bmi_exp_dat<-read_exposure_data(filename = bmi,sep = ",",snp_col ="SNP",beta...
write.csv(b, file="exposure.csv") (#此步骤是最不同于在线读取的地方。我们需将写出的文件,整理为下方代码读取的格式,然后放在TwoSampleMR包所在的位置处) bmi<-system.file("exposure.csv",package="TwoSampleMR") bmi_exp_dat<-read_exposure_data(filename = bmi,sep =",",snp_col ="SNP",beta_...
代码语言:javascript 复制 hyperten_tophits<-tophits(id="ukb-b-12493",clump=0)hyperten_gwas<-rename(hyperten_tophits,c("SNP"="rsid","effect_allele.exposure"="ea","other_allele.exposure"="nea","beta.exposure"="beta","se.exposure"="se","eaf.exposure"="eaf","pval.exposure"="p","N...