search for: stuetzl

Has a test for bimodality been implemented in R? Thanks, John NIWA is the trading name of the National Institute of Water & Atmospheric Research Ltd.

I have a Date x Stock (223 x 520) matrix of "trading volume". I can calculate a 5-day (past) average in about 1 second using: R> apply(vol, 1, filter, filter=c(0, rep(1/5,5)), sides=1) I would like to do the same with a 5-day median, e.g.: R> mymed <- function(x, n=5) { R> r <- rep(NA, length(x)) R> for (i in (n+1):length(x)) r[i] <- median(x[i-(1:n)]) R>

...CPP, CXXCPP and CXX98 and similar is deprecated. (CPP is found from the system make and may well not be set.) Use $CC -E and $CXX -E instead of CPP and CXXCPP. BUG FIXES: * runmed(x, *) when x contains missing values now works consistently for both algorithm="Stuetzle" and "Turlach", and no longer segfaults for "Turlach", as reported by Hilmar Berger. * apply(diag(3), 2:3, mean) now gives a helpful error message. * dgamma(x, shape, log=TRUE) now longer overflows to Inf for shape < 1 and very small x, fixing PR#1...

...CPP, CXXCPP and CXX98 and similar is deprecated. (CPP is found from the system make and may well not be set.) Use $CC -E and $CXX -E instead of CPP and CXXCPP. BUG FIXES: * runmed(x, *) when x contains missing values now works consistently for both algorithm="Stuetzle" and "Turlach", and no longer segfaults for "Turlach", as reported by Hilmar Berger. * apply(diag(3), 2:3, mean) now gives a helpful error message. * dgamma(x, shape, log=TRUE) now longer overflows to Inf for shape < 1 and very small x, fixing PR#1...

...CPP, CXXCPP and CXX98 and similar is deprecated. (CPP is found from the system make and may well not be set.) Use $CC -E and $CXX -E instead of CPP and CXXCPP. BUG FIXES: * runmed(x, *) when x contains missing values now works consistently for both algorithm="Stuetzle" and "Turlach", and no longer segfaults for "Turlach", as reported by Hilmar Berger. * apply(diag(3), 2:3, mean) now gives a helpful error message. * dgamma(x, shape, log=TRUE) now longer overflows to Inf for shape < 1 and very small x, fixing PR#1...

Dear all, I just found this issue: dd1 = c(rep(NaN,82), rep(-1, 144), rep(1, 74)) xx = runmed(dd1, 21) -> R crashes reproducibly in R 3.4.3, R3.4.4 (Ubuntu 14.04/Ubuntu 16.04) With GDB: Program received signal SIGSEGV, Segmentation fault. swap (l=53, r=86, window=window at entry=0xc59308, outlist=outlist at entry=0x12ea2e8, nrlist=nrlist at entry=0x114fdd8, print_level=print_level at

...is deprecated. (CPP is found from the system > make and may well not be set.) > > Use $CC -E and $CXX -E instead of CPP and CXXCPP. > > BUG FIXES: > > * runmed(x, *) when x contains missing values now works > consistently for both algorithm="Stuetzle" and "Turlach", and no > longer segfaults for "Turlach", as reported by Hilmar Berger. > > * apply(diag(3), 2:3, mean) now gives a helpful error message. > > * dgamma(x, shape, log=TRUE) now longer overflows to Inf for shape > < 1...

...lar is deprecated. (CPP is found from the system > make and may well not be set.) > > Use $CC -E and $CXX -E instead of CPP and CXXCPP. > > BUG FIXES: > > * runmed(x, *) when x contains missing values now works > consistently for both algorithm="Stuetzle" and "Turlach", and no > longer segfaults for "Turlach", as reported by Hilmar Berger. > > * apply(diag(3), 2:3, mean) now gives a helpful error message. > > * dgamma(x, shape, log=TRUE) now longer overflows to Inf for shape > < 1 a...

...he system > > make and may well not be set.) > > > > Use $CC -E and $CXX -E instead of CPP and CXXCPP. > > > > BUG FIXES: > > > > * runmed(x, *) when x contains missing values now works > > consistently for both algorithm="Stuetzle" and "Turlach", and no > > longer segfaults for "Turlach", as reported by Hilmar Berger. > > > > * apply(diag(3), 2:3, mean) now gives a helpful error message. > > > > * dgamma(x, shape, log=TRUE) now longer overflows to Inf for...

Hello, I apologize that this is off-topic. I am seeking information on perception of graphical data, in an effort to improve the plots I produce. Would anyone point me to literature reviews in this area? (Or keywords to try on google?) Is this located somewhere near cognitive science, psychology, human factors research? For example, some specific questions I have are: I recall as a child

...;In-Reply-To= To: r-help@r-project.org Message-ID: <25220627.post@talk.nabble.com> Content-Type: text/plain; charset=us-ascii Hi John, >> Has a test for bimodality been implemented in R? You may find the code at the URL below useful. It was written by Jeremy Tantrum (a PhD of Werner Stuetzle's). Amongst other things there is a function to plot the unimodal and bimodal Gaussian smoothers closest to the observed data. A dip-test statistic is also calculated. Regards, Mark. http://www.stat.washington.edu/wxs/Stat593-s03/Code/jeremy-unimodality.R John Sansom wrote: > > Has...