search for: angel_lul

...istian) that partly solves the problem, but I still would like to be able to restrict some of the independent variables in a nls model to be always >0, (is there a way to do it)?? Thanks, Angel >From: "Christian Ritz" <ritz at dina.kvl.dk> >To: "Angel -" <angel_lul at hotmail.com> >CC: <r-help at stat.math.ethz.ch> >Subject: Re: [R] nls: Missing value or an Infinity produced when evaluating >the model >Date: Wed, 23 Apr 2003 15:22:37 +0200 > >Hi Angel, > >I tried reparametrise your model, setting: > >BirthMass^0.25=u &g...

Which is the best editor for .R files? I currently use kate on my linux as it has R highlighting and allows me to split the window into two: in one I edit the .R file and in the other I have a shell so I run R and can easily copy and paste the code. There are some features that I don't like and I am having a look on some alternatives. I've heard wonders of emacs with ess but I am a

Hi, I am trying to fit a sigmoid curve to some data with nls but I am getting into some trouble. Seems that the optimization method is getting down to some parameter estimates that make the equation unsolvable. This is an example: >growth<-data.frame(Time=c(5,7,9,11,13,15,17,19,21,23,25,27),BodyMass=c(45,85,125,210,300,485,570,700,830,940,1030,1120))

Why do x<-b%*%ginv(A) and x<-solve(A,b) give different results?. It seems that I am missing some basic feature of matrix indexing. e.g.: A<-matrix(c(0,-4,4,0),nrow=2,ncol=2) b<-c(-16,0) x<-b%*%ginv(A);x x<-solve(A,b);x Thanks in advance, Angel

...ellmund/indhold.html For native Xemacs for Windows: www.xemacs.org Xemacs is easy to install, use and configure (more easy than the GNU-Version). Under Linux both comes with the most distributions (I use both under SuSE 8.2). HTH Thomas -----Urspr?ngliche Nachricht----- Von: Angel [mailto:angel_lul at hotmail.com] Gesendet: Donnerstag, 20. November 2003 22:27 An: r-help at stat.math.ethz.ch Betreff: [R] best editor for .R files Which is the best editor for .R files? I currently use kate on my linux as it has R highlighting and allows me to split the window into two: in one I edit the .R...

I want to 1) Select for each of the n points in a matrix A, those of the m points in B that lay within a given radius. 2) Of those points within the radius, select the k nearest ones. What I now do is 1) Create an n*m matrix C were I put the distances from all the points in B to the points in A and make NA those cells were the distance is larger than the radius. (The points are geographical

I needed a symbolic and algebraic computation software. The best solution I've found so far is maxima computer algebra system. I was wondering if there is any R package that has similar features as maxima's or somebody could point me to reference manuals on how to use R for this purpose. Thanks, Angel

Hi, Is there a function like rank but that solves the ties by randomly assigning a value (doesn't average ranks of ties). This is what I actually need: I want to make NA all elements of each column in an array that are ranked in a position larger that rankmax for each column. # Say I've got an array b: b<-cbind(c(1:5,5:1),c(1,12,14,2,5,4:8)) #> b # [,1] [,2] #[1,] 1 1

Is there a way to make a shading interpolation on an image plot? Something similar to matlab 'shading interp', I think it is called Gouraud shading. What I want is to make a image plot look nicer. with image() it looks very facetted, and I would like to make it look smoother. I've tried with interp.surface() in fields package but it (obviously) makes nan values at the borders and

### First, is there a way to access a univariate time series as a matrix instead of a vector? # For example: data(UKLungDeaths) # If I do apply(mdeaths,1,cumsum) # Gives an error as mdeaths is not a matrix but a vector, although when I look at it : mdeaths # the ts object has a matrix like "appearance" # The only way of doing it I've found is: mdeaths2<-as.matrix(mdeaths)

In plot(), when using option asp=1 the xlim and ylim have no effect because they are changed changed in order to fill the whole plot region. Is there a way to automatically set xlim and ylim when asp has been set to 1? For example: #This is a box of the plot ranges I want: boxxy=rbind(c(-1,2),c(-1,-1),c(1,-1),c(1,2),c(-1,2)) #Without asp=1 I get what I want (i.e. I can't see the box because it

I have always been intrigued by why ?for (or ?if,?while,etc) leave R wanting for more: > ?for + I know the help for these is in ?Control, but I sometimes make the mistake of typing ?for instead. What is R expecting me to say to finish the statement? Thanks, Angel [[alternative HTML version deleted]]

I have updated my R to the current R version (1.8.1) using apt-get of the R binaries for debian sid in cran. Now, when I do: png(file="myplot.png") I get: Error in X11(paste("png::", filename, sep = ""), width, height, pointsize, : unable to start device PNG In addition: Warning message: No png support in this version of R what do I have to do to get the

I have two dataframes A and B consisting of latitude longitude coordinates of points. For each point in A I want to find the k-nearest neighbours in B. Currently, I calculate the distance from each point in A to all the points in B (using rdist.earth() in fields package), sort the points of B by distance and select the k nearest ones. Is there a more efficient way to do this? Function knearneigh

I have a binary file representing a matrix with columns of different variable type/mode (i.e. the file was saved from C using double and int variables) I want to import it into R, I've been reading the R Import/Export and I am able to use readBin for a file that contains only one variable type (either all doubles or all ints) but I can not find the way to use it when you have a file with

Hi, I have a binary file with mixture of short integers (2 bytes) and floats (4 bytes). To get the data into R I use readBin reading each value one at a time (n=1) but it is giving me headaches, two questions: -In the "what" option is there any difference between "integer" and "int"? -To read the short int I use what="int" and size=2 and to read the

In function daystoyears in package pastecs, I get this (wrong?) result with 1995: > daystoyears(1,datemin="1/1/1997",dateformat="m/d/Y") 1997.001 > daystoyears(1,datemin="1/1/1995",dateformat="m/d/Y") 1994.999 Any insights? Thanks Angel

When I do a plot, e.g. plot(1:10) and resize the window so that the x-axis becomes too small to hold all the x-axis labels, R automatically makes some of the labels disappear so that the remaining fit in the available space. I would like to be able to tell R which labels should not be removed. I've tried plotting the axis with axis() but this behaviour is still there.

Does anybody know of any R functions to perform chronological clustering as explained in: Legendre, P., S. Dallot & L. Legendre. 1985. Succession of species within a community: chronological clustering, with applications to marine and freshwater zooplankton. American Naturalist 125: 257-288. http://www.fas.umontreal.ca/BIOL/legendre/reprints/succession_of_species.pdf Thanks, Angel

Dear all, Trying to solve a problem I had (see thread "putting NAs at the end" ) I've noticed a difference in system time requirements between using apply and sort (or order) to order each row or column of a matrix compared to a vectorized function I wrote. Using apply is much faster when the number of loops (number of rows or columns to order) is low BUT much slower when number of