«Оптимальный дизайн» с фиксированными точками

Question

Я использую библиотеку OptimalDesign в R, чтобы найти d-оптимальные проекты с функцией "od.AA. Мое физическое ограничение заключается в том, что у меня есть некоторые фиксированные точки, которые я должен использовать.

Есть ли:

• способ указать, что некоторое подмножество контрольных точек является фиксированным / обязательным.
• другая библиотека, которую я должен использовать, которая имеет эту функциональность.

Это пример кода:

# Create the regressors matrix for the third-degree trigonometric model 
# on a "partial circle".
F.trig <- F.cube(~I(cos(x1)) + I(sin(x1)) + 
          I(cos(2 * x1)) + I(sin(2 * x1)) + 
          I(cos(3 * x1)) + I(sin(3 * x1)), 
          -pi/2, pi/2, 200)

# Compute the D-, A-, and IV-optimal standardized approximate designs.
res.trig.D <- od.AA(F.trig, 1, alg = "doom", crit = "D", 
              graph=1:7, t.max=4)

Он выводит это в виде графика:

Он также дает их как подмножество из 200 кандидатовкоторые составляют точки:

> F.trig[which(res.trig.D$w.best>0),]
    (Intercept)   I(cos(x1))   I(sin(x1)) I(cos(2 * x1)) I(sin(2 * x1)) I(cos(3 * x1)) I(sin(3 * x1))
1             1 6.123032e-17 -1.000000000     -1.0000000  -1.224606e-16  -1.836910e-16     1.00000000
22            1 3.254852e-01 -0.945547126     -0.7881187  -6.155232e-01  -8.385272e-01     0.54485968
58            1 7.832356e-01 -0.621724976      0.2269161  -9.739143e-01  -4.277781e-01    -0.90388380
59            1 7.929527e-01 -0.609283155      0.2575481  -9.662655e-01  -3.845058e-01    -0.92312256
100           1 9.999688e-01 -0.007893367      0.9998754  -1.578624e-02   9.997196e-01    -0.02367813
101           1 9.999688e-01  0.007893367      0.9998754   1.578624e-02   9.997196e-01     0.02367813
142           1 7.929527e-01  0.609283155      0.2575481   9.662655e-01  -3.845058e-01     0.92312256
143           1 7.832356e-01  0.621724976      0.2269161   9.739143e-01  -4.277781e-01     0.90388380
179           1 3.254852e-01  0.945547126     -0.7881187   6.155232e-01  -8.385272e-01    -0.54485968
200           1 6.123032e-17  1.000000000     -1.0000000   1.224606e-16  -1.836910e-16    -1.00000000

Если бы я должен был выбрать какую-то случайную точку, такую ​​как индекс # 83, как «должен быть исправлен», как я могу повторить это с этим, и получить обновленный d-оптимальный набор, который включает в себя это место выборки?

> F.trig[83,]
   (Intercept)     I(cos(x1))     I(sin(x1)) I(cos(2 * x1)) I(sin(2 * x1)) I(cos(3 * x1)) I(sin(3 * x1)) 
     1.0000000      0.9620794     -0.2727697      0.8511934     -0.5248522      0.6757518     -0.7371292