EDISON: Network Reconstruction and Changepoint Detection

Package EDISON (Estimation of Directed Interactions from Sequences Of Non-homogeneous gene expression) runs an MCMC simulation to reconstruct networks from time series data, using a non-homogeneous, time-varying dynamic Bayesian network. Networks segments and changepoints are inferred concurrently, and information sharing priors provide a reduction of the inference uncertainty.

Install the latest version of this package by entering the following in R:
AuthorFrank Dondelinger, Sophie Lebre
Date of publication2016-03-30 21:04:12
MaintainerFrank Dondelinger <fdondelinger.work@gmail.com>

View on CRAN

Man pages

AcceptableMove: Check if move is acceptable.

addProposalNetworkInfo: Add the proposed new network to the new.nets list.

BinoHyperMove: Makes a binomial hyperparameter move.

BinoHyperRatio: Calculates the MH ratio of the binomial prior.

bp.computeAlpha: Computes the acceptance ratio of two changepoint...

buildXY: Builds response Y and predictor X.

CalculateChanges: Function to calculate the number of differences between...

calculateCPPGlobal: Calculated the global changepoint probabilities.

calculateCPProbabilities: Calculate the changepoint probabilities.

calculateEdgeProbabilities: Calculate the edge probabilities.

calculateEdgeProbabilitiesSegs: Calculate edge probabilities for fixed segments.

calculateEdgeProbabilitiesTimePoints: Calculate the edge posterior probabilities for each...

CalculateLikelihoodRatio: Calculates the ratio of two likelihoods in a structure move.

CalculatePriorRatio: Calculates the network prior ratio.

CollectNetworkInfo: Collects all the network information in one list.

computePx: Compute projection matrix.

computeRho4: Calculate proposal frequencies for changepoint moves.

convert_nets: Convert internal representation of networks.

cp.birth: Make changepoint birth move.

cp.death: Make changepoint death move.

cp.shift: Makes a changepoint shift move.

defaultOptions: Set the default options for the MCMC simulation.

dinvgamma: Calculate inverse gamma distribution.

EDISON-package: Allows for network reconstruction and changepoint detection.

EDISON.run: Wrapper function for starting an MCMC simulation

ExpHyperMove: Makes an exponential hyperparameter move.

ExpHyperRatioTarget: Calculates the ratio of an exponential hyperparameter move.

fix_eigenvalues: Modify network to ensure stationarity.

generateNetwork: Generate a random network.

HyperparameterMove: Make a hyperparameter move.

HyperParms: Sets up initial values of hyperparameters.

init: Initialise the MCMC simulation.

main: Main function of the MCMC simulation.

make_structure_move: Makes a structure move.

NetworkProbBino: Calculates the prior probability of the network segments...

NetworkProbExp: Calculates the prior probability of the network using the...

NetworkRatioBino: Calculates the ratio of binomial prior probabilites.

NetworkRatioExp: Calculates the ratio of exponential network prior...

output: Collects and saves output.

phase.update: Make a network structure or hyperparameter move.

PriorRatioPoisson: Calculate network prior ratio with Poisson prior.

proposalTuning: Tune the proposal width for betas.

proposeContinuous: Propose a new real hyperparameter value.

ProposeDiscrete: Propose a new discrete value.

psrf: Calculates the potential scale reduction factor.

psrf_check: Check the potential scale reduction factors for all...

psrf_check_hyper: Checks the potential scale reduction factor for the...

readDataTS: Read target data.

rinvgamma: Samples from the inverse gamma distribution.

runDBN: Setup and run the MCMC simulation.

sampleBinit: Sample initial regression coefficients.

sampleBxy: Sample regression coefficients.

sampleDelta2: Sample delta squared.

sampleK: Sample initial number of changepoints.

sampleParms: Sample initial parameters for the MCMC simulation.

sampleSig2: Sample initial sigma squared.

simulateNetwork: Generate network and simulate data.

updateSigMulti: Update sigma squared variances.

updateSigSolo: Sample new values for sigma squared.


