Hi @yollct ,
In the transcript-level workflow, I am not getting the same graph as the one provided in the tutorial.
Although the ascov table has a different order from yours, it seems to show the same results:
ascov[ascov["gene_symb"]=="CDK4"]
gene gene_symb major_transcript minor_transcript switch_prob corr diff event_importance exclusive_domains p_value adj_pval
71 1019 CDK4 ENST00000552862.1 ENST00000257904.11 1.000000 0.822139 0.066160 0.343667 [] 0.00298 0.002838
438 1019 CDK4 ENST00000312990.10 ENST00000257904.11 0.888889 0.990277 0.061661 0.321748 [] 0.00086 0.000851
spy.switch_plot("CDK4", dset, ascov)
spy.switch_plot("CDK4", dset, ascov, all_isoforms=True)
(Completely identical as the previous one)
the following graphs are identical to yours
Also, the read_gtf function is not defined.
Thanks for your answers.
Best,
Alex
Hi @yollct ,
In the transcript-level workflow, I am not getting the same graph as the one provided in the tutorial.
Although the ascov table has a different order from yours, it seems to show the same results:
ascov[ascov["gene_symb"]=="CDK4"]
gene gene_symb major_transcript minor_transcript switch_prob corr diff event_importance exclusive_domains p_value adj_pval
71 1019 CDK4 ENST00000552862.1 ENST00000257904.11 1.000000 0.822139 0.066160 0.343667 [] 0.00298 0.002838
438 1019 CDK4 ENST00000312990.10 ENST00000257904.11 0.888889 0.990277 0.061661 0.321748 [] 0.00086 0.000851
spy.switch_plot("CDK4", dset, ascov)
spy.switch_plot("CDK4", dset, ascov, all_isoforms=True)
(Completely identical as the previous one)
the following graphs are identical to yours
Also, the read_gtf function is not defined.
Thanks for your answers.
Best,
Alex