cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
Try the Materials Informatics Toolkit, which is designed to easily handle SMILES data. This and other helpful add-ins are available in the JMP® Marketplace
Choose Language Hide Translation Bar
vincem
Level II

Scaled Logitic Model

Hi- need to fit the following logistic model to data, where Y=binary (1=diseased, 0=no-disease) and X is continuous. The model is similar to a logistic but with an added parameter, lambda. Any advice on how to proceed in JMP will be appreciated.

 

vincem_0-1722487151613.png

vincem_1-1722487390178.png

 

4 REPLIES 4
MRB3855
Super User

Re: Scaled Logitic Model

Hi @vincem : That is the Logistic 3p model in the Nonlinear platform Model Library, where lambda (your model) = theta1, alpha ( your model) = ln(theta2), beta (your model) = theta3.

MRB3855_0-1722517255371.png

 

vincem
Level II

Re: Scaled Logitic Model

thank you very much. however my response is biinary (diseased =1, or not diseased=0) im unsure whether the non-linear platfrom can handle a binary response. appreciate your help. thanks again.

MRB3855
Super User

Re: Scaled Logitic Model

Hi @vincem  It can be done via the method described by @peng_liu  here. The log-likelihood will be different though. See here.

https://community.jmp.com/t5/Discussions/fitting-truncated-data-to-a-continuous-distribution/m-p/772...

And see section 12.2.1 here (you will need to adapt it slightly to include lambda).

https://www.stat.cmu.edu/~cshalizi/uADA/12/lectures/ch12.pdf

for the log-likelihood function (use the negative log-likelihood as @peng_liu explains)

Questions, come back.

vincem
Level II

Re: Scaled Logitic Model

Thanks again..I have fitted  the nonlinear model and fitted the below, where prop=n/N in below table. Parameter estimates noted below. i think its worked as its should..

 

vincem_0-1722767229388.png

 

vincem_3-1722767576833.png

 

 

vincem_1-1722767312218.png

Data is from the publication

https://pubmed.ncbi.nlm.nih.gov/21983214/

 

vincem_2-1722767366690.png

 

 

 

 

titresvaccinenNprop
5Non-Flu61190.05042
20Non-Flu020
40Non-Flu010
80Non-Flu040
160Non-Flu070
320Non-Flu0130
453Non-Flu010
640Non-Flu040
2560Non-Flu010
3620Non-Flu010
5TIV Control2250.08
10TIV Control2360.055556
20TIV Control4470.085106
28TIV Control010
40TIV Control2420.047619
57TIV Control010
80TIV Control4280.142857
113TIV Control010
160TIV Control0320
226TIV Control010
320TIV Control0170
640TIV Control0160
905TIV Control020
1280TIV Control0350
1810TIV Control040
2560TIV Control0190
3620TIV Control020
5120TIV Control040
5TIV adj040
40TIV adj010
80TIV adj010
160TIV adj0210
226TIV adj020
320TIV adj0630
453TIV adj070
640TIV adj1760.013158
905TIV adj040
1280TIV adj1630.015873
1810TIV adj040
2560TIV adj0460
3620TIV adj090
5120TIV adj0100