Re: [R] Accessing Results from cenmle function in NADA package

From: Tom La Bone <>
Date: Wed, 12 Nov 2008 05:48:51 -0800 (PST)

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]]


Tom La Bone wrote:
> The cenmle function is used to fit two sets of censored data and test if
> they are significantly different. I can print out the results of the
> analysis on the screen but can't seem to figure out how to access these
> results in R and assign them to new variables, e.g., assign the slope
> calculated with cenmle to the variable m. Any suggestions?
> Tom

View this message in context:
Sent from the R help mailing list archive at

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Wed 12 Nov 2008 - 13:52:07 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 12 Nov 2008 - 15:30:26 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive