(qc <- QcMetric())
qcdata(qc)
try(qcdata(qc, "x"))
x <- rnorm(10)
qcdata(qc, "qc1") <- x
qcdata(qc, "qc2") <- 1:10
qcdata(qc)
all.equal(qcdata(qc, "qc1"), x)
all.equal(qcdata(qc, "qc2"), 1:10)
name(qc) <- "My test QcMetric"
description(qc) <- "This qc metric describes bla bla bla, indicating possible issues in the third step of protocol bla bla bla."
status(qc) <- FALSE
qc
## or
e <- new.env()
e$qc1 <- rnorm(100)
e$qc2 <- 1:100
qcdata(qc) <- e
length(qcdata(qc, "qc1"))
head(qcdata(qc, "qc2"))
show(qc)
show(qc) <- function(object) cat("Updated show method
")
show(qc)
show(qc) <- qcshow
qc
plot(qc)
plot(qc) <-
function(object, ...)
plot(qcdata(object, "qc2"),
qcdata(object, "qc1"),
xlab = "qc1",
ylab = "qc2",
...)
plot(qc)
plot(qc, col = "red", pch = 19)
## generate a report
qcReport(qcm)Run the code above in your browser using DataLab