Я пытаюсь объединить некоторые значения p для тестов статистической значимости из ряда источников, но некоторые отсутствуют.
Похоже, что команда p.tpm из пакета TFisher R может помочь, но я сталкиваюсь с некоторыми проблемами при его реализации.
В принципе, если я помещаю все значения p, включая NA, в командную строку, я получаю результат NA, который, конечно, бесполезен:
install.packages("TFisher")
library(TFisher)
pval <- c(0.53, 0.01, 0.16, NA, 0.02, NA, 0.18, 0.8, NA, 0.25)
tpmstat <- stat.tpm(p=pval, tau1=0.05)
p.tpm(q=tpmstat, n=10, tau1=0.05)
Я даже пытался указать фактические значения p, исключая NA.Но он все равно не работает должным образом, так как дает итоговое значение p, намного превышающее значение каждого из них:
install.packages("TFisher")
library(TFisher)
pval <- c(0.53, 0.01, 0.16, 0.02, 0.18, 0.8, 0.25)
tpmstat <- stat.tpm(p=pval, tau1=0.05)
p.tpm(q=tpmstat, n=10, tau1=0.05)
Как эту команду следует применять в этом случае?