similar to: ylim in plot(corresp(,df=2))

Thanks to the many people who adviced me on this topic. xfig works very well, at least for what I have to do now. I'd like to answer the message by Paul Marrell as well as encourage others to do so: On Fri, 23 Nov 2001, Paul Murrell wrote: > I would be very interested to know the sorts of things you do to add things > interactively to plots. > > As one reply mentioned, there

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

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

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

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

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

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

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

Regarding user's home and working directories, is Rprofile just an old name of Renviron? or (can) should we have both files with different type of definitions? (talking about R.1.4.0 in linux). 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

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

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

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