Skip to content

crosstabulation

1 message · Victor Moreno

Original
Victor Moreno
# 
Hi,
I find difficult to read crosstabulated data without percentages, so I wrote 
this dirty function that may be useful to others. It only works up to 3 
dimensions since I am not very good at programming in R. I would appreciate 
if someone else has a better similar function or can improve this one.

As an example of use:
x     (%)
0    71 ( 71.0)
1    29 ( 29.0)
y
x   0     (%)     1     (%)
  0    52 ( 73.2)    19 ( 26.8)
  1    17 ( 58.6)    12 ( 41.4)
y
x   0     (%)     1     (%)
  0    52 ( 75.4)    19 ( 61.3)
  1    17 ( 24.6)    12 ( 38.7)
z = 0
   y
x   0     (%)     1     (%)
  0    38 ( 71.7)    15 ( 28.3)
  1    14 ( 63.6)     8 ( 36.4)
 
 z = 1
   y
x   0     (%)     1     (%)
  0    14 ( 77.8)     4 ( 22.2)
  1     3 ( 42.9)     4 ( 57.1)
#################################################
Table<-
function(..., margin = 1)
{
proportion.table<-
function(data, MARGIN = 1)
{
if(is.null(d <- dim(data)))
stop("data is not an array")
if(any(data < 0) || any(trunc(data) != data))
stop("data is not an array of counts")
if(length(MARGIN))
sweep(data, MARGIN = MARGIN, apply(data, MARGIN = MARGIN, sum), "/")
else data/sum(data)
}
tt <- table(...)
dd <- dim(tt)
dnam <- dimnames(tt)
if(length(dd) > 3)
stop("max 3 dimensions")
if(length(dd) == 1) {
dd <- c(dd, 1)
dim(tt) <- dd
margin <- 2
dnam<-list(dnam[[1]],names(dnam[1]))
}
dd[2] <- dd[2] * 2
rr <- array(dim = dd)
rr.lab <- rep(dnam[[2]], 2)
if(length(dd) == 3) {
for(j in 1:dd[3]) {
pp <- proportion.table(tt[,  , j], margin)
for(i in seq(1, dd[2], 2)) {
rr[, i, j] <- formatC(tt[, i/2 + 0.5, j], 0,5,format="f" )
rr[, i + 1, j] <- paste("(", formatC(100 * pp[, i/2 + 0.5], 1, 5,format="f"), 
")", sep = "")
rr.lab[i] <- dnam[[2]][i/2 + 0.5]
rr.lab[i + 1] <- "(%)"
} } }
else {
# dim<=2
pp <- proportion.table(tt, margin)
for(i in seq(1, dd[2], 2)) {
rr[, i] <- formatC(tt[, i/2 + 0.5], 0, 5,format="f")
rr[, i + 1] <- paste("(", formatC(100 * pp[, i/2 + 0.5], 1, 5,format="f"), 
")", sep = "")
rr.lab[i] <- dnam[[2]][i/2 + 0.5]
rr.lab[i + 1] <- "(%)"
}
}
dnam[[2]] <- rr.lab
dimnames(rr) <- dnam

if(length(dim(rr)) == 2) {
print(rr, quote = F)
}
else {
for(i in 1:dim(rr)[3]) {
cat("\n",names(dnam[3]),"=", dimnames(rr)[[3]][i], "\n")
print(rr[,  , i], quote = F)
}
}

invisible(rr)
}

#################################################