percent memisc 0.99.22

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
        54        100
percent(f)
 a   b   c   N
32  34  34 100
genTable(
 cbind(percent(x>0),
       percent(y>0),
       percent(z>0)) ~ f
 )
, , f = a

             percent(x > 0) percent(y > 0) percent(z > 0)
  Percentage         65.625         46.875          43.75
  N                  32.000         32.000          32.00

, , f = b

             percent(x > 0) percent(y > 0) percent(z > 0)
  Percentage       44.11765       47.05882       47.05882
  N                34.00000       34.00000       34.00000

, , f = c

             percent(x > 0) percent(y > 0) percent(z > 0)
  Percentage       52.94118       52.94118       58.82353
  N                34.00000       34.00000       34.00000
gt <- genTable(
 cbind("x > 0" = percent(x>0,ci=TRUE),
       "y > 0" = percent(y>0,ci=TRUE),
       "z > 0" = percent(z>0,ci=TRUE)) ~ f
 )
ftable(gt,row.vars=3:2,col.vars=1)
         Percentage    lower    upper
f
a x > 0    65.62500 46.80690 81.42809
  y > 0    46.87500 29.09398 65.25632
  z > 0    43.75000 26.36381 62.33743
b x > 0    44.11765 27.18502 62.11424
  y > 0    47.05882 29.77867 64.87072
  z > 0    47.05882 29.77867 64.87072
c x > 0    52.94118 35.12928 70.22133
  y > 0    52.94118 35.12928 70.22133
  z > 0    58.82353 40.69694 75.35293
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              53.6   0   0
  x >= 0 & x < 50    46.4 100  52
  x >= 50 & x < 100   0.0   0  48
  x >= 100            0.0   0   0
  N                 250.0 250 250

, , sd = 10

                   mean
                        0    25    50
  x < 0              51.2   0.0   0.0
  x >= 0 & x < 50    48.8  99.6  51.2
  x >= 50 & x < 100   0.0   0.4  48.8
  x >= 100            0.0   0.0   0.0
  N                 250.0 250.0 250.0

, , sd = 100

                   mean
                        0    25    50
  x < 0              44.0  38.8  28.4
  x >= 0 & x < 50    22.8  19.6  20.4
  x >= 50 & x < 100  11.6  16.8  17.2
  x >= 100           21.6  24.8  34.0
  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
2  Female    A 82.407407 17.59259 108
3    Male    B 63.035714 36.96429 560
4  Female    B 68.000000 32.00000  25
5    Male    C 36.923077 63.07692 325
6  Female    C 34.064081 65.93592 593
7    Male    D 33.093525 66.90647 417
8  Female    D 34.933333 65.06667 375
9    Male    E 27.748691 72.25131 191
10 Female    E 23.918575 76.08142 393
11   Male    F  5.898123 94.10188 373
12 Female    F  7.038123 92.96188 341