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; } } }

Use of Logistic 3 Parameter (L3P) model

CMSB

Hi,

We're looking at some biological data representing an amplification which we're using JMP to fit the L3P model to
- and are thereby obtaining the growth rate, asymptote and inflection point on the sigmoid curves which result.

We're then using the growth rate parameter specifically in an ANOVA test to detect whether the curves (ie the biological test) are adversely affected by challenges, under the expectation of an adverse challenge resulting in us rejecting the null hypotheis.

The problem that I'm experiencing is that the more adverse the challenge - the greater the likelihood of either some or all of the assays not working (curves not being generated).

When this occurs, I cannot obtain the 3 parameters which I'm using in ANOVA and so am somewhat stuck on how to handle this particularly informative data.

The pattern (with increasing challenge) is of increasing variability in the curve that is modelled (towards flat lining) and then increasing numbers of failed amplifications (actual flat lines) with severity of challenge - which culminate in all replicates actually flat lining and no curve fit being possible.

Thanks for the help.

Accepted answers

Mark_Bailey

You could analyze such a table (challenge, number of failures, total number of tests - 16) using a probit model.

Select Analyze > Fit Model. Select both the Failures and Tests data columns and click Y. (Be sure to enter them in that order.) Select the Challenge data column and click Add. Click Standard Least Squares and select Generalized Linear Model. Select the Binomial distribution and the Probit link function.

Click Run.

You can extend this analysis to include the experiment.

Add Experiment as two effects: the main effect and crossed with Challenge.

The estimates can be used to test the significance of the effects and the model can be used to predict.

The Experiment example here could be another categorical or continuous factor. There could be multiple factors in the probit model.

All comments

Mark_Bailey

Your continuous factor is challenge. Perhaps your factor range is too aggressive and leads to non-responses. The point of the experiment is to model the data and test the parameter estimates. Can you shift the factor range lower? For example, if challenge was 50-100, try 25-75? That way you might be able to get more usable data for your model.

CMSB

Hiya,

Many thanks for the response and I think I understand.

The problem is that we are required to test the assay to failure and so our data ranges from perfectly superimposed replicate sigmoidal curves in all cases (with little challenge) at one end
- through (at the other end) - to perfectly superimposed replicate flat lines.

Between these two extremes there's a pattern of increasing sporadic failures in replicate samples.

I was thinking that if I grouped all of our replicates (by challenge) together into generating just one curve - then as long as we've at least one successful sigmoidal curve, that we'll be able to obtain the 3 parameters on each challenge ie will be able to perform a valid curve fit.

The problem though that then arises - is that instead of having values for the growth rate, inflection point and asymptote in each replicate - that we then have just 1 value for all replicates - and I'm not too sure how to take this one value through to t testing, (M)ANOVAs and equivalence testing because it's the actual replicate values themselves that're required to perform this type of analysis.

My next line of thinking was to use the summary data which JMP delivers on mean, SE and confidence intervals on that single 'summary' sigmoidal curve and to attempt to perform an ANOVA from Data Summary Table ie by simply lifting mean and SE out of the 'By group' L3P fit and attempting this procedure (which you've been kind enough to supply):
https://community.jmp.com/t5/JMP-Scripts/ANOVA-from-Data-Summary/ta-p/21487

What worries me a little is whether this form of grouping of 'flat line' data with actual sigmoidal curve data will invalidate some underlying assumption ie will change the nature of the curve. The inflexion point is a constant in our curve fitting exercise, and so I'm wondering whether this 'by group' approach will alter the inflexion point away from where we expect it to be.

Genuinely appreciate the help!

Thanks. C.

CMSB

So - the data we might observe might look like (illustrated with 16 replicates):

challenge units - successful curve fit - unsuccessful curve fit*

-------1------------------16---------------------00-----

-------2------------------15---------------------01-----

-------3------------------12---------------------04-----

-------4------------------08---------------------08-----

-------5------------------04---------------------12-----

-------6------------------02---------------------14-----

-------7------------------00---------------------16-----

However the next experiment might look like:

challenge units - successful curve fit - unsuccessful curve fit*

-------1------------------16---------------------00-----

-------2------------------16---------------------00-----

-------3------------------15---------------------01-----

-------4------------------16---------------------00-----

-------5------------------15---------------------01-----

-------6------------------15---------------------01-----

-------7------------------14---------------------02-----

* no reaction occurring

So - the first experiment feels as though it lends itself to a binary outcome, and the second to quantitative outcome (increased variability in sigmoid curves with challenge as assessed by asymptote/growth rate) - however, I somehow need to roll the two together into one test.

Mark_Bailey

You could analyze such a table (challenge, number of failures, total number of tests - 16) using a probit model.

Click Run.

You can extend this analysis to include the experiment.

Add Experiment as two effects: the main effect and crossed with Challenge.

The estimates can be used to test the significance of the effects and the model can be used to predict.

The Experiment example here could be another categorical or continuous factor. There could be multiple factors in the probit model.