function (k=3,nvec=c(20,10,8),fac=1) 
{
nvec=fac*nvec
rho=c(1,2,3)
x=list()
response=NULL
treatment=NULL
for(i in 1:k){
x[[i]]=rnorm(nvec[i],0,rho[i])
response=c(response,x[[i]])
treatment=c(treatment,rep(i,nvec[i]))
}
treatment=as.factor(treatment)

par(mfrow=c(2,1))
boxplot(x)

stripchart(x,vertical=T)
out1=kruskal.test(x)
out2=adk.test(x)
out3=anova(lm(response~treatment))
list(p.value.KW=out1$p.value,p.value.adk=out2$adk[1,2],p.value.F=out3[1,5])
}