Description Usage Arguments Details Value Author(s) See Also Examples
Get minimum free energy and ensemble free energy for specified variants of RNA sequences and their differences from wild type (calculated by RNAfold); vectorized.
1 2 |
wt_sequences |
Character vector representing the wild type RNA sequences to (DNA also ok... will substitute U for T) |
variants |
Optional logical. Indicates whether variants are provided that should be analyzed (TRUE) or whether the wild type sequence should be analyzed (FALSE). Default is FALSE. |
start |
Optional integer vector same length as wt_sequences. Specifies the starting nucleotide positions of the variants. Must be provided if variants is TRUE. |
stop |
Optional integer vector same length as wt_sequences. Specifies the ending nucleotide positions of the variants. Must be provided if variants is TRUE. |
wt |
Optional character vector same length as wt_sequences. Not necessary to provide even if variant = TRUE. For error-checking purposes only. |
mut |
Optional character vector same length as wt_sequences. Specifies what the sequence between start and stop nucleotides should be changed to. Must be provided if variants is TRUE. |
MFEwt |
Optional numeric vector, the minimum free energies of the wild type sequences. This parameter is provided rather than calculated to avoid repetitive MFE calculations for vectors where all of the wild type sequences are identical. |
ensFEwt |
Optional numeric vector, the ensemble free energies of the wild type sequences. This parameter is provided rather than calculated to avoid repetitive MFE calculations for vectors where all of the wild type sequences are identical. |
A more detailed description of what the function does.
Data frame with 4 columns, MFE (minimum free energy), MFE_diff (minimum free energy difference, mutant minus wild type), ensFE (ensemble free energy), ensFE_diff (ensemble free energy difference, mutant minus wild type)
Rachel Silverstein (aut)
RNAfold
Called by this function
mutSeqs
Called by this function
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | # first get the wild type free energes using varfold without a variant
seq1 <- "AAAAAAGUUUGGGGGGCCCCCCTTTTTT"
seq2 <- "UUUUUUUUAAAAAAAUUUUUUUUAAAAAAAAUUUUUUUUAAAAAAAA"
wt_data <- varFold(wt_sequences = c(seq1, seq2), variants = FALSE)
MFEwt <- wt_data$MFE
ensFEwt <- wt_data$ensFE
# > wt_data
# MFE MFE_diff ensFE ensFE_diff
# 1 -7.4 No WT MFE data provided. -8.35 No WT ensFE data provided.
# 2 -10.0 No WT MFE data provided. -11.54 No WT ensFE data provided.
# now get the variant information with differences using the wild info just calculated
var_data <- varFold(wt_sequences = c(seq1, seq2), variants = TRUE, start = c(1, 2), stop = c(2, 3), wt = c("AA", "UU"), mut = c("GG", "CC"), MFEwt = MFEwt, ensFEwt = ensFEwt)
# > var_data
# MFE MFE_diff ensFE ensFE_diff
# MFE MFE_diff ensFE ensFE_diff
# 1 -7.6 -0.2 -8.49 -0.14
# 2 -8.2 1.8 -9.42 2.12
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.