[R] Print std. Error separately from mle-class object

Aziz, Muhammad Fayez aziz4 at illinois.edu
Thu Apr 5 02:43:08 CEST 2012

Thank you Jeff for your input. I did realize just after posting the question that the summary fn returns an object whose "slots" can be read and that the slot containing the coefficient information is a matrix. Here is how I got the desired results through the following code in "mypanel" function of xyplot:

	tmp <- power.law.fit(y, xmin=1) 
	coefMatrix <- slot(summary(tmp), "coef")
	print(coefMatrix["alpha", "Estimate"]) # alpha:
	print(coefMatrix["alpha", "Std. Error"]) # Std. Error
	print(slot(summary(tmp), "m2logL")) # -2 log L:

From: Jeff Newmiller [jdnewmil at dcn.davis.ca.us]
Sent: Wednesday, April 04, 2012 5:54 PM
To: Aziz, Muhammad Fayez; r-help at r-project.org
Subject: Re: [R] Print std. Error separately from mle-class object

Are you aware that the summary function normally returns a data value that you can extract values from and format to your hearts desire?
and read ?mle-class
(and try to provide a reproducible example next time)
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k
Sent from my phone. Please excuse my brevity.

"Aziz, Muhammad Fayez" <aziz4 at illinois.edu> wrote:

>I am using power.law.fit to get an mle-class object in tmp and print
>summary(tmp), coef(tmp) and logLik(tmp). I wanted to print the std.
>error for alpha separately as I want to show these values concisely in
>a graph legend. coef(summary(tmp)) displays the alpha and std. error
>jointly, while I need to print them separately on two lines.
>R-help at r-project.org mailing list
>PLEASE do read the posting guide
>and provide commented, minimal, self-contained, reproducible code.

More information about the R-help mailing list