Распределение Гаусса для цензурированной регрессии с cenreg

Я пытаюсь использовать функцию cenreg в пакете NADA в R для запуска цензурированной регрессии с распределением Гаусса. Функция работает с логнормальным распределением, но возвращает ошибку, когда я пытаюсь использовать распределение Гаусса.

Вот воспроизводимый пример с фиктивными данными:

## Simulate data
d <- data.frame(x = rep(1:10, times = 3))
d$conc <- d$x * 2 + rnorm(30, 0, sd = 0.25)
d$group <- rep(1:3, each = 10)
d$conc <- d$conc + d$group
d$group <- as.factor(d$group)
d$conc[ d$conc < 8] <- 0
d$conc <- d$conc + .01
d$cen <- FALSE
d$cen[ d$conc < 8] <- TRUE

library(NADA)
## Fit model with log-normal distribution 
m1 <- with(d, cenreg(Cen(conc, cen) ~ x + group - 1,
                     dist = 'lognormal'))
m1
## attempt to fit model with Gaussian distribution
m2 <- with(d, cenreg(Cen(conc, cen) ~ x + group - 1,
                     dist = 'gaussian'))

который дает мне это сообщение об ошибке:

Ошибка в eval(expr, envir, enclos): объект 'x' не найден

Я использую R версии 3.2.3 и NADA версии 1.5.6. Это самая последняя версия пакета NADA.

Как использовать функцию cenreg с распределением Гаусса?


person Richard Erickson    schedule 22.03.2016    source источник


Ответы (1)


Похоже на ошибку, но это работает:

cenreg(Cen(d$conc, d$cen) ~ d$x + d$group - 1, dist = 'gaussian')
person Julius Vainora    schedule 22.03.2016
comment
Да, этот код работает :-). Я тоже пробовал это, но не должен был. Кстати, я изначально использовал data.table и получал похожие сообщения об ошибках. - person Richard Erickson; 23.03.2016