search for: rje42

...,1L), c(1L,1L,1L,1L,1L), 1L, c(2L,2L,2L,2L,2L), 5L, ptr1, ptr2)[[7]] > x [1] -1 -1 1 1 1 1 -1 -1 1 1 -1 -1 -1 -1 1 1 1 1 -1 -1 -1 -1 1 1 -1 -1 1 1 1 1 -1 -1 -- Robin Evans Statistical Laboratory University of Cambridge blog: itsastatlife.blogspot.com web: www.statslab.cam.ac.uk/~rje42 Causal Inference Workshop July 15th: http://www.statslab.cam.ac.uk/~rje42/uai13/main.htm

..., 1, 0.5) system.time(sapply(x, function(x) isTRUE(all.equal(x, 0)))) system.time(abs(x) < .Machine$double.eps^0.5) The latter version is fast, but potentially dangerous. Any suggestions? Thanks, Robin -- Robin Evans Statistics Department University of Washington www.stat.washington.edu/~rje42

Hello, I'm trying to make a table like windfreq.dat in the rose diagram example of the climatol package. It looks like: N NNE NE ENE E ESE SE SSE S SSW SW WSW W WNW NW NNW 0-3 59 48 75 90 71 15 10 11 14 20 22 22 24 15 19 33 3-6 3 6 29 42 11 3 4 3 9 50 67 28

...al function with the same name (I presume a C function), and I'd like to see how it works. I've searched through the header files in R-2.x.x/include/ (Windows version), but to no avail. Thanks, Robin -- Robin Evans Statistics Department University of Washington www.stat.washington.edu/~rje42