Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- df = data.frame(x_1 = rnorm(200),
- x_2 = rnorm(200),
- group = c(rep('A',100),rep('B',100)),
- clus = c(rep('a_1',50),rep('a_2',50),rep('b_1',50),rep('b_2',50)))
- clusWilcox.test(x_1,paired = TRUE,cluster = "clus",data = df)
- temp_test <-
- df %>%
- group_by(group) %>%
- summarise_each(funs(clusWilcox.test(.,paired = TRUE,cluster = "clus")$p.value), vars = c('x_1','x_2'))
- Error in clusWilcox.test.default(x_1, paired = FALSE, cluster = "clus") :
- 'group' is required for the clustered rank sum test
- temp_test <-
- df %>%
- group_by(group) %>%
- summarise_each(funs(clusWilcox.test(.,paired = TRUE,cluster = "clus",data = df)$p.value), vars = c('x_1','x_2'))
- > temp_test
- # A tibble: 2 x 3
- group vars1 vars2
- <fct> <dbl> <dbl>
- 1 A 0.168 0.136
- 2 B 0.168 0.136
- temp_test <-
- df %>%
- group_by(group) %>%
- summarise_each(funs(t.test(.)$p.value), vars = c('x_1','x_2'))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement