Ошибка в функции R.Попытка выбрать менее одного элемента в intergerOneIndex

Question

Я создал некоторые случайные данные испытаний следующим образом:

Trial1 <- AminoAcidTrialGenerator2(10, 6, aacodons$MOD_Frequency); Trial1

Это сделано с помощью функции, которую я написал для генерации последовательностей аминокислот в 3-буквенном коде на основе распределения вероятностей.

Это выводит:

[[1]]
[1] cys  glu  glu  stop phe  asn 
21 Levels: ala arg asn asp cys gln glu gly his ile leu lys ... val

[[2]]
[1] cys pro ile lys cys arg
21 Levels: ala arg asn asp cys gln glu gly his ile leu lys ... val

[[3]]
[1] asp gln gln ala gln his
21 Levels: ala arg asn asp cys gln glu gly his ile leu lys ... val

[[4]]
[1] thr pro lys thr lys phe
21 Levels: ala arg asn asp cys gln glu gly his ile leu lys ... val

[[5]]
[1] ala  gly  his  arg  stop trp 
21 Levels: ala arg asn asp cys gln glu gly his ile leu lys ... val

Что моя функция emperical.codon.freq должна сделать, это создать новый список TrialData, который содержит подсчеты частоты, с которой каждый трехбуквенный символ встречается в каждом испытании.

emperical.Codon.Freq <- function(trialvec) {
  TrialData <- list()
  ecf.df <- data.frame(Codon = aacodons$Amino_Acid)
  for (i in 1:length(unlist(trialvec[i]))) {
   indexv <- is.element(ecf.df$Codon, unlist(trialvec[i]))
   name <- paste("Trial", i, sep = "_")
   currentT <- unlist(trialvec[i])
    for (f in 1:length(indexv)) {
      c[i] <- sum(is.element(currentT, 
                             aacodons$Amino_Acid[indexv][f]), 
                             na.rm = TRUE)
      TrialData[[i]] <- unlist(c)
    }
  }
  TrialData
}

Когда я пытаюсь запустить функцию emperical.Codon.frequency, я получаю сообщение об ошибке:

ОШИБКА в TrialData [[i]] <- unlist (c): попытка выделить менее одного элемента в intergerOneIndex. </p>

dput(Trial1)

list(structure(c(5L, 7L, 7L, 17L, 14L, 3L), .Label = c("ala", 
"arg", "asn", "asp", "cys", "gln", "glu", "gly", "his", "ile", 
"leu", "lys", "met", "phe", "pro", "ser", "stop", "thr", "trp", 
"tyr", "val"), class = "factor"), structure(c(5L, 15L, 10L, 12L, 
5L, 2L), .Label = c("ala", "arg", "asn", "asp", "cys", "gln", 
"glu", "gly", "his", "ile", "leu", "lys", "met", "phe", "pro", 
"ser", "stop", "thr", "trp", "tyr", "val"), class = "factor"), 
    structure(c(4L, 6L, 6L, 1L, 6L, 9L), .Label = c("ala", "arg", 
    "asn", "asp", "cys", "gln", "glu", "gly", "his", "ile", "leu", 
    "lys", "met", "phe", "pro", "ser", "stop", "thr", "trp", 
    "tyr", "val"), class = "factor"), structure(c(18L, 15L, 12L, 
    18L, 12L, 14L), .Label = c("ala", "arg", "asn", "asp", "cys", 
    "gln", "glu", "gly", "his", "ile", "leu", "lys", "met", "phe", 
    "pro", "ser", "stop", "thr", "trp", "tyr", "val"), class = "factor"), 
    structure(c(1L, 8L, 9L, 2L, 17L, 19L), .Label = c("ala", 
    "arg", "asn", "asp", "cys", "gln", "glu", "gly", "his", "ile", 
    "leu", "lys", "met", "phe", "pro", "ser", "stop", "thr", 
    "trp", "tyr", "val"), class = "factor"), structure(c(1L, 
    19L, 5L, 10L, 3L, 7L), .Label = c("ala", "arg", "asn", "asp", 
    "cys", "gln", "glu", "gly", "his", "ile", "leu", "lys", "met", 
    "phe", "pro", "ser", "stop", "thr", "trp", "tyr", "val"), class = "factor"), 
    structure(c(14L, 14L, 19L, 3L, 12L, 14L), .Label = c("ala", 
    "arg", "asn", "asp", "cys", "gln", "glu", "gly", "his", "ile", 
    "leu", "lys", "met", "phe", "pro", "ser", "stop", "thr", 
    "trp", "tyr", "val"), class = "factor"), structure(c(13L, 
    5L, 9L, 21L, 16L, 13L), .Label = c("ala", "arg", "asn", "asp", 
    "cys", "gln", "glu", "gly", "his", "ile", "leu", "lys", "met", 
    "phe", "pro", "ser", "stop", "thr", "trp", "tyr", "val"), class = "factor"), 
    structure(c(2L, 21L, 8L, 16L, 3L, 17L), .Label = c("ala", 
    "arg", "asn", "asp", "cys", "gln", "glu", "gly", "his", "ile", 
    "leu", "lys", "met", "phe", "pro", "ser", "stop", "thr", 
    "trp", "tyr", "val"), class = "factor"), structure(c(10L, 
    11L, 20L, 11L, 5L, 17L), .Label = c("ala", "arg", "asn", 
    "asp", "cys", "gln", "glu", "gly", "his", "ile", "leu", "lys", 
    "met", "phe", "pro", "ser", "stop", "thr", "trp", "tyr", 
    "val"), class = "factor"))