我有一个包含SNP(基因型)的大数据集,分别是501行和42844列。
数据是这样的:
OAR19_64675012.1 OAR19_64803054.1 OAR1_88143.1 s09912.1 s36301.1 s34880.1
AB AB BB BB AA AA
AB AB AB AA AB AB
AB AB BB AB BB AA
AA BB BB AB AA AA
我在质量控制过程中遇到了很大的问题。 我在R中使用QCGWAS库,当我点击命令以启动质量控制时:
headers<- names(snp) ## names of the columns, snp is the dataframe above
# make them as a dataframe with characters
h_data<- as.character(data.frame(headers , stringsAsFactors = F))
new_data <- QC_GWAS(filename = "data1.txt", # command for QC data1.txt is my snps data
header_translations = h_data,
save_final_dataset = TRUE)
我收到此错误:
Error: file.exists(paste(dir_references, header_translations, sep = "/")) is not TRUE
我不知道自己要做什么......我一直在寻找东西,但我一直都会遇到错误......我尝试按照/中的快速入门指南进行操作R / library / QCGWAS / doc,但我得到了错误。
如果有人有想法和时间来帮助我,我会想到这一点。为了继续这个项目,解决这个问题对我来说非常重要。我使用过其他一些图书馆,但我无法处理质量控制,因为他们没有给我干净的数据。