similar to: pch=NA (PR#7737)

I'd like to iterate my earlier request (#7737) to change the documentation for bxp(). The argument outpch=" " needs to be replaced with outpch=NA in two places. I actually wrote this part of the documentation myself at one point, but have now realized that pch=NA and pch=" " are not the same: x <- split(rlnorm(26e4), letters) ## NA generates small file

Dear r-bugs, I'd like to be able to draw boxplots with solid or dotted whisker lines. The lty argument is currently ignored by boxplot(), and indeed hardwired in bxp(). It's not very difficult to customize these functions to generate various types of boxplots, but generally it might be an idea to allow the user to pass a few more graphical arguments to boxplot(). The approach might as

The write.csv() function is currently implemented as function (..., col.names=NA, sep=",", qmethod="double") { write.table(..., col.names=NA, sep=",", qmethod="double") } Surely, it should be function (..., col.names=NA, sep=",", qmethod="double") { write.table(..., col.names=col.names, sep=sep,

Hi, I was just trying to figure out how to beautify the output of my bwplot-output. Altogether I figured most of the things out on my own. The one thing which puzzles me though are the symbols for the outliers. I can easily change the form of the median symbol by using "pch" but I don't know how to do this for outliers. Obviously the "outpch" of the

Full_Name: Arni Magnusson Version: 1.6.2 OS: Windows XP Submission from: (NULL) (210.48.49.68) Hi there. I'm experiencing unexpected behaviour from as.POSIXct: > as.POSIXct("1969-12-24") [1] "1969-12-23 23:00:00 New Zealand Standard Time" > as.POSIXlt("1969-12-24") [1] "1969-12-24" > as.POSIXlt("1969-12-24")+1 [1] "1969-12-23

One descriptive statistic that is conspicuously missing from core R is the statistical mode - the most frequent value in a discrete distribution. I would like to propose adding the attached 'statmode' (or a similar function) to the 'stats' package. Currently, it can be quite cumbersome to calculate the mode of a distribution in R, both for experts and beginners. The lack of a

Dear r-bugs, The 'see also' link in help(area,html=T) is broken, see line 91 in library/MASS/html/area.html: "../../integrate/html/integrate.html" should be "../../base/html/integrate.html" Regards, Arni OS: Windows XP R: 1.7.1 MASS: 7.1-8

In Windows XP, the matrix() function crashes the program when 'dimnames' is an empty list: matrix(1:4, nrow=2, dimnames=list()) # R has encountered a problem and needs to close ... This bug is specific to WinXP, as Linux64 handles this situation more gracefully: matrix(1:4, nrow=2, dimnames=list()) Error in matrix(1:4, nrow = 2, dimnames = list()) : invalid type

I present you with a function that solves a problem that has bugged me for many years. I think the problem may be general enough to at least consider adding this function, or a revamped version of it, to the 'stats' package, with the other multinomial functions reside. I'm using R to export data to text files, which are input data for an external model written in C++. Parts of the

The weighted.mean() function replaces NA values with 0.0 when the user specifies na.rm=TRUE: x <- c(101, 102, NA) mean(x, na.rm=TRUE) # 101.5, correct weighted.mean(x, na.rm=TRUE) # 67.66667, wrong weighted.mean(x, w=c(1,1,1), na.rm=TRUE) # 67.66667, wrong weighted.mean(x, w=c(1,1,1)/3, na.rm=TRUE) # 67.66667, wrong The weights are

The model frame shows the response and predictors in a data frame with nicely labelled columns: fm <- lm(wt~qsec+log(hp)+sqrt(disp), data=mtcars) model.frame(fm) # ok When the left hand side consists of more than one response, those response variables still look good, inside a matrix: fm <- lm(cbind(qsec,hp,disp)~wt, data=mtcars) model.frame(fm)[[1]] # ok A problem arises when

hello list, i'm a bit puzzled by the error message i get when i copy past this in R:

Dear r-devel, I have implemented a function to create color palettes with improved contrast and logical order, compared with the built-in ones: source("arni.colors.R") # code is given below barplot(rep(1,100), col=arni.colors(100), space=0, border=0, axes=F) par(mfrow=c(2,1)) # rainbow() has too much green... barplot(rep(1,50), col=rev(rainbow(50,end=0.7)), space=0, border=0,

I have been reading about autocorrelation in linear models over the last couple of days, and I have to say the more I read, the more confused I get. Beyond confusion lies enlightenment, so I'm tempted to ask R-Help for guidance. Most authors are mainly worried about autocorrelation in the residuals, but some authors are also worried about autocorrelation within Y and within X vectors

Dear Martin, Thank you for the well-reasoned response. I realized I was rather late to make this suggestion for 4.0.0, changing a somewhat low-level function that can indeed affect packages. I was just reviewing some R user scripts that were using type.convert(), mainly on data frames. In all cases, people were passing as.is=TRUE, so I was reminded that I would not be the only user who would

The audit trail indicates that PR#3663 is not reproducible and/or fixed, but I believe windows(record=T) is still not working as documented: windows(record=T) for(i in 1:5) plot(0, 0, cex=4, pch=as.character(i)) After cycling through the plots with PgUp and PgDn, the user finds out that plot 5 is lost. Indeed, .SavedPlots reports that only 4 plots were saved. I'm running the example in

Hi! I 've a problem to have a lognorm distribution with mean=1 and var (or sigma)=1. rlnorm(1000,0,0) rlnorm(1000,1,1) rlnorm(1000,0,1) .... ? Can you help me?

Hey, Here is a snippet that generated a boxplot and separates points so that they do not overlap. I have a problem to avoir text overlapping. Any help would be helpful. > attach(InsectSprays) boxplot(count ~ spray, data = InsectSprays, outpch = NA) stripchart(count ~ spray, data = InsectSprays, vertical = TRUE, method = "jitter", pch = 21, col =

Hi Simon, thank you for the concise reply. Do you mean the reported behavior of drop() is not a bug? It looks like a borderline bug to me (see below), but I'm not the judge of that. If this is the intended behavior and serves an actual purpose, then that could be explicitly documented in a \note{} on the help page. Such a note would slightly reduce the surprise of users running into this

dear all, could somebody tell me how I can plot mild outliers as a circle(?) and extreme outliers as an asterisk(*) in a box-whisker plot? Thanks very much in advance -- View this message in context: http://www.nabble.com/mild-and-extreme-outliers-in-boxplot-tp25040545p25040545.html Sent from the R help mailing list archive at Nabble.com.