如何用R来检验Power Law Distribution

1.假如不是power law,而是weibull会比较方便,比较简单.

如果你是要估计weibull参数,在matlab只需wblfit

如果是做Kolmogorov-Smirnov goodness-of-fit hypothesis test

可用 kstest- single sample

kstest2- two-sample

简单范例:

%generate random numbers from a Weibull distribution

x = wblrnd(1, 2, 100, 1);

kstest(x, [x wblcdf(x, 1, 2)])

2.因为估power law 的参数,要用到zeta function所以应用

matlab zeta(1,x) zeta(x)很容易算出power-law exponent.

很感谢楼上的诸位的帮忙，呵呵，假如我需要估计{1-(1-q)x/kappa}^(1/1-q)『物理学中的Tsllie统计』，需要自定义函数这种类型的，matlab中没的，那如何啊？，

还有bootstrap,它是怎么样的思想啊，如果是样本的重复抽样那有什么意义？

谢谢各位，谢谢yiyo900

[此贴子已经被作者于2007-6-6 12:26:52编辑过]

1.不太清楚什么是『物理学中的Tsllie统计』,你的单一方程中

有三个variables:q,x,kappa,你要解的是x?q?or kappa?

在matlab解非线性单一方程的解,是用fzero

很容易使用,可参考Optimization Toolbox.pdf page-8-116

内有范例,及敎你如何自订函数.

2.至于bootstrap的观念,如何抽样,如何....

可参考matlab statistics.pdf page-13-60

及s-plus STATMAN2 chap 35

有了初步观念,方便讨论.

不好意思，说不清楚，x是观测值，kappa和q是待估计两参数；基本明白我再研究下

谢谢；

以下是引用yiyo900在2007-4-17 9:58:00的发言：

1.MLE估计power-law exponent

可参考"Problems with Fitting to the Power-Law Distribution"

公式(3) the derivative of the log-likelihood function

http://arxiv.org/PS_cache/cond-mat/pdf/0402/0402322v3.pdf

公式(3)的右半部,根据数据就可算出.左半部查表,或是利用

Matlab,maple,mathematica 都可算出.

文献中公式右半部算出是0.2739,利用matlab可算出

power-law exponent : gamma(MLE)等于2.544

2.S(x), F*(x)

X S(x) F*(x) |F*(x)-S(x)|

1 0.8302 0.83025

2 0.93794 0.93669

3 …… …….

4 …… …….

5 …… …….

6 …… …….

空白的部份,请自行运算.

3.在s-plus运行bootstrap很简单,请先参考STATMAN2.pdf

Chap35,page-538.内有范例

请问STATMAN2.pdf Chap35,page-538哪里有啊？谢谢啊

splus-statman2-chap35 124078.pdf (137.64 KB)

2007-6-8 08:13:00 上传

如何用R来检验Power Law Distribution

太感谢了，很有用！

请问各位啊，能做伽马分布的检验吗？主要是想知道用r能生成伽马分布函数吗？

几位大侠给的材料很有用，多谢啦！