hahmmr1.0.0 package

Haplotype-Aware Hidden Markov Model for RNA

analyze_allele

Analyze allele profile

analyze_joint

Analyze allele and expression profile

annot_cm

Annotate genetic distance between markers

annot_segs

Annotate copy number segments after HMM decoding

approx_phi_post

Laplace approximation of the posterior of expression fold change phi

approx_theta_post_s2

Find optimal theta of a CNV segment using forward-backward; uses a 2-s...

approx_theta_post_s3

Find optimal theta of a chromosome using forward-backward; uses a 3-st...

calc_allele_lik_s2

Calculate allele likelihoods for 2-state allele HMM

calc_allele_lik_s3

Calculate allele likelihoods for 3-state allele HMM

calc_allele_LLR

Calculate LLR for an allele HMM

calc_exp_LLR

Calculate LLR for an expression HMM

calc_trans_mat_s15

Calculate the transition matrix for 15-state joint HMM

calc_trans_mat_s7

Calculate the transition matrix for 7-state joint HMM

check_cols

Check if columns are present in dataframe

check_matrix

Check the format of a count matrix

classify_alleles

classify alleles using viterbi and forward-backward

combine_bulk

Combine allele and expression pseudobulks

dbbinom

Beta-binomial distribution density function A distribution is beta-bin...

dpoilog

Returns the density for the Poisson lognormal distribution with parame...

filter_genes

filter for mutually expressed genes

fit_gamma

fit gamma maximum likelihood

fit_lnpois_cpp

Fit MLE of log-normal Poisson model

fit_ref_sse

Fit a reference profile from multiple references using constrained lea...

forward_back_allele

Forward-backward algorithm for allele HMM

generate_postfix

Generate alphabetical postfixes

get_allele_bulk

Aggregate into pseudobulk alelle profile

get_allele_hmm_s2

Make a 2-state allele HMM - no transitions to netural state

get_allele_hmm_s3

Make a 3-state allele HMM - allow transitions to netural state

get_bulk

Produce combined bulk expression and allele profile

get_exp_bulk

Aggregate into bulk expression profile

get_trans_probs_s15

Helper function to calculate transition porbabilities for 15-state joi...

get_trans_probs_s7

Helper function to calculate transition porbabilities for 7-state join...

l_bbinom

calculate joint likelihood of allele data

l_lnpois

calculate joint likelihood of a PLN model

likelihood_allele

Only compute total log likelihood from an allele HMM

logSumExp

logSumExp function

plot_bulks

Plot a group of pseudobulk HMM profiles

plot_psbulk

Plot a pseudobulk HMM profile

run_allele_hmm_s3

Run a 3-state allele HMM

run_allele_hmm_s5

Run a 5-state allele-only HMM - two theta levels

run_joint_hmm_s15

Run 15-state joint HMM on a pseudobulk profile

run_joint_hmm_s7

Run 7-state joint HMM on a pseudobulk profile

smooth_segs

Smooth the segments after HMM decoding

switch_prob_cm

predict phase switch probablity as a function of genetic distance

theta_hat_seg

Estimate of imbalance level theta in a segment

viterbi_allele

Viterbi algorithm for allele HMM

viterbi_joint

Viterbi algorithm for joint HMM, can only handle one set of observatio...

viterbi_joint_mat

Generalized viterbi algorithm for joint HMM, can handle multiple sets ...

Haplotype-aware Hidden Markov Model for RNA (HaHMMR) is a method for detecting copy number variations (CNVs) from bulk RNA-seq data. Additional examples, documentations, and details on the method are available at <https://github.com/kharchenkolab/hahmmr/>.

  • Maintainer: Teng Gao
  • License: MIT + file LICENSE
  • Last published: 2023-10-25