notebook.community

Edit and run

Construct an ER graph





In [1]:



    


using Gadfly, Graphons



    


















    






Warning: New definition 
    *(S<:Union(GenericQuadExpr{Float64,Variable},GenericAffExpr{Float64,Variable},Variable),T<:GenericQuadExpr{Float64,Variable}) at /Users/madeleine/.julia/v0.3/JuMP/src/operators.jl:266
is ambiguous with: 
    *(T<:GenericQuadExpr{Float64,Variable},S<:Union(GenericQuadExpr{Float64,Variable},GenericAffExpr{Float64,Variable},Variable)) at /Users/madeleine/.julia/v0.3/JuMP/src/operators.jl:264.
To fix, define 
    *(T<:GenericQuadExpr{Float64,Variable},T<:GenericQuadExpr{Float64,Variable})
before the new definition.

Sample from SBM





In [2]:



    


k = 5
A = SBM(fill(1/k,k))
sample(A, 5)



    


















    
Out[2]:








5x5 Array{Int64,2}:
 0  0  0  0  1
 0  0  1  0  0
 0  1  0  0  0
 0  0  0  0  1
 1  0  0  1  0

Compute graph statistics





In [3]:



    


erdos = SBM([1], fill(.3,1,1))
edge_density(erdos,100)
plot(x=degree_distribution(erdos,1000), Geom.histogram)



    


















    
Out[3]:









  
    x
  
  
    0.05
    0.10
    0.15
    0.20
    0.25
    0.30
    0.35
    0.40
    0.45
    0.50
    0.55
    0.60
    0.00
    0.01
    0.02
    0.03
    0.04
    0.05
    0.06
    0.07
    0.08
    0.09
    0.10
    0.11
    0.12
    0.13
    0.14
    0.15
    0.16
    0.17
    0.18
    0.19
    0.20
    0.21
    0.22
    0.23
    0.24
    0.25
    0.26
    0.27
    0.28
    0.29
    0.30
    0.31
    0.32
    0.33
    0.34
    0.35
    0.36
    0.37
    0.38
    0.39
    0.40
    0.41
    0.42
    0.43
    0.44
    0.45
    0.46
    0.47
    0.48
    0.49
    0.50
    0.51
    0.52
    0.53
    0.54
    0.55
    0.0
    0.2
    0.4
    0.6
    0.09
    0.10
    0.11
    0.12
    0.13
    0.14
    0.15
    0.16
    0.17
    0.18
    0.19
    0.20
    0.21
    0.22
    0.23
    0.24
    0.25
    0.26
    0.27
    0.28
    0.29
    0.30
    0.31
    0.32
    0.33
    0.34
    0.35
    0.36
    0.37
    0.38
    0.39
    0.40
    0.41
    0.42
    0.43
    0.44
    0.45
    0.46
    0.47
    0.48
    0.49
    0.50
    0.51
    0.52
    0.53
    0.54
    0.55
    0.56
  
  
    
      
    
    
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
    
    
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
    
    
      
        
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
        
      
    
    
      
        
        
      
      
        
      
      
        
      
      
        
        
      
    
  
  
    -50
    -40
    -30
    -20
    -10
    0
    10
    20
    30
    40
    50
    60
    70
    80
    90
    -40
    -38
    -36
    -34
    -32
    -30
    -28
    -26
    -24
    -22
    -20
    -18
    -16
    -14
    -12
    -10
    -8
    -6
    -4
    -2
    0
    2
    4
    6
    8
    10
    12
    14
    16
    18
    20
    22
    24
    26
    28
    30
    32
    34
    36
    38
    40
    42
    44
    46
    48
    50
    52
    54
    56
    58
    60
    62
    64
    66
    68
    70
    72
    74
    76
    78
    80
    -50
    0
    50
    100
    -40
    -35
    -30
    -25
    -20
    -15
    -10
    -5
    0
    5
    10
    15
    20
    25
    30
    35
    40
    45
    50
    55
    60
    65
    70
    75
    80
  



  





















In [4]:



    


k = 3
sbm2 = SBM(fill(1/k,k))
plot(x=degree_distribution(sbm2,1000), Geom.histogram)



    


















    
Out[4]:









  
    x
  
  
    -0.4
    -0.3
    -0.2
    -0.1
    0.0
    0.1
    0.2
    0.3
    0.4
    0.5
    0.6
    0.7
    0.8
    0.9
    1.0
    1.1
    1.2
    1.3
    -0.30
    -0.28
    -0.26
    -0.24
    -0.22
    -0.20
    -0.18
    -0.16
    -0.14
    -0.12
    -0.10
    -0.08
    -0.06
    -0.04
    -0.02
    0.00
    0.02
    0.04
    0.06
    0.08
    0.10
    0.12
    0.14
    0.16
    0.18
    0.20
    0.22
    0.24
    0.26
    0.28
    0.30
    0.32
    0.34
    0.36
    0.38
    0.40
    0.42
    0.44
    0.46
    0.48
    0.50
    0.52
    0.54
    0.56
    0.58
    0.60
    0.62
    0.64
    0.66
    0.68
    0.70
    0.72
    0.74
    0.76
    0.78
    0.80
    0.82
    0.84
    0.86
    0.88
    0.90
    0.92
    0.94
    0.96
    0.98
    1.00
    1.02
    1.04
    1.06
    1.08
    1.10
    1.12
    1.14
    1.16
    1.18
    1.20
    1.22
    -0.5
    0.0
    0.5
    1.0
    1.5
    -0.30
    -0.25
    -0.20
    -0.15
    -0.10
    -0.05
    0.00
    0.05
    0.10
    0.15
    0.20
    0.25
    0.30
    0.35
    0.40
    0.45
    0.50
    0.55
    0.60
    0.65
    0.70
    0.75
    0.80
    0.85
    0.90
    0.95
    1.00
    1.05
    1.10
    1.15
    1.20
    1.25
  
  
    
      
    
    
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
    
    
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
    
    
      
        
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
          
        
      
    
    
      
        
        
      
      
        
      
      
        
      
      
        
        
      
    
  
  
    -50
    -40
    -30
    -20
    -10
    0
    10
    20
    30
    40
    50
    60
    70
    80
    90
    -40
    -38
    -36
    -34
    -32
    -30
    -28
    -26
    -24
    -22
    -20
    -18
    -16
    -14
    -12
    -10
    -8
    -6
    -4
    -2
    0
    2
    4
    6
    8
    10
    12
    14
    16
    18
    20
    22
    24
    26
    28
    30
    32
    34
    36
    38
    40
    42
    44
    46
    48
    50
    52
    54
    56
    58
    60
    62
    64
    66
    68
    70
    72
    74
    76
    78
    80
    -50
    0
    50
    100
    -40
    -35
    -30
    -25
    -20
    -15
    -10
    -5
    0
    5
    10
    15
    20
    25
    30
    35
    40
    45
    50
    55
    60
    65
    70
    75
    80

Content source: madeleineudell/Graphons.jl

Similar notebooks:

notebook.community | gallery | about