tests/test_final_deg.R

library("mcPAFit")


alpha <- 1
dat <- simple_net(time_step = 10000, alpha = alpha, p_max = 0.1, p_min = 0.01, mode_p = 1)

result <- empirical_estimate(dat$deg_vec)
plot(result$k + 1, result$A, log = "xy")
lines(result$k + 1,(result$k + 1)^alpha)
polygon(c(result$k + 1,rev(result$k + 1)) , 
        c(result$lower_A[as.character(result$k)],
          rev(result$upper_A[as.character(result$k)])), 
        col = rgb(1,0,0,0.1),border = NA)

#prepare the custom_PA function
custom_PA <- rep(0,max(result$k) + 1)
custom_PA[result$k + 1] <- result$A

zero_entry <- which(custom_PA == 0)
for (i in 1:length(zero_entry)) {
   if (zero_entry[i] != zero_entry[i + 1] - 1) 
      custom_PA   
}
# 
thongphamthe/mcPAFit documentation built on May 20, 2019, 10:23 p.m.