我有一个实验设计。我想计算它的D效率。 我认为R包AlgDesign可以提供帮助。我找到了生成设计的函数optFederov，并且 - 如果用户需要 - 返回其D效率。但是，我不想使用optFederov来生成设计 - 我已经有了我的设计！ 我试过eval.design（〜。，mydesign）。但它给我的唯一指标是：行列式，A，对角性和gmean.variances。也许有一种方法可以从行列式或A到D效率（我不是数学家，所以我不确定）。或者也许是其他一些方法来计算D效率＆＃34;手动＆＃34;那么说呢？
非常感谢任何暗示！
答案 0 :(得分：2)
我正在做一个类似的项目。我在这个link中找到了这个公式Deff =（| X'X | ^（1 / p））/ ND。其中X是模型矩阵，p是线性模型中的beta数，ND是实验运行的数量。你可以制作这样的代码，它就可以了。
det(t(X)%*%X)^(1/beta)/(numRuns)
我使用JMP为我的项目测试了结果，所以我相信这是正确的公式