R script for edgeR allelic imbalance scores