search for: gematrix

...ik.uni-frankfurt.de) # Reference for Escoufier's equivalents vectors (test III.5): # Escoufier Y., 1970. Echantillonnage dans une population de variables # aleatoires r??les. Publ. Inst. Statis. Univ. Paris 19 Fasc 4, 1-47. ... Error in getClass(Class, where = topenv(parent.frame())) : "geMatrix" is not a defined class Calls: new -> getClass Execution halted

...ow fast OLS coefficients can be obtained. ### including tests to exclude missing observations where necessary. ### ### for a more informed article (and person), see Bates, 'Least Squares ### Calculations in R', Rnews 2004-1. the code here does not test his sparse ### matrix examples, or his geMatrix/poMatrix examples. ### ### Basic Results: for the examples that I tried, typical relative time ### factors of the algorithms were about ### ### lm lmfit solve crossprod cholesky (special-case 2vars) ### 1.0 0.5 0.3 0.15 0.17 0.1 ### ### there was no advantage to chol...

I think one needs an LU decomposition rather than QR. However, I couldn't find anything off the shelf to do an LU, other than learning that determinant() now uses LU instead of QR or SVD, so the code to do it must be in there for those that want it. You'll probably need to divide rows of U by the first entry if you insist on the unique reduced REF. However, I can't see any reason

Hello, A collegue of mine has compared the runtime of a linear model + anova in SAS and S+. He got the same results, but SAS took a bit more than a minute whereas S+ took 17 minutes. I've tried it in R (1.9.0) and it took 15 min. Neither machine run out of memory, and I assume that all machines have similar hardware, but the S+ and SAS machines are on windows whereas the R machine is Redhat