I have results from using the DEXSeqDataSetFromHTSeq() function in DEXSeq. I want to calculate exon PSI (given each transcript has its own groupID and each exon has its own featureID), can I just convert the exon usage coefficients to counts by taking 2^x then do (ExonCount/SumOfAllExonCounts)?
I am aware that most approaches just use the raw counts, but I am avoiding it because that doesn't seem to take advantage of DEXSeq's statistical model, I think?
Also, would it be more accurate to calculate a delta PSI value using the log2fold value? Is there any easy/quick way I can achieve this?