library(hdrm)

# bcTCGA
Tab4.1()
Fig4.4()

attachData(bcTCGA)
A <- cor(X)
rho <- A[upper.tri(A)]
hist(rho, breaks=seq(-1,1,by=.05), col='gray', border='white', main='', yaxt='n', ylab='', xlab=expression(rho))
mean(abs(rho)<0.4)

# Scheetz2006
Tab4.2()

attachData(Scheetz2006)
B <- cor(X)
rho <- B[upper.tri(B)]
hist(rho, breaks=seq(-1,1,by=.05), col='gray', border='white', main='', yaxt='n', ylab='', xlab=expression(rho))
mean(abs(rho)<0.4)