R语言Tobit模型的分组回归

R语言Tobit模型的分组回归

用R语言进行了2步的subset取子集

原始数据：cfps2014（样本量：37147） —— 去除缺省值后新数据为cfps14_all（样本量为：10353）

cfps14_all<-subset(cfps2014, age!="NA" & marrige!="NA"）

想进一步分组进行Tobit回归（分别按照 城市男、城市女、农村男、农村女）进行回归

现在2个思路：

1、subset取子集后分别进行回归

cfps_man_urban<-subset(cfps14_all, gender == "男" & urban == "城市")

结果发现新生成的数据集cfps_man_urban的样本数为**：7951！**（实际应该为2000多）

如果加$限定数据集，则样本量变为0

cfps_man_urban<-subset(cfps14_all, cfps14_all$gender == "男" & cfps14_all$urban == "城市")

2、直接在tobit模型中进行分组回归

m_u <- censReg( hw ~ age + edu + gapincome, data = cfps14_all[gender=="男" & curban=="城市"])
summary(m_u)

跑出来的结果样本是按照cfps14_all的，没有实现分组回归的目的。

想问一下如何实现Tobit模型的分组回归？