Dear all,
I am trying to fit a straight line to the degree distribution of my network (which is scale-free), using following commands:
g <- barabasi.game(1000, directed =F)
d <- degree(g)
dd <- degree.distribution(g,cumulative=TRUE)
plot(dd,log="xy")
alpha <- power.law.fit(d, xmin=10)
lines(5:50, 10*(5:50)^(-coef(alpha)+1))
but if I change the value xmin, the value of alpha changes and hence different straight line are generated! Is there a way to know which xmin is best to use? and which alpha and fitted line is more accurate?
Thanks in advance,
Narsis