stephen sefick wrote:
> I am using the describe function in prettyR. I would like to add the
> 25% 75% quartiles to the summary table
>
> how do I do this
>
> I have tried
>
> describe(x.f, num.desc=c("mean", "median", "sd", "min", "max",
> "skewness", "quantile(x.f, na.rm=T, probs=seq(0.25, 0.75))",
> "valid.n"))
>
Hi Stephen,
If you want to run your own summary function within "describe", you have
to write a "wrapper" function that can be called by its name (unless I
radically rewrite "describe.numeric" so that a list of arguments can be
included). Additionally, the function should return a single value
(otherwise the "pretty" columnar display is messed up). So:
q25<-function(x,na.rm) {
return(quantile(x,probs=0.25,na.rm=na.rm))
}
q75<-function(x,na.rm) {
return(quantile(x,probs=0.75,na.rm=na.rm))
}
describe(x.f, num.desc=c("mean", "median", "sd", "min", "max",
"skewness","q25","q75","valid.n"))
You have anticipated the help file of the latest version (1.3) that
shows how to write a little wrapper for describe.numeric. Very good.
Jim
______________________________________________
R-help@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.
