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

From: Rolf Turner <>
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(! ?)


         Rolf Turner mailing list PLEASE do read the posting guide 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 Please read the posting guide before posting to the list.

list of date sections of archive