How to make script run for each sample in a list if there are multiple entries per sample in data table.

J-Leonard

I have a data table where the same samples were measured multiple times (see example below). The samples are separated by their identifier.  

 

I want my script to be able to perform the same set of calculations on each sample so I can determine control limits and make a control chart for  A, B, C, etc. I do not know the identifier of each sample. I already have a script that can do my calculations and make the charts for me, it just cannot do separate calculations for each identifier. 

 

My initial through was to use some kind of loop to do this. I have a script that returns each unique identifier in a list and gives me the number of unique identifiers as a set up for the loop. 

Identifiers = Column( "Identifier") << GetAsMatrix();
Identifiers set = Associative Array( Identifiers);
Unique Identifiers = Identifiers set << Get Keys;

nUnique Identifiers = N Items (Unique Identifiers);

I just can't seem to get a loop to work though. Is there a better way to do this? Or could someone show me how to do the loop?

 

I am using JMP 10.

 

Thank you!

ih

If you do not want to modify your script you could run it against a subset of each level of the identifier, like this:

 

Names default to here(1);

//Open some sample data
dt = Open( "$Sample_data/big class.jmp" );

//Your code
Identifiers = Column( "age") << GetAsMatrix();
Identifiers set = Associative Array( Identifiers);
Unique Identifiers = Identifiers set << Get Keys;

nUnique Identifiers = N Items (Unique Identifiers);

//A window to hold a bunch of graphs
win = New Window("Graphs", graphbox = v list box());

//Run an analysis on a subset of the data table for each level of a variable
for( i=1, i<=nUnique Identifiers, i++,
	//Make a subset of the table with only values for this level
	dt << Select Where( :age == UniqueIdentifiers[i]);
	dtSub = dt << Subset("Selected Rows");
	
	//run your existing script on the subset
	//For demonstration, taking a picture of the graph builder platform
	// and getting the R2 value from a linear model and adding those to a window
	pic = (gb = dtSub << Graph Builder(
			Size( 534, 456 ),
			Show Control Panel( 0 ),
			Variables( X( :height ), Y( :weight ) ),
			Elements( Points( X, Y, Legend( 2 ) ) ),
			Local Data Filter(
				Close Outline( 1 ),
				Add Filter(
					columns( :age ),
					Where( :age == Unique Identifiers[i] ),
					Display( :age, N Items( 6 ) )
				)
			)
		)) << Get Picture;
	graphbox << Append( pic );
	gb << Close Window;
	
	fm = dtSub << Fit Model(
		Y( :height ),
		Effects( :weight ),
		Personality( "Standard Least Squares" ),
		Emphasis( "Effect Leverage" ),
		Run(
			:height << {Summary of Fit( 1 ), Analysis of Variance( 1 ),
			Parameter Estimates( 1 ), Scaled Estimates( 0 ),
			Plot Actual by Predicted( 1 ), Plot Residual by Predicted( 1 ),
			Plot Studentized Residuals( 0 ), Plot Effect Leverage( 1 ),
			Plot Residual by Normal Quantiles( 0 ), Box Cox Y Transformation( 0 )}
		)
	);
	Wait( 0 );
	dtSum = Report( fm )[Outline Box( "Response height" )][Outline Box( "Whole Model" )][
		Outline Box( "Summary of Fit" )][Table Box( 1 )] << Make Into Data Table;
	Report( fm ) << Close Window;
	Rsq = dtSum:Column 2[1];
	
	graphbox << Append( Text Box( "RSquared was " || char(Rsq) ) );
	
	dtSum << Close Window;
	
	//close the subset
	dtSub << Close Window;
	
);