R语言画图，幂律分布

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问一下如何画幂律分布的拟合曲线，我用power.law.fit 只能求出幂指数，不知道怎么画出那条曲线

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：R语言 Power 拟合曲线 不知道 law 如何

用curve函数可以实现，我用的例子是power.law.fit函数帮助文档下面的example:

1. library(igraph)
   
2. g <- barabasi.game(1000) # increase this number to have a better estimate
   
3. d <- degree(g, mode="in")
   
4. fit1 <- power.law.fit(d+1, 10)
   
5. ##　这里可以作出拟合的图出来
   
6. curve(x^fit1$alpha,from=1,to=10)

复制代码

ps：以后再提问时，最好能够把数据什么的一并贴出来，这样可以方便别人写程序。

jiangbeilu 发表于 2014-11-24 22:59
用curve函数可以实现，我用的例子是power.law.fit函数帮助文档下面的example:

老师，您好！



我也想画幂律分布图，我的数据是K（如下）：想画类似的图1的幂律分布图，我的数据分布式如图2，我的应该是像图1的那样的两条直线拟合。
感谢您的指导！
                           
                                                         图1                                                                图2
> k<-data.frame(e,c)
>k   

e          c1   20 300.0038262   30 125.0035873   40  77.7817464   50  56.2544845   60  44.0050516   70  36.1167547   80  30.6184938   90  26.5693629  100  23.46427210 110  21.00810511 120  19.01699612 130  17.37047113 140  15.98632214 150  14.80654115 160  13.78903316 170  12.90251417 180  12.12325118 190  11.43291319 200  10.81712020 210  10.26443321 220   9.76563722 230   9.31322923 240   8.90103524 250   8.52393125 260   8.17762626 270   7.85850327 280   7.56348828 290   7.28995629 300   7.03564730 310   6.79860931 320   6.57714532 330   6.36977333 340   6.17519434 350   5.99226335 360   5.81996836 370   5.65741137 380   5.50379238 390   5.35839639 400   5.22058440 410   5.08978041 420   4.96546542 430   4.84717043 440   4.73447144 450   4.62698245 460   4.52435246 470   4.42626047 480   4.33241548 490   4.24254849 500   4.15641350 510   4.07378451 520   3.99445452 530   3.91822953 540   3.84493354 550   3.77440255 560   3.70648256 570   3.64103557 580   3.57792758 590   3.51703859 600   3.45825360 610   3.40146761 620   3.34658062 630   3.29350063 640   3.24214064 650   3.19242065 660   3.14426266 670   3.09759667 680   3.05235368 690   3.00847269 700   2.96589270 710   2.92455771 720   2.88441572 730   2.84541573 740   2.80751074 750   2.77065675 760   2.73481076 770   2.69993277 780   2.66598578 790   2.63293279 800   2.60073980 810   2.56937581 820   2.53880782 830   2.50900783 840   2.47994784 850   2.451601

jiangbeilu 发表于 2014-11-24 22:59
用curve函数可以实现，我用的例子是power.law.fit函数帮助文档下面的example:

您好，想请一下如果要用自己的数据的话，您公式中的参数我要怎么改啊...