lm运行，但是censReg没有相同的数据和变量

问题描述

我遇到了与here所解释的问题类似的问题。

示例代码也来自这个问题。

Happiness <- c(60,80,39,40,69,90,100,30)

City <- as.factor(c("New York","Chicago","New York","New York"))
Gender <- as.factor(c(0,1,1)) # 0 = man,1 = woman.
Employment <- c(0,1) # 0 = unemployed,1 = employed.
Worktype <- as.factor(c("Unemployed","Whitecolor","Unemployed","Bluecolor","Whitecolor"))
Holiday <- as.factor(c(0,2,1))
levels(Holiday) <- c("Unemployed","1 day a week","2 day a week")

data <- data.frame(Happiness,City,Gender,Employment,Worktype,Holiday)

head(data,8)
str(data)

reg <- lm(Happiness ~ City + Gender + Employment:(Worktype + Holiday))
summary(reg)

正如问题的答案所解释的那样，使用lm时这些NA不会造成任何问题。

但是，我没有使用lm，而是使用censReg，所以我这样做了：

reg <- censReg(Happiness ~ City + Gender + Employment:(Worktype + Holiday),left=0,right=100)
summary(reg)

现在有一个问题，因为censReg无法正确处理此问题并给出错误：

Error in printCoefmat(coef(x,logSigma = logSigma),digits = digits) : 
  'x' must be coefficient matrix/data frame

我第一次遇到这个问题时，我对虚拟对象进行了重新编码，但是当（在我的实际数据中）虚拟对象具有7级并且交互是连续的时，这变得有些烦人。

有人可以帮助我解释为什么censReg遇到问题以及我应该怎么做吗？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

lm r r