wsci <- function(mat,meth=0,interval=95,...)

{n <- nrow(mat);  k <- ncol(mat)

 cilimits <- c((1-interval/100)/2,.5+interval/200)

 forgraph <- {};  adjvalue <- {}

 for (i in 1:n)

  {adjvalue <- c(adjvalue, mean(mat) - mean(mat[i,]))}

 for (j in 1:k)

  {cols <- {}

   for (i in 1:n)

    {cols <- c(cols, mat[i,j] + adjvalue[i])}

   if (meth==0)

    {x <- c(mean(cols),mean(cols)+qt(cilimits,n-1)*stdev(cols)/sqrt(n))}

   if (meth==1)

    {b <- bootstrap(cols,mean)

     x <- c(mean(cols),limits.bca(b, probs=cilimits))}

   forgraph <- cbind(forgraph,x)

  }

dimnames(forgraph) <- list(c("mean", "lower", "upper"),names(mat))

y <- cbind(names(mat),as.data.frame(t(forgraph)))

guiPlot("High Low", DataSetValues = data.frame(y))

guiModify("YAxisTitle", Name=guiGetAxisTitleName(), Title=paste(interval,"% WSCI"))

guiModify("XAxisTitle", Name=guiGetAxisTitleName(), Title="Repeated variables")

return(forgraph)

}