pacf in devel seems by default to return a different number of lags than 1.7.1 for $pacf. I don't see any mention of this in the NEWS file, or any change in the documentation, so I suspect it is and error, though it may be an undocumented improvement. (Newbie question: How is the simplest way to display a function like pacf.default that is not exported from a namespace?) Paul Gilbert
Try (in R-devel) getAnywhere("pacf.default") or even ts:::pacf.default -roger Paul Gilbert wrote:> pacf in devel seems by default to return a different number of lags > than 1.7.1 for $pacf. I don't see any mention of this in the NEWS > file, or any change in the documentation, so I suspect it is and > error, though it may be an undocumented improvement. > > (Newbie question: How is the simplest way to display a function like > pacf.default that is not exported from a namespace?) > > Paul Gilbert > > ______________________________________________ > R-devel@stat.math.ethz.ch mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-devel >-- Together, we can stop attaching Word documents http://www.fsf.org/philosophy/no-word-attachments.html
Looking even more carefully, it seems both sampleT and the lag.max correction for nser need to be fixed: sampleT <- if (is.matrix(x)) nrow(x) else length(x) if (is.null(lag.max)) lag.max <- floor(10 * (log10(sampleT) - if (is.matrix(x)) log10(ncol(x)) else 0)) Paul Gilbert Paul Gilbert wrote:> On further examination, it looks like > sampleT <- length(x) > in pacf.default is not handling the matrix case properly and probably > needs to be replaced by something like > sampleT <- if (is.matrix(x)) nrow(x) else length(x) > > Paul Gilbert > Paul Gilbert wrote: > >> pacf in devel seems by default to return a different number of lags >> than 1.7.1 for $pacf. I don't see any mention of this in the NEWS >> file, or any change in the documentation, so I suspect it is and >> error, though it may be an undocumented improvement. >> >> (Newbie question: How is the simplest way to display a function like >> pacf.default that is not exported from a namespace?) >> >> Paul Gilbert >> >> ______________________________________________ >> R-devel@stat.math.ethz.ch mailing list >> https://www.stat.math.ethz.ch/mailman/listinfo/r-devel >> >