Table of Percentages with Percentage Base

Description

percent returns a table of percentages along with the percentage base. It will be useful in conjunction with Aggregate or genTable.

Usage

  percent(x,...)
  ## S4 method for signature 'default'
percent(x,weights=NULL,total=!(se || ci),
      se=FALSE,ci=FALSE,ci.level=.95,
      total.name="N",perc.label="Percentage",...)
  ## S4 method for signature 'logical'
percent(x,weights=NULL,total=!(se || ci),
      se=FALSE,ci=FALSE,ci.level=.95,
      total.name="N",perc.label="Percentage",...)

Arguments

x

a numeric vector or factor.

weights

a optional numeric vector of weights of the same length as x.

total

logical; should the total sum of counts from which the percentages are computed be included into the output?

se

logical; should standard errors of the percentages be included?

ci

logical; should confidence intervals of the percentages be included?

ci.level

numeric; nominal coverage of confidence intervals

total.name

character; name given for the total sum of counts

perc.label

character; label given for the percentages if the table has more than one dimensions, e.g. if se or ci is TRUE.

...

for percent.mresp: one or several 1-0 vectors or matrices otherwise, further arguments, currently ignored.

Value

A table of percentages.

Examples

x <- rnorm(100)
y <- rnorm(100)
z <- rnorm(100)
f <- sample(1:3,100,replace=TRUE)
f <- factor(f,labels=c("a","b","c"))
percent(x>0)
Percentage          N
        51        100
percent(f)
 a   b   c   N
36  26  38 100
genTable(
 cbind(percent(x>0),
       percent(y>0),
       percent(z>0)) ~ f
 )
, , f = a

                    1        2        3
  Percentage 47.22222 44.44444 44.44444
  N          36.00000 36.00000 36.00000

, , f = b

                    1        2        3
  Percentage 57.69231 57.69231 57.69231
  N          26.00000 26.00000 26.00000

, , f = c

                    1        2        3
  Percentage 50.00000 47.36842 34.21053
  N          38.00000 38.00000 38.00000
gt <- genTable(
 cbind(percent(x>0,ci=TRUE),
       percent(y>0,ci=TRUE),
       percent(z>0,ci=TRUE)) ~ f
 )
print(gt)
, , f = a

                    1        2        3
  Percentage 47.22222 44.44444 44.44444
  lower      30.40506 27.93542 27.93542
  upper      64.51362 61.90232 61.90232

, , f = b

                    1        2        3
  Percentage 57.69231 57.69231 57.69231
  lower      36.91804 36.91804 36.91804
  upper      76.64780 76.64780 76.64780

, , f = c

                    1        2        3
  Percentage 50.00000 47.36842 34.21053
  lower      33.37890 30.98073 19.63291
  upper      66.62110 64.18166 51.35268
ftable(gt,row.vars=2,col.vars=c(3,1))
f a b c
Percentage lower upper Percentage lower upper Percentage lower upper

1 47.22222 30.40506 64.51362 57.69231 36.91804 76.64780 50.00000 33.37890
  66.62110
2 44.44444 27.93542 61.90232 57.69231 36.91804 76.64780 47.36842 30.98073
  64.18166
3 44.44444 27.93542 61.90232 57.69231 36.91804 76.64780 34.21053 19.63291
  51.35268
ex.data <- expand.grid(mean=c(0,25,50),sd=c(1,10,100))[rep(1:9,rep(250,9)),]
ex.data <- within(ex.data,x <- rnorm(n=nrow(ex.data),mean=ex.data$mean,sd=ex.data$sd))
ex.data <- within(ex.data,x.grp <- cases( x < 0,
                                           x >= 0 & x < 50,
                                           x >= 50 & x < 100,
                                           x >= 100
                                         ))
genTable(percent(x.grp)~mean+sd,data=ex.data)
, , sd = 1

                   mean
                        0    25    50
  x < 0              52.4   0.0   0.0
  x >= 0 & x < 50    47.6 100.0  52.0
  x >= 50 & x < 100   0.0   0.0  48.0
  x >= 100            0.0   0.0   0.0
  N                 250.0 250.0 250.0

, , sd = 10

                   mean
                        0    25    50
  x < 0              46.8   1.6   0.0
  x >= 0 & x < 50    53.2  98.0  45.6
  x >= 50 & x < 100   0.0   0.4  54.4
  x >= 100            0.0   0.0   0.0
  N                 250.0 250.0 250.0

, , sd = 100

                   mean
                        0    25    50
  x < 0              49.6  36.8  26.4
  x >= 0 & x < 50    20.4  20.4  20.4
  x >= 50 & x < 100  14.8  21.2  20.0
  x >= 100           15.2  21.6  33.2
  N                 250.0 250.0 250.0
Aggregate(percent(Admit,weight=Freq)~Gender+Dept,data=UCBAdmissions)
   Gender Dept  Admitted Rejected   N
1    Male    A 62.060606 37.93939 825
3  Female    A 82.407407 17.59259 108
5    Male    B 63.035714 36.96429 560
7  Female    B 68.000000 32.00000  25
9    Male    C 36.923077 63.07692 325
11 Female    C 34.064081 65.93592 593
13   Male    D 33.093525 66.90647 417
15 Female    D 34.933333 65.06667 375
17   Male    E 27.748691 72.25131 191
19 Female    E 23.918575 76.08142 393
21   Male    F  5.898123 94.10188 373
23 Female    F  7.038123 92.96188 341