rna_majiq.Heterogen

class rna_majiq.Heterogen(psi1, psi2, min_experiments_f=0.5, name1='grp1', name2='grp2', show_progress=False)

Compare Psi between two groups of PsiCoverage (independence assumption)

Compare Psi between two groups of PsiCoverage under the assumption that each experiment is independent. Perform null-hypothesis testing on posterior means and/or samples under assumption that underlying values of PSI per experiment are independent and identically distributed (i.i.d.).

Parameters:

psi1, psi2 (Union[PsiCoverage, PsiGroup]) – PsiCoverage or PsiGroup for the two groups to compare
min_experiments_f (float) – Number or proportion of experiments required to pass in each group for quantification
name1, name2 (str) – Names to indicate group identity
show_progress (bool) – Show progress of intermediate computations required for setting up this object.

__init__(psi1, psi2, min_experiments_f=0.5, name1='grp1', name2='grp2', show_progress=False)

Methods

`__init__`(psi1, psi2[, min_experiments_f, ...])
`approximate_stats`([ec_idx, quantiles, ...])	Statistics on means, samples from approximate posteriors
`dataset`([pvalue_quantiles, use_stats, ...])
`events_to_zarr`(path, mode[, consolidated])	Save events information to specified path
`get_events`(introns, junctions)	Construct `Events` using saved dataset and introns, junctions
`raw_stats`([ec_idx, use_stats])	Statistics on means, samples from raw posteriors

Attributes

`EVENTS_EXPECTED_VARIABLES`
`event_size`	array(ec_idx) total number of connections from same event
`labels`	labels associated with prefixes in self.psi_concatenated
`lsv_idx`	array(ec_idx) index identifying event it belongs to
`lsv_offsets`	array(e_offsets_idx) offsets for events into ec_idx
`num_connections`	Total number of connections over all events
`num_events`	Total number of events
`psi_concatenated`	Concatenated values from self.psi1 and self.psi2 for statistics