R 曼惠特尼U检验

我想创建一个函数，该函数可以计算R中函数的U值，并使用10000次复制的精确分布计算p值。 我的R函数

U<-function(x,y){
 n1<-length(x)
n2<-length(y)
d<-n1+n2
}

---

U为什么不使用R的专用函数：
wilcox.test（）
？我可以使用，但我被要求先组合向量x和y，然后使用秩来计算x的U值，形成一个函数。这是家庭作业问题吗？如果是，请添加自学标签