0

对于我的生活,我不知道这里出了什么问题,并认为我可能做错了什么:

就像我问过的所有其他问题一样,这与遗传学有关

我正在尝试使用 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 中的问题。

4

0 回答 0