Hello, I am trying to fit my data to the default autoregressive models in R. I'm trying to apply the effectiveSize function from the coda package to a list of data frames using;
sapply(splits, function(x) coda::effectiveSize(x["V5"])) However when I do, I get the error; Error in ar.yw.default(x, aic = aic, order.max = order.max, na.action = na.action, : 'order.max' must be >= 1 For some reason the command above works fine if my list is smaller (otherwise, not properly split). And it works perfectly fine on one of my data sets. If I apply it to anything else (data of a similar set up) I get the error. I've tried looking for solutions on the internet but so far I've come up with nothing. Has this happened to anyone before? Can I bypass this order. max call? [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.