void average_rows()
{
	// Bring up file dialog and get user to select all OPJs to be processed
	int iNumFiles;
	StringArray saFilePaths;
	iNumFiles = GetMultiOpenBox( saFilePaths, "*.OPJ");

	// If user selected less than 2 files, quit
	if( iNumFiles < 2 )
	{
		out_str("Less than 2 OPJs selected!");
		return;
	}

	// Declare vector to hold data and also counts
	// Initialize size to be large enough and set to zero
	vector vecEData(100);
	vector<int> vecECount(100);
	vecEData = 0;
	vecECount = 0;

	// Loop over all files
	for(int ii = 0; ii < iNumFiles; ii++)
	{
		// If successful in opening project...
		if( Project.Open(saFilePaths[ii]) )
		{
			printf("Processing: %s\n", saFilePaths[ii]);
			// Declare dataset for data1 wks col E
			// Change this as desired to point to your specific wks and col name
			Dataset dsE("data1_e");
			// If dataset is valid...
			if( dsE.IsValid() )
			{
				// Loop over all elements and add to data vector and count vector
				for(int jj = 0; jj < dsE.GetSize(); jj++)
				{
					vecEData[jj] += dsE[jj];
					vecECount[jj] += 1;
				}
			}
			else
				out_str("Data1_E column not found in current project!")
		}	
	}

	// Done with all files - output the results to script window
	// One could instead create a new OPJ and place results in a worksheet
	printf("Row #\tCount\tAverage\n");
	for(ii = 0; ii < 100; ii++)
	{
		if( 0 == vecECount[ii]) break;
		printf("%d\t%d\t%f\n", ii, vecECount[ii], vecEData[ii] / vecECount[ii]);
	}
}

#define MAX_ROWS 100

void average_rows(string strFolderPath)
{
	// Add a \ to the path if necessary
	if( 0 != strFolderPath.Right(1).Compare("\\") )
		strFolderPath += "\\";
		
	// Get list of all OPJ files in specified folder
	StringArray saOPJFiles;
	bool bRet = FindFiles(saOPJFiles, strFolderPath, "OPJ");
	if( !bRet )
	{
		out_str("Failed to find OPJ files in specified path!");
		return;
	}

	// Get file count
	int nFileCount = saOPJFiles.GetSize();
	if( nFileCount < 2 )
	{
		out_str("Less than 2 OPJs found!");
		return;
	}

	// Declare vector to hold X values
	vector vecX();
	// Declare vector to hold summed E column data and counts
	// Initialize size to be large enough and set to zero
	vector vecEData(MAX_ROWS);
	vector<int> vecECount(MAX_ROWS);
	vecEData = 0;
	vecECount = 0;

	// Loop over all files
	for(int ii = 0; ii < nFileCount; ii++)
	{
		string strOPJFileName = strFolderPath + saOPJFiles[ii];
		// If successful in opening project...
		if( Project.Open(strOPJFileName) )
		{
			printf("Processing: %s\n", strOPJFileName);
			// Assume X data is in col A
			// and there is a col E
			Dataset dsA("data1_a");
			Dataset dsE("data1_e");
			if( dsA && dsE )
			{
				// If dsA is larger than current X vector, then update X vector
				if( dsA.GetSize() > vecX.GetSize() )
					vecX = dsA;
				// Loop over all elements of dsE and add to E data vector and count vector
				for(int jj = 0; jj < dsE.GetSize(); jj++)
				{
					vecEData[jj] += dsE[jj];
					vecECount[jj] += 1;
				}
			}
			else
				out_str("   Dataset Data1_A or Data1_E not found in this project!");
		}
		else
			printf("Failed to open %s\n", strOPJFileName);
	}

	// Done with all files
	
	// Create a new OPJ in same folder for storing result
	Project.Open();
	// Create a new worksheet and add cols for output
	WorksheetPage wpg;
	wpg.Create("Origin");
	Worksheet wks = wpg.Layers(0);
	while( wks.DeleteCol(0) );
	wks.AddCol("XValue");
	wks.Columns(0).SetType(OKDATAOBJ_DESIGNATION_X);
	wks.AddCol("AverageE");
	Dataset dsXValue(wks, 0);
	Dataset dsEAverage(wks, 1);
	for(ii = 0; ii < MAX_ROWS; ii++)
	{
		if( 0 == vecECount[ii]) break;
		dsXValue.Add(vecX[ii]);
		dsEAverage.Add( vecEData[ii] / vecECount[ii] );
	}	
	// Save this result OPJ
	string strResultOPJ = strFolderPath + "AverageResult.OPJ";
	if( !Project.Save(strResultOPJ) )
		out_str("Failed to save result OPJ!");
}

quote:

---

Curve cuvFit1("Data2_XValue[X]", "Data2_AverageE[Y}");
int nPlot = gl.AddPlot(cuvFit1, IDM_PLOT_FLOWVECTOR);

---

quote:

---

vecEData[jj] += dsE[jj];
vecECount[jj] += 1;

printf("%d\t%d\t%f\n", ii, vecECount[ii], vecEData[ii] / vecECount[ii]);

---

quote:

---

You can skip the missing values in your computations like this...

if( vecEData[ii]!=NANUM || vecECount[ii]!=NANUM )
{
// do math if values are not missing
}
else
{
// set result to some default value
}

Mike Buess
Origin WebRing Member

---

quote:

---

Simple case: only col B contains missing values...

Worksheet wks = Project.ActiveLayer();
Dataset ds1(wks,0);
Dataset ds2(wks,1);
Dataset ds3(wks,2);
ds3.SetSize(ds1.GetSize());
for(int i=0;i<ds1.GetSize();i++)
{
if( ds2[i]==NANUM )
{
ds3[i] = ds1[i];
}
else
{
ds3[i] = (ds1[i] + ds2[i]) / 2.0;
}
}

Mike Buess
Origin WebRing Member

Edited by - Mike Buess on 04/11/2007 10:52:30 AM

---

quote:

---

Statistics > Descriptive Statistics > Statistics on Rows.

---