This should work!! rmse<-function (x){ dquared<-x^2 sum1<-sum(x^2,na.rm=TRUE) rmse<-sqrt((1/length(x))*sum1) rmse} -- View this message in context: http://r.789695.n4.nabble.com/Integrate-na-rm-in-own-defined-functions-tp3462492p3462615.html Sent from the R help mailing list archive at Nabble.com.