I would try e.g. RcppArmadillo if I wanted to bypass those @gappy3000 I am concerned about overhead to check for NAs #rstats