R help - May 2007 - error using boxplot.stats (but boxplot works¿?)

Hi

The answer to this may be obvious, but it's got me floored.

I'm unable to get boxplot.stats to work for me!

My session looks something like this:

	> ia=read.table('/tmp/prueba.csv', header=TRUE, sep=",")
	> attach(ia)

	> boxplot.stats(X8weeks~Orden)

	Error in sort (na.last, decreasing, ...) : 
        	argument 1 is not a vector

	In addition: Warning messages:
	1: is.na() applied to non-(list or vector) in: is.na(x) 
	2: is.na() applied to non-(list or vector) in: is.na(x) 
	3: is.na() applied to non-(list or vector) in: FUN(X[[1]], ...) 

	[the error message in R 2.3.1 is slightly different: it says it is not 
	atomic and shows only the 1st two warning messages]

However, boxplot(X8weeks~Orden) prints the boxplot properly (no error messages,
data accurately represented).
> is.numeric(X8weeks)=TRUE
> is.factor(Orden)=TRUE
"X8weeks" has no NA's 


So, the question is: how can I access the data used to generate the boxplot (and
why does boxplot.stats fail)?


Thanks!

Hugo Fern?ndez

-------------
	> sessionInfo()
	R version 2.4.1 (2006-12-18) 
	i486-pc-linux-gnu 
	
	locale:
	LC_CTYPE=es_ES.UTF-8 at euro;LC_NUMERIC=C;LC_TIME=es_ES.UTF-8 at euro;
	LC_COLLATE=es_ES.UTF-8 at euro;LC_MONETARY=es_ES.UTF-8 at euro;
	LC_MESSAGES=es_ES.UTF-8 at euro;LC_PAPER=es_ES.UTF-8 at euro;
	LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;
	LC_MEASUREMENT=es_ES.UTF-8 at euro;LC_IDENTIFICATION=C

	attached base packages:
	[1] "stats"     "graphics"  "grDevices"
"utils"     "datasets"  "methods"
	[7] "base"

On 06/05/2007 4:34 AM, Hugo M.Fernandez Bellon wrote:
> Hi
> 
> The answer to this may be obvious, but it's got me floored.
> 
> I'm unable to get boxplot.stats to work for me!
> 
> My session looks something like this:
> 
> 	> ia=read.table('/tmp/prueba.csv', header=TRUE,
sep=",")
> 	> attach(ia)
> 
> 	> boxplot.stats(X8weeks~Orden)
> 
> 	Error in sort (na.last, decreasing, ...) : 
>         	argument 1 is not a vector
See ?boxplot.stats.  It can't handle a formula argument, it wants a 
numeric vector.
> 
> 	In addition: Warning messages:
> 	1: is.na() applied to non-(list or vector) in: is.na(x) 
> 	2: is.na() applied to non-(list or vector) in: is.na(x) 
> 	3: is.na() applied to non-(list or vector) in: FUN(X[[1]], ...) 
> 
> 	[the error message in R 2.3.1 is slightly different: it says it is not 
> 	atomic and shows only the 1st two warning messages]
> 
> However, boxplot(X8weeks~Orden) prints the boxplot properly (no error
messages, data accurately represented).
You can use boxplot(X8weeks ~ Orden, plot=FALSE).

Duncan Murdoch
> 
>> is.numeric(X8weeks)=TRUE
>> is.factor(Orden)=TRUE
> "X8weeks" has no NA's 
> 
> 
> So, the question is: how can I access the data used to generate the boxplot
(and why does boxplot.stats fail)?
> 
> 
> Thanks!
> 
> Hugo Fern?ndez
> 
> -------------
> 	> sessionInfo()
> 	R version 2.4.1 (2006-12-18) 
> 	i486-pc-linux-gnu 
> 	
> 	locale:
> 	LC_CTYPE=es_ES.UTF-8 at euro;LC_NUMERIC=C;LC_TIME=es_ES.UTF-8 at euro;
> 	LC_COLLATE=es_ES.UTF-8 at euro;LC_MONETARY=es_ES.UTF-8 at euro;
> 	LC_MESSAGES=es_ES.UTF-8 at euro;LC_PAPER=es_ES.UTF-8 at euro;
> 	LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;
> 	LC_MEASUREMENT=es_ES.UTF-8 at euro;LC_IDENTIFICATION=C
> 
> 	attached base packages:
> 	[1] "stats"     "graphics"  "grDevices"
"utils"     "datasets"  "methods"
> 	[7] "base"
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.