/* * Copyright [1999-2015] Wellcome Trust Sanger Institute and the EMBL-European Bioinformatics Institute * Copyright [2016-2017] EMBL-European Bioinformatics Institute * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. */ /** * File: InterproFeatureTest.java * Created by: dstaines * Created on: May 6, 2009 * CVS: $$ */ package org.ensembl.healthcheck.testcase.eg_core; import java.sql.ResultSet; import java.sql.SQLException; import java.util.Map; import java.util.Map.Entry; import org.ensembl.healthcheck.DatabaseRegistryEntry; import org.ensembl.healthcheck.ReportManager; import org.ensembl.healthcheck.util.CollectionUtils; import org.ensembl.healthcheck.util.MapRowMapper; import org.ensembl.healthcheck.util.SqlTemplate; import org.ensembl.healthcheck.util.TemplateBuilder; import org.ensembl.healthcheck.util.TestCaseUtils; /** * Test to see if interpro xrefs and features agree * @author dstaines * */ public class InterproFeature extends AbstractEgCoreTestCase { private final static String GET_XREFS = "select t.translation_id,count(*) from translation t " + "join object_xref ox on (ox.ensembl_object_type='Translation' and " + "ox.ensembl_id=t.translation_id) join xref x using (xref_id) where x.external_db_id=1200 group by t.translation_id"; private final static String GET_FEATURES = "select t.translation_id,count(*) from translation t " + "join protein_feature pf using (translation_id) join analysis a using (analysis_id) " + " where a.logic_name in " + "($inlist$)" + " group by t.translation_id"; MapRowMapper<Integer, Integer> featureMapper = new MapRowMapper<Integer, Integer>() { public void existingObject(Integer currentValue, ResultSet resultSet, int position) throws SQLException { throw new UnsupportedOperationException(); } public Integer getKey(ResultSet resultSet) throws SQLException { return resultSet.getInt(1); } public Map<Integer, Integer> getMap() { return CollectionUtils.createHashMap(); } public Integer mapRow(ResultSet resultSet, int position) throws SQLException { return resultSet.getInt(2); } }; private final String feature_query; public InterproFeature() { super(); feature_query = TemplateBuilder.template(GET_FEATURES, "inlist", TestCaseUtils.resourceToInList("/org/ensembl/healthcheck/testcase/eg_core/interpro_names.txt")); } /* * (non-Javadoc) * * @see org.ensembl.healthcheck.testcase.SingleDatabaseTestCase#run(org.ensembl.healthcheck.DatabaseRegistryEntry) */ protected boolean runTest(DatabaseRegistryEntry dbre) { SqlTemplate srv = getTemplate(dbre); boolean passes = true; // 1. get features for each translation into a map Map<Integer, Integer> features = srv.queryForMap(feature_query, featureMapper); // 2. get xrefs into a map Map<Integer, Integer> xrefs = srv.queryForMap(GET_XREFS, featureMapper); for (Entry<Integer, Integer> e : xrefs.entrySet()) { if(features.containsKey(e.getKey())) { int nFeatures = features.get(e.getKey()); int nXrefs = e.getValue(); if (nFeatures < nXrefs) { passes = false; // note failure ReportManager.warning(this, dbre.getConnection(), nXrefs + " InterPro xrefs but only " + nFeatures + " interpro features found for translation " + e.getKey()); } } } if (passes) { ReportManager.correct(this, dbre.getConnection(), "Sufficient InterPro features found for translations with InterPro xrefs"); } return passes; } /* (non-Javadoc) * @see org.ensembl.healthcheck.testcase.eg_core.AbstractEgCoreTestCase#getEgDescription() */ @Override protected String getEgDescription() { return "Test to see if interpro xrefs and features agree"; } }