Many thanks for supporting the igraph package.
It continues to be indispensable for my research on digitally networked social movements.
I have 3 questions about the function "fit_power_law".
- The igraph documentation suggests that the bfgs function is used to estimate the power law alpha, but I think the C implementation relies on the Broyden-Fletcher-Goldfarb-Shanno optimization function of the lbfgs library instead. Is that correct?
- The fit_power_law function relies on the MLE function of the stat4 package. I am curious why this was deprecated, given the availability of plfit and MLE parameters. Is this simply a memory issue?
- How to interpret the p-value of the Kolmogorov-Smirnov test? The igraph help file states: "Small p-values (less than 0.05) indicate that the test rejected the hypothesis that the original data could have been drawn from the fitted power-law distribution" . The C implementation of the KS test in igraph uses the Hurwitz Zeta function. Shouldn't this mean that high p-values indicate a good model fit, as suggested by Clauset et al (2009:678)?
Apologies if these questions are obtuse. I am not a trained programmer nor statistician, and could not solve these questions from any documentation or other source online.
Thank you for your time and help.
Your work allows me to learn much about R (and networks!).