def test_add_track(self): new_track_name, new_track_file = self.new_track # Open new track genome = Genome(self.gdfilepath, mode="r+") with genome: self.assertEqual(genome.num_tracks_continuous, 0) genome.add_track_continuous(new_track_name) # Load data for new track load_data(self.gdfilepath, new_track_name, test_filename(new_track_file), verbose=self.verbose) # Close data with new track close_data(self.gdfilepath, verbose=self.verbose) # Make sure addition was successful genome = Genome(self.gdfilepath) with genome: # Track ordering should now end with dnase self.assertEqual(genome.tracknames_continuous, [new_track_name]) # Given track ordering, check single track data retrieval self.assertArraysEqual(genome["chr1"][305:310, new_track_name], [-2.65300012, 0.37200001, 0.37200001, 0.37200001, 0.37099999])
def test_add_track(self): new_track_name, new_track_file = self.new_track # Open new track genome = Genome(self.gdfilepath, mode="r+") with warnings.catch_warnings(): warnings.simplefilter("ignore", GenomedataDirtyWarning) with genome: self.assertEqual(genome.num_tracks_continuous, 0) genome.add_track_continuous(new_track_name) # Load data for new track load_data(self.gdfilepath, new_track_name, test_filename(new_track_file), verbose=self.verbose) # Close data with new track close_data(self.gdfilepath, verbose=self.verbose) # Make sure addition was successful genome = Genome(self.gdfilepath) with genome: # Track ordering should now end with dnase self.assertEqual(genome.tracknames_continuous, [new_track_name]) # Given track ordering, check single track data retrieval self.assertArraysEqual( genome["chr1"][305:310, new_track_name], [-2.65300012, 0.37200001, 0.37200001, 0.37200001, 0.37099999])