RNAlib-2.4.4
|
|
Partition function for two RNA sequences. More...
Go to the source code of this file.
Functions | |
vrna_dimer_pf_t | vrna_pf_co_fold (const char *seq, char *structure, vrna_ep_t **pl) |
Calculate partition function and base pair probabilities of nucleic acid/nucleic acid dimers. More... | |
vrna_dimer_pf_t | co_pf_fold (char *sequence, char *structure) |
Calculate partition function and base pair probabilities. More... | |
vrna_dimer_pf_t | co_pf_fold_par (char *sequence, char *structure, vrna_exp_param_t *parameters, int calculate_bppm, int is_constrained) |
Calculate partition function and base pair probabilities. More... | |
vrna_ep_t * | get_plist (vrna_ep_t *pl, int length, double cut_off) |
void | compute_probabilities (double FAB, double FEA, double FEB, vrna_ep_t *prAB, vrna_ep_t *prA, vrna_ep_t *prB, int Alength) |
Compute Boltzmann probabilities of dimerization without homodimers. More... | |
void | init_co_pf_fold (int length) |
FLT_OR_DBL * | export_co_bppm (void) |
Get a pointer to the base pair probability array. More... | |
void | free_co_pf_arrays (void) |
Free the memory occupied by co_pf_fold() More... | |
void | update_co_pf_params (int length) |
Recalculate energy parameters. More... | |
void | update_co_pf_params_par (int length, vrna_exp_param_t *parameters) |
Recalculate energy parameters. More... | |
Variables | |
int | mirnatog |
Toggles no intrabp in 2nd mol. | |
double | F_monomer [2] |
Free energies of the two monomers. | |
Partition function for two RNA sequences.
vrna_dimer_pf_t co_pf_fold | ( | char * | sequence, |
char * | structure | ||
) |
Calculate partition function and base pair probabilities.
This is the cofold partition function folding. The second molecule starts at the cut_point nucleotide.
sequence | Concatenated RNA sequences |
structure | Will hold the structure or constraints |
vrna_dimer_pf_t co_pf_fold_par | ( | char * | sequence, |
char * | structure, | ||
vrna_exp_param_t * | parameters, | ||
int | calculate_bppm, | ||
int | is_constrained | ||
) |
Calculate partition function and base pair probabilities.
This is the cofold partition function folding. The second molecule starts at the cut_point nucleotide.
sequence | Concatenated RNA sequences |
structure | Pointer to the structure constraint |
parameters | Data structure containing the precalculated Boltzmann factors |
calculate_bppm | Switch to turn Base pair probability calculations on/off (0==off) |
is_constrained | Switch to indicate that a structure contraint is passed via the structure argument (0==off) |
DO NOT USE THIS FUNCTION ANYMORE
void compute_probabilities | ( | double | FAB, |
double | FEA, | ||
double | FEB, | ||
vrna_ep_t * | prAB, | ||
vrna_ep_t * | prA, | ||
vrna_ep_t * | prB, | ||
int | Alength | ||
) |
Compute Boltzmann probabilities of dimerization without homodimers.
Given the pair probabilities and free energies (in the null model) for a dimer AB and the two constituent monomers A and B, compute the conditional pair probabilities given that a dimer AB actually forms. Null model pair probabilities are given as a list as produced by assign_plist_from_pr(), the dimer probabilities 'prAB' are modified in place.
FAB | free energy of dimer AB |
FEA | free energy of monomer A |
FEB | free energy of monomer B |
prAB | pair probabilities for dimer |
prA | pair probabilities monomer |
prB | pair probabilities monomer |
Alength | Length of molecule A |
void init_co_pf_fold | ( | int | length | ) |
DO NOT USE THIS FUNCTION ANYMORE
FLT_OR_DBL* export_co_bppm | ( | void | ) |
Get a pointer to the base pair probability array.
Accessing the base pair probabilities for a pair (i,j) is achieved by
FLT_OR_DBL *pr = export_bppm(); pr_ij = pr[iindx[i]-j];
void free_co_pf_arrays | ( | void | ) |
Free the memory occupied by co_pf_fold()
void update_co_pf_params | ( | int | length | ) |
Recalculate energy parameters.
This function recalculates all energy parameters given the current model settings.
length | Length of the current RNA sequence |
void update_co_pf_params_par | ( | int | length, |
vrna_exp_param_t * | parameters | ||
) |
Recalculate energy parameters.
This function recalculates all energy parameters given the current model settings. It's second argument can either be NULL or a data structure containing the precomputed Boltzmann factors. In the first scenario, the necessary data structure will be created automatically according to the current global model settings, i.e. this mode might not be threadsafe. However, if the provided data structure is not NULL, threadsafety for the model parameters dangles, pf_scale and temperature is regained, since their values are taken from this data structure during subsequent calculations.
length | Length of the current RNA sequence |
parameters | data structure containing the precomputed Boltzmann factors |