similar to: Getting rid of loops?

Dear listers, I'm having a bad R day. I just can't think of the vectorized equivalent of: for (ii in 1:n) aa[ii] = bb[ii,cc[ii]] Any suggestion received with embarrassment and gratitude Simon Gatehouse CSIRO Exploration and Mining, Newbigin Close off Julius Ave North Ryde, NSW Mail: PO Box 136, North Ryde NSW 1670, Australia

Dear R users, I am using R 5.1 on Windows 2000. Is there a system call or somesuch way of identifying the graphic window which currently has focus and using that information to make it the active device for futher plotting etc. I have done a quick seaqrch in R the archives using 'focus active device' but came up with nothing useful. Any information would be most welcome. Thanks in

Dear Any, I am trying to initialise Rjava with the objective of running ORCA within R. I am running R1.5.1 on Windows2000 and get the same problem with Sun Java run time jre-1.2.2_013 Java j2re-1.4.1_01 Microsoft Java Virtual Machine Build 3805 .JavaInit() is failing at : .C("s_start_VM", as.character(classPath), as.character(classPath),

The following might be of benefit to those who find themselves in a similar predicament with an apparently corrupted R in Windows 2000. The problem was traced back to a corrupted register in my user profile in Windows 2000 apparently caused by a DLL in Visual Basic installation not self registering correctly. The fix was to copy a new uncorrupted register from the default user. Though some of

All: Are there any functions out there for selecting all the points in a region of a plot. I envision something like the identify() function except one could circle a cloud of points (and perhaps a vector would be returned of the same length as the points plotted indicating logical membership in the circled cloud). Perhaps someone has done something with the locator() function that would

Dear R users, I run R on windows 2000 and windows 98 I am seeking opinion, or better still a code snippet, on what might be the quickest, not necessarily the best or preferred way, of making the Windows color selection common dialog available in R. Would it be an executable called with system() or shell() or perhaps a tcltk module or have I missed something obvious? I don't believe that I

Dear All, Has anyone worked with linking R to both ArcView and MapInfo GIS systems? My desire is to capture object selection from GIS and pass to R and to return spatially based statistical model results from R to GIS. If both could simultaneously access the same data file this would be a bonus. My inclination is to put together some routines based on ODBC linking but thought I would seek like

Microsoft Windows R users, I am operating Windows 2000 (build2195) with R1.7. It may be pertinent that I am using a dual head screen with the initial RGui filing the entire area of both screens. When starting R I use the .First() function to add menu items to the RGui interface using winMenuAdd() and winMenuAddItem(). The menus do not display until the RGui window is physically manipulated in

Greets, I'm trying to use nlminb() to estimate the parameters of a bivariate normal sample and during one of the iterations it passes a parameter vector to the likelihood function resulting in an invalid covariance matrix that causes dmvnorm() to throw an error. Thus, it seems I need to somehow communicate to nlminb() that the final three parameters in my parameter vector are used to

Hi, I've noticed I get different results fitting a function to some data on my laptop to when I do it on my computer at work. Here's a code snippet of what I do: ##------------------------------------------------------------------ require(circular) ## for Bessel function I.0 ## Data: dd <- c(0.9975948929787, 0.9093316197395, 0.7838819026947, 0.9096108675003, 0.8901804089546,

Hi, I'm trying to write a simple function like case1 <- function (m, cov, Q, R) { theta <- (acos(R/sqrt(Q^3))) beta <- (-2)*sqrt(Q)*cos(theta/3)+m[1]/3 rho1 <- (-2)*sqrt(Q)*cos((theta+2*pi)/3)+m[1]/3 rho2 <- (-2)*sqrt(Q)*cos((theta-2*pi)/3)+m[1]/3 stderrb <- deltamethod( ~(-2)*sqrt(Q)*cos(theta/3)+x1/3,m,cov) stderrr1 <- deltamethod(

Hi everyones, my function like; e <- rnorm(n=50, mean=0, sd=sqrt(0.5625)) x0 <- c(rep(1,50)) x1 <- rnorm(n=50,mean=2,sd=1) x2 <- rnorm(n=50,mean=2,sd=1) x3 <- rnorm(n=50,mean=2,sd=1) x4 <- rnorm(n=50,mean=2,sd=1) y <- 1+ 2*x1+4*x2+3*x3+2*x4+e x2[1] = 10 #influential observarion y[1] = 10 #influential observarion data.x <- matrix(c(x0,x1,x2,x3,x4),ncol=5) data.y

Dear anybody dealing with TCLTK on R I'm trying, without any success, to get a simple TCLTK window to remain floating on top of the RGui and RConsol windows when the later have focus. I cannot find reference in the TCKTK documentation and fear that I'm missing something quite fundamental. My objective is to make a simple button populated menu to invoke R functions. I wish the TCLTK menu

Hi Dear All, This is a function which contains Covariance Ratio and Likelihood Distance values (CVRi, LDi). i want to compute the all row's values, that is run this function for nrow(X) times. The X and Y matrices are;

Using the heatmap function in mva, it seems to be hard to use different colors in the edges leading to different groups of objects, as commonly done in many heatmaps in the microarray graphics. Any suggestions? Thanks. max

Hello! I have a problem with integrate() in my function nctspa(). Integrate produces an error message "evaluation of function gave a result of wrong length". I don't know what that means. Could anyone suggest me what is wrong with my function? These are the examples of function calls that work OK: nctspa(a=1:10,n=5) nctspa(a=1:10, n=5, mu=2, theta=3, renorm=0) This does not work:

Hello, I'm relatively new to R. I've read the intro guide, but I can't quite figure out how to: I have a function: Jcost <- function (theta, in, out) { a <- output - input %*% theta 1/2 * t(a) %*% a } where "theta" is a 2x1 matrix "in" is a 20x2 matrix "out" is a 20x1 matrix return value is a scaler This works well when I only want to

Dear All, I need to plot some points on the surface of a sphere, but I am not sure about how to proceed to achieve this in R (or if it is suitable for this at all). In any case, I am not looking for really fancy visualizations; for instance you can consider the images between formulae 5 and 6 at http://bit.ly/hOgK9h Any suggestion is appreciated. Cheers Lorenzo

Not an R question directly, but has anyone got a method for placing a moderately large number of (near) equi-spaced points on a sphere? I have a nasty feeling platonic solids are needed for exact solutions and I'm thinking of samplings involving around 200 - 1000 regularly-spaced points, Thanks, Richard Rowe Richard Rowe Senior Lecturer Department of Zoology and Tropical Ecology, James

If possible I would like to add two sub-menus to the R Console under Windows. For example, I would like to add: winMenuAddItem("File", "Load CSV...", "loadCSV()") winMenuAddItem("File", "Save CSV...", "saveCSV()") and have them appear under the initial 'File' item rather than add a new 'File' menu item. I seem to