Из модели выпали избыточные предикторы? - PullRequest
Новая
       31

Из модели выпали избыточные предикторы?

0 голосов
/ 01 октября 2019

Я выполняю метаанализ и выполняю двойное дугообразное преобразование в моей модели. В настоящее время я делаю анализ подгруппы по переменной модератора (лабораторный метод). Переменная модератора имеет 2 группы (лаборатория 1 и лаборатория 2). В lab1 есть 6 объектов и 9 в lab2. Мой фрейм данных (файл CSV) содержит 15 строк данных.

Я использую следующий код: 

subganal.lab_method=rma(yi, vi, data=ies.da, mods=~lab_method)
pes.da.lab1=rma(yi, vi, data=ies.da, mods=~lab_method=="lab1")
pes.da.lab2=rma(yi, vi, data=ies.da, mods=~lab_method=="lab2")
pes.subg.lab_method=predict(subganal.lab_method,
                           transf=transf.ipft.hm,
                           targ=list(ni=dat$samples_total))
dat.samevar=data.frame(estimate=c((pes.da.lab1$b)[1], (pes.da.lab2$b)[1]),
                       stderror=c((pes.da.lab1$se)[1], (pes.da.lab2$se)[1]),
                       tau2=subganal.lab_method$tau2)
pes.da.lab_method=rma(estimate, sei=stderror,
                      method="DL",
                      data=dat.samevar)
pes.lab_method=predict(pes.da.lab_method, transf=transf.ipft.hm, targ=list(ni=dat$samples_total))
print(pes.subg.lab_method ["lab1"]) #display subgroup 1 summary effect size
print(pes.subg.lab_method["lab2"]) #display subgroup 2 summary effect size
print(subganal.lab_method) #display subgroup analysis results
print(pes.lab_method) #display recomputed summary effect size

Я получаю сообщение об ошибке "Избыточные предикторы выпали из модели" дважды. Один раз для lab1 и один раз для lab2. Кроме того, сводные эффекты подгрупп для lab1 и 2 не печатаются. Он печатает «NULL», где должно быть напечатано значение.

Общий анализ подгруппы печатается.

"dat" - это имя моего фрейма данных

...