cenregによる打ち切り回帰のガウス分布

Question

RにあるNADAパッケージのcenreg関数を使用して、ガウス分布を持つ検閲済み回帰を実行しようとしています。関数は対数正規分布で動作しますが、ガウス分布を使用しようとするとエラーを返します。ここで

は、ダミーデータで再現可能な例である：evalの中

エラー（exprの、ENVIR、enclos）：私は、このエラーメッセージが表示できます

## 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')) 

オブジェクトが

見つからない 'X'

Rバージョン3.2.3とNADAバージョン1.5.6を使用しています。これは、NADAパッケージの最新バージョンです。

cenreg関数をガウス分布で使用するにはどうすればよいですか？

Answer 1

はバグのように見えるが、これは動作します：

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