0

我找到了一种方法来执行我想做的程序,但我一直坚持如何重复 100 或 1000 次。

基本上,我想从数据集(花)中的每个组(group_ID)中抽取一个随机观察值(indvd_ID)。然后,我想使用 Mann-Whitney U 检验测试 A 组或 B 组的个人的平均命中数 (No_hits) 是否更高。

我为此编写的代码是:

library(dplyr)

random_indvdIDs <- flowers %>% group_by(group_ID) %>% sample_n(1)
wilcox.test(random_indvdIDs$No_hits ~ random_indvdIDs$A_or_B)

我想我需要编写一个 for 循环,然后让它每次打印 wilcox.test 的 p 值,但我正在努力这样做!

我有点菜鸟,所以非常感谢任何知道如何提供帮助的人。:)

4

0 回答 0