对于我的生活,我不知道这里出了什么问题,并认为我可能做错了什么:
就像我问过的所有其他问题一样,这与遗传学有关
我正在尝试使用 A/A 等字符串拆分遗传变量,并使用以下 R 代码:
gt$A1 <- sapply(strsplit(as.character(gt[c(4)]), "/"), function(x) x[1])
gt$A2 <- sapply(strsplit(as.character(gt[c(4)]), "/"), function(x) x[2])
然而,出来的是
A1 = ' c("G '
A2 = ' G", "G '
即使基因型不是 G/G,也适用于每一个变体。
我的文件示例如下所示:
ID MGT FGT CGT
001 A/A A/G G/A
002 T/C T/C C/C
003 T/C C/C T/C
有没有理由为什么这不能干净地分裂?我假设字符串的长度可能会搞砸 - 但不确定这是否是 R 中的问题。