Вычислить D-эффективность экспериментального desgin в R

У меня есть экспериментальный дизайн. Я хочу рассчитать его D-эффективность. Я думал, что R-пакет AlgDesign может помочь. Я нашел функцию optFederov, которая генерирует дизайн, и - если пользователь хочет - возвращает свою эффективность D. Тем не менее, я не хочу использовать optFederov для создания дизайна - у меня уже есть мой дизайн! Я пробовал eval.design (~., Mydesign). Но единственными метриками, которые она дает мне, являются: детерминант, А, диагональность и gmean.variances. Может быть, есть способ получить от детерминанта или от А до D-эффективности (я не математик, поэтому я не уверен). Или, может быть, какой-то другой способ рассчитать D-эффективность «вручную», так сказать?Вычислить D-эффективность экспериментального desgin в R

Большое спасибо за любой намек!

источник

2015-03-20 user2323534

Я никогда не использовал его, но 'eval.design' говорит, что выводит«Dea - Нижняя граница по эффективности D для приближенных теоретических конструкций. Он равен exp (1 - 1/Ge). " В нем также говорится: «I, Ge и Dea рассчитываются только при вводе X». Так что, возможно, вам нужно пройти в 'X'. – Gregor

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

Я работал над подобным проектом. Я нашел эту формулу Deff = (| X'X |^(1/p))/ND в этом link. Где X - модельная матрица, p - количество бета в вашей линейной модели, а ND - количество прогонов вашего эксперимента. Вы могли бы просто сделать такой код, и он будет делать трюк. ?

det(t(X)%*%X)^(1/beta)/(numRuns)

Я проверил результаты, используя JMP для моего проекта, поэтому я считаю, что это правильная формула

источник

2016-06-24 19:24:16 JoseBatAbi

Вычислить D-эффективность экспериментального desgin в R

ответ

Смежные вопросы