在R中执行区间回归

我试图运行区间回归，其中我的因变量y由14个区间组成，代表收入。我有5000次观察。我试图用六个自变量来预测我的y

我尝试按照此处执行的步骤进行操作：

所以，我实际上得到了y的精确值，但我正试图从中学习如何进行区间回归。首先，我把y转换成一个区间

---

Income[Income<10000]10000&Income<20001]20000&Income<30001]300000]同时，我将尝试对程序包intreg（请参见此处）进行详细的问题描述，但没有实际数据。对不起，我应该如何附加数据？可以是指向可下载站点的链接，也可以使用dput创建文本“图像”在你的问题中放入代码块。好的，谢谢。我从使用intReg软件包列出的链接中得到了结果，它们看起来不错，所以我可能不再需要帮助了。虽然我也会尝试执行多项式概率回归，但我需要先自己尝试一下。不过我很感激有兴趣帮助我。
 lIncome <- rep(0,5000)#lower income bound
uIncome <- rep(0,5000)#upper income bound

for (i in 1:5000){
    if(Income[i] == 1){
        lIncome[i] = 0
        uIncome[i] = 10000      
    }
    if(Income[i] == 2){
        lIncome[i] = 10001
        uIncome[i] = 20000      
    }
    ...
    if(Income[i] == 14){
        lIncome[i] = 300001
        uIncome[i] = Inf    
    }
}

TEST <- Surv(lIncome, uIncome, event = rep(3,5000))

    Warning message:
In Surv(lIncome, uIncome, event = rep(3, 5000)) :
  Invalid status value, converted to NA

TEST <- Surv(lIncome, uIncome, event = rep(2,5000))

m <- survreg(TEST ~ Age + AgeSq + ... , dist="gaussian")

    Error in survreg(TEST ~ Age + AgeSq + NoDegree, dist = "gaussian") : 
  Invalid survival type