search for: powerlaw

Dear All, I'm getting confused with the concept R uses to do regression using lm. I'm afmiliar with gnuplot and the build-in fit command, but couldn't get R's lm to work on my data. I know that my data follows a powerlaw or maybe an exponential function, and I'd like to determine the best fitting factors for the formula: a*x^b where b < 0. I've tried thge follwoing: s <- lm(y ~ x) > summary(s) Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -18.454 -7.577 -2.861...

p { margin-bottom: 0.08in; } Dear All, I will appreciate some suggestions of R packages for "ESTIMATION OF THE EXPONENT OF POWER-LAW FREQUENCY DISTRIBUTIONS". I have been searching at the R-help list several keywords for this subject and I did not find a very specific package, except the useful normalp package. I believe there are others but I was not able to identify it. I have

...s. This provided me with a visual cue of a potential power-law distribution (at least, that's what I am boldly thinking). So I went ahead and performed the procedure as suggested by Clauset, Shalizi, and Newman (2007/09), by reusing the essentials found in http://tuvalu.santafe.edu/~aaronc/powerlaws/plfit.r. I obtained and tested the following fitting estimates based on dat1: xmin <- 0.01715686 alpha <- 2.381581 I then simply wanted to print the observed and fitted dists in one plot, so I ran: library(ggplot2) library(VGAM) dat2 <- rpareto(length(dat1), location=xmin, shape=alp...