n = 80; x = rep(1, n)
dat = rARCens(n=n, beta=2, phi=.6, sig2=.3, x=x, cens='right', pcens=.05,
innov='t', nu=4)
fit = ARtCensReg(dat$data$cc, dat$data$lcl, dat$data$ucl, dat$data$y, x,
M=5, pc=.12, tol=0.001)
summary(fit)
Run the code above in your browser using DataLab