src/hg/instinct/hgGeneset/hgGenesets.c 1.6
1.6 2010/01/28 23:05:25 jsanborn
added option to grab first N genes in dataset
Index: src/hg/instinct/hgGeneset/hgGenesets.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/instinct/hgGeneset/hgGenesets.c,v
retrieving revision 1.5
retrieving revision 1.6
diff -b -B -U 4 -r1.5 -r1.6
--- src/hg/instinct/hgGeneset/hgGenesets.c 28 Jan 2010 22:59:07 -0000 1.5
+++ src/hg/instinct/hgGeneset/hgGenesets.c 28 Jan 2010 23:05:25 -0000 1.6
@@ -856,8 +856,10 @@
char *sampleIds = cartOptionalString(cart, hghSampleIds);
char *featureIds = cartOptionalString(cart, hghFeatureIds);
int dataset_id = cartUsualInt(cart, hghDatasetId, -1);
+int getFirst = cartUsualInt(cart, hghGetFirst, -1);
+
char *metric = cartOptionalString(cart, hghClusterMetric);
char *method = cartOptionalString(cart, hghClusterMethod);
struct datasets *da = getDatasetById(conn, dataset_id);
@@ -869,17 +871,19 @@
sampleIds = getAllIdsInDataset(conn, da, "sample_id");
//errAbort("%s is required\n", hghSampleIds);
}
+if (!featureIds && getFirst > 0)
+ featureIds = getNumAnalysisFeatureIdsInDataset(conn, da, getFirst);
+
if (!featureIds)
errAbort("%s is required\n", hghFeatureIds);
struct samples *samples = getSamplesByIds(conn, sampleIds);
if (!samples)
errAbort("No samples matching ids = %s\n", sampleIds);
-// To test LOTS of features
-//featureIds = getNumAnalysisFeatureIdsInDataset(conn, da, 1000);
+
struct analysisFeatures *afList = getAnalysisFeaturesByIds(conn, featureIds);
if (!afList)
errAbort("No features matching ids = %s\n", featureIds);