## Load the input object from Marioni's data as returned by readData()
data(myCounts)
## Computing differential expression probability on RPKM-normalized data by NOISeq-real using factor "Tissue"
mynoiseq = noiseq(mydata, k = 0.5, norm = "rpkm", replicates = "technical", factor="Tissue",
pnr = 0.2, nss = 5, v = 0.02, lc = 1)
## Computing differential expression probability on Upper Quartile normalized data by NOISeq-real
## using factor "TissueRun" and comparing samples in Run 1 (levels "Kidney_1" and "Liver_1")
mynoiseq.uqua = noiseq(mydata, k = 0.5, norm = "uqua", replicates = "technical", factor="TissueRun",
conditions = c("Kidney_1","Liver_1"), pnr = 0.2, nss = 5, v = 0.02, lc = 1)
Run the code above in your browser using DataLab