diff --git a/NEWS b/NEWS index b26be73..0a7616c 100644 --- a/NEWS +++ b/NEWS @@ -1,3 +1,6 @@ +1.7.1: +* fix array dimensionality consistency for summary statistics on bigWig files + 1.7.0: * adapted existing python interface to open bigWig files diff --git a/genomedata/_bigwig.py b/genomedata/_bigwig.py index 3bc335d..c3ce680 100644 --- a/genomedata/_bigwig.py +++ b/genomedata/_bigwig.py @@ -69,23 +69,23 @@ def __getitem__(self, name): # NB: One implicit trackname across chromosomes for stat retrieval @property def mins(self): - return array(self.bw_file_header['minVal']) + return array([self.bw_file_header['minVal']]) @property def maxs(self): - return array(self.bw_file_header['maxVal']) + return array([self.bw_file_header['maxVal']]) @property def sums(self): - return array(self.bw_file_header['sumData']) + return array([self.bw_file_header['sumData']]) @property def sums_squares(self): - return array(self.bw_file_header['sumSquared']) + return array([self.bw_file_header['sumSquared']]) @property def num_datapoints(self): - return array(self.bw_file_header['nBasesCovered']) + return array([self.bw_file_header['nBasesCovered']]) def tracknames_continuous(self): # Return filepath to bigWig as implicit trackname diff --git a/test/run_tests.py b/test/run_tests.py index 5942701..4c5883a 100755 --- a/test/run_tests.py +++ b/test/run_tests.py @@ -304,7 +304,7 @@ def test_interface(self): self.assertEqual(genome.maxs, [93473]) self.assertEqual(genome.sums, [32720078]) self.assertEqual(genome.sums_squares, [1280737190372]) - self.assertAlmostEqual(genome.vars, 414615.8372, places=4) + self.assertAlmostEqual(genome.vars[0], 414615.8372, places=4) # Test chromosome retrieval chr1 = genome["chr1"] # memoization