# NOT RUN {
# An invalid distribution (it should have "init =" on line 2)
# surveg would give an error message
mycauchy <- list(name='Cauchy',
init<- function(x, weights, ...)
c(median(x), mad(x)),
density= function(x, parms) {
temp <- 1/(1 + x^2)
cbind(.5 + atan(temp)/pi, .5+ atan(-temp)/pi,
temp/pi, -2 *x*temp, 2*temp^2*(4*x^2*temp -1))
},
quantile= function(p, parms) tan((p-.5)*pi),
deviance= function(...) stop('deviance residuals not defined')
)
survregDtest(mycauchy, TRUE)
# }
Run the code above in your browser using DataLab