運行預測（）封裝內託比特（）AER

Question

後，我做的一個數據集，其中因變量（可以稱之爲Y）一託比特分析留在截尾0，所以這是我做的：

library(AER) 
fit <- tobit(data=mydata,formula=y ~ a + b + c) 

這很好。現在我想運行「預測」函數來獲得擬合的值。理想情況下，我對未觀測的潛變量「y *」和觀察到的刪失變量「y」的預測值感興趣[參見參考文獻1]。

我檢查了predict.survreg [參考文獻2]的文檔，我不認爲我理解哪個選項給了我預測的截尾變量（或潛在變量）。

大多數例子我在網上找到告知：

predict(fit,type="response"). 

再次，它不清楚這些是什麼樣的預測。

我的猜測是預測函數中的「類型」選項在這裏是關鍵，type =「response」用於檢查變量預測，type =「linear」用於潛在變量預測。

有人在這裏有一些經驗，請給我一些亮光嗎？

非常感謝！

參考文獻：

1. http://en.wikipedia.org/wiki/Tobit_model

2. http://astrostatistics.psu.edu/datasets/2006tutorial/html/survival/html/predict.survreg.html

Answer 1

一般predict- 「響應」 結果已經背轉化爲數據的原始規模無論從任何模型變換爲用於迴歸，而「線性」預測是鏈接變換尺度上的線性預測。在具有身份鏈接的導軌的情況下，它們應該是相同的。

你可以很容易地檢查我的元預測。我只是用例子檢查了它的?tobit頁：

plot(predict(fm.tobit2, type="response"), predict(fm.tobit2,type="linear")) 

Answer 2

我貼在stats.stackexchange了類似的問題，我得到的回答可能對您有用：

https://stats.stackexchange.com/questions/149091/censored-regression-in-r

有包的作者之一顯示瞭如何計算$ Y $的平均值（即預測值），其中$ Y = max（Y^*，0）$。使用包裝AER這必須有點「手工」完成。