Running ‘test_CV.R’
  Running ‘test_CV_2.R’ [9s/10s]
  Running ‘test_GetStatistics_compress_2.R’
  Running ‘test_GetStatistics_undirected.R’
  Running ‘test_PAFit.R’
  Running ‘test_PAFit_2.R’
  Running ‘test_overall.R’ [12s/14s]
Running the tests in ‘tests/test_PAFit.R’ failed.
Complete output:
  > #testing CV-related functions
  > library(PAFit)
  > prob_m <- "TRUE"
  > inc    <- "TRUE"
  > log    <-  c("TRUE")
  > mode_f_value <- c("Constant_PA", "Log_linear")[1]  
  > i <- 1
  > net  <- GenerateNet(N = 50, m = 10,prob_m = prob_m, increase = inc, log = log,
  +                     mode = i, shape = 100, rate = 100)
  > 
  > net  <- GenerateNet(N = 50, m = 10,prob_m = prob_m, increase = inc, log = log,
  +                     mode = i, shape = 100, rate = 100, specific_start = 10, 
  +                     mode_f = "log_normal")
  > 
  > net  <- GenerateNet(N = 50, m = 10,prob_m = prob_m, increase = inc, log = log,
  +                     mode = i, shape = 100, rate = 100, specific_start = 10, 
  +                     mode_f = "power_law")
  > 
  > net  <- GenerateNet(N = 50, m = 10,prob_m = prob_m, increase = inc, log = log,
  +                     mode = i, shape = 100, rate = 100, custom_PA = 1:1000,
  +                     specific_start = 10, 
  +                     mode_f = "power_law")
  > 
  > 
  > net_stats <- GetStatistics(net$graph,deg_threshold = 1, 
  +                              net_type = "directed",
  +                              Binning = TRUE, G = 10) 
  > 
  > net_stats <- GetStatistics(net$graph,deg_threshold = 0, 
  +                            net_type = "directed",
  +                            Binning = TRUE, G = 10, CompressMode = 1,
  +                            CompressRatio = 1, only_true_deg_matrix = TRUE)
  > 
  > net_stats <- GetStatistics(net$graph,deg_threshold = 0, 
  +                            net_type = "directed",
  +                            Binning = TRUE, G = 10, CompressMode = 1,
  +                            CompressRatio = 1)
  > 
  > 
  > net_stats <- GetStatistics(net$graph,deg_threshold = 0, 
  +                            net_type = "directed",
  +                            Binning = TRUE, G = 10, CompressMode = 1,
  +                            CompressRatio = 2)
  > 
  > for (q in 1:3)
  +  for (j in 0:2)
  +    for (uu in 0:1)   
  +       result <- PAFit(net_stats, r = 0.01, mode_reg_A = j, weight_PA_mode = uu,
  +                     stop_cond = 10^-2, q = q, normalized_f = TRUE,debug = TRUE)
  [1] -860.5583
  [1] -661.4044
  [1] -630.1815
  [1] -623.1953
  [1] -620.1298
  [1] -860.7449
  [1] -661.7311
  [1] -629.7881
  [1] -622.9338
  [1] -619.9441
  [1] -859.1164
  [1] -659.3852
  [1] -630.1926
  [1] -623.3654
  [1] -620.2491
  Error in upper_aa[ok_log] <- upper_aa[ok_log] + lambda * hessian_of_regularization_log(theta[non_zero_theta[ok_log]]) : 
    replacement has length zero
  Calls: PAFit
  Execution halted