Download .zip (8kB), or download .tar.bz2 (6kB).
This code is mainly written in Cython.
apt-get install cython should be all you need.
It was originally written in plain Python, and could be converted back
fairly easily if necessary.
Many of the algorithms used scale horribly with the number of nodes in
the GDAG. One useful optimization would be to find a way to enumerate
non-isomorphic GDAGs directly, rather than generating and then removing
isomorphic ones as is currently done in
To list small irreducibly interesting GDAGs, as shown in Appendix E, edit the settings in list_interesting.py and then run it.
Each line of output corresponds to a GDAG. The first interesting one is of size 4:
(['', '', 'AB', 'AC'], 'BCD')
The bit in square brackets describes the underlying DAG:
The final part lists the observable nodes B, C, and D, which means that node A is unobservable.
To counts the total number of small GDAGs, and how many we can identify as having C = I as shown in Table 1, edit the settings in count_interesting.py and then run it.
The total number minus the boring number will generally be higher than the number output by list_interesting.py since no account is taken of reducibility.
Each line of output corresponds to a size n. The first line will be:
(1, 2, 2)
This means: considering GDAGs of size 1, there are 2 in total (a single node of either type), and of those 2 are boring (C = I).
Matthew F. Pusey, November 2014