Ok thanks Erik
Indeed, "exclude" in "xtabs" seems working whith character vectors:
> x <- c(rep(c("A","B","C"), 2))
> x
[1] "A" "B" "C" "A" "B" "C"
> xtabs(~ x, exclude = "B")
x
A C
2 2
and not directly with factors:
> x <- factor(rep(c("A","B","C"), 2))
> x
[1] A B C A B C
Levels: A B C
> levels(x)
[1] "A" "B" "C"
> xtabs(~ x, exclude = "B")
x
A B C
2 2 2
while function "table" does it:
> x <- factor(rep(c("A","B","C"), 2))
> x
[1] A B C A B C
Levels: A B C
> table(x, exclude = "B")
x
A C
2 2
However, one point remains confusing for me in function "table". When I add
a NA in the vector, it still works:
> x <- factor(c(rep(c("A","B","C"), 2), NA))
> x
[1] A B C A B C <NA>
Levels: A B C
> table(x, exclude = "B")
x
A C
2 2
but not anymore when I use argument "exclude = NULL" in function "factor":
> x <- factor(c(rep(c("A","B","C"), 2), NA), exclude = NULL)
> x
[1] A B C A B C <NA>
Levels: A B C <NA>
> table(x, exclude = "B")
x
A B C <NA>
2 2 2 1
Finally, if I remove NA, it works again:
> x <- factor(rep(c("A","B","C"), 2), exclude = NULL)
> x
[1] A B C A B C
Levels: A B C
> table(x, exclude = "B")
x
A C
2 2
it looks like "exclude" does not work when NA is a level
Regards
Matthieu
> -----Message d'origine-----
> De : Erik Iverson [mailto:iver...@biostat.wisc.edu]
> Envoyé : 24 April 2009 19:13
> À : Matthieu Lesnoff
> Cc : r-help@r-project.org
> Objet : Re: [R] argument 'exclude' in function xtabs
>
>
> > I was willing to use argument 'exclude' in function xtabs to remove
> > some levels of factors (xtabs help page says '"exclude: a vector of
> > values to be excluded when forming the set of levels of the
> classifying factors").
>
>
> I think I see what's happening, and it's a little confusing
> to me, too.
> If your classifying factor to xtabs is not actually a factor,
> but say, a character vector, it will be converted to a factor
> while using the "exclude" argument. However, if it already
> is a factor, this does not happen.
> So, for example, contrast:
>
> trt.char <- c("A","B","C","A","B","C")
> xtabs(~trt.char, exclude = "B")
>
> with
>
> trt.fac <- factor(c("A","B","C","A","B","C"))
> xtabs(~trt.fac, exclude = "B")
>
> Hope that helps,
> Erik Iverson
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
