Hola @angela ,
Este tema parece similar a este, que trata sobre "Mezclas de mezclas": Cómo tener un diseño de mezcla de mezclas - Comunidad de Usuarios JMP.Pueden ser posibles dos opciones:
- Trabaje con componentes individuales de sus "grupos" de formulación: descripción clara del impacto de cada componente individual, pero puede representar una gran cantidad de experimentos, o
- Trabaje con mezclas de grupos previamente probadas/predefinidas: se necesitan muchos menos experimentos, pero puede perder cierta información en el proceso, por ejemplo, interacciones entre componentes individuales.
Como está interesado en examinar los factores del experimento, tengo dos preguntas:
- ¿Es un diseño de mezcla la única opción (lo que quiero decir aquí es que los diferentes componentes SIEMPRE se introducen en la mezcla y sus cantidades en relación con otras partes de la formulación)? Intentar seleccionar factores con un diseño de mezcla puede resultar complicado y dar lugar a una gran cantidad de experimentos, de ahí la pregunta.
- ¿Son necesarias todas las restricciones? Al reducir el espacio experimental, es posible que se pierda alguna información.
Intenté probar un diseño D-óptimo con todas sus restricciones, pero parece que no puedo introducir la restricción 4 cuando se usan ingredientes del grupo G como factores de mezcla. El diseño sin esta restricción representaría 105 experimentos y se vería así (script aquí para probarlo, tabla de datos incluida si es más conveniente para usted):
DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ), Add Factor( Mixture, 0, 1, "S1", 0 ),
Add Factor( Mixture, 0, 1, "S2", 0 ), Add Factor( Mixture, 0, 1, "F1", 0 ),
Add Factor( Mixture, 0, 1, "F2", 0 ), Add Factor( Mixture, 0, 1, "T1", 0 ),
Add Factor( Mixture, 0, 1, "T2", 0 ), Add Factor( Mixture, 0, 1, "T3", 0 ),
Add Factor( Mixture, 0, 1, "T4", 0 ), Add Factor( Mixture, 0, 1, "H1", 0 ),
Add Factor( Mixture, 0, 1, "H2", 0 ), Add Factor( Mixture, 0, 1, "H3", 0 ),
Add Factor( Mixture, 0, 1, "G1", 0 ), Add Factor( Mixture, 0, 1, "G2", 0 ),
Add Factor( Mixture, 0, 1, "G3", 0 ), Add Factor(
Continuous, 1.5, 3, "Water", 0
), Set Random Seed( 800692047 ), Number of Starts( 2 ),
Add Constraint(
[-1 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.35,
1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0.55,
0 0 -1 -1 0 0 0 0 0 0 0 0 0 0 0 -0.24,
0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0.34,
0 0 0 0 -1 -1 -1 -1 0 0 0 0 0 0 0 -0.07,
0 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0.14,
0 0 0 0 0 0 0 0 0 0 0 -1 -1 -1 0 -0.001,
0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0.015]
), Add Term( {1, 1} ), Add Term( {2, 1} ), Add Term( {3, 1} ),
Add Term( {4, 1} ), Add Term( {5, 1} ), Add Term( {6, 1} ), Add Term( {7, 1} ),
Add Term( {8, 1} ), Add Term( {9, 1} ), Add Term( {10, 1} ), Add Term( {11, 1} ),
Add Term( {12, 1} ), Add Term( {13, 1} ), Add Term( {14, 1} ),
Add Term( {15, 1} ), Add Alias Term( {1, 1}, {2, 1} ),
Add Alias Term( {1, 1}, {3, 1} ), Add Alias Term( {1, 1}, {4, 1} ),
Add Alias Term( {1, 1}, {5, 1} ), Add Alias Term( {1, 1}, {6, 1} ),
Add Alias Term( {1, 1}, {7, 1} ), Add Alias Term( {1, 1}, {8, 1} ),
Add Alias Term( {1, 1}, {9, 1} ), Add Alias Term( {1, 1}, {10, 1} ),
Add Alias Term( {1, 1}, {11, 1} ), Add Alias Term( {1, 1}, {12, 1} ),
Add Alias Term( {1, 1}, {13, 1} ), Add Alias Term( {1, 1}, {14, 1} ),
Add Alias Term( {1, 1}, {15, 1} ), Add Alias Term( {2, 1}, {3, 1} ),
Add Alias Term( {2, 1}, {4, 1} ), Add Alias Term( {2, 1}, {5, 1} ),
Add Alias Term( {2, 1}, {6, 1} ), Add Alias Term( {2, 1}, {7, 1} ),
Add Alias Term( {2, 1}, {8, 1} ), Add Alias Term( {2, 1}, {9, 1} ),
Add Alias Term( {2, 1}, {10, 1} ), Add Alias Term( {2, 1}, {11, 1} ),
Add Alias Term( {2, 1}, {12, 1} ), Add Alias Term( {2, 1}, {13, 1} ),
Add Alias Term( {2, 1}, {14, 1} ), Add Alias Term( {2, 1}, {15, 1} ),
Add Alias Term( {3, 1}, {4, 1} ), Add Alias Term( {3, 1}, {5, 1} ),
Add Alias Term( {3, 1}, {6, 1} ), Add Alias Term( {3, 1}, {7, 1} ),
Add Alias Term( {3, 1}, {8, 1} ), Add Alias Term( {3, 1}, {9, 1} ),
Add Alias Term( {3, 1}, {10, 1} ), Add Alias Term( {3, 1}, {11, 1} ),
Add Alias Term( {3, 1}, {12, 1} ), Add Alias Term( {3, 1}, {13, 1} ),
Add Alias Term( {3, 1}, {14, 1} ), Add Alias Term( {3, 1}, {15, 1} ),
Add Alias Term( {4, 1}, {5, 1} ), Add Alias Term( {4, 1}, {6, 1} ),
Add Alias Term( {4, 1}, {7, 1} ), Add Alias Term( {4, 1}, {8, 1} ),
Add Alias Term( {4, 1}, {9, 1} ), Add Alias Term( {4, 1}, {10, 1} ),
Add Alias Term( {4, 1}, {11, 1} ), Add Alias Term( {4, 1}, {12, 1} ),
Add Alias Term( {4, 1}, {13, 1} ), Add Alias Term( {4, 1}, {14, 1} ),
Add Alias Term( {4, 1}, {15, 1} ), Add Alias Term( {5, 1}, {6, 1} ),
Add Alias Term( {5, 1}, {7, 1} ), Add Alias Term( {5, 1}, {8, 1} ),
Add Alias Term( {5, 1}, {9, 1} ), Add Alias Term( {5, 1}, {10, 1} ),
Add Alias Term( {5, 1}, {11, 1} ), Add Alias Term( {5, 1}, {12, 1} ),
Add Alias Term( {5, 1}, {13, 1} ), Add Alias Term( {5, 1}, {14, 1} ),
Add Alias Term( {5, 1}, {15, 1} ), Add Alias Term( {6, 1}, {7, 1} ),
Add Alias Term( {6, 1}, {8, 1} ), Add Alias Term( {6, 1}, {9, 1} ),
Add Alias Term( {6, 1}, {10, 1} ), Add Alias Term( {6, 1}, {11, 1} ),
Add Alias Term( {6, 1}, {12, 1} ), Add Alias Term( {6, 1}, {13, 1} ),
Add Alias Term( {6, 1}, {14, 1} ), Add Alias Term( {6, 1}, {15, 1} ),
Add Alias Term( {7, 1}, {8, 1} ), Add Alias Term( {7, 1}, {9, 1} ),
Add Alias Term( {7, 1}, {10, 1} ), Add Alias Term( {7, 1}, {11, 1} ),
Add Alias Term( {7, 1}, {12, 1} ), Add Alias Term( {7, 1}, {13, 1} ),
Add Alias Term( {7, 1}, {14, 1} ), Add Alias Term( {7, 1}, {15, 1} ),
Add Alias Term( {8, 1}, {9, 1} ), Add Alias Term( {8, 1}, {10, 1} ),
Add Alias Term( {8, 1}, {11, 1} ), Add Alias Term( {8, 1}, {12, 1} ),
Add Alias Term( {8, 1}, {13, 1} ), Add Alias Term( {8, 1}, {14, 1} ),
Add Alias Term( {8, 1}, {15, 1} ), Add Alias Term( {9, 1}, {10, 1} ),
Add Alias Term( {9, 1}, {11, 1} ), Add Alias Term( {9, 1}, {12, 1} ),
Add Alias Term( {9, 1}, {13, 1} ), Add Alias Term( {9, 1}, {14, 1} ),
Add Alias Term( {9, 1}, {15, 1} ), Add Alias Term( {10, 1}, {11, 1} ),
Add Alias Term( {10, 1}, {12, 1} ), Add Alias Term( {10, 1}, {13, 1} ),
Add Alias Term( {10, 1}, {14, 1} ), Add Alias Term( {10, 1}, {15, 1} ),
Add Alias Term( {11, 1}, {12, 1} ), Add Alias Term( {11, 1}, {13, 1} ),
Add Alias Term( {11, 1}, {14, 1} ), Add Alias Term( {11, 1}, {15, 1} ),
Add Alias Term( {12, 1}, {13, 1} ), Add Alias Term( {12, 1}, {14, 1} ),
Add Alias Term( {12, 1}, {15, 1} ), Add Alias Term( {13, 1}, {14, 1} ),
Add Alias Term( {13, 1}, {15, 1} ), Add Alias Term( {14, 1}, {15, 1} ),
Set Sample Size( 105 ), Simulate Responses( 0 ), Save X Matrix( 0 ), Make Design
}
);
Dado que los ingredientes del grupo G tienen una concentración muy baja en comparación con otros grupos de ingredientes, una solución alternativa puede ser considerar estos ingredientes como factores continuos para introducir la restricción 4 (aún no probada, lo intentaré más adelante cuando tenga más tiempo).
Estos son sólo algunos aportes para conocer mejor su tema, de modo que otros miembros puedan unirse a las discusiones y brindar ideas relevantes.