Я пытаюсь выполнить мета-анализ с помощью метода вероятности профиля в пакете metaLik в R. В качестве примера я использую набор данных, приведенный ниже из этой статьи о пакете metaLik: "https://www.jstatsoft.org/article/view/v050i07 "
Набор данных диуретиков
Study Y Sigma2
1 0.04 0.16
2 -0.92 0.12
3 -1.12 0.18
4 -1.47 0.30
5 -1.39 0.11
6 -0.30 0.01
7 -0.26 0.12
8 1.09 0.69
9 0.14 0.07
Я могу выполнить мета-анализ с помощью следующего кода:
m <- metaLik(Y~1, data=diuretics, sigma2=Sigma2)
summary(m)
Но я не могу вычислить доверительные интервалы. В статье выше говорится «Доверительный интервал для параметра модели β на основе статистики Вальда c получается стандартным методом confint (), в то время как доверительные интервалы вероятности профиля предоставляются функцией profile ()»
Однако, когда я пытаюсь это сделать, я получаю следующие результаты:
confint(m)
2.5 % 97.5 %
(Intercept) NA NA
И
profile(m)
Error in UseMethod("profile") :
no applicable method for 'profile' applied to an object of class "metaLik"
Я использую точный набор данных и методы, как в примере. Кто-нибудь знает, что я делаю не так?