Являются ли веса в пакете truncreg сломанным?

Question

Я что-то упустил?

library(truncreg) 

n <- 10^4 
lambda <- 0.3 # Proba y is taken from component 0 

df <- data.frame(x=rnorm(n)) 
df$y0 <- pmax(rnorm(n, 10 + df$x, 5), 0) 
df$y1 <- pmax(rnorm(n, 2 - 5*df$x, 2), 0) 
df$component <- ifelse(runif(n) < lambda, 0, 1) 
df$y <- ifelse(df$component == 0, df$y0, df$y1) # Mixture of censored regressions 

plot(df$x, df$y) 

model <- truncreg(y ~ x, data=df) # All data 
model.w <- truncreg(y ~ x, data=df, weights=component) # Only component 1, using weights 
model.subset <- truncreg(y ~ x, data=subset(df, component == 1)) # Only component 1, using subset 

identical(coefficients(model), coefficients(model.w)) # True -- I expected this to be false 
identical(coefficients(model.w), coefficients(model.subset)) # False -- I expected this to be true 

## For comparison, here is the same using lm: 
model <- lm(y ~ x, data=df) 
model.w <- lm(y ~ x, data=df, weights=component) 
model.subset <- lm(y ~ x, data=subset(df, component == 1)) 

identical(coefficients(model), coefficients(model.w)) # False as expected 
identical(coefficients(model.w), coefficients(model.subset)) # True as expected 

Answer 1

Да, я могу воспроизвести вашу проблему.
Затем я попытался установить method="model.frame" в lm пробегах и получить те же «неожиданные» результаты, которые вы получили, т. Е. Те же коэффициенты с применением или без весов. Я заглянул в источник truncreg и не увидел никакого очевидного места, которое «выбирает» не для использования method="model.frame"; затем я ворвался в источник truncreg.fit и снова не видел ссылки на значения весов. Мне непонятно, что делается , поэтому весы могут быть переданы в код соответствия, но я могу начать с более тщательного копания через этот код.