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])
