Медиация для цензурированного результата с использованием функции vglm в R - PullRequest
Новая
       83

Медиация для цензурированного результата с использованием функции vglm в R

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

У меня есть переменная лечения, медиатора и результата. Переменная результата остается подвергнутой цензуре.

Как указано в этой ссылке (http://lists.r -forge.r-project.org / pipermail / mediation-information / 2011-July / 000002.html ), используя vglm (пакет vgam) для результата, функция-посредник должна работать.

Когда я использую lm для результата, функция работает. Однако, с vglm это не работает, и показывает эту ошибку:

Ошибка в eval (predvars, data, env): объект 'NBGPercent' не найден

Если кто-нибудь может мне помочь, я очень ценю это. Вот код: 

library(mediation)
library(VGAM)

model.m<-lm(NS34 ~ CRTscore, data=stock_binary)
model.y<-vglm(NBGPercent ~ CRTscore + NS34 + CRTscore:NS34 ,tobit(Lower =0),data=stock_binary)
out.1 <- mediate (model.m, model.y,  sims = 1000, treat = "CRTscore", mediator = "NS34")
out.2 <- mediate (model.m, model.y, sims = 1000, treat = "CRTscore", mediator = "NS34", data=stock_binary)
summary(out.2)
...