n = 50; x = rep(1, n)
dat = rARCens(n=n, beta=2, phi=.5, sig2=.3, x=x, cens='left', pcens=.1)
fit = ARCensReg(dat$data$cc, dat$data$lcl, dat$data$ucl, dat$data$y, x,
M=5, pc=.12, tol=0.001, show_se=FALSE)
plot(fit)
Run the code above in your browser using DataLab