R 如何在wilcoxon测试中重复代码块10000次以找到V值

所以我从10个均匀分布的RV中找到了Wilcoxon和排序检验的V值。 我可以这样做一次

uniform=runif(10,min=0,max=1)
diff <- c(uniform - .5) 
 diff <- diff[ diff!=0 ] 
 diff.rank <- rank(abs(diff)) 
 diff.rank.sign <- diff.rank * sign(diff) 
 ranks.pos <- sum(diff.rank.sign[diff.rank.sign > 0])

uniform=runif（10，最小值=0，最大值=1）
diff尝试设置一个数据帧来保存结果，并使用这样的循环。绘图将使用
ggplot2
进行设计。代码如下：

library(dplyr)
library(ggplot2)
#Dataframe for storing results
N <- 10000
df <- data.frame(Var=rep(NA,N))
#Loop
for(i in 1:N)
{
  #Code
  uniform=runif(10,min=0,max=1)
  diff <- c(uniform - .5) 
  diff <- diff[ diff!=0 ] 
  diff.rank <- rank(abs(diff)) 
  diff.rank.sign <- diff.rank * sign(diff) 
  ranks.pos <- sum(diff.rank.sign[diff.rank.sign > 0])
  #Save
  df$Var[i] <- ranks.pos
}
#Plot
ggplot(df,aes(x=Var))+
  geom_histogram(fill='cyan3',color='black')

库（dplyr）
图书馆（GG2）
#用于存储结果的数据帧
N
hist（replicate（10000，{your code}））
Awesome看起来对我也很有用！现在我只需要找到5%sig lvl的排斥区域，谢谢@尼克：太好了！如果您有兴趣帮助，请检查此项