anovabinom <- function(nj,pj,alpha=0.05)
{
  a <- length(nj)
  mYj <- pj/nj  # prumery ve skupinach
  mY <- sum(pj)/sum(nj)   # celkovy prumer
  Q <- 1/(mY*(1-mY))*sum(nj*(mYj-mY)^2)
  kvantil <- qchisq(1-alpha,a-1)
  
  if (Q<kvantil)
  {
    writeLines("We can not reject H0, samples have the same mean.")
    anovabinom <- NULL
  }
  else
  {   writeLines("We reject H0, samples have not the same mean.")
      cat("p-value: ")
      pval<-1-pchisq(Q,df=a-1)
      print(pval,digits=4,quote=F)
      cat("\n")
      # mnohonasobne porovnavani
      Aj <- asin(sqrt(mYj))
      pom1<- matrix(Aj,a,1)%*%rep(1,times=4)
      pom2<- rep(1,times=4)%*%matrix(Aj,1,a)
      ML <- abs(pom1-pom2)   # matice levych stran
      pom1<- matrix(1/nj,a,1)%*%rep(1,times=4)
      pom2<- rep(1,times=4)%*%matrix(1/nj,1,a)
      MR <- sqrt((pom1+pom2)/8)*qtukey(1-alpha,a,Inf)   # matice pravych stran
      detect <- ML>MR
      detect <- detect*upper.tri(detect)
      anovabinom <- which(detect==1,arr.ind=T)
      n <- dim(anovabinom)[1]
      for (k in 1:n)
      {
        cat(c("different groups: ",as.character(anovabinom[k,]),"\n"))
      }
  }
  return(anovabinom)
}