কোনও বিতরণ একটি পাওয়ার আইন অনুসরণ করে কিনা তা কীভাবে পরীক্ষা করবেন?

Question

কতজন ব্যবহারকারী কতগুলি প্রশ্ন পোস্ট করে সে সম্পর্কে আমার কাছে ডেটা রয়েছে। উদাহরণ স্বরূপ,

[UserCount, QuestionCount] 
[2, 100] 
[9, 10] 
[3, 80] 
... ...

এর অর্থ হ'ল 2 জন ব্যবহারকারী প্রতিটি 100 টি প্রশ্ন পোস্ট করেছেন, 9 জন ব্যবহারকারী প্রতিটি 10 ​​টি প্রশ্ন পোস্ট করেছেন। সুতরাং, UserCount, QuestionCountবিতরণটি কোনও পাওয়ার আইন অনুসরণ করে কিনা তা আমি কীভাবে নির্ধারণ করতে পারি ?

আমি পাওয়ারআলা প্যাকেজটি পেয়েছি । যাইহোক, আমি মূল্যায়ন করতে কেবলমাত্র একটি গ্রুপের পাস করতে পারি। (এই প্যাকেজে প্রদত্ত উদাহরণটি শব্দটির ফ্রিকোয়েন্সি)) সুতরাং আমি কীভাবে এই প্যাকেজটি ব্যবহার করব? নাকি আমার কিছু ভুল আছে? আমার কাছে প্রতিটি ব্যবহারকারীর প্রশ্ন গণনার ডেটাও রয়েছে, যেমন [100, 100, 10, 10, 10 ... ]। আমি যদি এই ডেটা প্যাকেজে পাস করি তবে কী পাব?

Comments

এই নিবন্ধটি পাওয়ার আইন বিতরণের পাশাপাশি আর কোডের জন্য কীভাবে পরীক্ষা করতে হয় তা গাণিতিক বিবরণ সরবরাহ করে। ক্লজেট এট আল। "এমিরিকাল ডেটাতে পাওয়ার আইন বিতরণ।"

— সাইকোরাক্স বলছেন

ধন্যবাদ। তবে আমি কাগজে সমস্ত কন্টেন্ট ধরতে পারি না। আমি বৈধতা করতে কিছু পন্থা চাই। 'একক গ্রুপের ডেটা প্যাকেজে পাস করার অর্থ' সম্পর্কে আপনার কী ধারণা আছে? উদাহরণস্বরূপ, ডেটা সেট হ'ল শব্দ ফ্রিকোয়েন্সি। আবার ধন্যবাদ.

— তৃতীয় দিন

দুঃখিত, আমি এই প্যাকেজটির সাথে পরিচিত নই।

— সাইকোরাক্স বলছেন

Answer 1

ক্লাউসেট এট আল। এর মতে, আপনি poweRlawপ্যাকেজ সহ পাওয়ার আইন লেজটি এইভাবে পরীক্ষা করেন :

1. শক্তি আইন বিতরণ অবজেক্টটি তৈরি করুন। এই ক্ষেত্রে, আপনার ডেটা পৃথক, সুতরাং শ্রেণীর বিচ্ছিন্ন সংস্করণটি ব্যবহার করুন

data <- c(100, 100, 10, 10, 10 ...)
data_pl <- displ$new(data)

2. $x_{min}$$\alpha$

est <- estimate_xmin(data_pl)
data_pl$xmin <- est$xmin
data_pl$pars <- est$pars

শেষ দুটি রেখাটি একটি লাইন হিসাবে আবারও লেখা যেতে পারে

data_pl$xmin <- est

এছাড়াও, এই সময়ে, আপনি কেএস পরিসংখ্যান দেখতে পারেন:

est$KS

3. $p$

bs <- bootstrap_p(data_pl)
bs$p

এটি কিছু সময় নিতে পারে, তাই গিয়ে এক কাপ চা পান করুন ...

4. $p$poweRlaw

data_alt <- dislnorm$new(data)
data_alt$xmin <- est$xmin
data_alt$pars <- estimate_pars(data_alt)
comp <- compare_distributions(data_pl, data_alt)

$x_{min}$compare_distributions$x_{min}$compcomp$test_statisticdata_plcomp$p_two_side

সঙ্গে এই পদক্ষেপ পুনরাবৃত্তি করুন disexp, dispoisঐ বিকল্পসহ ক্ষমতা আইন তুলনা ক্লাস।

Comments

দয়া করে @ ক্লিভেল, আপনি কি কিছু উপায়ে আপনার ডেটা উদাহরণের অন্তর্ভুক্ত করতে পারেন? আমার ডেটা পরীক্ষা করতে আমার সমস্যা আছে তাই আমি একটি চলমান উদাহরণকে প্রশংসা করব ... ধন্যবাদ, আমি আপনার পোস্টটি খুব অনুপ্রেরণামূলক পেয়েছি!

— maycca

আপনি কমান্ডটি টাইপ করতে data("moby")পারেন এবং আমার উদাহরণের mobyপরিবর্তে ব্যবহার করতে পারেন data, যেমন পরিচয়

— জেব্রা প্রপালশন ল্যাব