Re: [R] Multidimensional contingency tables

From: hadley wickham <h.wickham_at_gmail.com>
Date: Tue, 22 Apr 2008 08:54:02 -0500

On Tue, Apr 22, 2008 at 8:24 AM, hadley wickham <h.wickham_at_gmail.com> wrote:
> On Mon, Apr 21, 2008 at 9:46 PM, Robert A. LaBudde <ral@lcfltd.com> wrote:
> > How does one ideally handle and display multidimenstional contingency
> > tables in R v. 2.6.2?
> >
> > E.g.:
> >
> > > prob1<- data.frame(victim=c(rep('white',4),rep('black',4)),
> > + perp=c(rep('white',2),rep('black',2),rep('white',2),rep('black',2)),
> > + death=rep(c('yes','no'),4), count=c(19,132,11,52,0,9,6,97))
> > > prob1
> > victim perp death count
> > 1 white white yes 19
> > 2 white white no 132
> > 3 white black yes 11
> > 4 white black no 52
> > 5 black white yes 0
> > 6 black white no 9
> > 7 black black yes 6
> > 8 black black no 97
> >
> > The xtabs() function doesn't seem appropriate, as it has no means of
> > using 'count'.
> >
> > This must be a common problem.
>
> You can also use the reshape package (http://had.co.nz/reshape)
>
> cast(prob1, victim ~ perp, sum, value="count")
> cast(prob1, victim ~ perp ~ death, sum, value="count")
> cast(prob1, death + victim ~ perp, sum, value="count")

Ooops, that's using my development version. For the currently available version you need:

prob1 <- rename(prob1, c("count" = "value"))

cast(prob1, victim ~ perp, sum)
cast(prob1, victim ~ perp ~ death, sum)
cast(prob1, death + victim ~ perp, sum)

Hadley

-- 
http://had.co.nz/

______________________________________________
R-help_at_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.
Received on Tue 22 Apr 2008 - 14:02:02 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 22 Apr 2008 - 14:30:31 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive