X-Git-Url: https://git.auder.net/%24%7BgetWhatsApp%28link%29%7D?a=blobdiff_plain;f=pkg%2FR%2Fplot.R;h=372abf0a99e086cceb9fec1dc5f9da9c5907da90;hb=360a5d41bee019d4c0e39f0f9b26be4f35b243ad;hp=95dc0f77abf6aba1a925ae409a99152b47a31d6f;hpb=98e958cab563866f8e00886b54336018a2e8bc97;p=talweg.git diff --git a/pkg/R/plot.R b/pkg/R/plot.R index 95dc0f7..372abf0 100644 --- a/pkg/R/plot.R +++ b/pkg/R/plot.R @@ -9,7 +9,7 @@ plotCurves <- function(data, indices=seq_len(data$getSize())) { series = data$getSeries(indices) - yrange = quantile(series, probs=c(0.05,0.95), na.rm=TRUE) + yrange = quantile(series, probs=c(0.025,0.975), na.rm=TRUE) par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5) for (i in seq_along(indices)) { @@ -182,7 +182,7 @@ computeFilaments <- function(data, index, limit=60, plot=TRUE) # Complete series with (past and present) tomorrows ref_serie = c(ref_serie,data$getCenteredSerie(index+1)) centered_series = rbind( centered_series, data$getCenteredSeries(fdays+1) ) - yrange = quantile(cbind(ref_serie,centered_series), probs=c(0.05,0.95), na.rm=TRUE) + yrange = quantile(cbind(ref_serie,centered_series), probs=c(0.025,0.975), na.rm=TRUE) par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5, lwd=2) for (i in nn:1) { @@ -192,7 +192,7 @@ computeFilaments <- function(data, index, limit=60, plot=TRUE) } # Also plot ref curve, in red plot(ref_serie, ylim=yrange, type="l", col="#FF0000", xlab="", ylab="") - abline(v=24, lty=2, col=colors()[56]) + abline(v=24, lty=2, col=colors()[56], lwd=1) } list("index"=index,"neighb_indices"=fdays[sorted_distances$ix[1:nn]],"colors"=colors) @@ -222,8 +222,9 @@ plotFilamentsBox = function(data, fil, ...) usr <- par("usr") yr <- (usr[4] - usr[3]) / 27 par(new=TRUE) - plot(data$getSerie(fil$index), type="l", lwd=2, lty=2, + plot(c(data$getSerie(fil$index),data$getSerie(fil$index+1)), type="l", lwd=2, lty=2, ylim=c(usr[3] + yr, usr[4] - yr), xlab="", ylab="") + abline(v=24, lty=2, col=colors()[56]) } #' Plot relative conditional variability / absolute variability