def test_triadic_census(self): cens = triadic_census.triadic_census(self.random_graph) f = open("unit_test_data/random_triadic_census.pickle", "r") shouldBe = cPickle.load(f) f.close() self.assertEqual(cens, shouldBe)
def test_triadic_census(self): cens = triadic_census.triadic_census(self.random_graph) f = open("unit_test_data/random_triadic_census.pickle","r") shouldBe = cPickle.load(f) f.close() self.assertEqual(cens,shouldBe)
triadGraph.add_edge(nodes[id(neuron2)], nodes[id(neuron1)]) for innervation in network.innervations(): triadGraph.add_edge(nodes[id(innervation.neurite.rootObject())], nodes[id(innervation.muscle)]) for pathway in network.pathways(): if pathway.region1Projects: triadGraph.add_edge(nodes[id(pathway.region1)], nodes[id(pathway.region2)]) if pathway.region2Projects: triadGraph.add_edge(nodes[id(pathway.region2)], nodes[id(pathway.region1)]) for stimulus in network.stimuli(): triadGraph.add_edge(nodes[id(stimulus)], nodes[id(stimulus.target.rootObject())]) for synapse in network.synapses(): for postPartner in synapse.postSynapticPartners: triadGraph.add_edge(nodes[id(synapse.preSynapticNeurite.rootObject())], nodes[id(postPartner.rootObject())]) updateProgress(gettext('Finding motifs...')) counts = triadic_census.triadic_census(triadGraph) # Report the results in the console in case the user wants to copy/paste print('Triad motif occurences:') print('Motif\tCount') for triad in triadOrder: print(triad + '\t' + str(counts[triad])) # Open a new display that shows the results graphically network = createNetwork() display = displayNetwork(network) display.setDefaultFlowSpacing(0.15) edgeLength = 0.5 xDiff = edgeLength / 2.0 yDiff = sqrt(edgeLength ** 2 - xDiff ** 2) / 2.0 nodeSize = edgeLength / 5.0