[Dplr-commits] r688 - branches/redfit/R
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Thu Sep 12 18:37:03 CEST 2013
Author: mvkorpel
Date: 2013-09-12 18:37:03 +0200 (Thu, 12 Sep 2013)
New Revision: 688
Modified:
branches/redfit/R/redfit.R
Log:
Edited a comment
Modified: branches/redfit/R/redfit.R
===================================================================
--- branches/redfit/R/redfit.R 2013-09-06 19:43:00 UTC (rev 687)
+++ branches/redfit/R/redfit.R 2013-09-12 16:37:03 UTC (rev 688)
@@ -229,7 +229,7 @@
ww[, i] <- redfitWinwgt(twk, iwin2)
}
## determine autospectrum of input data
- lmfitfun <- tryCatch(match.fun(".lm.fit"),
+ lmfitfun <- tryCatch(match.fun("lm.fit"),
error = function(...) match.fun("lm.fit"))
gxx <- .Call(dplR.spectr, t2, x2, np, ww, tr[[1]], tr[[2]], tr[[3]],
nseg, nfreq, avgdt, freq, dn50, segskip, lmfitfun)
@@ -330,8 +330,8 @@
## can assume that real(nout/2) was supposed to be real(nout)/2.
## sqrtHalfNfreq <- sqrt(nfreq %/% 2)
sqrtHalfNfreq <- sqrt(nfreq / 2)
- ## dplR: NOTE! Is round() the right function to use? Maybe floor()
- ## for the lower limit and ceiling for the higher limit?
+ ## dplR: NOTE! Is round() the right function to use? Maybe ceiling
+ ## for the lower limit and floor for the higher limit?
rcritlo <- round((-0.79557086 + 1.0088719 * sqrtHalfNfreq)^2)
rcrithi <- round(( 0.75751462 + 0.9955133 * sqrtHalfNfreq)^2)
} else {
More information about the Dplr-commits
mailing list