我找到了一种方法来执行我想做的程序,但我一直坚持如何重复 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 值,但我正在努力这样做!
我有点菜鸟,所以非常感谢任何知道如何提供帮助的人。:)