nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Fit YbyX Non Parametric Quantile Density Contours

Byron_JMP

In Fit Y by X, the Non Parametric Quantile Density Contours are super useful for visualizing high density data.

In HLS people use software like FlowJo to make simple bivariate plots that have the data points colored by density. JMP gets close enough, plus there is an option to select rows based on a density range.

It would be super convenient if there was a slider for interactively selecting the range, rather than guessing at ranges until the most appropriate one is discovered iteratively.

Screen Shot 2022-03-25 at 12.39.38 PM.png

Status: In Review

Comments

Sarah-Sylvestre

@Byron_JMP -Thank you for your suggestion! We have captured your request and will take it under consideration.

SamGardner

Thanks @Byron_JMP . We will look into this.

Here is a prototype script that will do this, I'm sure this could be adapted into an add-in.

Names Default To Here( 1 );
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
obj = dt << Bivariate( Y( :Weight ), X( :Height ) );
robj=report(obj);
obj << Nonpar Density( {"5% Contours"n( 0 )} );
obj << (curve[1] << "5% Contours"n( 1 ));

sliderLowerValue=0;
sliderUpperValue=1;
robj["Quantile Density Contours"]<< Append(
Panel Box( "Select Points by Non-Parametric Density",
		tb1 = Text Box( "Lower Density: " || Char( sliderLowerValue ) ),
		tb2 = Text Box( "Upper Density: " || Char( sliderUpperValue ) ),
		sb = Range Slider Box(
			0,
			1,
			sliderLowerValue,
			sliderUpperValue,
			tb1 << Set Text( "Lower Density: " || Char( sliderLowerValue ) );
			tb2 << Set Text( "Lower Density: " || Char( sliderUpperValue ) );
			obj << (curve[1] << Select Points by Density( sliderLowerValue, sliderUpperValue ));
		)
	)
	);

Byron_JMP

//exactly what I had in mind
//works really well with high density data 

Names Default To Here( 1 );
dt = Open("$SAMPLE_DATA/Cytometry.jmp");
obj = dt << Bivariate( Y( :CD3 ), X( :CD8 ) );
robj=report(obj);
obj << Nonpar Density( {"5% Contours"n( 0 )} );
obj << (curve[1] << "5% Contours"n( 1 ));

sliderLowerValue=0;
sliderUpperValue=1;
robj["Quantile Density Contours"]<< Append(
Panel Box( "Select Points by Non-Parametric Density",
		tb1 = Text Box( "Lower Density: " || Char( sliderLowerValue ) ),
		tb2 = Text Box( "Upper Density: " || Char( sliderUpperValue ) ),
		sb = Range Slider Box(
			0,
			1,
			sliderLowerValue,
			sliderUpperValue,
			tb1 << Set Text( "Lower Density: " || Char( sliderLowerValue ) );
			tb2 << Set Text( "Lower Density: " || Char( sliderUpperValue ) );
			obj << (curve[1] << Select Points by Density( sliderLowerValue, sliderUpperValue ));
		)
	)
	);

mia_stephens

Thanks for the suggestion @Byron_JMP, and for developing the script @SamGardner! The developer has agreed to work on this in JMP 18.