我仍然发现 R 中的 ifelse 结构有点令人困惑,我有以下数据框:
df <- structure(list(snp = structure(1:11, .Label = c("AL0009", "AL00014", "AL0021", "AL00046", "AL0047", "AS0005", "AS0014", "AS00021", "AS0047", "AS0071", "DR0001" ), class = "factor"), CHROMOSOME = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), COUNT_ALLELE = structure(c(1L, 1L, 1L, 3L, 1L, 1L, 1L, 2L, 3L, 3L, 1L), .Label = c("A", "C", "G"), class = "factor"), OTHER_ALLELE = structure(c(3L, 3L, 2L, 1L, 3L, 2L, 2L, 1L, 1L, 1L, 3L), .Label = c("A", "C", "G"), class = "factor"), `116601888` = c(0L, 0L, 0L, 2L, 2L, 0L, 0L, 0L, 0L, 0L, 2L ), `116621563` = c(0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 1L), `117253533` = c(0L, 0L, 0L, 2L, 2L, 0L, 0L, 0L, 1L, 0L, 2L), `117423827` = c(1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 2L)), .Names = c("snp", "CHROMOSOME", "COUNT_ALLELE", "OTHER_ALLELE", "11688", "11663", "11533", "13827" ), row.names = c(NA, 11L), class = "data.frame")
# snp CHROMOSOME COUNT_ALLELE OTHER_ALLELE 11688 11663 11533 13827
# 1 AL0009 1 A G 0 0 0 1
# 2 AL00014 1 A G 0 0 0 1
# 3 AL0021 1 A C 0 0 0 1
# 4 AL00046 1 G A 2 1 2 1
# 5 AL0047 1 A G 2 1 2 1
# 6 AS0005 1 A C 0 0 0 0
# 7 AS0014 1 A C 0 0 0 0
# 8 AS00021 1 C A 0 1 0 0
# 9 AS0047 1 G A 0 0 1 1
# 10 AS0071 1 G A 0 0 0 1
# 11 DR0001 1 A G 2 1 2 2
使用TranslateAllele
函数 我想用相应的两个字母代码替换从第 5 列开始的列中的数字:
TranslateAllele <- function(COUNT_ALLELE, OTHER_ALLELE, genotype){
if(genotype==0){
print(paste(OTHER_ALLELE, OTHER_ALLELE, sep=""))
} else if(genotype==1){
print(paste(COUNT_ALLELE, OTHER_ALLELE, sep=""))
} else if(genotype==2){
print(paste(COUNT_ALLELE, COUNT_ALLELE, sep=""))
}
}
所以期望的输出如下:
# snp CHROMOSOME COUNT_ALLELE OTHER_ALLELE 11688 11663 11533 13827
# 1 AL0009 1 A G GG GG GG AG
# 2 AL00014 1 A G GG GG GG AG
# 3 AL0021 1 A C CC CC CC AC
# 4 AL00046 1 G A GG GA GG GA
# 5 AL0047 1 A G AA AG AA AG
# 6 AS0005 1 A C CC CC CC CC
# 7 AS0014 1 A C CC CC CC CC
# 8 AS00021 1 C A AA CA AA AA
# 9 AS0047 1 G A AA AA GA GA
# 10 AS0071 1 G A AA AA AA GA
# 11 DR0001 1 A G AA AG AA AA
最终我需要对 160 万行乘 1M 列执行此操作,因此我将无法简单地使用 for 循环:(