Code:
>A <- readcufflinks(dir=system.file(“data”, package=”cummerbund”)) >A CuffSet instance with: 6 samples 23284 genes 30073 isoforms 25872 TSS 24748 CDS 348555 promoters 388080 splicing 305730 relCDS > sigMatrix(A, level=’genes’, alpha =0.05) >length(getSig(A, alpha=0.05, x= ‘sampleX’, y= ‘sampleY’)) [1] 140