Nothing
###############################################################################
# run the simulations for boosting of functional data
# code based on code by Fabian Scheipl
# author: Sarah Brockhaus
###############################################################################
##### test gain in time for computation as array model compared to normal computation
print(R.Version()$version.string)
#library(refund)
library(FDboost)
library(splines)
pathResults <- NULL
library(pryr)
source("simfuns.R")
##################################### M=50
set.seed(18102012)
settingsTime <- makeSettings(
M=c(50, 100, 500),
ni=c(1),
Gy=c(30, 100),
Gx=c(30),
snrEps=c(2),
snrE=c(0),
snrB=c(2),
scenario=3,
balanced=c(TRUE),
nuisance=c(0),
rep=1:10) # 1:10
length(settingsTime)
usecores <- 10
options(cores=usecores)
res <- try(doSimFDboost(settings=settingsTime, oneRepTime))
save(res, file=paste(pathResults, "resTime.Rdata", sep=""))
res$M <- factor(res$M)
res$G <- factor(res$Gy)
res$mod <- factor(res$model-2, levels=c(0,1), labels = c("FLAM", "long"))
#### generate boxplots of errors
library(ggplot2)
pdf("arrayLongTime.pdf", width=9, height=7)
p <- ggplot(res, aes(y=time.elapsed, x=M, fill=mod))
p + geom_boxplot(aes(colour = mod), outlier.size=.6) +
facet_grid( ~ G, labeller="label_both") +
scale_colour_manual('mod', values = c('FLAM' = 'black', 'long' = 'grey30')) +
scale_fill_manual(values=c("grey80", "grey95")) +
scale_y_continuous(breaks=c(1, 2, 5, 10, 20, 60, 120, 300, 600, 1200, 2700,
5400, 10800, 21600, 43200),
trans="log10",
labels=c("1s", "2s", "5s", "10s", "20s", "1 min", "2 min", "5 min",
"10 min", "20 min", "45 min", "90 min", "3h", "6h", "12h")) +
labs(title="", x="N", y="time") +
theme(legend.position="top", legend.direction="horizontal",
text=element_text(size = 25), legend.title=element_blank()) +
theme(text=element_text(size = 20)) #+ theme_bw()
dev.off()
pdf("arrayLongMemory.pdf", width=9, height=7)
p <- ggplot(res, aes(y=memory, x=M, fill=mod))
p + geom_boxplot(aes(colour = mod), outlier.size=.6) +
facet_grid( ~ G, labeller="label_both") +
scale_colour_manual('mod', values = c('FLAM' = 'black', 'long' = 'grey30')) +
scale_fill_manual(values=c("grey80", "grey95")) +
scale_y_continuous(trans="log10") +
labs(title="", x="N", y="memory") +
theme(legend.position="top", legend.direction="horizontal",
text=element_text(size = 25), legend.title=element_blank()) +
theme(text=element_text(size = 20)) #+ theme_bw()
dev.off()
######################################
print(sessionInfo())
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.