power.law.fit
fits a power-law distribution to a
data set.power.law.fit(x, xmin = NULL, start = 2, ...)
x
will be used. This argument makes it
possible to fit only the tail of the distribution.mle
.mle
mle-class
for details.The additional arguments are passed to the mle function, so it is possible to change the optimization method and/or its parameters.
mle
# This should approximately yield the correct exponent 3
g <- barabasi.game(1000) # increase this number to have a better estimation
d <- degree(g, mode="in")
power.law.fit(d+1, 20)
Run the code above in your browser using DataLab