Jak sprawdzić, czy dystrybucja jest zgodna z prawem mocy?

Mam dane o tym, ilu użytkowników publikuje ile pytań. Na przykład,

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

Oznacza to, że 2 użytkowników opublikowało 100 pytań, 9 użytkowników opublikowało 10 pytań i tak dalej. Jak więc ustalić, czy UserCount, QuestionCountrozkład jest zgodny z prawem mocy?

Znalazłem pakiet poweRlaw . Mogę jednak przekazać tylko jedną grupę liczb, aby dokonać oceny. (Przykładem podanym w tym pakiecie jest częstotliwość słów.) Jak więc korzystać z tego pakietu? Czy mam coś nie tak? Mam również dane o liczbie pytań każdego użytkownika, tj [100, 100, 10, 10, 10 ... ]. Jeśli przekażę te dane do paczki, co dostanę?

Według Clauset i wsp. W ten sposób testujesz ogon prawa mocy za pomocą poweRlawpakietu:

1. Skonstruuj obiekt dystrybucji prawa mocy. W takim przypadku dane są dyskretne, więc użyj dyskretnej wersji klasy

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

dwie ostatnie linie można przepisać jako jedną linię

data_pl$xmin <- est

Ponadto w tym momencie możesz zobaczyć statystyki KS:

est$KS

3. $p$

bs <- bootstrap_p(data_pl)
bs$p

Może to zająć trochę czasu, więc idź i napij się herbaty ...

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

Powtórz tę czynność z disexp, dispoiszajęciach porównać Prawo energetyczne z tych alternatyw.