Index: trunk/psModules/src/detrend/pmFringeStats.c
===================================================================
--- trunk/psModules/src/detrend/pmFringeStats.c	(revision 31067)
+++ trunk/psModules/src/detrend/pmFringeStats.c	(revision 31069)
@@ -1029,19 +1029,24 @@
     psVector *science_values   = psVectorAllocEmpty(4000,PS_TYPE_F32);
     psVector *science_errors   = psVectorAllocEmpty(4000,PS_TYPE_F32);
-
-    psStats *binStats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
+    psVector *science_counts   = psVectorAllocEmpty(4000,PS_TYPE_S32);
+
+    psStats *binStats = psStatsAlloc(PS_STAT_CLIPPED_MEAN | PS_STAT_CLIPPED_STDEV);
     for (int i = 0; i < 4000; i++) {
       psVector *bin = bins->data[i];
-      if (bin->n > 2) {
+      if (bin->n > 10) {
 	psStatsInit(binStats);
 
 	psVectorStats(binStats,bin,NULL,NULL,1);
 	
-	if (isfinite(binStats->robustStdev) &&
-	    isfinite(binStats->robustMedian) &&
-	    binStats->robustStdev > 0) {
+	if (isfinite(binStats->clippedStdev) &&
+	    isfinite(binStats->clippedMean) &&
+	    (binStats->clippedStdev > 0) &&
+	    (binStats->clippedNvalues > 10) &&
+	    (binStats->clippedNvalues > 0.5 * bin->n)
+	    ) {
 	  psVectorAppend(fringe_positions,-0.1 + i * 5e-5);
-	  psVectorAppend(science_values, binStats->robustMedian);
-	  psVectorAppend(science_errors, binStats->robustStdev);
+	  psVectorAppend(science_values, binStats->clippedMean);
+	  psVectorAppend(science_errors, binStats->clippedStdev);
+	  psVectorAppend(science_counts, bin->n);
 	}
       }
@@ -1052,8 +1057,9 @@
 
     for (int i = 0; i < fringe_positions->n; i++) {
-      psTrace("psModules.detrend",7,"FITDATA: %f %f %f\n",
+      psTrace("psModules.detrend",7,"FITDATA: %f %f %f %d\n",
 	      fringe_positions->data.F32[i],
 	      science_values->data.F32[i],
-	      science_errors->data.F32[i]);
+	      science_errors->data.F32[i],
+	      science_counts->data.S32[i]);
     }
 /*     // Begin switch from old outlier removal and fitting code. */
@@ -1071,4 +1077,5 @@
     psFree(science_values);
     psFree(science_errors);
+    psFree(science_counts);
     
     for (int i = 0; i <= poly->nX; i++) {