Revision 948ad5967298b4dc174f4d96511af60f38e9279e authored by Roger Koenker on 27 July 2019, 09:38 UTC, committed by cran-robot on 27 July 2019, 09:38 UTC
1 parent 2b3a85a
Raw File
stack.R
# log likelihood for stackloss fit
require(quantreg)
data(stackloss)
logLik.rq.process <- function(fit){
    y <- model.response(model.frame(fit))
    fhat <- predict(fit, type = "fhat")
    fy <- mapply(function(f,y) f(y), fhat, y)
    sum(log(fy))
}
# First try with full process estimates
f0 <- rq(stack.loss ~ 1, tau=-1)
f1 <- rq(stack.loss ~ stack.x, tau=-1)
l0 <- logLik(f0)
l1 <- logLik(f1)
# Now try with discrete process estimates
f0 <- rq(stack.loss ~ 1, tau=1:19/20)
f1 <- rq(stack.loss ~ stack.x, tau=1:19/20)
l0 <- logLik(f0)
l1 <- logLik(f1)

back to top