AcceptableMove Man page
addProposalNetworkInfo Man page
BinoHyperMove Man page
BinoHyperRatio Man page
bp.computeAlpha Man page
buildXY Man page
CalculateChanges Man page
calculateCPPGlobal Man page
calculateCPProbabilities Man page
calculateEdgeProbabilities Man page
calculateEdgeProbabilitiesSegs Man page
calculateEdgeProbabilitiesTimePoints Man page
CalculateLikelihoodRatio Man page
CalculatePriorRatio Man page
CollectNetworkInfo Man page
computePx Man page
computeRho4 Man page
convert_nets Man page
cp.birth Man page
cp.death Man page
cp.shift Man page
defaultOptions Man page
dinvgamma Man page
EDISON Man page
EDISON-package Man page
EDISON.run Man page
ExpHyperMove Man page
ExpHyperRatioTarget Man page
fix_eigenvalues Man page
generateNetwork Man page
HyperparameterMove Man page
HyperParms Man page
init Man page
main Man page
make_structure_move Man page
NetworkProbBino Man page
NetworkProbExp Man page
NetworkRatioBino Man page
NetworkRatioExp Man page
output Man page
phase.update Man page
PriorRatioPoisson Man page
proposalTuning Man page
proposeContinuous Man page
ProposeDiscrete Man page
psrf Man page
psrf_check Man page
psrf_check_hyper Man page
readDataTS Man page
rinvgamma Man page
runDBN Man page
sampleBinit Man page
sampleBxy Man page
sampleDelta2 Man page
sampleK Man page
sampleParms Man page
sampleSig2 Man page
simulateNetwork Man page
updateSigMulti Man page
updateSigSolo Man page


tests/testthat/test-structure-moves.r tests/testthat/test-fixed-edges.r tests/test-all.R
R/sampleK.R R/calculateCPProbabilities.R R/proposalTuning.R R/EDISON-package.R R/cp.shift.R R/computeRho4.R R/ProposeDiscrete.R R/output.R R/CalculateChanges.R R/psrf_check_hyper.R R/calculateEdgeProbabilities.R R/BinoHyperMove.R R/updateSigMulti.R R/calculateEdgeProbabilitiesSegs.R R/fix_eigenvalues.R R/CollectNetworkInfo.R R/PriorRatioPoisson.R R/psrf_check.R R/HyperparameterMove.R R/sampleBxy.R R/updateNetworkInfo.R R/sampleBinit.R R/bp.computeAlpha.R R/rinvgamma.R R/cp.death.R R/readDataTS.R R/NetworkRatioExp.R R/NetworkProbExp.R R/proposeContinuous.R R/simulateNetwork.R R/EDISON.run.R R/defaultOptions.R R/NetworkProbBino.R R/CalculateLikelihoodRatio.R R/BinoHyperRatio.R R/make_structure_move.R R/init.R R/CalculatePriorRatio.R R/generateNetwork.R R/sampleDelta2.R R/sampleParms.R R/NetworkRatioBino.R R/sampleSig2.R R/ExpHyperMove.R R/dinvgamma.R R/convert_nets.R R/buildXY.R R/phase.update.R R/runDBN.R R/calculateCPPGlobal.R R/updateSigSolo.R R/computePx.R R/calculateEdgeProbabilitiesTimePoints.R R/ExpHyperRatioTarget.R R/psrf.R R/addProposalNetworkInfo.R R/AcceptableMove.R R/main.R R/cp.birth.R R/HyperParms.R
man/psrf_check.Rd man/NetworkProbBino.Rd man/fix_eigenvalues.Rd man/defaultOptions.Rd man/readDataTS.Rd man/proposalTuning.Rd man/ExpHyperMove.Rd man/init.Rd man/NetworkRatioExp.Rd man/convert_nets.Rd man/NetworkProbExp.Rd man/calculateEdgeProbabilitiesTimePoints.Rd man/addProposalNetworkInfo.Rd man/EDISON.run.Rd man/cp.shift.Rd man/calculateEdgeProbabilitiesSegs.Rd man/main.Rd man/sampleK.Rd man/calculateCPProbabilities.Rd man/proposeContinuous.Rd man/calculateEdgeProbabilities.Rd man/NetworkRatioBino.Rd man/PriorRatioPoisson.Rd man/CollectNetworkInfo.Rd man/ExpHyperRatioTarget.Rd man/psrf_check_hyper.Rd man/cp.birth.Rd man/updateSigMulti.Rd man/BinoHyperRatio.Rd man/computeRho4.Rd man/EDISON-package.Rd man/buildXY.Rd man/CalculateChanges.Rd man/output.Rd man/sampleSig2.Rd man/runDBN.Rd man/AcceptableMove.Rd man/HyperparameterMove.Rd man/HyperParms.Rd man/phase.update.Rd man/BinoHyperMove.Rd man/CalculatePriorRatio.Rd man/psrf.Rd man/generateNetwork.Rd man/sampleBinit.Rd man/ProposeDiscrete.Rd man/CalculateLikelihoodRatio.Rd man/calculateCPPGlobal.Rd man/rinvgamma.Rd man/sampleBxy.Rd man/cp.death.Rd man/computePx.Rd man/make_structure_move.Rd man/sampleParms.Rd man/bp.computeAlpha.Rd man/dinvgamma.Rd man/sampleDelta2.Rd man/updateSigSolo.Rd man/simulateNetwork.Rd

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.