[R] A Harder Score Test Question

[Charles C. Berry]

On Sat, 25 Jul 2009, Sean wrote:

> Does anyone know how get the score and information under the null from
> coxph? I know that I can get the chi-square value of the score from coxph,
> but I need the two components separately.  I have a function that computes
> the two components when I do not have ties but I would like to leverage the
> options(ties and strata components) already available in the coxph function.
> The function  coxph.detail has the score under the mle, but not the null. I
> looked through the coxph code but I don't see where the score test gets
> calculated. If there is no easy way to calculate these values from coxph I
> can always program it myself.

I'd try

  fit0 <- coxph( ... , iter.max = 0 )
  fit1 <- coxph( ... , iter.max = 1 )

and see if I have the pieces I need to work it out from vcov(fit0) and coefficients(fit1)

HTH,

Chuck

>
>
>
> Thanks in advance.
>
>
>
> Sean
>
>
>
> Sean Brummel M.S.
>
> Doctoral Candidate
>
> Dept. Statistics
>
> University of California, Irvine
>
> Irvine, CA 92797-1250
>
>
>
> Ph: 949-307-1811
>
> Fax: 949-824-9863
>
>
>
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

Charles C. Berry                            (858) 534-2098
                                             Dept of Family/Preventive Medicine
E mailto:cberry at tajo.ucsd.edu	            UC San Diego
http://famprevmed.ucsd.edu/faculty/cberry/  La Jolla, San Diego 92093-0901