percent memisc 0.99.25.4

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
31  41  28 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        51.6129       35.48387       54.83871
  N                 31.0000       31.00000       31.00000

, , f = b

             percent(x > 0) percent(y > 0) percent(z > 0)
  Percentage       46.34146       56.09756       46.34146
  N                41.00000       41.00000       41.00000

, , f = c

             percent(x > 0) percent(y > 0) percent(z > 0)
  Percentage       57.14286             50             50
  N                28.00000             28             28
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    51.61290 33.06060 69.84543
  y > 0    35.48387 19.22674 54.63044
  z > 0    54.83871 36.03423 72.68350
b x > 0    46.34146 30.65594 62.57516
  y > 0    56.09756 39.74984 71.53129
  z > 0    46.34146 30.65594 62.57516
c x > 0    57.14286 37.17936 75.53761
  y > 0    50.00000 30.64710 69.35290
  z > 0    50.00000 30.64710 69.35290
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              49.6   0   0.0
  x >= 0 & x < 50    50.4 100  50.4
  x >= 50 & x < 100   0.0   0  49.6
  x >= 100            0.0   0   0.0
  N                 250.0 250 250.0

, , sd = 10

                   mean
                        0    25    50
  x < 0              49.2   1.2   0.0
  x >= 0 & x < 50    50.8  98.8  48.8
  x >= 50 & x < 100   0.0   0.0  51.2
  x >= 100            0.0   0.0   0.0
  N                 250.0 250.0 250.0

, , sd = 100

                   mean
                        0    25    50
  x < 0              48.8  37.2  34.0
  x >= 0 & x < 50    21.2  22.8  17.6
  x >= 50 & x < 100  14.4  13.6  17.2
  x >= 100           15.6  26.4  31.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
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