library(fabia)
setwd("c:/sepp/work/fabia/experiments/benchmark/data_R")
for (i in 1:10) {
file <- paste("exp_",as.character(i),".RData",sep="")
load(file)
n <- nrow(dat$X)
l <- ncol(dat$X)
p <- 13
lll <- 12
yl <- c(1,round(n*(1:lll)/lll))
yll <- rep("",n)
yll[yl] <- as.character(yl)
x11()
x11()
myImagePlot(dat$X,title=paste("Experiment ",as.character(i),": Noisy Data",sep=""),yLabels= yll)
dev.print(file=paste("exp_",as.character(i),"_noisy.pdf",sep=""),device=pdf)
x11()
myImagePlot(dat$Y,title=paste("Experiment ",as.character(i),": Noise Free Data",sep=""),yLabels= yll)
dev.print(file=paste("exp_",as.character(i),"_noise_free.pdf",sep=""),device=pdf)
}