Доверительные интервалы в пакете metaLik в R - PullRequest
0 голосов
/ 28 апреля 2020

Я пытаюсь выполнить мета-анализ с помощью метода вероятности профиля в пакете 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"

Я использую точный набор данных и методы, как в примере. Кто-нибудь знает, что я делаю не так?

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...