Re: [R] Integrate na.rm in own defined functions

From: Rolf Turner <rolf.turner_at_xtra.co.nz>
Date: Thu, 21 Apr 2011 10:16:41 +1200

On 20/04/11 22:25, vioravis wrote:
> This should work!!
>
> rmse<-function (x){
> dquared<-x^2
> sum1<-sum(x^2,na.rm=TRUE)
> rmse<-sqrt((1/length(x))*sum1)
> rmse}

Shouldn't the divisor be the number of non-missing values in x? Rather than the length of x? (Like, e.g. sum(!is.na(x)) ?)

     cheers,

         Rolf Turner



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. Received on Wed 20 Apr 2011 - 22:20:25 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 20 Apr 2011 - 22:40:32 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive