Examples for 'igraph::fit_power_law'


Fitting a power-law distribution function to discrete data

Aliases: fit_power_law power.law.fit

Keywords: graphs

### ** Examples


# This should approximately yield the correct exponent 3
g <- barabasi.game(1000) # increase this number to have a better estimate
d <- degree(g, mode="in")
fit1 <- fit_power_law(d+1, 10)
fit2 <- fit_power_law(d+1, 10, implementation="R.mle")

fit1$alpha
[1] 3.192652
stats4::coef(fit2)
   alpha 
3.192646 
fit1$logLik
[1] -73.17276
stats4::logLik(fit2)
'log Lik.' -73.17275 (df=1)

[Package igraph version 1.3.1 Index]