получить все "NA" для p-значения при использовании tryCatch для t.test в R

Question

У меня есть данные о лечении и контроле, каждый из 3 повторов.Столбцы с 1 по 3 в данных относятся к лечению, а столбцы с 4 по 6 - соответственно.

Я выполнил t.test в R следующим образом:

df$p.value <- apply (df, 1, function(x) t.test(x[1:3], x[4:6],alternative = "two.sided", paired = TRUE)$p.value)

Когда я выполняю t.test с моими данными.Я получил ошибку:

Error in t.test.default(x[1:3], x[4:6], alternative = "two.sided", paired = TRUE) :data are essentially constant

Поэтому я хочу превратить ошибку в «NA», когда есть ошибка для определенных строк, используя tryCatch, как показано ниже:

df$p.value <- tryCatch({
  apply(df, 1, function(x) t.test(x[1:3], x[4:6],alternative = "two.sided", paired = TRUE)$p.value)
}, error = function(e) {NA})

Однако, он возвращает все "NA", несмотря на то, что у меня есть некоторые сопоставимые данные.Как я могу исправить этот скрипт?

Заранее спасибо.

Answer 1

Можно попробовать

apply (df, 1, 
   function(x) tryCatch({
     t.test(x[1:3], x[4:6],alternative = "two.sided", paired = TRUE)$p.value
   }, error = function(e) {NA}))