如何在R中获取3个可变表的行百分比?

I have following code:

我有如下代码:

> library(MASS)
> with(bacteria, table(trt, y, week))
, , week = 0

         y
trt        n  y
  placebo  2 19
  drug     2 12
  drug+    1 14

, , week = 2

         y
trt        n  y
  placebo  1 19
  drug     2 11
  drug+    1 10

, , week = 4

         y
trt        n  y
  placebo  4 14
  drug     5  7
  drug+    2 10

, , week = 6

         y
trt        n  y
  placebo  1 16
  drug     5  6
  drug+    5  7

, , week = 11

         y
trt        n  y
  placebo  4 16
  drug     4  8
  drug+    4  8

How can I get row percentages for each row in above tables. I tried prop.table() which works well for 2 variable table but for above it does not produce correct row percentages. Thanks for your help.

如何获取表中每一行的行百分比。我尝试了prop.table()，它对于两个可变表来说效果很好，但是对于上面的表，它并不能产生正确的行百分比。谢谢你的帮助。

Edit: I tried:

编辑:我试着:

prop.table(with(bacteria, table(trt, y, week)))
prop.table(with(bacteria, table(trt, y, week)),1)
prop.table(with(bacteria, table(trt, y, week)),2)
prop.table(with(bacteria, table(trt, y, week)),3)

None of above give correct row percentages, which would be following for first table:

以上没有一个给出正确的行百分比，这将用于第一个表:

         y
trt        n  y
  placebo  0.095 0.905
  drug     0.143 0.857
  drug+    0.067 0.933

Alternatively, getting these percentages in data.frame obtained by above table would also be perfect:

或者，将这些百分比作为数据。上表得到的框架也将是完美的:

> data.frame(with(bacteria, table(trt, y, week)))
       trt y week Freq
1  placebo n    0    2
2     drug n    0    2
3    drug+ n    0    1
4  placebo y    0   19
5     drug y    0   12
6    drug+ y    0   14
7  placebo n    2    1
8     drug n    2    2
9    drug+ n    2    1
10 placebo y    2   19
11    drug y    2   11
12   drug+ y    2   10
13 placebo n    4    4
14    drug n    4    5
15   drug+ n    4    2
16 placebo y    4   14
17    drug y    4    7
18   drug+ y    4   10
19 placebo n    6    1
20    drug n    6    5
21   drug+ n    6    5
22 placebo y    6   16
23    drug y    6    6
24   drug+ y    6    7
25 placebo n   11    4
26    drug n   11    4
27   drug+ n   11    4
28 placebo y   11   16
29    drug y   11    8
30   drug+ y   11    8

1 个解决方案

#1

Here are a couple of options (h.t. @user20650 for Option 1)

这里有几个选项(h.t. @user20650选项1)

library(MASS)
x <- with(bacteria, table(trt, y, week))

## Option 1
out1 <- prop.table(x, c(1,3))

## Option 2
## apply() to perform prop.table row-wise;
## aperm() to permute the output back to same form as in input:
out2 <- aperm(apply(x,c(1,3),prop.table), c(2,1,3))

out2
# , , week = 0
# 
#          
# trt                n         y
#   placebo 0.09523810 0.9047619
#   drug    0.14285714 0.8571429
#   drug+   0.06666667 0.9333333
# 
# , , week = 2
# 
#          
# trt                n         y
#   placebo 0.05000000 0.9500000
#   drug    0.15384615 0.8461538
#   drug+   0.09090909 0.9090909
# 
# , , week = 4
# 
#          
# trt               n         y
#   placebo 0.2222222 0.7777778
#   drug    0.4166667 0.5833333
#   drug+   0.1666667 0.8333333
# 
# , , week = 6
# 
#          
# trt                n         y
#   placebo 0.05882353 0.9411765
#   drug    0.45454545 0.5454545
#   drug+   0.41666667 0.5833333
# 
# , , week = 11
# 
#          
# trt               n         y
#   placebo 0.2000000 0.8000000
#   drug    0.3333333 0.6666667
#   drug+   0.3333333 0.6666667