如何解决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 举报,一经查实,本站将立刻删除。