Contrainte de disponibilité des matières premières dans un DOESalut,
Et si vous faites un DOE et à chaque fois que vous faites face à une contrainte de disponibilité de matière première. Bien entendu, la quantité de réactifs est limitée.
Vous pouvez sûrement résoudre le problème en exigeant le nombre total d'exécutions, mais je me demande si l'on pourrait utiliser des contraintes. Idéalement, je voudrais stipuler comme contrainte que la somme totale d'un réactif ne dépasse pas une certaine valeur. nous pourrions également vouloir utiliser exactement cette quantité afin de ne pas perdre de matière première.
Donc avoir des contraintes le long des essais.
Quelqu'un a-t-il déjà rencontré ce problème ?
Ma proposition serait d'exécuter un ensemble de DOE et de comparer la somme des réactifs limités afin de sélectionner le meilleur DOE. Quelqu'un at-il essayé cela?
Merci pour vos suggestions.
Florence
Re : Contrainte de disponibilité des matières premières dans un DOEPour contourner le problème, j'ai développé un petit script. Ce script construit un DOE avec les facteurs de la table associée (changez le chemin à la ligne 2).
J'ai considéré que le facteur Réactif n'était disponible qu'en 240 unités (ligne 4).
Ensuite, je lance le DOE de sa taille minimale (j'ai récupéré le script DOE de la plateforme) à une taille maximale en nombre d'essais (ici, de 12 à 25 essais) et je teste la somme de chaque DOE sur la variable réactif et vérifier qu'il est inférieur ou égal à 240. Il vous laisse les tableaux possibles (titre du tableau = la somme)
Cette approche n'est pas la plus idéale car je ne prends pas en compte cette contrainte lors de l'optimisation du critère d'optimalité comme on le ferait avec une contrainte inter-facteurs mais elle permet de garder les plans acceptables par cette contrainte sur un ensemble de plans.
FactorsTable= open("Reactor Factors.jmp");
ReactantMaxSum = 240;
For( i = 12, i < 25, i++,
DOeInit = DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ), Set Random Seed( 4547922 ),
Optimality Criterion( Name( "Make D-Optimal Design" ) ), Simulate Responses( 0 ),
Save X Matrix( 0 )}
);
DOeInit << Load Factors( FactorsTable );
DOeInit << Set Sample Size( i );
DOeInit << Make Table();
temp = Current Data Table();
tempdata = temp << Summary( Sum( :Reactant ), Freq( "Aucun(e)" ), Weight( "Aucun(e)" ) );
mySum = Column( 2 )[1];
Show( mySum );
If( mySum > ReactantMaxSum,
Close( temp, NoSave );
Close( tempdata, NoSave );
DOeInit<<CLOSE window=""></CLOSE>Re : Contrainte de disponibilité des matières premières dans un DOEQuelqu'un m'a posé une question similaire récemment. Je pense qu'une autre approche consisterait à traiter le réactif comme un facteur covariable. Si vous souhaitez simplement tester l'effet principal du facteur, l'idéal est que le facteur ne prenne que les limites inférieure et supérieure de la plage de facteurs (10 et 20 dans cet exemple).
Pour 16 exécutions, vous pouvez équilibrer correctement avec 8 exécutions à Reactant = 10 et 8 exécutions à Reactant = 20. Je m'attendrais donc à ce qu'une conception optimale pour l'effet principal du réactif réponde à la contrainte de réactif total <= 240.
Pour plus d'exécutions, vous pouvez calculer toutes les valeurs de Reactant qui donnent le plus d'exécutions à 20, et donc le meilleur équilibre, tout en respectant la contrainte de Reactant <= 240. Utilisez-les ensuite comme tables d'ensembles candidats pour une conception optimale.
J'ai déterminé les niveaux optimaux de réactif pour 16, 18 et 20 cycles (fichier .xlsx joint). Tableau de covariables JMP à 20 exécutions également joint.
À titre d'exemple, j'ai utilisé ceci pour créer la conception ci-jointe pour 20 exécutions.Je pense que c'est la conception D-optimale pour les effets principaux, compte tenu de la contrainte de réactif total <= 240.
C'est un problème intéressant. J'espère que ma solution sera utile à quelqu'un !
Phil
