[R] Accessing Results from cenmle function in NADA package
Tom La Bone
booboo at gforcecable.com
Wed Nov 12 15:30:37 CET 2008
Oh, I like your answer better than mine! Thanks.
Tom
Richard Cotton wrote:
>
>> I figured it out. In case anyone else ever has this question -- given
> the
>> following output from cenmle:
>>
>> >fit.cen <- cenmle(obs, censored, groups)
>> > fit.cen
>> Value Std. Error z p
>> (Intercept) 1.19473 0.0772 15.4695 5.58e-54
>> groups1 0.00208 0.0789 0.0263 9.79e-01
>> Log(scale) -0.40221 0.0168 -23.9034 2.82e-126
>>
>> Scale= 0.669
>>
>> Log Normal distribution
>> Loglik(model)= -3908.9 Loglik(intercept only)= -3908.9
>> Loglik-r: 0.0006265534
>>
>> Chisq= 0 on 1 degrees of freedom, p= 0.98
>> Number of Newton-Raphson Iterations: 1
>> n = 1766
>>
>> The p-value is (for example):
>>
>> p.value <- summary(fit.cen)[[9]][[11]]
>
> Or slightly more transparently
> summary(fit.cen)$table[,'p']
>
> Regards,
> Richie.
>
> Mathematical Sciences Unit
> HSL
>
>
> ------------------------------------------------------------------------
> ATTENTION:
>
> This message contains privileged and confidential inform...{{dropped:20}}
>
> ______________________________________________
> 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.
>
>
> -----
> Regards,
> Richie.
>
> Mathematical Sciences Unit
> HSL
>
--
View this message in context: http://www.nabble.com/Accessing-Results-from-cenmle-function-in-NADA-package-tp20437420p20461407.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list