Corrected conidence bands in plot.pca.fd

Giles Hooker · Giles Hooker · commit 8658fc1b47ee · 2022-04-06T09:14:03.000-07:00
diff --git a/R/plot.pca.fd.R b/R/plot.pca.fd.R
@@ -55,7 +55,7 @@ plot.pca.fd <- function(x, nx = 128, pointplot = TRUE, harm = 0,
 #        
       iharm <- harm[jharm] 
       if(expand == 0) {
-        fac <- sqrt(pcafd$values[iharm])
+        fac <- 2*sqrt(pcafd$values[iharm])
       } else {
         fac <- expand
       }
@@ -123,7 +123,7 @@ plot.pca.fd <- function(x, nx = 128, pointplot = TRUE, harm = 0,
         iharm <- harm[jharm]
 #        
         fac <- {
-          if (expand == 0) sqrt(pcafd$values[iharm]) 
+          if (expand == 0) 2*sqrt(pcafd$values[iharm]) 
           else expand
         }
         

Original file line number	Diff line number	Diff line change
`@@ -55,7 +55,7 @@ plot.pca.fd <- function(x, nx = 128, pointplot = TRUE, harm = 0,`
`55`	`55`	`#`
`56`	`56`	`iharm <- harm[jharm]`
`57`	`57`	`if(expand == 0) {`
`58`		`- fac <- sqrt(pcafd$values[iharm])`
	`58`	`+ fac <- 2*sqrt(pcafd$values[iharm])`
`59`	`59`	`} else {`
`60`	`60`	`fac <- expand`
`61`	`61`	`}`
`@@ -123,7 +123,7 @@ plot.pca.fd <- function(x, nx = 128, pointplot = TRUE, harm = 0,`
`123`	`123`	`iharm <- harm[jharm]`
`124`	`124`	`#`
`125`	`125`	`fac <- {`
`126`		`- if (expand == 0) sqrt(pcafd$values[iharm])`
	`126`	`+ if (expand == 0) 2*sqrt(pcafd$values[iharm])`
`127`	`127`	`else expand`
`128`	`128`	`}`
`129`	`129`