cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
Check out the JMP® Marketplace featured Capability Explorer add-in
Choose Language Hide Translation Bar
TRR21
Level III

script for creating a table of p values

I have a large table of 40 rows (4 groups of 10 samples) and 10000 columns of values—

1. Is there a way for me to quickly run t-tests/ANOVA comparing particular group against another group, for each of those columns and list the p-values in a table to quickly sort through to understand the ones that are significantly different? 

2. Also, is it possible to set the p-value is such test to correct for multiple correction?

 

Thanks

1 ACCEPTED SOLUTION

Accepted Solutions
txnelson
Super User

Re: script for creating a table of p values

Run the Response Screening Platform and it will do that for you

     Analyze==>Screening==>Response Screening

Jim

View solution in original post

11 REPLIES 11
txnelson
Super User

Re: script for creating a table of p values

Run the Response Screening Platform and it will do that for you

     Analyze==>Screening==>Response Screening

Jim
TRR21
Level III

Re: script for creating a table of p values

Thanksf or the fast response - wouldn't have guessed this is what it would do.

 

Is there a way to trigger this from a script, selecting only particular rows for this analysis, based on data in couple other columns (ie. some way to subset/filter), without having to make that table first?

txnelson
Super User

Re: script for creating a table of p values

You can use the ability to set the Row State for given rows to be Excluded from the subsequent analysis.

Is this what you are trying to do?

Jim
TRR21
Level III

Re: script for creating a table of p values

Thanks, excluding all the rows I DON'T want does get me where I need to go.
TRR21
Level III

Re: script for creating a table of p values

Is there a way to choose between parametric vs. non-parametric/rank test to be used in this? None of the options available in the dialog seem to apply. Thanks.
TRR21
Level III

Re: script for creating a table of p values

Is there a way to force the system to used a paired test to generate the p-vals? Thanks.

txnelson
Super User

Re: script for creating a table of p values

JMP 13 and JMP 14

     Analyze==>Specialized Modeling==>Matched Pairs

JMP 12 and earlier

     Analyze==>Matched Pairs

 

For future reference, this information can be found in the Statistics Index, searching on "paired"

     Help==>Statistics Index

Jim
TRR21
Level III

Re: script for creating a table of p values

Thanks - but am not sure if this is the one I want. What I want is more like, what I get if I use:
Analyze==>Fit Y by X; and from red triangle pull down, use "Matching column"
This gives me the pre-vs-post paired comparison I am looking for, but would like to run this for the 20K columns as mentioned above.
txnelson
Super User

Re: script for creating a table of p values

Here is a sample script that might be close to what you want.  I had to guess a bit at what matching you were going to do.  But, it creates a table of output for all 127 columns in the data table.

Names Default To Here( 1 );

// Create a data table somewhat similar to your description
dtorig = Open( "$SAMPLE_DATA\semiconductor capability.jmp" );

dt = dtorig << Subset( Sample Size( 10 ), Selected columns only( 0 ), Stratify( :SITE ) );

Close( dtorig, nosave );

// Get all of the tests
colNamesList = dt << get column names( continuous );

// In my example, I am using the first column as the matching column
// so remove it from the list
colNamesList = Remove( colNamesList, 1, 1 );

// Run the Fit Y by X(Oneway)
fg = Fit Group(
	Oneway(
		Y( Eval( colNamesList ) ),
		X( :SITE ),
		Matching Column( :NPN1 ),
		X Axis Proportional( 0 ),
		Matching Lines( 1 ),
		SendToReport(
			Dispatch( {}, "Oneway Plot", FrameBox, {Grid Line Order( 2 ), Reference Line Order( 3 )} ),
			Dispatch( {}, "Matching Fit", OutlineBox, {Close( 0 )} )
		)
	)
);

// Create the output table
Report( fg )[Table Box( 3 )] << make combined data table;

// Close the report output
Report( fg ) << close window;
Jim