Как добавить первый и последний элемент в каждый вектор (разной длины) списка в R, не замедляя код?

Question

Я новичок в R и StackOverFlow. Так что, если я что-то упустил, пожалуйста, дайте мне знать. Я работаю над моделированием неоднородных процессов, чтобы лучше понять поведение нейронов. Мой код работает так, как я, скажем 20 длительные испытания 2 секунд каждая (каждое испытание представляет шипованный поезд). Тогда список SpikeTimes дает мне список 20 векторы, где каждый вектор соответствует отметкам времени, когда всплески произошли в этом конкретном испытании. [Например.SpikeTimes[1] который выглядит так, 0.002250802 0.053934034...1.971574170 2.023090384 означает, что в первом Спайк-поке 0.002250802, 0.053934034 и так далее. Я не знаю, почему это также вызывает отметки времени за пределами моего лимита времени 2 секунд, но я буду работать над этим позже. Мой код выглядит так:

nTrials=20
t_max=2.0000000


LambdaInv<- function(x){ifelse( x< 15, x/30,
                    ifelse(x >= 15 & x < 38, ((x-15)/46)+0.5,
                           ifelse(x >= 38 & x <53, ((x-38)/30)+1.0,
                                  ifelse(x >= 53 & x< 67.4, ((x-53)/72)+1.5, 
                                         ((x-67.4)/30)+1.7))))}

 t = 0 
 s = 0
 X = numeric(0)


 NonHomoSpikes <- function(t_max){
  while(t <= t_max){
    u <- runif(1)
    s <- s-log(u)
    t <- LambdaInv(s)
    X <- c(X,t)
       }
    return(X)
   }


SpikeTimes <- lapply(1:nTrials, function(x) NonHomoSpikes(2))

Моя проблема в том, что для каждого вектора в списке SpikeTimes; который дает отметки времени шипов, я также хочу включить начало (то есть 0и конец (то есть 2) шипового поезда. Поэтому я хочу добавить этот список, чтобы каждый вектор включал первую запись как 0 и последняя запись как 2,

мой SpikeTimes[1] будет выглядеть 0 0.002250802 0.053934034...1.971574170 2 и другие SpikeTimes[i] будет выглядеть похоже. Я старался SpikeTimes <- c(0, SpikeTimes) для входа 0 в начале, но это только сделало список имеет 21 вектор вместо 20 с вектором 0 в качестве первого элемента (я имею в виду, я понимаю, почему это произошло). Как я могу сделать это так, чтобы мой код не замедлялся? Я новичок в R, и чтение в Интернете не помогло с этой конкретной проблемой. Буду признателен за любой вклад.

Answer 1

Решение: изменить NonHomoSpikesreturn заявление к return(c(0, X, 2)),

"Вручную" остановить if t < 2 (как описано в комментариях):

NonHomoSpikes <- function(t_max){
  while(t <= t_max){
    u <- runif(1)
    s <- s-log(u)
    t <- LambdaInv(s)
    if(t > 2) break
    X <- c(X,t)
       }
    return(c(0,X,0))
   }