/* * The MIT License (MIT) * * Copyright (c) 2007-2015 Broad Institute * * Permission is hereby granted, free of charge, to any person obtaining a copy * of this software and associated documentation files (the "Software"), to deal * in the Software without restriction, including without limitation the rights * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell * copies of the Software, and to permit persons to whom the Software is * furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included in * all copies or substantial portions of the Software. * * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN * THE SOFTWARE. */ /* * To change this template, choose Tools | Templates * and open the template in the editor. */ package org.broad.igv.data; import org.broad.igv.AbstractHeadlessTest; import org.broad.igv.util.ResourceLocator; import org.broad.igv.util.TestUtils; import org.junit.Test; import java.util.ArrayList; import java.util.List; import java.util.Set; import java.util.TreeSet; import static org.junit.Assert.assertEquals; import static org.junit.Assert.assertTrue; /** * @author jrobinso */ public class IGVDatasetTest extends AbstractHeadlessTest { String cnFile = TestUtils.DATA_DIR + "igv/MIP_44.cn"; /** * Test loading a dataset and accessing some random data cells */ @Test public void testDataset() { IGVDataset ds = new IGVDataset(new ResourceLocator(cnFile), genome); // Get the start locations and data from sample yw280-4_44 on chr 1 int[] startLocations = ds.getStartLocations("chr1"); float[] data = ds.getData("yw280-4_44", "chr1"); assertEquals(startLocations.length, data.length); // Check the first data point assertEquals(794055, startLocations[0]); assertEquals(1.86971, data[0], 1.0e-4); // Check the last data point int end = startLocations.length - 1; assertEquals(245171540, startLocations[end]); assertEquals(3.7871208, data[end], 1.0e-4); } /** * Test of the parser getHeadings method */ @Test public void testGetHeadings() { String firstHeading = "yw280-4_44"; String secondHeading = "yw280-4_48"; String lastHeading = "DN30_3_03"; String headingsLine = "SNP Chromosome PhysicalPosition yw280-4_44 yw280-4_48 yw280-5_40 yw280-5_41 yw280-5_42 yw280-5_43 yw280-5_44 yw280-5_46 yw280-5_47 yw280-5_48 yw280-5_45 DN30_01 DN30_02 DN30_03 DN30_04 DN30_06 DN30_07 DN30_09 DN30_10 DN30_11 DN30_13 DN30_14 DN30_15 DN30_16 DN30_17 DN30_18 DN30_19 DN30_20 DN30_21 DN30_32 DN30_33 DN30_34 DN30_35 DN30_36 DN30_41 DN30_42 DN30_43 DN30_44 DN30_45 DN30_46 DN30_47 DN30_48 DN30_3_02 DN30_3_03"; String[] tokens = headingsLine.split("\t"); IGVDataset ds = new IGVDataset(new ResourceLocator(cnFile), genome); IGVDatasetParser parser = new IGVDatasetParser(new ResourceLocator(cnFile), genome); //parser.scan necessary to set values needed in getHeadings List<ChromosomeSummary> summaries = parser.scan(ds); String[] headings = parser.getHeadings(tokens, 1); assertEquals(firstHeading, headings[0]); assertEquals(secondHeading, headings[1]); assertEquals(lastHeading, headings[headings.length - 1]); } /** * Test of scan method. Test that we find * all chromosomes we expect */ @Test public void testScanDataset() { IGVDataset ds = new IGVDataset(new ResourceLocator(cnFile), genome); IGVDatasetParser parser = new IGVDatasetParser(new ResourceLocator(cnFile), genome); List<ChromosomeSummary> summaries = parser.scan(ds); assertEquals(24, summaries.size()); Set<String> chromos = new TreeSet<String>(); ArrayList<String> chromos_L = new ArrayList<String>(24); for (int ii = 1; ii <= 22; ii++) { chromos_L.add("chr" + ii); } chromos_L.add("chrX"); chromos_L.add("chrY"); chromos.addAll(chromos_L); for (ChromosomeSummary cSum : summaries) { assertTrue(chromos.contains(cSum.getName())); } } }