# NOT RUN {
infile <- system.file("extdata", "example1_hst.hst", package = "scbursts")
open = hst.read(infile, extract="open")
closed = hst.read(infile, extract="closed")
header = hst.extract_header(infile)
### Do stuff
hst.write(open, closed, file=file.path(tempdir(), "new_histogram.hst"), header=header)
# }
Run the code above in your browser using DataLab