cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
  • Due to inclement weather, JMP support response times may be slower than usual during the week of January 26.
    To submit a request for support, please send email to support@jmp.com.
    We appreciate your patience at this time.
  • Register to see how to import and prepare Excel data on Jan. 30 from 2 to 3 p.m. ET.

JMP Wish List

We want to hear your ideas for improving JMP. Share them here.
Choose Language Hide Translation Bar
0 Kudos

Box Tidwell for logistic regression

Hi JMP community,

One of the assumptions of logistic regression is linearity of variables and log odds. However, I still have not figured out how to test this in JMP. The specific test I have seen elsewhere is the Box Tidwell test. Is there a way to get this test on JMP? Failing that, is there another way to test linearity of a logistic regression model?

Thanks!

3 Comments
Status changed to: Acknowledged

Hi @aande2, thank you for your suggestion! We have captured your request and will take it under consideration.

SamGardner
Level VII
Status changed to: Investigating
 
SamGardner
Level VII
Status changed to: Not Planned For Now

Hello @aande2 after some research we have decided to not pursue this.  You can do this in JMP manually by adding in the interactions or each continuous term with its log transform. That is, if Logit(Y) = aX, then fit the model Logit(Y) = a X + b X*Ln(X).  If the parameter estimate for the X*Ln(X) effect is significant, then Logit(Y) is deemed not linearly dependent on X. 

 

Run this script to see how to do this:

 

dt=Open( "$SAMPLE_DATA/Ingots.jmp" );

// Logistic regression model
dt << Fit Model(
	Freq( :count ),
	Y( :ready ),
	Effects( :heat, :soak ),
	Personality( "Nominal Logistic" ),
	Run( Likelihood Ratio Tests( 1 ), Wald Tests( 0 ) )
);

// model with Box Tidwell effects
dt <<Fit Model(
	Y( :ready ),
	Effects( :heat, :soak, :heat * :"Log[heat]"n, :soak * :"Log[soak]"n ),
	Target Level( "Ready" ),
	Personality( "Nominal Logistic" ),
	Run( Likelihood Ratio Tests( 1 ), Wald Tests( 0 ) )
);