[R] extractAIC

2011-08-04 Thread Roby
Use extractAIC in frailty cox model (estimated with coxph function, gaussian random effect) i obtaided > extractAIC(fit.cox.f) [1] 11.84563 8649.11736 but I don't know why I can't use the classic formulation of the AIC where the degree of freedom are the number of the parameter (in my case 3).

Re: [R] coxme and random factors

2011-07-30 Thread Roby
But, if I use coxph() function with gaussian frailty effect and I extract the residuals from this model, it's correct? I noticed that I obtained different estimates for same coefficients (very small). -- View this message in context: http://r.789695.n4.nabble.com/coxme-and-random-factors-tp32366

Re: [R] coxme frailty model standard errors?

2011-07-26 Thread Roby
Hi, but why we do the difference : ltemp <- 2 * diff(tfit$loglik[1:2]) ?? Where I can find information about Integrate Likelihooh and null like lihood?? Thank you very much, Roby -- View this message in context: http://r.789695.n4.nabble.com/coxme-frailty-model-standard-err

[R] Applying a function to a list of arguments ...

2008-11-07 Thread roby . brunelli
How can I apply function f, that I get as an argument as in func <- function(f, ...) { . . . } to a list of arguments list(a, b, c) (eg the ... argument of func above) in order to obtain f(a, b, c) Thanks a lot, Roberto [[alternative HTML version deleted]] ___