Hello, Reading the code I have found that in the line 653 (Pnew[i + r * j] = tmp;) of http://svn.r-project.org/R/trunk/src/library/stats/src/arima.c (latest as of now) the multiplier for j is r instead of rd (which is the dimension of the matrix). Was it intentional? That seems at least worth a comment why it is the case instead of Pnew[i + rd * j] = tmp; Thanks and Kind Regards, gabor
______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
