similar to: FYI: APL in R

I wrote some functions for multiway CANDECOMP, i.e. for least squares fitting of a_{i_1\cdots i_m}\approx\sum_{s=1}^p x^1_{i_1s}x^1_{i_1s}\cdots x^m_{i_ms} with arrays of arbitrary dimension. Reminded me of the good old APL days. I could not find this in the archives, but if it's already there, I would appreciate if someone let me know.

Hi, is anybody out there who knows APL and would help me to translate 52 lines of APL code into propper R? best, Torsten

anacor-0.9.0 is on CRAN (by De Leeuw and Mair) anacor does correspondence analysis and canonical correspondence analysis. It can make row plots, column plots, joint plots, Benz?cri plots, regression plots, and transformation plots. Where appropriate, plots can be in 3d using either rgl or scatterplot3d. Row and column points can be in standard scaling, Benz?cri scaling, Goodman scaling,

anacor-0.9.0 is on CRAN (by De Leeuw and Mair) anacor does correspondence analysis and canonical correspondence analysis. It can make row plots, column plots, joint plots, Benz?cri plots, regression plots, and transformation plots. Where appropriate, plots can be in 3d using either rgl or scatterplot3d. Row and column points can be in standard scaling, Benz?cri scaling, Goodman scaling,

http://idisk.mac.com/jdeleeuw-Public/jacobi This is paper/software for various techniques based on Jacobi plane rotations. There is R code for -- classical cyclical Jacobi Eigen diagonalization -- Jacobi-based SVD diagonalization -- approximate simultaneous diagonalization of symmetric matrices (De Leeuw/Pruzansky 1978) -- approximate simultaneous diagonalization of rectangular matrices

Hello! The Array Programming Languages Conference will be from 24th - 27th July 2000 Please look at: http://stat.cs.tu-berlin.de/APL-Berlin-2000/ Chairman's Address APL Berlin 2000 will take place in a city which is worth a visit even without array processing languages. There is a fortunate coincidence: Much activity in scientific and commercial institutions is based on such languages.

Hello! I've been trying 3.0.8 on OS X 10.6.7 today and have been getting some errors. The command I'm using is: /usr/local/bin/rsync -aNHAXx --protect-args --fileflags --force-change --progress /Users/davidallie/Documents/ /Volumes/david.allie/Documents/ The target Documents folder is on an iDisk mounted as david.allie . The errors are happening at the beginning and end. BEGINNING:

I'm very new to R and utterly blown away by not only the language but the unbelievable set of packages and the documentation and the documentation standards and... I was an early APL user and never lost my love for it and in R I find most of the essential things I loved about APL except for one thing. At this early stage of my learning I can't yet determine if there is a way to

Hi, I have mentioned before that we have added some sdt dtrace probes in SAX, our APL interpreter. Encouraged by Angelo and Jignesh, I have created a small document (5 pages) describing our experience with it, together with some problems we have encountered and some scripts we use for pretty-printing dtrace outputs. The said document can be found at

nowhere does it say that centos is approved for use in DoD. it is not on the APL, only RedHat and SuSE

x<-seq(-1,1,length=10) y<-seq(-1,1,length=10) a<-matrix(c(1,2,2,1),2,2) b<-matrix(c(2,1,1,2),2,2) fv<-function(x,y) { m<-x*a+y*b t<-m[1,1]+m[2,2]; d<-m[1,1]*m[2,2]-m[1,2]^2 return((t-sqrt(t^2-4*d))/2) } gv<-function(x,y) { t<-x*(a[1,1]+a[2,2])+y*(b[1,1]+b[2,2]) d<-(x*a[1,1]+y*b[1,1])*(x*a[2,2]+y*b[2,2])-(x*a[1,2]+y*b[1,2])^2 return((t-sqrt(t^2-4*d))/2) }

Hello people, Is there some simple way of spinning and/or flipping arrays in R? Here's what I mean. Suppose that foo is a 2x3x4 array with the following contents: (I know this is different than typing 'foo' at and R prompt, but I'm so used to row major order from using APL, I have a hard time with R's output) > foo[1,,] [,1] [,2] [,3] [,4] [1,] 1 2 3 4

I'd like my compiler to emit proper line number information. The docs talk about Instruction::setDebugLoc(), but that method doesn't actually have to be in my 2.7 LLVM Debian package. What's the correct way of doing this? In addition, can anyone point me at an example of how to emit a comment attached to an instruction (or function)? -- ┌─── ｄｇ＠ｃｏｗｌａｒｋ．ｃｏｍ ─────

Dear r-users, Suppose I have three datasets: Dataset-1: Date x y Jan-1,2005 120 230 Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21 Jan-7,2005 11 299 Mar-5,2005 200 311 Dataset-2: Date x y Jan-2,2005 123 -125 Jan-3,2005 -110 300 Jan-4,2005 114 -21

Given a structure like this (using C syntax rather than LLVM because I'm still not fluent with LLVM assembly): struct Object { int i1; int i2; int i3; }; Then, if I have an int* pointer which I know is pointing to the i3 element, what's the best way of recovering a pointer to the structure as a whole? My fallback option is to cast the pointer to an int64, use getelementptr to

Given an nxm matrix A I want to compute the nxm matrix B whose ij-th element is the sum of the elements of A lying on the diagonal that ends with element ij, i.e., b_ij = a_ij + a_(i-1)(j-1) + a_(i-2)(j-2) + ... In APL (which I no longer use), I would use the 'rotate' operator to derive an array whose columns are diagonals of the given array and then cumulate down columns. Is

I bogged down about half way through reading the Green Book, in part because it became increasingly difficult to understand how some of the ideas related to R, as opposed to S (which I have not used). Does any reader know whether there is a document that points out differences between S and R that would be helpful in reading the Green Book? Ideally, perhaps, I need a "crib sheet" to

Dear developeRs (and other abuseRs ;-), I would like to contribute a patch against functions head() and tail() of package utils to allow for a negative 'n' argument. This allows to extract all but the first/last 'n' elements/rows/lines of an object, similar to the "drop" operator of APL. [1] I put the patched head.R and head.Rd files, along with diff files in

OpenSSH-2.5.2.p1 won't connect to OpenSSH-2.5.1p2 using version 2 protocol, quitting with the error message: [dunlap at tesla dunlap]$ ssh -2 kraken 7a 90 3f 39 37 67 0d 9e ac 43 74 c3 83 83 f5 a2 Disconnecting: Bad packet length 2056273721. tesla is Linux tesla.apl.washington.edu 2.2.16-3 #1 Mon Jun 19 19:11:44 EDT 2000 i686 unknown Intel RHL6.2 with OpenSSH-2.5.2.p1 compiled from sources

I see Thomas has already nailed this one, so it becomes a non-issue. Nevertheless I feel moved to say I think the idea would have been a step in the wrong direction in the first place. It comes from a desire to make R behave "a bit more like matlab" and that is ultimately unhelpful. Having tried to teach generations of students how to use the system (S-PLUS, but it could equally well