Displaying 1 result from an estimated 1 matches for "behivour".
Did you mean:
behaivour
2001 Apr 05
1
PR#896
...able(dr)
dz <- (sum(dr[seq(along = x)]) - n.x * (n.x + 1)/2 - n.x * n.y/2)
CORRECTION.CI<-0
SIGMA.CI <- sqrt((n.x * n.y/12) * ((n.x + n.y + 1) -
sum(NTIES.CI^3 - NTIES.CI)/((n.x + n.y) * (n.x + n.y - 1))))
dz <- (dz - CORRECTION.CI)/SIGMA.CI
abs(dz - zq)
}
To examine the behivour I plotted the course of the wdiff for a lot of
d and three zq (0.05, 0.5, 0.95) and I let optimize to compute the
minimums. Then I plotted it.
mumin<-min(x)-max(y)
mumax<-max(x)-min(y)
lll<-seq(mumin,mumax,by=0.01)
wdl<-apply(cbind(lll),1,wdiff,zq=qnorm(0.05))
wdm<-apply(cbind(...