search for: rabj

...$ouest ansari.test(classe1, classe2, alternative = "greater") # AB statistic : AB1.calc<-sum(subset(obs4[,1],obs4[,3]=="sud")) AB1.calc # Z statistic : n1<-length(classe1) n2<-length(classe2) n<-n1+n2 n # n is even mAB1<-n1*(n+2)/4 tj<- table(obs4$xr) rABj<- obs4$rAB[c(which(diff(obs4$xr)!=0),n)] s2.AB<-(n1*n2*(16*sum(tj*rABj^2)-n*(n+2)^2))/(16*n*(n-1)) Z.calc<-(AB1.calc-mAB1)/ (s2.AB)^0.5 Z.calc # p value : pnorm(Z.calc, lower.tail = FALSE) pnorm(Z.calc) Z is negative. So at the right of the normal law, the p value should be over 0.5 but...