Displaying 2 results from an estimated 2 matches for "newcalc".
Did you mean:
newcall
2006 Oct 05
1
unexpected behavior of boxplot(x, notch=TRUE, log="y")
A function I've been using for a while returned a surprising [to me,
given the data] error recently:
Error in plot.window(xlim, ylim, log, asp, ...) :
Logarithmic axis must have positive limits
After some digging I realized what was going on:
x <- c(10460.97, 10808.67, 29499.98, 1, 35818.62, 48535.59, 1, 1,
42512.1, 1627.39, 1, 7571.06, 21479.69, 25, 1, 16143.85, 12736.96,
2006 Oct 09
1
boxplot, notches, etc.
...,sd=1)
f <- factor(rep(1:2,each=npts))
par(mfrow=c(1,2))
boxplot(X~f,notch=TRUE)
A possible fix is to truncate the notches
(and issue a warning) when this happens,
in src/library/grDevices/R/calc.R:
[WATCH OUT FOR LINE WRAPPING BELOW!]
*** calc.R 2006-10-07 17:44:49.000000000 -0400
--- newcalc.R 2006-10-07 19:25:38.000000000 -0400
***************
*** 16,21 ****
--- 16,26 ----
if(any(out[nna])) stats[c(1, 5)] <- range(x[!out], na.rm = TRUE)
}
conf <- if(do.conf) stats[3] + c(-1.58, 1.58) * iqr / sqrt(n)
+ if (do.conf) {
+ if (conf[1]<stats[2] || co...