data(damage)
# Step 1
# Get the probabilities of differential expression
# for the knockout of ATM in the healthy cells
probs.healthy.ATM= differential.probs(data.healthy[,"ATM",FALSE], NULL)
# Get the probabilities of differential expression
# for the knockout of ATM in the damaged cells
probs.damage.ATM= differential.probs(data.damage[,"ATM",FALSE], NULL)
# Step 2
# Regulation scores for a dataset with only one regulator
# equal the signed probabilities
# Step 3
# Get the deregulation scores
deregulation.ATM= deregulation.scores(probs.healthy.ATM, probs.damage.ATM, TRUE)
## Not run:
# # Step 1
# probs.healthy= differential.probs(data.healthy, beliefs.healthy)
# probs.damage= differential.probs(data.damage, beliefs.damage)
#
# # Step 2
# regulation.healthy= regulation.scores(probs.healthy, model.healthy)
# regulation.damage= regulation.scores(probs.damage, model.damage)
#
# # Step 3
# deregulation= deregulation.scores(regulation.healthy, regulation.damage, TRUE)
# ## End(Not run)
Run the code above in your browser using DataLab