search for: plfit

...s and predicted values, blah... ... then... lines( sort( $predicted ) ~ sort( ht ) ) which results in a line that isn't smooth (which I knew would happen). I've checked the FAQ,docs and archives and I'm not sure if there's function that will so what Heut et. al (2004) do with their plfit(). So, is there already an R function, or process to do this, or will I have to write one? Thanks, Jeff. --- Jeff D. Hamann Forest Informatics, Inc. PO Box 1421 Corvallis, Oregon USA 97339-1421 (office) 541-754-1428 (cell) 541-740-5988 jeff.hamann at forestinformatics.com www.forestinformatics.co...

...rovided 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=alpha) hi...