ht=c(1, 1, 8, 10, NA))
> result = fitdistcens(data, 'chisq', control=list(reltol=1e-14))
> result
Fitting of the distribution ' chisq ' on censored data by maximum
likelihood
Parameters:
         estimate
df 5.060329
r-