data(usz_13c)
if (FALSE) {
str(usz_13c)
# Plot all records; this needs some time
pdf(file.path(tempdir(), "usz_13c.pdf"), height= 30)
# null_fit makes data plotable without fitting a model
plot(null_fit(usz_13c))
dev.off()
}
# Plot a subset
suppressPackageStartupMessages(library(dplyr))
usz_part = usz_13c %>%
filter(patient_id %in% c("norm_001","norm_002", "pat_001", "pat_002"))
plot(null_fit(usz_part))
Run the code above in your browser using DataLab