希望能得到各位的提示。
我目前的模拟代码如下,采用的是Helmut Strasser(1999)的General Independence Tests,但是power of test实在太小,不知道用什么检验方法才好?
<br />
library(coin)<br />
<br />
n.size=1000<br />
n.exp=1000<br />
res=rep(list(matrix(rep(NA,2*n.exp),ncol=2)),4)<br />
<br />
for(i.exp in 1:n.exp){<br />
for(i.dist in 1:4){<br />
dat=switch(i.dist<br />
,dat=rnorm(n.size)<br />
,dat=runif(n.size,1,10)#sample(c(runif(n.size,-3,-2),runif(n.size,111,112)),n.size)<br />
,dat=rexp(n.size)<br />
,dat=rgamma(n.size,500,.5)<br />
)<br />
res[[i.dist]][i.exp,]=c(mean(dat),var(dat))<br />
<br />
}<br />
}<br />
<br />
i=2#均匀分布情况<br />
x=res[[i]][,1]<br />
y=res[[i]][,2]<br />
<br />
independence_test(x~y)<br />
<br />
plot(x~y)<br />
<br />
summary(lm(x~y))<br />