plotPredReal <- function(data, pred, index)
{
prediction = pred$getForecast(index)
plotPredReal <- function(data, pred, index)
{
prediction = pred$getForecast(index)
yrange = range(measure, prediction)
par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5, lwd=3)
plot(measure, type="l", ylim=yrange, xlab="Time (hours)", ylab="PM10")
yrange = range(measure, prediction)
par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5, lwd=3)
plot(measure, type="l", ylim=yrange, xlab="Time (hours)", ylab="PM10")
par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5)
rainbow::fboxplot(series_fds, "functional", "hdr", xlab="Time (hours)", ylab="PM10",
plotlegend=FALSE, lwd=2)
par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5)
rainbow::fboxplot(series_fds, "functional", "hdr", xlab="Time (hours)", ylab="PM10",
plotlegend=FALSE, lwd=2)
- ref_var = c( apply(data$getSeries(fil$neighb_indices),1,sd),
- apply(data$getSeries(fil$neighb_indices+1),1,sd) )
+ ref_var = c( apply(data$getSeries(fil$neighb_indices-1),1,sd),
+ apply(data$getSeries(fil$neighb_indices),1,sd) )
tdays = .getNoNA2(data, 2, fil$index)
global_var = c(
apply(data$getSeries(tdays-1),1,sd),
tdays = .getNoNA2(data, 2, fil$index)
global_var = c(
apply(data$getSeries(tdays-1),1,sd),