Вставка NA после теста

Question

У меня есть список, полный данных. Кадры с двумя столбцами, временем и сигналом. Данные представляют собой результаты хроматографического анализа на основе процесса, который периодически отбирали.

Я хочу сравнить данные gc, которые я собрал.

Я написал функцию преобразования времен и площадей пиков в процентные области (исключая пик растворителя) и относительное время удерживания.

В связи с характером процесса разные GC имеют различное количество пиков, и поэтому сравнение не является простым. Примеси появляются в разных частях моего процесса и, следовательно, дают дополнительные пики.

Я хочу перечислить свой список и найти самый длинный вектор относительного времени удерживания (без проблем). Я хочу использовать самый длинный вектор в качестве компаратора и устанавливать значения NA в относительное время хранения, которые появляются одновременно с компаратором, но не отображаются в других кадрах данных.

Отсюда результаты следующего списка относительного времени удерживания,

prac <- list(a=c(0.203,0.305,0.444,0.780,1.000,1.101,1.403), 
      b=c(0.201,0.306,0.442,0.778,1.000,1.101,1.208,1.401)) 

где b- вектор компаратор, должен выглядеть

0.203 0.305 0.444 0.780 1.000 1.101 NA 1.403 
0.201 0.306 0.442 0.778 1.000 1.101 1.208 1.401 

Может кто-нибудь подскажет, как я мог бы быть в состоянии начать ?

Моя первая мысль была для цикла, но я не думаю, что это сработает. Обратите внимание, что иногда требуется более 1 значения NA.

(Я планирую сопоставить процентные области с относительным временем удерживания компаратора для всех хроматограмм, если только я могу выйти за рамки этой проблемы).

Answer 1

Это одно решение (перебор) только с одним отсутствующего значением:

prac <- list(a=c(0.203,0.305,0.444,0.780,1.000,1.101,1.403), 
       b=c(0.201,0.306,0.442,0.778,1.000,1.101,1.208,1.401)) 


NA.index <- which(abs(prac$b[1:length(prac$a)] - prac$a) > 0.05) 
newlist.a <- c(prac$a[1:NA.index-1], NA, prac$a[NA.index]) 

это здесь должно быть generizeable (в зависимости от того, как данные на самом деле структурированы):

prac <- list(a=c(0.203,0.305,0.444,0.780,1.000,1.101,1.403), 
       b=c(0.201,0.306,0.442,0.778,1.000,1.101,1.208,1.401)) 

for(i in seq_along(prac$a)) { 
    if(abs(prac$b[i] - prac$a[i]) < 0.05) { 
     prac$a[i] <- prac$a[i] 
    } else { 
     prac$a[i+1] <- prac$a[i] 
     prac$a[i] <- NA 
    } 
} 

любопытных трудно чтобы рассказать, как обобщить это на примерах с множественными NA, если вы не дадите другой воспроизводимый пример, потому что прямо сейчас я просто нахожусь в темноте о том, как структурированы ваши данные.