我试图为一个被审查的data.frame估计一个Tobit和一个截断的回归模型,但是遇到一些问题。首先,如果它使用“ tobit”函数,我将无法获得任何结果。代码如下,其中YY是我的依存者,XX是我的独立矩阵。
tobit1<-tobit(YY~XX,left=0,right=Inf,data=consumo_familiare_2016,iterlim=200)
我收到的错误消息是
轨道误差(YY〜XX,左= 0,右= Inf,数据= consumo_familiare_2016,: 未使用的参数(left = 0,right = Inf,data = consumo_familiare_2016,iterlim = 200)
但是,如果我使用censReg包,则可以得到一些结果。我使用的代码是
tobit2<-censReg(YY~XX, left=0,right=Inf,data=consumo_familiare_2016,iterlim=200)
我得到的结果是回归输出
Newton-Raphson最大化,迭代28次
返回码2:公差范围内的连续功能值 对数似然:14 Df上的-40229.48
返回码2是否有问题?只有一种算法可行? 此外,如果我尝试使用相同的数据估计两步截断模型,则会出现奇异性问题。我正在使用的代码是
> probit1<-glm(I(YY>0)~XX,data=consumo_familiare_2016,family=binomial(link="probit"))
> summary(probit1)
> data.truncated<-subset(consumo_familiare_2016,YY>0)
> truncreg1<-truncreg(YY~XX,point=0,direction="left",data=data.truncated)
> summary(truncreg1)
尝试总结truncregression的输出后得到的错误是
警告信息: 在sqrt(diag(vcov(object)))中:产生的NaNs
与所得的sigma一起显示st的NA。错误,t统计量和p值。
我相信这两个问题是相关的,并且我的数据或代码有一些问题。这是否与我检查过许多观测结果(10615中的7082)有关?有人可以帮我吗?
非常感谢