similar to: Programming using formulas

Dear R experts, I wrote already a couple of weeks ago but did not get any reply. Therefore I'll rephrase my question and hope some of you can give me some answers or some advice where better to look for an answer. In general I'm looking for information on how to program functions using formulas in order to write a nonparametric multivariate regression function with possible

Dear R experts, we are preparing an R-package to compute the Oja Median which contains some C++ code in which random numbers are needed. To generate the random numbers we use the following Mersenne-Twister implementation: // MersenneTwister.h // Mersenne Twister random number generator -- a C++ class MTRand // Based on code by Makoto Matsumoto, Takuji Nishimura, and Shawn Cokus // Richard J.

Dear R useRs, The new contributed packages 'ICS' and 'ICSNP' are available on CRAN. Descriptions: The 'ICS' package implements the 2 scatter matrix transformation to obtain an invariant coordinate system or independent components, depending on the underlying assumptions. The result of the transformation is an object of the S4 class ics which is provided by this

Dear R useRs, The new contributed packages 'ICS' and 'ICSNP' are available on CRAN. Descriptions: The 'ICS' package implements the 2 scatter matrix transformation to obtain an invariant coordinate system or independent components, depending on the underlying assumptions. The result of the transformation is an object of the S4 class ics which is provided by this

Dear list: since it took me a little while to figure out how to make use of naresid, I thought that the below R code might be useful as an example on the help page. Regards, Markus # generate some data x1 <- runif(20) y <- 10 + 5*x1 + rnorm(20) summary(lm.0 <- lm(y ~ x1)) # append some NA's to y y <- c(y, rep(NA, 5)) # generate some further x1s x1 <- c(x1, runif(5)) #

Dear list: since it took me a little while to figure out how to make use of naresid, I thought that the below R code might be useful as an example on the help page. Regards, Markus # generate some data x1 <- runif(20) y <- 10 + 5*x1 + rnorm(20) summary(lm.0 <- lm(y ~ x1)) # append some NA's to y y <- c(y, rep(NA, 5)) # generate some further x1s x1 <- c(x1, runif(5)) #

I discovered recently that the phrase "Oja median" produces no hits in Jonathan Baron's very valuable R search engine. I found this surprising since I've long regarded this idea as one of the more interesting notions in the multivariate robustness literature. To begin to remedy this oversight I wrote a bivariate version and then decided that writing a general p-variate version

Hi, Can we get the code for calculating Oja median for multivariate data Thanks and Regards Rahul Agarwal Analyst Equities Quantitative Research UBS_ISC, Hyderabad On Net: 19 533 6363 [[alternative HTML version deleted]]

Hi Roger, As we know that The Oja median has (finite) breakdown point 2/n, i.e., is not robust in any reasonable sense, and is quite expensive to compute, so do we have some better methodology to compute multivariate median Rahul Agarwal Analyst Equities Quantitative Research UBS_ISC, Hyderabad On Net: 19 533 6363 [[alternative HTML version deleted]]

L?hetit postia @uta.fi osoitteille: suvi.virtanen Postiasi ei voi toimittaa perille, koska yliopiston s?hk?posti- j?rjestelm?ss? ei ole t?llaista yksik?sitteist? osoitetta. Voit etsi? tarkoittamasi henkil?n s?hk?postiosoitetta osoitehaun avulla http://www.uta.fi/hakemisto/sahkopostiosoitteet.html Huomaa kuitenkin, ett? s?hk?postiosoitteiden listassa ovat mukana vain ne yliopistoyhteis?n

Hello, A little informal input on the TCP/UDP checksum deferral/offload: it fails for me. I have one Xen host, running a snapshot of xen-unstable at 2006-03-24 (9435:11fee62328cc). Each of its domUs are allocated their own /30, with the help of a hacked-up vif script. Prior to working around the problem, Domain-0 could establish TCP connections to other domains. In the case of SSH, it would

Gee, it's been a while since I posted a new version of winetricks. This is rather a big update. As always, online at http://kegel.com/wine/winetricks Please submit bug reports to http://code.google.com/p/winezeug/issues/list svn log since version 20081127: r264 | daniel.r.kegel | 2009-01-16 Don't set winver to win2K when installing .net 1.1. r263 | daniel.r.kegel | 2009-01-16 Add

Dear R-users: The following code produces two cones in two panels. What I would like to have is to have them in one, and to meet in the origin. Does anyone have any good ideas how to do this? Thanks for your help Jaakko library(lattice) A<-matrix(ncol=2, nrow=64) for(i in 0:63) { A[i+1,1]<-sin(i/10) A[i+1,2]<-cos(i/10) }

I use kmeans to classify spectral events in high and low 1/3 octave bands: #Do cluster analysis CyclA<-data.frame(LlowA,LhghA) CntrA<-matrix(c(0.9,0.8,0.8,0.75,0.65,0.65), nrow = 3, ncol=2, byrow=TRUE) ClstA<-kmeans(CyclA,centers=CntrA,nstart=50,algorithm="MacQueen") This works well when the actual data shows 1,2 or 3 groups that are not "too close" in a cross plot.

Hi, I am trying to install the rJava package in a Linux environment (Platform: x86_64-suse-linux-gnu (64-bit)), where I am facing an issue which I cannot solve so far and where I could require some help: The installation seems to go fine, however in the last step of the installation where loading the package is tested, the process fails with the following error: ** testing if installed package

Dear useRs, version 1.0-0 of the "Formula" package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y ~ x1 + x2 | z1 + z2 where there are two groups of explanatory variables, e.g.,

Dear useRs, version 1.0-0 of the "Formula" package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y ~ x1 + x2 | z1 + z2 where there are two groups of explanatory variables, e.g.,

Dear All, I want to compare formulas resulting from using stepAIC() on glm() objects. I have 2 different formulas, say F1 and F2, which consist of about 10 factors (Y~A+B+C+...) but one has an interaction term (A:B) while the other doesn't. But when I perform F1==F2 it comes back as TRUE. Is there a proper way to compare formulas? ------------------------------------------- Jeffrey M

Does the latest winelib compile and work on the Sparc? Are there any special patches/instructions for this? Are there any mirrors that have the precompiled binary winelib for Sparc? Thanks in advance.

Hello! I hope, i'm doing right, writing here. I have a question about the R statistical formulas: were can I look, what are they? E. g. I'd like to know which exactly mathematical formula is used in Wilcoxon test, and others. Is there somewhere on the Internet information about that? Thanks a lot! Ilya. PS Sorry for any mistakes in my English, it is not my mother-tongue.