Calculate PSI from bam files from one sample(with 2 duplicates)

Question

[pxd]

Dear all, I just want to focus on the exon skipping events and I have two bam files for the 2 replicates, and both of them are pair-ended. I want to calculate the PSI value for every exon and so that to give a label for which exon are more likely spliced out, and set a threshold for the psi value for downstream analysis. I checked many softwares both most of them need more than 2 samples so that to analysis the differential patterns for alternative splicing, I am not sure how to do that and I want some methods or softwares to solve the problem.