similar to: data.restore: warnings

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: >

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

Thanks to all who answered my question (appended to this message). Two solutions: 1. If the user answers with the 4 lim values (i.e., 10,20,1,100), limStr <- readline("Enter xlim and ylim values: "); lim <- as.numeric(unlist(strsplit(limStr, ","))); 2. If the user answers in R vector notation (i.e., c(10,20,1,100) ) limStr <- readline("Enter xlim and ylim

Does anybody know if there are functions to: 1. Define a convex hull on a space with more than 2 dimensions? chull just works for a plane. 2. Numerically select elements within a given complex hull. 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

Agustin Lobo asks (and I interpolate in his question) > -----Original Message----- > From: Agustin Lobo [mailto:alobo at ija.csic.es] > Sent: Monday, 23 July 2001 7:45 AM > To: r-help > Subject: [R] lapply and for > > > Given a list as such: > > > milista > $"1": > [1] 23 25 11 > > $"2": > [1] 34 2 > >

-----Forwarded Message----- > From: Ernesto Jardim <ernesto at ipimar.pt> > To: Agustin Lobo <alobo at ija.csic.es> > Subject: Re: [R-gui] NEdit Highligth patterns for R > Date: 10 Oct 2002 14:36:24 +0100 > > Hi > > You have to install the *.pats file. > > Use > > nedit -import R-5.1.pats > > then "save defaults". > > If

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

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 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'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

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! 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'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

Dear R makers and users. I summarize answers to my question. Thanks for your help. The conclusion seems to be that, in terms of performance with large datasets, Matlab=Octave > R > Splus5.1 > Splus5.0 but Matlab/Octave lack many of the stat functions of R/S, and, R does not have as many time-series functions as S+ (but see last message to r-announce by adrian.trapletti at

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

I would like to fit a model of the type Y = Xf + e where f are cover fractions (or percent cover). Therefore, I must constrain the fit to 1. sum(f) = 1 2. 0<= f <=1 How can I introduce these 2 constraints in a ML fit? Currently, I'm using nlregb and I'm using the "upper and low" keywords to set the bounds and use the following to ensure that sum(f) = 1: let's say