pythonalgorithmrgraphnetworkx# Generating a graph with certain degree distribution?

I am trying to generate a random graph that has small-world properties (exhibits a power law distribution). I just started using the networkx package and discovered that it offers a variety of random graph generation. Can someone tell me if it possible to generate a graph where a given node's degree follows a gamma distribution (either in R or using python's networkx package)?

Solution

If you want to use the configuration model something like this should work in NetworkX:

```
import random
import networkx as nx
z=[int(random.gammavariate(alpha=9.0,beta=2.0)) for i in range(100)]
G=nx.configuration_model(z)
```

You might need to adjust the mean of the sequence z depending on parameters in the gamma distribution. Also z doesn't need to be graphical (you'll get a multigraph), but it does need an even sum so you might have to try a few random sequences (or add 1)...

The NetworkX documentation notes for configuration_model give another example, a reference and how to remove parallel edges and self loops:

```
Notes
-----
As described by Newman [1]_.
A non-graphical degree sequence (not realizable by some simple
graph) is allowed since this function returns graphs with self
loops and parallel edges. An exception is raised if the degree
sequence does not have an even sum.
This configuration model construction process can lead to
duplicate edges and loops. You can remove the self-loops and
parallel edges (see below) which will likely result in a graph
that doesn't have the exact degree sequence specified. This
"finite-size effect" decreases as the size of the graph increases.
References
----------
.. [1] M.E.J. Newman, "The structure and function
of complex networks", SIAM REVIEW 45-2, pp 167-256, 2003.
Examples
--------
>>> from networkx.utils import powerlaw_sequence
>>> z=nx.create_degree_sequence(100,powerlaw_sequence)
>>> G=nx.configuration_model(z)
To remove parallel edges:
>>> G=nx.Graph(G)
To remove self loops:
>>> G.remove_edges_from(G.selfloop_edges())
```

Here is an example similar to the one at http://networkx.lanl.gov/examples/drawing/degree_histogram.html that makes a drawing including a graph layout of the largest connected component:

```
#!/usr/bin/env python
import random
import matplotlib.pyplot as plt
import networkx as nx
def seq(n):
return [random.gammavariate(alpha=2.0,beta=1.0) for i in range(100)]
z=nx.create_degree_sequence(100,seq)
nx.is_valid_degree_sequence(z)
G=nx.configuration_model(z) # configuration model
degree_sequence=sorted(nx.degree(G).values(),reverse=True) # degree sequence
print "Degree sequence", degree_sequence
dmax=max(degree_sequence)
plt.hist(degree_sequence,bins=dmax)
plt.title("Degree histogram")
plt.ylabel("count")
plt.xlabel("degree")
# draw graph in inset
plt.axes([0.45,0.45,0.45,0.45])
Gcc=nx.connected_component_subgraphs(G)[0]
pos=nx.spring_layout(Gcc)
plt.axis('off')
nx.draw_networkx_nodes(Gcc,pos,node_size=20)
nx.draw_networkx_edges(Gcc,pos,alpha=0.4)
plt.savefig("degree_histogram.png")
plt.show()
```

- Pandas: new column, sampled from list, based on column value
- How to correctly install rpy2 in python?
- How do I change the dataset when I Fine-Tuning the Whisper model?
- Customizing legend in Seaborn histplot subplots
- python issues with KeyError: skip over errors when data is not in the json?
- Print string to text file
- Explanation of reportAttributeAccessIssue of pylance
- Why is a test (100/100) leading to an unexpected output?
- Arrays Merging Instead of Appending in Python
- Why did I got an error ModuleNotFoundError: No module named 'distutils'?
- Python / Boto3: Float types are not supported. Use Decimal types instead
- Why is this regex not capturing the last input?
- How to fully disassemble Python source
- Following Links in HTML using BeautifulSoup
- Floating-point efficiency of Python library function implementations
- Convert set to a hashable type
- Why is my python code so slow (leetcode)?
- Getting individual colors from a color map in matplotlib
- How do self-assessing sites like Hackerrank and Leetcode check whether an answer is correct?
- Combine list of dictionaries in python
- Numpy, from a list of points, extract ones that have the same range of coordinates without for loop
- PostgreSQL - how to run VACUUM from code outside transaction block?
- Why is heap slower than sort for K Closest Points to Origin?
- How to get specific unique combinations of a dataframe using only dataframe operations?
- Make this data cleaning loop more Python
- How to set the working directory for debugging a Python program in VS Code?
- How to plot multiple slices of a 3D brain image with nilearn and matplotlib subplots
- Using nibabel to save nifti in "SPM" style
- I want to apply an algorithm to different images from a folder and save the new ones in another folder
- installing Brain Scaffold Builder with broken rtree dependency