Hi
Keith Ponting wrote: > 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) Both of these problems were bugs, which have now been fixed in the development version of R. Thanks for the report! Paul > 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. -- Dr Paul Murrell Department of Statistics The University of Auckland Private Bag 92019 Auckland New Zealand 64 9 3737599 x85392 [EMAIL PROTECTED] http://www.stat.auckland.ac.nz/~paul/ ______________________________________________ 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.
