similar to: Fit with constraints (fwd)

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

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

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

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

Thanks for your help on moving files from S+ to R and on attaching directories. If attaching directories is NOT possible in R, is there any other way to separate objects into different folders or something similar? For example, to keep user functions for multivariate analysis in a different "place" than user functions on time series, or to keep data objects of project_1 in a different

Thanks for your help on moving files from S+ to R and on attaching directories. If attaching directories is NOT possible in R, is there any other way to separate objects into different folders or something similar? For example, to keep user functions for multivariate analysis in a different "place" than user functions on time series, or to keep data objects of project_1 in a different

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

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

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

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

Saving lda.default as text files from within gui="none" and from within gui="gnome" and comparing afterwards with diff, I get: alobo at humboldt:> diff lda.default.gnome.txt lda.default.nognome.txt 1c1 < function (x, grouping, prior = proportions, tol = 1, method = c("moment", --- > function (x, grouping, prior = proportions, tol = 1e-04, method =

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