projects
/
talweg.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
'update'
[talweg.git]
/
pkg
/
R
/
plot.R
diff --git
a/pkg/R/plot.R
b/pkg/R/plot.R
index
8675c62
..
b5a8e4b
100644
(file)
--- a/
pkg/R/plot.R
+++ b/
pkg/R/plot.R
@@
-30,24
+30,24
@@
plotError <- function(err, cols=seq_along(err))
{
if (!is.null(err$abs))
err = list(err)
{
if (!is.null(err$abs))
err = list(err)
- par(mfrow=c(2,2), mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5
, lty=1
)
+ par(mfrow=c(2,2), mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5)
L = length(err)
yrange = range( sapply(1:L, function(i) ( err[[i]]$abs$day ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$abs$day ), type="l",
L = length(err)
yrange = range( sapply(1:L, function(i) ( err[[i]]$abs$day ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$abs$day ), type="l",
- xlab="Time (hours)", ylab="Mean |y - y_hat|", ylim=yrange, col=cols, lwd=2 )
+ xlab="Time (hours)", ylab="Mean |y - y_hat|", ylim=yrange, col=cols, lwd=2
, lty=1
)
yrange = range( sapply(1:L, function(i) ( err[[i]]$abs$indices ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$abs$indices ), type="l",
yrange = range( sapply(1:L, function(i) ( err[[i]]$abs$indices ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$abs$indices ), type="l",
- xlab="Time (days)", ylab="Mean |y - y_hat|", ylim=yrange, col=cols, lwd=2 )
+ xlab="Time (days)", ylab="Mean |y - y_hat|", ylim=yrange, col=cols, lwd=2
, lty=1
)
yrange = range( sapply(1:L, function(i) ( err[[i]]$MAPE$day ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$MAPE$day ), type="l",
yrange = range( sapply(1:L, function(i) ( err[[i]]$MAPE$day ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$MAPE$day ), type="l",
- xlab="Time (hours)", ylab="Mean MAPE", ylim=yrange, col=cols, lwd=2 )
+ xlab="Time (hours)", ylab="Mean MAPE", ylim=yrange, col=cols, lwd=2
, lty=1
)
yrange = range( sapply(1:L, function(i) ( err[[i]]$MAPE$indices ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$MAPE$indices ), type="l",
yrange = range( sapply(1:L, function(i) ( err[[i]]$MAPE$indices ) ), na.rm=TRUE )
matplot( sapply( seq_len(L), function(i) err[[i]]$MAPE$indices ), type="l",
- xlab="Time (days)", ylab="Mean MAPE", ylim=yrange, col=cols, lwd=2 )
+ xlab="Time (days)", ylab="Mean MAPE", ylim=yrange, col=cols, lwd=2
, lty=1
)
}
#' Plot measured / predicted
}
#' Plot measured / predicted