r - 分层 Fisher 和 Wilcoxon 检验

Question

我想进行分层费希尔测试。

我试过tabulate没有成功。

这些是我的数据：

> db$site
    [1] 2 2 2 3 1 1 1 1 1 1 1 1 1 1 2 3 1 1 1 1 1 1 1 1 2 2 3 3 1 1 3 1 2 1 1 2 1 1 1 1 1 3 1 1 1 1 1 3 1
 [50] 2 1 1 3 3 3 1 1 1 1 1 1 1 1 1 1 1 1 1 3 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 1 3 3 1
 [99] 3 1 1 1 1 1 1 1 3 3 1 1 1 1 1 3 2 1 1 1 1 1 3 3 1 3 3 3 3 1 3 1 3 3 1 3 1 1 3 3 3 2 3 3 3 3 1 3 3
[148] 3 2 3 3 1 3 1 3 3 3 3 3 3 1 3 3 3 1 3 3 1 3 1 1 1 1 1 1 1 3 2 1 3 2 2 2 3 2 3 2 2 2 2 2 2 2 2 3 2
[197] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 3 1 3 2 1 2 3 1 3 3 1 1 1 1 3 3 3 1 3 2 2 1 3 1 2 3 1
[246] 1 1 1 1 1 2 3 2 1 2 3 3 3 1 1 2 3 2 3 3 3 2 2 1 3 2 3 1 1 3 3 2 1 1 1 1 1 2 1 1 2 2 1 2 3 3 1 1 1
[295] 3 2 2 3 1 1 2 2 2 3 3 2 1 3 1 2 1 3 1 1 3 1 1 3 2 2 2 2 2 1 3 1 1 2 3 3 3 1 3 1 3 2 3 1 1 1 3 3 3
[344] 3 1 2 2 2 3 1 3 1 1 3 1 3 2 1 3 2 2 2 2 2 2 2
Levels: 1 2 3

> db$phq_cat
  [1] 1 2 2 3 1 2 2 2 1 1 1 2 1 1 1 2 2 1 2 1 3 2 1 1 1 5 1 2 3 2 3 1 2 4 2 1 1 2 2 1 1 1 1 2 1 2 2 2 2
 [50] 2 1 3 1 2 3 2 2 2 3 2 1 1 3 2 2 2 2 2 3 1 1 2 3 2 2 5 3 1 3 1 2 3 2 2 3 3 1 3 1 1 2 2 1 2 2 1 2 4
 [99] 1 1 2 2 2 2 2 1 3 2 2 1 1 3 2 1 2 2 2 1 3 2 2 3 2 1 1 1 2 2 2 2 1 3 1 3 2 2 1 2 2 2 2 1 1 3 2 2 1
[148] 3 2 4 1 2 1 2 3 1 3 2 2 2 2 2 2 1 2 2 2 2 2 2 1 1 1 1 1 3 3 2 1 1 3 2 2 3 3 2 4 2 2 2 3 2 1 4 1 2
[197] 3 1 2 2 2 2 3 2 3 3 3 3 5 1 3 2 3 1 3 3 3 2 2 1 2 2 1 3 4 4 2 1 2 2 3 4 3 3 2 1 2 4 1 1 2 1 3 2 3
[246] 3 1 1 2 2 3 3 1 1 4 3 2 1 1 2 2 2 2 1 2 2 2 3 1 1 2 4 1 1 2 3 3 2 1 1 1 2 5 1 2 3 2 2 3 1 1 3 3 1
[295] 3 5 2 1 1 1 2 2 1 1 1 4 4 2 1 1 2 1 2 3 1 2 1 2 4 1 1 1 3 1 4 1 1 1 1 4 3 1 2 1 3 3 3 1 2 1 3 1 1
[344] 2 1 1 2 2 2 1 4 2 2 1 4 1 1 4 2 1 4 2 3 2 1 1
Levels: 1 2 3 4 5

>  db$area
  [1] 3    1    1    0    3    2    5    3    1    3    3    3    1    5    4    0    3    5    5   
 [20] 1    3    2    3    1    3    3    0    0    3    3    0    1    3    1    3    3    3    3   
 [39] 3    3    3    0    1    2    2    2    2    0    2    1    1    1    0    0    0    3    3   
 [58] 3    3    4    3    3    3    3    3    1    3    3    0    3    3    5    3    2    3    5   
 [77] 5    3    3    3    3    5    3    2    5    3    3    3    2    0    0    0    0    0    5   
 [96] 0    0    3    0    3    5    3    3    3    1    3    0    0    3    3    2    1    3    0   
[115] 3    2    5    2    5    1    0    0    5    0    0    0    0    1    0    1    0    0    3   
[134] 0    3    3    0    0    0    3    0    0    0    0    3    0    0    0    3    0    0    2   
[153] 0    3    0    0    0    0    0    0    3    0    0    0    3    0    0    3    0    5    3   
[172] 5    3    3    3    3    0    2    3    0    3    2    3    0    3    0    3    2    5    3   
[191] 2    3    5    5    0    3    5    3    2    3    3    3    3    2    2    3    1    3    3   
[210] 3    3    3    5    3    3    3    3    3    0    3    0    3    2    1    0    3    0    0   
[229] 5    3    3    1    0    0    0    3    0    5    1    3    0    3    3    0    1    5    5   
[248] 3    1    3    5    0    2    3    2    0    0    0    3    3    5    0    3    0    0    0   
[267] 3    3    3    0    3    0    3    4    0    0    3    3    3    3    5    5    3    3    1   
[286] 3    1    2    3    0    0    5    3    1    0    5    3    0    3    3    3    2    3    0   
[305] 0    <NA> 3    0    3    5    5    0    5    3    0    3    1    0    5    3    3    3    2   
[324] <NA> 0    3    5    1    0    0    0    5    0    5    0    5    0    3    3    3    0    0   
[343] 0    0    2    3    3    3    0    5    0    3    3    0    3    0    5    1    0    3    3   
[362] 2    3    3    3    3   
Levels: 0 1 2 3 4 5

library(survival)
library(Exact)
library(plyr)

b<-tabulate(db$site, db$phq_cat, db$area,  tests=c("fisher"))

我收到此错误消息：

制表错误(db$site, db$phq_cat, db$area, tests = c("fisher")) : 未使用的参数 (db$AREEDISCIPL, tests = c("fisher"))

我该如何处理？

我也想进行分层 wilcoxon 秩和检验。有办法吗？

谢谢！

标签： rplyrsurvival