Does a bad powerlaw fit not even be closer to the data points??

19 views
Skip to first unread message

Nixon Raj

unread,
May 25, 2021, 12:27:16 PM5/25/21
to powerlaw-general
I am trying to fit a data with power law, to show it is a scale free network.

When i use powerlaw, i could able to get a better fit with alpha=2.03 and xmin 3.52 (please check the attached figure 1). 

So extracted the values less than 3.5 as separate dataset and wanted to know the alpha value for that dataset. This is just to differentiate how scale free the higher values that had a power law fit is from the lower values. (In other words, the better fit should have alpha between 2 to 3 and the lower values with bad fit should have alpha value very low.)

So i tried to power law fit the lower values, ie dataset with values less than 3.5, i get a weird fit, which is away from the data points. (please check the attached figure 2).

Figure 3, shows the overall power law fit (ignore the lomax fit in it)
Questions:

1. Why is the power_law.plot_pdf is very away from the data points?
2. Why is it even not in the trendline of the data? Any reasons.
3. The lower value data points visually looks very linear, so atleast i should have got a line trend through the data points even if it is a bad poer law fit right!! May be i'm wrong. Please help me sort this out.

Figure 1: The tail region (alpha:2.03, amin:3.52)
Figure1.png
Figure 2: Values < 3.52
Figure2.png
Figure 3: (All data)
Figure3.png


Jeff Alstott

unread,
May 25, 2021, 6:50:10 PM5/25/21
to powerlaw...@googlegroups.com
Thanks for using powerlaw!

You've created an xmax in the data. Do you get the same effect if when you fit with powerlaw you add xmax=3.5?

I dimly recall there *might* be some issue with values smaller than 1, but I doubt it. Worth trying multiplying all the data by the smallest value, just in case. But I apologize that I haven't maintained powerlaw in years.


--
You received this message because you are subscribed to the Google Groups "powerlaw-general" group.
To unsubscribe from this group and stop receiving emails from it, send an email to powerlaw-gener...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/powerlaw-general/e17c4b01-4e56-475c-b67b-ab1a0fcb559cn%40googlegroups.com.

Nixon Raj

unread,
May 25, 2021, 10:01:20 PM5/25/21
to powerlaw...@googlegroups.com
Adding 3.5 as xmax does not help either. And multiplying the smallest value with the data also didn't help. 

I have attached the data along with this mail, can you check it. I just need a better fit and the alpha value of it. Please help me.

You received this message because you are subscribed to a topic in the Google Groups "powerlaw-general" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/powerlaw-general/RjGtD2mk-gY/unsubscribe.
To unsubscribe from this group and all its topics, send an email to powerlaw-gener...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/powerlaw-general/CAGoRfgE7QJD31A23REs9dARDoML%2BHu3DSJR01YCV0LkAUWxBvg%40mail.gmail.com.


--
Regards 

Nixon Raj N
Department of Biological Science and Technology
Institute of Bioinformatics and Systems Biology
National Chiao Tung University
208 Lab Building 1, 75 Bo-Ai St.
Dong District, Hsinchu, Taiwan 30062
(R.O.C.)
Mob:+886-975426494 
0ffice ext: 56997
data.txt
Reply all
Reply to author
Forward
0 new messages