Deal with par issue

Author: JamesRamsay5

DESCRIPTION

@@ -2,9 +2,15 @@ Package: fda
 Version: 6.0.5
 Date: 2022-07-2
 Title: Functional Data Analysis
-Author: J. O. Ramsay <ramsay@psych.mcgill.ca> [aut,cre], 
- Spencer Graves <spencer.graves@effectivedefense.org> [ctb], 
- Giles Hooker <gjh27@cornell.edu>  [ctb]
+Authors@R: c(person("James", "Ramsay", 
+          role=c("aut","cre"),
+          email="james.ramsay@mcgill.ca"), 
+          person("Giles", "Hooker",)
+          role="ctb",
+          email="gjh27@cornell.edu"),
+          person("Spencer, "Graves", 
+          role="ctb",
+          email="spencer.graves@effectivedefense.org"))
 Maintainer: J. O. Ramsay <ramsay@psych.mcgill.ca>
 Depends: R (>= 3.5), splines, fds, deSolve
 Suggests: 

R/smooth.morph.R

@@ -43,7 +43,7 @@ smooth.morph <- function(x, y, ylim, WfdPar,
   #  last modified 3 June 2022 by Jim Ramsay
 
   nobs <- length(x)        #  number of observations
-  wt=matrix(1,nobs,1)
+  wt   <- matrix(1,nobs,1)
   wt[nobs] <- 10
 
   # check consistency of x and y and convert to column matrices
@@ -103,8 +103,11 @@ smooth.morph <- function(x, y, ylim, WfdPar,
 
   #  check range of x
 
-  if (x[1] < xlim[1] || x[nobs] > xlim[2]) {
-    stop("Values in x are out of bounds.")
+  if (abs(x[1] - xlim[1]) > 1e-7 || abs(x[nobs] - xlim[2]) > 1e-7) {
+    stop("Values in x are out of bounds by more than 1e-7")
+  }  else {
+    x[1   ] <- xlim[1]
+    x[nobs] <- xlim[2]
   }
 
   #  initialize matrix Kmat defining penalty term

man/Data2fd.Rd

@@ -269,6 +269,7 @@ b1.1 <- create.bspline.basis(nbasis=1, norder=1)
 y12 <- 1:2
 # smooth data, giving a constant function with value 1.5
 fd1.1 <- Data2fd(y12, basisobj=b1.1)
+oldpar <- par(no.readonly=TRUE)
 plot(fd1.1)
 # now repeat the analysis with some smoothing, which moves the
 # toward 0.
@@ -401,5 +402,6 @@ AmRevLin.ct <- as.numeric(AmRevYrs.ct-AmRev.ct[1])
 AmRevYrs.ct <- as.numeric(AmRevYrs.ct)/rescale
 AmRevLin.ct <- as.numeric(AmRevLin.ct)/rescale
 fitLin.ct   <- Data2fd(AmRevYrs.ct, AmRevLin.ct, BspRev.ct)
+par(oldpar)
 }
 \keyword{smooth}

man/Fperm.fd.Rd

@@ -134,6 +134,7 @@ Fperm.fd(yfdPar, xfdlist, betalist, wt=NULL, nperm=200,
   \code{\link{Fstat.fd}}
 }
 \examples{
+oldpar <- par(no.readonly=TRUE)
 ##
 ## 1.  yfdPar = vector
 ##
@@ -216,6 +217,7 @@ betalist = list(fdPar(hgtbasis,2,1e-6),fdPar(hgtbasis,2,1e-6))
 # Should use nperm = 200 or so,
 # but use 10 to save test time
 Fres = Fperm.fd(hgtfd,xfdlist,betalist,nperm=10)
+par(oldpar)
 }
 
 }

man/Lfd.Rd

@@ -73,5 +73,4 @@ bfdobj  <- fd(Lcoef,Lbasis)
 bwtlist <- fd2list(bfdobj)
 harmaccelLfd <- Lfd(3, bwtlist)
 }
-% docclass is function
 \keyword{smooth}

man/StatSciChinese.Rd

@@ -33,6 +33,7 @@
   95, 9-15.
 }
 \examples{
+oldpar <- par(no.readonly=TRUE)
 data(StatSciChinese)
 
 i <- 3
@@ -48,5 +49,6 @@ plot(StatSci1[, 1:2], type='l')
 
 mark <- seq(1, 601, 12)
 points(StatSci1[mark, 1], StatSci1[mark, 2])
+par(oldpar)
 }
 \keyword{datasets}

man/arithmetic.fd.Rd

@@ -68,6 +68,7 @@ times.fd(e1, e2, basisobj=NULL)
   \code{\link{exponentiate.fd}}
 }
 \examples{
+oldpar <- par(no.readonly=TRUE)
 ##
 ## add a parabola to itself
 ##
@@ -165,6 +166,9 @@ all.equal(coef(-parab4.5), coefneg)
 plot(parab4.5-1)
 
 plot(1-parab4.5)
+
+par(oldpar)
+
 }
 \keyword{smooth}
 

man/as.fd.Rd

@@ -92,6 +92,7 @@ fdx. <- eval.fd(x., fd.)
 
 rfdx <- range(fdx.)
 
+oldpar <- par(no.readonly= TRUE)
 plot(range(x2), range(y2, fx., rfdx), type='n')
 points(x2, y2)
 lines(x., sin((x.-0.5)*pi), lty='dashed')
@@ -130,6 +131,7 @@ lines(cars.spl)
 sp. <- with(cars, seq(min(speed), max(speed), len=101))
 d. <- eval.fd(sp., cars.fd)
 lines(sp., d., lty=2, col='red', lwd=3)
+par(oldpar)
 }
 \keyword{smooth}
 \keyword{manip}

man/axisIntervals.Rd

@@ -82,6 +82,7 @@ axesIntervals(side=1:2, atTick1=fda::monthBegin.5, atTick2=fda::monthEnd.5,
   \code{\link{monthLetters}}
 }
 \examples{
+oldpar <- par(no.readonly= TRUE)
 daybasis65 <- create.fourier.basis(c(0, 365), 65)
 
 daytempfd <- with(CanadianWeather, smooth.basis(
@@ -109,7 +110,7 @@ axesIntervals()
 
 axisIntervals(3, labels=monthLetters, cex.lab=1.2, line=-0.5)
 # 'line' argument here is passed to 'axis' via '...'
-
+par(oldpar)
 }
 % docclass is function
 \keyword{smooth}

man/cca.fd.Rd

@@ -88,6 +88,9 @@ round(ccafd0$ccacorr[1:6],3)
 #  compute a VARIMAX rotation of the canonical variables
 ccafd <- varmx.cca.fd(ccafd0)
 #  plot the canonical weight functions
+oldpar <- par(no.readonly= TRUE)
 plot.cca.fd(ccafd)
+par(oldpar)
+
 }
 \keyword{smooth}