similar to: Logic Error (PR#8677)

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) }

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.

This is just a note that R would get a lot more citations if the recommended citation was an article in a recognised journal or from a recognised publisher. I use R in work leading to publications often, and I strongly want to give the R core team credit for their work. However I find that I can't persuade my biological collaborators to include the current R citation (below) in their

Many new things at http://www.cuddyvalley.org/psychoR/ The scalassoc package, which fits exponential distance association models to indicator matrices, it now at version 1.0.0. It seems to be robust and can analyze large examples easily. It is a major improvement (in speed and robustness) over the distassoc package, which is at the same site. scalassoc does something neat (if you like

http://www.cuddyvalley.org/psychoR/code now contains version 1.4 of ca.R There are already many versions of CA and CCA in various R packages, but this one has some unique features, so maybe it is useful. The function ca() can do simple CA, but it also allows for linear restrictions on the row and column scores (they can be restricted to be in the span of a number of covariates). Thus the

Over the last couple of years I have written quite a few R programs for various "psychometric" techniques, and I am regularly updating and expanding what is there. I now have (wholly or partially), or have planned -- gifi package (update to homals on CRAN). Code for multiple correspondence analysis, nonlinear principal component analysis, nonlinear multiset canonical

Get articles from http://www.jstatsoft.org Firth Bradley-Terry Models in R Volume 12, Issue 01 Sturtz, Ligges, and Gelman R2WinBUGS: A Package for Running WinBUGS from R Volume 12, Issue 03 Mineo and Ruggieri A Software Tool for the Exponential Power Distribution: The normalp Package Volume 12, Issue 04 Ritz and Streibig Bioassay Analysis Using R Volume 12, Issue 05 Baddeley and Turner

Get articles from http://www.jstatsoft.org Firth Bradley-Terry Models in R Volume 12, Issue 01 Sturtz, Ligges, and Gelman R2WinBUGS: A Package for Running WinBUGS from R Volume 12, Issue 03 Mineo and Ruggieri A Software Tool for the Exponential Power Distribution: The normalp Package Volume 12, Issue 04 Ritz and Streibig Bioassay Analysis Using R Volume 12, Issue 05 Baddeley and Turner

JSS has been trying to develop a "Code Snippet" section (Hornik and Koenker, editors), but so far there have not been many submissions, possibly because this is not widely known. if you have some code that is not a full-blown package (just one or two functions) and you think it will be of general use, think about submitting. There is an excellent example in Volume 11 by Duncan

JSS is trying to develop it's Code Snippets section. We have some snippets lined up and one published in the latest volume http://www.jstatsoft.org/index.php?vol=16 If you have small chunks of code of obvious relevance to statistical computing (need not be in R) consider submitting it -- why keep it to yourself ? Just a matter of pasting a minimal amount of TeX into our templates and

The Journal of Statistical Software was founded by Jan de Leeuw in 1996. Currently, there are 26 volumes with 260 aricles and the journal publishes about 4 new volumes per year. The web address is www.jstatsoft.org. All articles are pdf files, since 2005 all produced from JSS LaTeX templates. JSS is edited by Jan de Leeuw, since 2005 jointly with Achim Zeileis of the Wirtschaftsuniversität Wien.

This is OS X 10.6. I can build llvm-gcc and llvm-g++ just fine, but the gfortran front end cannot find a symbol. Where would it normally be defined ? ==================================================================== /Developer/usr/bin/g++-4.2 -g -O2 -DIN_GCC -W -Wall -Wwrite-strings -Wstrict-prototypes -Wmissing-prototypes -pedantic -Wno-long-long -Wno-variadic-macros

http://idisk.mac.com/jdeleeuw-Public/utilities/apl/apl.R Dedicated to the IBM 2741. Implemented for general multidimensional arrays: drop, take, reshape, shape, rank, select, generalized inner product, generalized outer product, representation, base value, join, expand, reduce, scan, member of, ravel, compress, tranpose, rotate Basically, the APL-I part is complete, and after some testing

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

but maybe of use to some: http://gifi.stat.ucla.edu/psychoR/derivatives Computes generalized eigenvalue solutions Ax=\lambda Bx and generalized singular value solutions Rz=\gamma Px and R'x=\gamma Qy for matrices that are differentiable functions of a vector of parameters. Along with the decomposition the code returns arrays with all first-order partial derivatives of the values/vector wrt

At www.jstatsoft.org we now have three special volumes Volume 22, Ecology and Ecological Modelling in R (Thomas Kneib and Thomas Petzoldt, eds.) Volume 20, Psychometrics in R (Jan de Leeuw and Patrick Mair, eds.) Volume 18, Spectroscopy and Chemometrics in R (Kate Mullen and Ivo van Stokkum, eds.) Coming soon -- Econometrics in R -- Political Methodology in R Gleam in my eye --

in 2007 we will publish seven volumes of JSS (four of them special volumes). Volume 18 Spectroscopy and Chemometrics in R Volume 19 regular Volume 20 Psychometrics in R Volume 21 regular Volume 22 Ecology and Ecological Modeling in R Volume 23 is being started, it is a regular volume, and may not fill up completely. Volume 24 is 'Statistical Modeling of Social Networks with

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,