Danes sem se odločil pogledati, kako je za dano porazdelitev z ustrezno kumulativno funkcijo.
Naj bo dana porazdelitev (pk). Tedaj definiramo zaporedje
qk = ∑(i=k,∞) pi
če za pk velja potenčni zakon pk = c k-α velja tudi za (qk)
qk ~ ∫(k,∞) c x-α dx ~ C k-α + 1
Torej moramo dobiti na sliki dve (skoraj) premici (npr. kot na sf.pdf).
setwd("C:/Users/batagelj/Documents/manuscripts/monika/zb/R")
F <- read.table("keywords_simple.vec", skip=1)$V1
# F <- read.table("scaleFree.vec", skip=1)$V1
T <- table(F)
x <- as.numeric(names(T)); t <- as.vector(T)
n <- max(x); zx <- 1:n
X <- numeric(n)
X[x] <- t Z <- X[n:1]
z <- cumsum(Z)[n:1]
plot(zx,z,pch=16,cex=0.5,col='red',log="xy")
points(x,t,pch=16,cex=0.5)
Če stvar narišemo za keywords , je zaporedje p precej premično, zaporedje q pa močno odstopa od tega (slika keywords.pdf). Torej je brezlestvičnost zaporedja p le navidezna.