# Load some data
data("recoup_test_data",package="recoup")
# Calculate coverages
test.tss <- recoup(
test.input,
design=NULL,
region="tss",
type="chipseq",
genome=test.genome,
flank=c(2000,2000),
selector=NULL,
plotParams=list(plot=FALSE,profile=TRUE,
heatmap=TRUE,device="x11"),
rc=0.5
)
# Plot coverage profiles
# Get the curve profile plot
pp <- getr(test.tss,"profile")
# Change some ggplot parameter
pp <- pp +
theme(axis.title.x=element_text(size=14))
# Store the new plot
test.tss <- setr(test.tss,"profile",pp)
## or even better
# test.tss <- setr(test.tss,list(profile=pp))
Run the code above in your browser using DataLab