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
}
#
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.