cenreg ile sansürlü regresyon için Gaussian dağılımı

Question

Bir Gauss dağılımı ile sansürlü bir regresyon gerçekleştirmek için NADA paketindeki cenreg işlevini kullanmaya çalışıyorum. Bu işlev, log normal dağıtma ile çalışır, ancak Gauss dağılımını kullanmaya çalıştığımda bir hata döndürür. Bana bu hata mesajını verir

## 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 içinde

Hatası (ifade, envir, enclos): Nesne

bulunamadı 'x' İşte

taklit veriler tekrarlanabilir örnektir

R sürüm 3.2.3 ve NADA sürüm 1.5.6 kullanıyorum. Bu, NADA paketinin en son sürümüdür.

Gauss dağılımı ile cenreg işlevini nasıl kullanırım?

Answer 1

bir böcek gibi görünüyor, ama bu işleri:

cenreg(Cen(d$conc, d$cen) ~ d$x + d$group - 1, dist = 'gaussian')