Newman | R Documentation |
This function implements a correction proposed in [1] of the original Newman's method in [2] to estimate the preferential attachment function.
Newman(net_object ,
net_stat = get_statistics(net_object),
start = 1 ,
interpolate = FALSE)
net_object |
an object of class |
net_stat |
An object of class |
start |
Positive integer. The starting time from which the method is applied. Default value is |
interpolate |
Logical. If |
Outputs an PA_result
object which contains the estimated attachment function. In particular, it contains the following field:
k
and A
: a degree vector and the estimated PA function.
center_k
and theta
: when we perform binning, these are the centers of the bins and the estimated PA values for those bins.
g
: the number of bins used.
alpha
and ci
: alpha
is the estimated attachment exponenet \alpha
(when assume A_k = k^\alpha
), while ci
is the mean plus/minus two-standard-deviation interval.
loglinear_fit
: this is the fitting result when we estimate \alpha
.
Thong Pham thongphamthe@gmail.com
1. Pham, T., Sheridan, P. & Shimodaira, H. (2015). PAFit: A Statistical Method for Measuring Preferential Attachment in Temporal Complex Networks. PLoS ONE 10(9): e0137796. (\Sexpr[results=rd]{tools:::Rd_expr_doi("10.1371/journal.pone.0137796")}).
2. Newman, M.. Clustering and preferential attachment in growing networks. Physical Review E. 2001;64(2):025102 (\Sexpr[results=rd]{tools:::Rd_expr_doi("10.1103/PhysRevE.64.025102")}).
See get_statistics
for how to create summerized statistics needed in this function.
See Jeong
, only_A_estimate
for other methods to estimate the attachment function in isolation.
library("PAFit")
net <- generate_net(N = 1000 , m = 1 , mode = 1 , alpha = 1 , s = 0)
net_stats <- get_statistics(net)
result <- Newman(net, net_stats)
summary(result)
# true function
true_A <- result$center_k
#plot the estimated attachment function
plot(result , net_stats)
lines(result$center_k, true_A, col = "red") # true line
legend("topleft" , legend = "True function" , col = "red" , lty = 1 , bty = "n")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.