On Sat, 2011-02-05 at 23:39 -0600, Paul Johnson wrote: > On Sat, Feb 5, 2011 at 9:19 AM, David Winsemius <dwinsem...@comcast.net> > wrote: <snip /> > > cbind(scalermca[,1] * 0.827094, scalermca[,2] * -0.7644828) > [,1] [,2] > 1 1.06070017 -0.81549999 > 2 0.77057891 0.63456780 > 3 1.07031764 -1.30675217 > 4 1.07031764 -1.30675217 > 5 0.23075886 0.90002547 > 6 0.69488883 0.60993995 > 7 0.10530240 0.78445402 > 8 -0.27026650 0.44225049 > 9 0.13426089 1.15670532 > 10 0.11861965 0.64778456 > 11 0.23807570 1.21775202 > 12 1.01156703 -0.01927226 > 13 0.28051938 -0.59805897 > 14 -1.17343686 -0.27122981 > 15 -0.83838041 -0.64003061 > 16 -0.05453708 -0.22925816 > 17 -0.91732401 -0.49899374 > 18 -0.92694148 -0.00774156 > 19 -1.30251038 -0.34994509 > 20 -1.30251038 -0.34994509 > > So, that does reproduce SAS exactly. And I'm a little frustrated I > can't remember the matrix command to get that multiplication done > without cbinding the 2 columns together that way.
You might have been thinking of sweep(): sweep(scalermca[,1:2], 2, c(0.827094,-0.7644828), "*") <snip/> HTH G -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/ UK. WC1E 6BT. [w] http://www.freshwaters.org.uk %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% ______________________________________________ R-help@r-project.org mailing list 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.
