Returns the centroid and dispersion of a genes or list of genes specified by argument 'genin'. When 'ignore_log' is True, it treat expression values as being in natural scale, even if self.log2 is True (used internally).
(self, genin, ignore_log=False)
| 1711 | return values |
| 1712 | |
| 1713 | def centroid(self, genin, ignore_log=False): |
| 1714 | """ |
| 1715 | Returns the centroid and dispersion of a genes or list of genes specified by argument 'genin'. |
| 1716 | When 'ignore_log' is True, it treat expression values as being in natural scale, even if self.log2 is |
| 1717 | True (used internally). |
| 1718 | """ |
| 1719 | dicge = self.get_gene(genin, ignore_log=ignore_log)[0] |
| 1720 | pel1 = 0.0 |
| 1721 | pel2 = 0.0 |
| 1722 | pel3 = 0.0 |
| 1723 | for node in self.pl: |
| 1724 | pel1 += self.dicdis[node]*dicge[node] |
| 1725 | pel2 += dicge[node] |
| 1726 | if pel2 > 0.0: |
| 1727 | cen = float(pel1)/float(pel2) |
| 1728 | for node in self.pl: |
| 1729 | pel3 += numpy.power(self.dicdis[node]-cen, 2)*dicge[node] |
| 1730 | return [(cen-self.po[1])/self.po[0], (numpy.sqrt(pel3/float(pel2))-self.po[1])/self.po[0]] |
| 1731 | else: |
| 1732 | return [None, None] |
| 1733 | |
| 1734 | def get_gene(self, genin, ignore_log=False, con=True): |
| 1735 | """ |