Singed rank plot

signedrank.r

signTable <- function(mu,x){
  x2 <- x-mu
  x2 <- x2[x2 != 0]
  data.frame(table(sign(x2)), loc=mu)
}

wilcoxTable <- function(mu,x){
  x2 <- x-mu
  x2 <- x2[x2 != 0]
  data.frame(w=sign(x2)*rank(abs(x2)), loc=mu)
}

x <- rt(30,5) + 1

muv <- seq(min(x),max(x),by=0.1)
dfw <- bind_rows(lapply(muv, wilcoxTable, x=x))

ggplot(dfw,aes(x=loc, y=w))+
  geom_tile(height=1, width=0.1, colour="lightgrey")+
  geom_vline(xintercept = 1, lty=1, colour="royalblue", alpha=0.7)+
  theme_bw()+
  labs(y="signed-rank", x = "location parameter (null)")

ggsave("wilcox.png")