x=rnorm(10000,0,1)
x.sample1=sample(x,100)
x.sample2=sample(x,100)
diff.median.sample=median(x.sample1)-median(x.sample2)
permutation.diff.median=rep(0,1000)
for(i in 1:1000){
ID=sample(1:200,100)
x.sample.all=c(x.sample1,x.sample2)
x.sample1.new=x.sample.all[ID]
x.sample2.new=x.sample.all[-ID]
permutation.diff.median[i]=median(x.sample1.new)-median(x.sample2.new)
}
P_values=sum(permutation.diff.median>diff.median.sample)/1000
你没有输出语句啊
在R语言中用的是 print() 函数来打印要输出的内容