# efg, Stowers Institute, 25 April 2006

#source("0-LoadData.R")

ShowZeroStats <- function()
{
  pdf("ZeroStats.pdf")
  par(oma=c(1,0,1,0))
  zero <- d == 0

  # 13.15555 %0s (177 arrays) or 13.19176 %0s (172 arrays)

  cat(100.0 * sum(zero) / prod(dim(d)), '%0s\n')

  hist( apply(zero, By.Row<-1, sum), breaks=50,  main="Zeros By Patient",
    xlab=paste("Number of zero intensities for", ncol(d), "probes"))
  hist( apply(zero, By.COL<-2, sum), breaks=100, main="Zeros By Probe",
    xlab=paste("Number of zero intensties for", nrow(d), "patients"))
  PlotFooter("1-Zeros.R: ShowZeroStats")

  dev.off()
}

sink("2-Zeros.txt")
ShowZeroStats()
sink()