[R] How to pass na.rm=T to a user defined function

[Jun Shen]

Because in reality the NA may appear in one variable but not others. For
example for ID=1, CL may be NA but not for others, For ID=2, V1 may be NA
etc. To keep all the IDs and all the variables in one data frame, it's
inevitable to see some NA

On Thu, Jul 28, 2016 at 10:22 PM, Jeff Newmiller <jdnewmil at dcn.davis.ca.us>
wrote:

> Why not remove it yourself before passing it to those functions?
> --
> Sent from my phone. Please excuse my brevity.
>
> On July 28, 2016 5:51:47 PM PDT, Jun Shen <jun.shen.ut at gmail.com> wrote:
> >Dear list,
> >
> >I write a small function to calculate multiple stats on multiple
> >variables
> >and export in a format exactly the way I want. Everything seems fine
> >until
> >NA appears in the data.
> >
> >Here is my function:
> >
> >do.stats <- function(data, stats.func, summary.var)
> >            as.data.frame(signif(sapply(stats.func,function(func)
> >mapply(func,data[summary.var])),3))
> >
> >A test dataset:
> >test <-
>
> >data.frame(ID=1:100,CL=rnorm(100),V1=rnorm(100),V2=rnorm(100),ALPHA=rnorm(100))
> >
> >a command like the following
> >do.stats(test, stats.func=c('mean','sd','median','min','max'),
> >summary.var=c('CL','V1', 'V2','ALPHA'))
> >
> >gives me
> >
> >         mean    sd  median   min  max
> >CL     0.1030 0.917  0.0363 -2.32 2.47
> >V1    -0.0545 1.070 -0.2120 -2.21 2.70
> >V2     0.0600 1.000  0.0621 -2.80 2.62
> >ALPHA -0.0113 0.919  0.0284 -2.35 2.31
> >
> >
> >However if I have a NA in the data
> >test$CL[1] <- NA
> >
> >The same command run gives me
> >         mean    sd  median   min  max
> >CL        * NA    NA      NA    NA   NA*
> >V1    -0.0545 1.070 -0.2120 -2.21 2.70
> >V2     0.0600 1.000  0.0621 -2.80 2.62
> >ALPHA -0.0113 0.919  0.0284 -2.35 2.31
> >
> >I know this is because those functions (mean, sd etc.) all have
> >na.rm=F by default. How can I
> >
> >pass na.rm=T to all these functions without manually redefining those
> >stats functions
> >
> >Appreciate any comment.
> >
> >Thanks for your help.
> >
> >
> >Jun
> >
> >       [[alternative HTML version deleted]]
> >
> >______________________________________________
> >R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> >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.
>
>

	[[alternative HTML version deleted]]