if (FALSE) {
#setwd("tests")
rm(list = ls())
library(PAFit)
set.seed(1)
prob_m <- "FALSE"
inc <- "FALSE"
log <- c("FALSE")
M <- 5
alpha_vec <- rep(0,M)
s_vec <- rep(0,M)
alpha_optimal <- rep(0,M)
r_optimal <- rep(0,M)
for (i in 1:M) {
#set.seed(1)
net <- generate_net(N = 1000, m = 50, num_seed = 500,
multiple_node = 50, alpha = 1,
s = 5)
#shape = 1, rate = 1)
net_stats <- get_statistics(net)
print(result <- joint_estimate(net, net_stats))
#result
alpha_vec[i] <- result$estimate_result$alpha
s_vec[i] <- result$estimate_result$shape
alpha_optimal[i] <- result$cv_result$alpha_optimal
r_optimal[i] <- result$cv_result$r_optimal
}
print(alpha_vec)
print(mean(alpha_vec))
print(s_vec)
print(mean(s_vec))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.