Note

Screenshots may differ slightly depending on software version.

Logistic Regression (Mixture)

Introduction

The process of separating mixtures of compounds is called chromatography. Ultra high performance liquid chromatography (UHPLC) is the gold standard for commercially available chromatography techniques, but it uses very expensive equipment to provide ultra high pressure.

You are a member of a group of researchers who are working to create a polymeric column that will have a pore size and skeleton thickness to achieve the same efficiency as UHPLC techniques but with lower pressure and therefore much lower cost.

Problems:
  • A formulation does not always create a polymer that is completely solid, which we will call a “homogenous” column.

  • The pores are note always interconnected, meaning the compound to be separated will not flow through the column.

First Step:
  • Design and experiment to model how the blending properties of water, monomer and surfactant effect the probability of creating an acceptable column, meaning homogenous and capable of flow.

Design

The emulsion used to create the column is a mixutre of a water, a monomer, surfactant and an initiator.

  • The monomer and water are two immiscible fluids, and the surfactant is a chemical agent that reduces the surface tension of the two fluids.

  • The monomer is what eventually becomes the solid structure of the monolith, and the water droplets in the mixture are replaced during the polymerization process by pores. The idea is to create a mixture where the water is well dispersed throughout the monomer in tiny droplets.

The researchers determined the constraints on the mixture components based on literature and experience, but are still somewhat arbitrary:

Component

Lower Bound

Upper Bound

A

Water

70

86

B

Monomer

10

25

C

Surfactant

2

8

D

Initiator*

0.02

0.25

E

Salt

0.8

0.8

Total =

100 wt %

*Initiator varies from 0.2% of monomer to 1.0% of monomer; i.e. 0.002 \(\leq\) D/B \(\leq\) 0.01

Having a binary response changes how we build the design:

  • Estimating a model for a response with two complementary outcomes requires a larger design than is required for modeling a continuous numeric response.

  • Logistic reqgression depends on having independent input factors. The collinearity in this constrained mixture requires having a larger design.

  • For a binary response it is good to have a space filling design.

  • Replicates aren’t necessary when the response is binary.

Because of the binary response we need to change the optimal design defaults.

../../_images/logmix-4.PNG

Defaults

We will:

  • Add “Additional model points” to overcome collinearity.

  • Add “Lack-of-fit points” to fill space.

  • Reduce “Replicate points” to zero.

Build the Design

1. Choose a 5-component Optimal Mixture design and enter the components and levels:

../../_images/logmix-5.PNG

Entering components and levels

2. Click on Edit constraints… and enter the constraint for “Initiator/Monomer” ratio:

../../_images/logmix-6.PNG

Entering the constraint

Click OK and then Next >>. Press Yes to accept the adjustments to the ranges.

3. Enter 10 Additional model points, 40 Lack-of-fit points, and 0 Replicate points:

../../_images/logmix-7.PNG

Specifying the design

Click Next >>.

  1. Enter two responses:

../../_images/logmix-8.PNG

Entering responses

  1. Click Finish.

Loading the Data

Rather than entering response data into this new design, you can load a design with pre-existing data by clicking on the Help, Tutorial Data menu and selecting Polymer Column.

Analyze the Design

Both homogeneous and flow are binary responses: homogeneous (1) or not (0); flow (1) or none (0). Analyze using Special Models specified for Logistic (Binary) regression.

Note

Try choosing a quadratic model and reducing it. Use backward selection with the AICc criterion.

../../_images/logmix-10.PNG

homogeneous - ANOVA

../../_images/logmix-11.PNG

homogeneous - Diagnostics

../../_images/logmix-12.PNG

homogeneous - Model Graphs

../../_images/logmix-13.PNG

flow - ANOVA

../../_images/logmix-14.PNG

flow- Diagnostics

../../_images/logmix-15.PNG

flow - Model Graphs

Optimize

For stretch target on Numerical goals for maximization, set the upper threshold limits on both responses to 0.999. This will generate the most desirable formulation for the polymeric column for homogeneity and flow.

../../_images/logmix-18.PNG

Maximize both homogeneous and flow