Hello,
I am having two difficulties with xspline:
1) Passing lty=2 (via ...) appears to have no effect
2) Apparently depending on the current plot axes, I may or may not
get the failure: add_point - reached MAXNUMPTS (25200)
This happens both on Linux/R 2.8.0 (details below) and Windows Vista/R
2.8.0 patched (2008-10-21 r46766)
Simplest case:
plot(c(0,1),c(0,1))
xspline(c(0,1),c(0,1)) # failure MAXNUMPTS
xspline(c(0,0.5,1),c(0,1,0),border=2,lty=2) # OK-ish, red, but no dashes
More complex case:
x <- c(-50.25, -50, -49.75, -0.75, -0.5, -0.25, 0.25, 0.5, 0.75, 49.75,
50, 50.25)
y <- c(0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1)
s <- c(1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1)
plot(c(0,1),c(0,1))
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2) # fails MAXNUMPTS
points <-
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2,draw=FALSE) #fails
MAXNUMPTS
plot(c(-50,50),c(0,1))
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2) # OK-ish: red, but
no dashes
points <-
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2,draw=FALSE) # works
I note from the help page that "Device coordinates are used (with an
anisotropy correction if needed.)", so I assume the difference between
the more complex cases which work and fail is to do with that, but I am
stumped as to how to use that information! (The only help message
mentioning MAXNUMPTS I can find is to do with grid.xspline, so I cannot
use the solution suggested there without recoding everything into grid.)
Thankyou,
Keith Ponting
R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> plot(c(0,1),c(0,1))
> xspline(c(0,1),c(0,1)) # failure MAXNUMPTS
Error in xspline(c(0, 1), c(0, 1)) :
add_point - reached MAXNUMPTS (25200)
> xspline(c(0,0.5,1),c(0,1,0),border=2,lty=2) # OK-ish, red, but no
dashes
> x <- c(-50.25, -50, -49.75, -0.75, -0.5, -0.25, 0.25, 0.5, 0.75,
49.75, 50, 50.25)
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2) # fails MAXNUMPTS
points <-
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2,draw=FALSE) #fails
MAXNUMPTS
> y <- c(0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1)
> s <- c(1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1)
> plot(c(0,1),c(0,1))
> xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2) # fails
MAXNUMPTS
Error in xspline(x, y, s, open = TRUE, repEnds = FALSE, border = 2, lty
= 2) :
add_point - reached MAXNUMPTS (25200)
> points <-
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2,draw=FALSE) #fails
MAXNUMPTS
Error in xspline(x, y, s, open = TRUE, repEnds = FALSE, border = 2, lty
= 2, :
add_point - reached MAXNUMPTS (25200)
> plot(c(-50,50),c(0,1))
> xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2) # OK-ish: red,
but no dashes
> points <-
xspline(x,y,s,open=TRUE,repEnds=FALSE,border=2,lty=2,draw=FALSE) # works
> sessionInfo()
R version 2.8.0 (2008-10-20)
x86_64-unknown-linux-gnu
locale:
LC_CTYPE=en_GB.UTF-8;LC_NUMERIC=C;LC_TIME=en_GB.UTF-8;LC_COLLATE=en_GB.U
TF-8;LC_MONETARY=C;LC_MESSAGES=en_GB.UTF-8;LC_PAPER=en_GB.UTF-8;LC_NAME=
C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_GB.UTF-8;LC_IDENTIFICATI
ON=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
> Sys.info()[c(1:3,5)]
sysname
release
"Linux"
"2.6.18.2-34-default"
version
machine
"#1 SMP Mon Nov 27 11:46:27 UTC 2006"
"x86_64"
Keith Ponting
Aurix Ltd, Malvern WR14 3SZ UK
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
