각각 2개 수준(0,1)이 있는 약 10개의 범주형 변수가 있는 디자인을 만들려고 합니다.내 디자인이 주요 효과와 모든 바이너리 상호 작용을 모델링하고 싶습니다.JMP 기본값을 사용하면 "1" 수준의 평균 범주형 변수 수가 약 5개인 D-최적 설계를 얻습니다. 최대 3개만 허용하도록 설계를 제한하는 방법이 있습니까 변수가 "1" 수준에 있습니까?(더 많은 실행이 필요할 수 있음을 이해하지만 괜찮습니다.)
감사합니다.
같지 않다고 말하면 알고리즘을 100회 실행만으로 너무 많이 제한하고 있는 것입니다. 모형 항을 모두 추정하지 못할 수도 있습니다.
예를 들어, 3개의 요인만 높은 설정으로 제한하면 120개의 가능한 조합이 있습니다. 120개의 조합을 모두 선택해도 주효과 모형을 추정할 수 없습니다. 왜? 내장된 특이점이 있기 때문입니다. 특히 절편(설계 행렬에서 1의 열)은 모든 요소의 선형 조합과 같습니다. 1/3*(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10 ). 따라서 분석을 완료할 수 없으며 해당 모델에 대한 설계를 찾을 수 없습니다.
3 또는 4 또는 원하는 수준에 대해 "더 이상"이라고 말하지 않는 것이 더 나을 수 있습니다.
Dan, 매우 유용한 응답에 다시 한 번 감사드립니다.나는 좋은 디자인을 찾아 행복하게 사냥하고 있지만, 또 다시 난관에 부딪히게 됩니다...
높은 수준에서 범주형 요인 수의 함수로 설계 진단(특히 관심 있는 예측의 평균 분산)을 표시하는 그래프를 만들려고 합니다.모든 경우에 10개의 2-수준 요인이 있고 정해진 런 수인 100만 허용합니다.높은 수준에서 "너무 많음" 또는 "너무 적음"이 허용되면 디자인이 저하될 것으로 예상합니다.
따라서 높은 수준에서 허용되는 요소 수에 대한 일련의 다른 제약 조건을 사용하여 DOE 스크립트를 반복해서 실행하고 있습니다.
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 2
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 삼
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 4
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 5
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 6
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 7
허용되지 않는 조합( X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 != 8
(이것을 루프에 포함하는 것이 좋겠지만 스크립트를 복사하여 붙여넣기만 합니다.)
내가 찾은 것은 높은 수준에서 2개 또는 8개의 요소가 있는 요소만 수렴하여 디자인을 제공한다는 것입니다!. 다른 경우에는 "최적 설계자가 수렴하지 못했습니다"라는 오류 메시지가 나타납니다.이 문제를 해결할 방법이 있습니까?
디자인을 만들기 위해 스크립트를 사용할 필요가 없습니다. 스크립트는 모든 작업을 설정하는 편리한 방법입니다.
제약 조건에 대해 게시한 것처럼 스크립트가 필요합니다. 다양한 시나리오를 시도하고 싶다면 GUI를 사용하는 것이 좋지만 선형 제약 조건을 "허용되지 않는 조합 스크립트"로 지정하고 제가 제공한 스크립트에 붙여넣으십시오. 이렇게 하면 무언가를 시도할 수 있으며 마음에 들지 않으면 뒤로 버튼을 클릭하고 변경한 다음 다시 시도하십시오.
런 수의 경우 최소 런 수는 추정하려는 모형 항의 총 수입니다. 10개의 2수준 요인에 대해 절편(1) + 주효과(10) + 양방향 상호작용(45)이 있습니다. 이는 총 56개 항을 제공하므로 이것이 최소 실행 수입니다.
기본 런 수에는 오류에 대한 자유도(일반적으로 4)를 제공하고 직교성(포인트까지)을 위해 노력하는 몇 가지 기본 제공 기준이 있습니다. 똑바른 공식이 아닙니다. 이와 동일한 상황에서 기본 런 수는 60입니다(Ryan이 해당 설계를 생성한 이유이며 오류에 대한 자유도가 4임을 알 수 있음). 기본 실행 수는 제약 조건을 고려하지 않습니다(제가 아는 한!). 따라서 제약 조건을 포함하거나 포함하지 않고 이 동일한 설계를 생성하는 경우 기본값은 60으로 유지됩니다. 기본값은 요구 사항을 충족하는 디자인을 보장하지 않습니다! 언급한 바와 같이 제약 조건을 추가하면 좋은 설계를 찾는 것이 더 어려울 수 있으며 종종 더 많은 실행이 필요합니다. 60은 끔찍한 결과를 보였으므로 더 높은 숫자로 이동하십시오. 평균 분산 예측에 대한 특정 수준을 얻는 데 필요한 실행 수를 알려주는 쉬운 계산은 없습니다(그러나 예측 분산에 대한 분산을 잊지 마십시오! 즉, Design Space Plot의 비율을 잊지 마십시오).
이것이 도움이 되기를 바라며 좋은 디자인을 위한 행복한 사냥을 하세요!
Ryan의 접근 방식은 내가 제안하려는 것과 정확히 일치합니다. 그러나 요인을 범주형으로 유지하려는 경우(이름이 x1 - x10이라고 가정) 다음과 같은 허용되지 않는 조합 스크립트를 지정할 수 있습니다.
(x1 == 2) + (x2 == 2) + (x3 == 2) + (x4 == 2) + (x5 == 2) + (x6 == 2) + (x7 == 2)
+(x8 == 2) + (x9 == 2) + (x10 == 2) > 3
요인은 범주형이므로 각 X# == 2는 실제로 해당 요인이 두 번째 수준에 있는지 여부를 나타냅니다.
이 문제가 발생하면 일반적으로 범주형을 0/1 이산 숫자 요소로 변환한 다음 허용되지 않는 조합을 사용하여 합계가 원하는 1 수를 초과하지 않도록 합니다. 10개의 입력(주 효과 모델이지만 아이디어를 제공하기 위한 것임)의 경우 다음과 같습니다.
DOE(
Custom Design,
{Add Response( Maximize, "Y", ., ., . ), Add Factor( Discrete Numeric, {0, 1}, "X1", 0 ),
Add Factor( Discrete Numeric, {0, 1}, "X2", 0 ), Add Factor( Discrete Numeric, {0, 1}, "X3", 0 ),
Add Factor( Discrete Numeric, {0, 1}, "X4", 0 ), Add Factor( Discrete Numeric, {0, 1}, "X5", 0 ),
Add Factor( Discrete Numeric, {0, 1}, "X6", 0 ), Add Factor( Discrete Numeric, {0, 1}, "X7", 0 ),
Add Factor( Discrete Numeric, {0, 1}, "X8", 0 ), Add Factor( Discrete Numeric, {0, 1}, "X9", 0 ),
Add Factor( Discrete Numeric, {0, 1}, "X10", 0 ), Set Random Seed( 1550400750 ), Number of Starts( 3017 ),
Add Term( {1, 0} ), 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} ), Set Sample Size( 16 ), Disallowed Combinations(
X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 > 3
), Make Design}
);
