similar to: Selecting transitive couples

Hi! I have a 3d array: > dim(ib5km15.dbc) [1] 190 241 19 and a set of positions to extract: > ib5km.lincol.random[1:3,] [,1] [,2] [1,] 78 70 [2,] 29 213 [3,] 180 22 Geting the values of a 2D array for that set of positions would be: > ima <- ib5km15.dbc[,,1] > ima[ib5km.lincol.random[1:10,]] but don't find the way for the case of the 3D array: >

Dear R-users and R-developpers, I've not been able to find 3d graphics in R, like Splus spin() or brush(). Are there any that I've missed and/or is anybody working on implementing these type of functions? Thanks Dr. Agustin Lobo Instituto de Ciencias de la Tierra (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona SPAIN tel 34 93409 5410 fax 34 93411 0012 alobo at ija.csic.es

I'd appreciate advice on the following: I've written an R function that uses 3 vectors (temperature, precipitation and potential evapotranspiration) for a given site, calculates a water budget (which implies few (<4) iterations), and, from thresults of this water budget, calculates a number of bioclimatic indices. Now I want to calculate these indices for a large number of sites

I'm experiencing segmentation faults from time to time that abruptely end my R sesion. In such cases, is there any way to recover the work done? (i.e., is there any temporal file that I could use?) Thanks Agus Dr. Agustin Lobo Instituto de Ciencias de la Tierra (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona SPAIN tel 34 93409 5410 fax 34 93411 0012 alobo at ija.csic.es

Hi! I'm trying to use readline() to ask for xlim() and ylim() values for a subsequent plot. I'm doing: lim <- readline("Enter xlim and ylim values: ") then would like to use the values in lim for a plot: plot(etc, xlim=lim[1:2],ylim=lim[3:4]) The problem is that lim is a character,i.e., "10,20,-10,50" or "c(10,20,-10,50)" depending on the input to

Hi there, Given 2 lists of integer vectors, i.e.: > lista1 $"1" [1] 1 34 5 $"2" [1] 2 1 $"3" [1] 3 10 15 > lista2 $"1" [1] 1 5 $"2" [1] 2 1 $"3" [1] 3 10 29 I want to obtain the union of both, defined as the union of the vectors, that is lista.union[[1]] <- union(lista1[[1]],lista2[[1]]): > lista.union

I've tried to move objects from S to R. In S+, I use data.dump() and data.restore(). I've made a file all.dum using data.dump(ls()) in S+, but the R command load("all.dum") gives an error: > load("/jaz/all.dum") Error: restore file corrupted -- no data loaded Is there any way to pass my objects from S+ to R? Thanks Agus Dr. Agustin Lobo Instituto de Ciencias de

I've found the following and I'm kind of confused: > summary(delme) Min. 1st Qu. Median Mean 3rd Qu. Max. 1 2950 5699 5756 8572 11680 > range(delme) [1] 1 11675 summary() and range() give different Max. value for the same vector! Agus Dr. Agustin Lobo Instituto de Ciencias de la Tierra (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona SPAIN tel

Hi! I'm writing down the methods that I used to migrate from Splus to R and expand a bit the R docs for this subject. Prior to send these notes to the FAQ curator I would like to make sure on the following: Given an ascii text file with several functions (i.e., myfuncs.txt), is there any equivalent to load(), so that the functions become usable in R? Or is it necessary to install it as a

Hi! I've been able to configure and make R (1,4,0 linux) except for the gnome support. I've checked all packages and libs mentioned in R-admin.pdf (and by Luke Tierny) to make sure that I have them and even using the following: ./configure --with-gnome --with-gnome-includes=/opt/gnome/include --with-gnome-libs=/opt/gnome/lib --with-libglade-config=opt/gnome/lib I still get: R is now

I have the following problem: I have a vector x of data (0<x<=1 ) with a U-shaped histogram and try to fit a beta distribution using fitdistr. In fact, hist(rbeta(100,0.1,0.1)) looks a lot like my data. The equivalent to the example in the manual sometimes work: > a <- rbeta(100,0.1,0.1) > fitdistr(x=a, "beta", start=list(shape1=0.1,shape2=0.1))1) > shape1

Hi! I run a function in which I fit a linear model, split the graphic window with layout, plot the data and the linear fit to the left, and the equation to the right using text() in this way (probably there's a better way): par(pty="m") plot(c(0,1),c(0,1), type="n",axes=F,xlab="",ylab="") text(0.25,0.5, paste("y =

Hi! Is it possible to display a 3d array as an RGB image? For example, given the following array: > matriz3d , , 1 [,1] [,2] [,3] [,4] [,5] [1,] 170 174 173 172 161 [2,] 171 178 174 166 149 [3,] 168 174 173 166 156 [4,] 171 170 173 166 164 [5,] 167 170 170 171 169 , , 2 [,1] [,2] [,3] [,4] [,5] [1,] 138 131 128 128 125 [2,] 138 127 129 122 134

I'm trying to export an R (1.4.0 on linux) matrix to python (+Pyclimate) via a binary file. The problem is that while R writes by columns, python reads by rows. For a 2D matrix, there is no problem: writeBin writes vectors, so writeBin(t(x)[1:length(x)],...) works fine. But the objects I want to export are 3D arrays (slice rasters of lat,lon for each time t). One way is to use a loop over

Hi! I'm upgrading to 1.6.2. (linux suse7.3). I see that there is a patch in R-release.diff.gz. 1. Is this patch stable? In other words, is it adviced that I install R-release.diff.gz ? 2. When I do: zcat R-release.diff.gz | patch -p1 -E I get many lines such as: patching file AUTHORS Reversed (or previously applied) patch detected! Assume -R? [n] Should I just accept the default

I'm geting ready to compile and install R-1.5.1 (R-patched_2002-07-15.tar.gz) on Linux (Suse 8.0, 2 AMD proc., ATLAS libs installed). My current version is R.1.4.0 I'd like to get advice on two doubts: 1. Is it better that I remove my current /usr/lib/R directory? 2. Regarding ATLAS, the R-1.5.1 R-admin manual states : "This is currently not supported. The problem is that SIGINT

Hi, I'm using plot.lda (R1.4.0 on linux) and, if dimen > 2, I get the error: Error in text.default(x, y, txt, cex = cex, font = font) : X11 font at size 16 could not be loaded The traceback indicates that the problem comes from pairs(). I've tried several cex and font values, always with the same error. The same command works fine with pdf() If I'm missing fonts, how

As far as I can see by searching the archieves, R could not take advantage of two processors in the same machine but there were plans for implementing multithreading in R. Is this work going on? On the other hand, if two processors are present in the same linux machine is it possible to focus one processor on tasks of the operative system and the other on R and/or other real computing? Could

Hi all, I have problems transforming a matrix into a data.frame: If I do: > is.matrix(parcelas.cobtot.conti) [1] TRUE > dim(parcelas.cobtot.conti) [1] 25 64 > a <- data.frame(parcelas.cobtot.conti) > is.data.frame(a) [1] TRUE > > dim(a) [1] 1600 3 Why a does not have the same dimensions than parcelas.cobtot.conti? I've tested with caith and the behaviour there

Hi! I'm trying to use R. In S+ I use to have different projects in different directories. To make the functions and other objects from a directory available and to use that directory as default, I just do: attach(directory, pos = 1) where directory is, for example, "/home/alobo/S/bolivia1" Acutually, I have a simpe function (work) to do this: > work function(directorio) {