论坛上面收索，可以找到的

已找到 谢谢

謝謝樓主

R> data("Affairs")
R> aff_tob <- tobit(affairs ~ age + yearsmarried +
+ religiousness + occupation + rating, data = Affairs)
R> summary(aff_tob)
Call:
tobit(formula = affairs ~ age + yearsmarried +
religiousness + occupation + rating, data = Affairs)
Observations:
Total Left-censored Uncensored Right-censored
601 451 150 0
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 8.1742 2.7414 2.98 0.0029
age -0.1793 0.0791 -2.27 0.0234
yearsmarried 0.5541 0.1345 4.12 3.8e-05
religiousness -1.6862 0.4038 -4.18 3.0e-05
occupation 0.3261 0.2544 1.28 0.2000
rating -2.2850 0.4078 -5.60 2.1e-08
Log(scale) 2.1099 0.0671 31.44 < 2e-16
Scale: 8.25
Gaussian distribution
Number of Newton-Raphson Iterations: 4
Log-likelihood: -706 on 7 Df
Wald-statistic: 67.7 on 5 Df, p-value: 3.1e-13

出处：Applied Econometrics with R，论坛就有的

install.packages("AER")
library(AER)