Re: [R] Weibull distribution

2006-07-20 Thread Leaf Sun
: William Asquith, [EMAIL PROTECTED] Sent: 2006-07-17, 16:18:31 To: Leaf Sun, [EMAIL PROTECTED] Subject: Re: [R] Weibull distribution Do not have answer per se, but if you are seeking some comparisons-- try three parameter Weibull as implemented by the lmomco package

[R] (no subject)

2006-07-20 Thread Leaf Sun
please? Any suggestion would be appreciated! Leaf - Original Message - From: William Asquith, [EMAIL PROTECTED] Sent: 2006-07-17, 16:18:31 To: Leaf Sun, [EMAIL PROTECTED] Subject: Re: [R] Weibull distribution Do not have answer per se, but if you are seeking some

Re: [R] Weibull distribution

2006-07-21 Thread Leaf Sun
Thanks for the suggestion! I switched to optimize(), al - optimize(f.fn, lower = 0.1, upper =100,tol=0.001); the warnings were gone and it works stably. But when I tried al - uniroot(f.fn, lower = 0.1, upper =100,tol=0.001); error occured: f() values at end points not of opposite sign. The

[R] number of iteration s exceeded maximum of 50

2006-06-14 Thread Leaf Sun
Hi all, I found r-site-research not work for me these days. When I was doing nls( ) , there was an error number of iterations exceeded maximum of 50. I set number in nls.control which is supposed to control the number of iterations but it didn't work well. Could anybody with this experience

Re: [R] number of iteration s exceeded maximum of 50

2006-06-16 Thread Leaf Sun
Sorry, I thought it was a straightforward question inside which I was stuck . I used nls( ) to estimate a and b in this function. nls(y~ a*x^b,start=list(a=a1,b=b1) seems the start list I gave was not able to reach convergence and it gave notes: number of iteration s exceeded maximum of 50.

Re: [R] number of iteration s exceeded maximum of 50

2006-06-19 Thread Leaf Sun
Thanks to Douglas and all others who responded. I applied nls(y ~ a*x^b, start = list(a = a1, b = b1), control = list(maxiter = 500), trace=TRUE) to increase the number of iterations, found it successful. The suggestion Douglas raised in plotting the data and then tracing the

[R] Weibull distribution

2006-07-17 Thread Leaf Sun
Hi all, By its definition, the mean and variance of two-par. Weibull distribution are: (www.wikipedia.org) I was wondering, if given mean and sd. could we parameterize the distribution? I tried this in R. gamma.fun - function(mu,sd,start=100) { f.fn - function(alpha)

[R] how to get the row name?

2006-01-30 Thread Leaf Sun
Hi R-listers, I have a simple question about a data frame. I sorted a data set by one of the variable in some condition (eg. X=0), the followed is part of the achieved. I was wondering how can I get the row name, i. e. (1202, 2077 , 2328, 3341,... ) and save them as a vector.

[R] Sorting a data frame by one of the variables

2005-10-15 Thread Leaf Sun
Dear all, I have a date frame like this: X Y Z 22 24 4.3 2.3 3.4 5.3 . 57.223.434 What my purpose is: to sort the data frame by either X, Y or Z. sample output is (sorted by X) : X Y Z 2.3 3.4 5.3 . .. 22 24 4.3 ... 57.2 23.4 34 I have

[R] Finding the neighbors of the point

2005-10-21 Thread Leaf Sun
Dear all, I got point data of trees. I was wondering if anybody has experience in searching the neighbors within a specified distance efficiently. XY Z 99 34 65 98 35 29 98 34 28 99 33 33 98 32 23 99 33 21 99 33

[R] Errors occured

2005-10-22 Thread Leaf Sun
Hi all, Has anybody have the experience in the errors: Error in data.frame(..., check.names=FALSE): arguments imply differing number of rows: 343,15 This is the error occured in the middle of the program. I don't think the data frame has any problem, if there is problem with the program, why

Re: [R] Finding the neighbors of the point

2005-10-25 Thread Leaf Sun
, 04:31:54 you wrote: === On Mon, 24 Oct 2005, Leaf Sun wrote: Running R 2.2.0 on winXP. Computer P4 CPU 3.2G and 1G of RAM. Please try the attached Windows binary package. Look at the help page for ann.dist(). It returns a list of three elements, the first, lnn, gives the index numbers

[R] breaks in hist()

2005-11-02 Thread Leaf Sun
Dear listers, A quick question about breaks in hist(). The histogram is highly screwed to the right, say, the range of the vector is [0, 2], but 95% of the value is squeezed in the interval (0.01, 0.2). My question is : how to set the breaks then make the histogram look even? Thanks in

Re: [R] Visualizing a Data Distribution -- Was: breaks in hist()

2005-11-03 Thread Leaf Sun
=== At 2005-11-02, 12:07:12 you wrote: === Leaf Sun wrote: The histogram is highly screwed to the right, say, the range of the vector is [0, 2], but 95% of the value is squeezed in the interval (0.01, 0.2). I guess the histogram is as you wrote. See http://web.maths.unsw.edu.au

[R] OLS variables

2005-11-06 Thread Leaf Sun
Dear all, Is there any simple way in R that can I put the all the interactions of the variables in the OLS model? e.g. I have a bunch of variables, x1,x2, x20... I expect then to have interaction (e.g. x1*x2, x3*x4*x5... ) with some combinations(2 way or higher dimensions). Is there

Re: [R] OLS variables

2005-11-06 Thread Leaf Sun
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Leaf Sun Sent: Sunday, November 06, 2005 3:11 AM To: r-help@stat.math.ethz.ch Subject: [R] OLS variables Dear all, Is there any simple way in R that can I put the all

Re: [R] Variogram

2005-11-09 Thread Leaf Sun
Sorry about this, I didn't know. I guess I have posted too much garbage on here. Thanks to Edzer for your answers! Leaf === At 2005-11-09, 02:14:21 you wrote: === Leaf, please note that r-help is not the appropriate place to ask package-specific questions. We have r-sig-geo for

[R] Point pattern to grid

2005-11-17 Thread Leaf Sun
Dear all, I'd like to change a point pattern to a grid of cells and use one of the variables as the output. e.g. The point pattern is of a window of (500*500) and several features such as pH, SoilType etc. I like to divide it into a grid with cell size 5*5, and use the mean of the point

Re: [R] Point pattern to grid

2005-11-18 Thread Leaf Sun
2005, Leaf Sun wrote: Dear all, I'd like to change a point pattern to a grid of cells and use one of the variables as the output. e.g. The point pattern is of a window of (500*500) and several features such as pH, SoilType etc. I like to divide it into a grid with cell size 5*5, and use

[R] read.list()

2005-11-24 Thread Leaf Sun
Hi all, I need to write and read a list in R. I did r.site.search, found there is a package rmutil doing this, unfortunately it is not on the list of package. In another words, I can't install it from any CRAN mirror. Anybody has idea about this? or any suggestion about the list? Thanks!