[1] "Created: Wed Jun 25 09:34:12 2014"

Load in the new version of the original data. Load code from previous work.

load("oldtranrev.rda")
source("funcs.R")
GPR <- TRUE

Set the kernel type to squared exponential and set the inverse width:

kerntype <- "exponential"
wvec <- 2e-4

Use the same class labels as previously.

levels(lcdbi$type) <- c("agn","blazar","cv","flare","sn")

Investigate the range of observation. Compute the date of the first observation:

fd <- sapply(lcdb, function(x) min(x$mjd))
plot(fd, ylab="MJD")