我正在尝试根据6(有时是7)数字代码对样本进行排序。当我有7位数字的代码时,我遇到了一个问题。
我的命名系统示例:
jT1-1A
jT1-1B
iT2-1A
iT2-1B
iT2-2A
iT2-2B
rT2-10A
rT2-10B
前四个数字代表“位置”,下一个数字(有时是2个数字)代表“个人”,最后,A或B代表样本(大多数人同时拥有A和B)。最终,我尝试对具有两个样本的样本进行排序,然后取平均值。
library(dplR)
bb_pixx <- read.tucson('ButteLake_Burn_PIXX_FINAL_DATED.TXT', header = NULL, long = FALSE, encoding = getOption("encoding"), edge.zeros = TRUE)
##Calulate mean across cores in a tree
bb_pixx.ids <- read.ids(bb_pixx, stc = c(4, 1, 1), ignore.site.case = TRUE)
我该如何编写代码以解决偶尔出现的2位数字?
希望这很清楚。