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

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

0 голосов
/ 11 мая 2018

Я новичок в 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, и чтение в Интернете не помогло с этой конкретной проблемой. Буду признателен за любой вклад.

1 Ответ

0 голосов
/ 11 мая 2018

Решение: замените оператор NonHomoSpikes return на 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))
   }
...