WilcoxonSignedRankTest.R

#Fonksiyonu yaratalim
WilcoxonSignedRankTest <- function(x,y){
  diff <- c(x - y)  #her gozlem ikilisi arasindaki farki tutan bir vektor hesaplayalim
  diff.rank <- rank(abs(diff)) #farklarin mutlak degerini alarak siralayalim ve siralarini kaydedelim.
  diff.rank.sign <- diff.rank * sign(diff) #her sira numarasini, karsilik gelen isaret ile etiketleyelim.
  ranks.pos <- sum(diff.rank.sign[diff.rank.sign > 0]) #pozitif degerli siralamalarin sira degerlerini toplayalim.
  ranks.neg <- -sum(diff.rank.sign[diff.rank.sign < 0]) #negatif degerli siralamalarin sira degerlerini toplayalim.
  result<- c(ranks.pos,ranks.neg)
  print(result)
  }

Agustos <- c(18.3,13.3,16.5,12.6,9.5,13.6,8.1,8.9,10,8.3,7.9,8.1,13.4)
Kasim <-c(12.7,11.1,15.3,12.7,10.5,15.6,11.2,14.2,16.3,15.5,19.9,20.4,36.8)

V<-WilcoxonSignedRankTest(Agustos, Kasim)
#[1] 16 75
#Bu test ile hesaplanan aralik (16, 75)'dir.

c(min(c(Agustos-Kasim)[order(c(Agustos-Kasim))]),max(c(Agustos-Kasim)[order(c(Agustos-Kasim))]))
#[1] -23.4   5.6
#Ikililer arasindaki farklara ait tablolanmis aralik ise (-23.4, 5.6)'dir.

#Tablolanmis aralik, hesaplanmis araligi kapsamadigi icin H0 bos hipotezi reddedilir.
#Yani Agustos ve Kasim aylari arasindaki medyan degisimi 0'dan farklidir.

#######################################################################
#VEYA R'in yerlesik fonksiyonu (built-in function) kullanilabilir.
wilcox.test(Agustos, Kasim, paired=TRUE)

#        Wilcoxon signed rank test
#data:  Agustos and Kasim
#V = 16, p-value = 0.03979
#alternative hypothesis: true location shift is not equal to 0

#p-degeri 0.05'ten kucuk oldugu icin H0 bos hipotezi reddedilir.
#Yani Agustos ve Kasim aylari arasindaki medyan degisimi 0'dan farklidir.