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