sampleTREuncertainty.Rd
Function to calculate the uncertainty in the TRE values by sampling the posterior distribution of the Rex parameters. Typically used to obtain uncertainty in the conformational signatures.
sampleTREuncertainty(
HdxData,
RexParamsList,
quantity = "TRE",
states,
whichChain = 1,
whichSamples = seq.int(1, 50),
whichTimepoint = 600,
num_montecarlo = 5000
)
The HdxData object
A list of RexParams objects
The quantity to use for the analysis. Default is "TRE"
The state name to use for the analysis. e.g. ligand used in differential analysis
The chain to use for the analysis. Default is 1
The samples to use for the analysis. Default is seq.int(1, 50)
The timepoint to use for the analysis. Default is 600.
The number of montecarlo samples to use for the analysis
A data frame containing the TRE values and the residue number over monte carlo interations
if (FALSE) {
library("RexMS")
data("LXRalpha_processed")
data("rex_lxr")
data("LXRalpha_compounds")
states <- names(LXRalpha_compounds)
TRE_dist <- sampleTREuncertainty(HdxData = LXRalpha_processed,
RexParamsList = rex_lxr,
quantity = "TRE",
states = states,
whichChain = 1,
whichSamples = seq.int(1, 50),
whichTimepoint = 600,
num_montecarlo = 5000)
}