我试图做一个普通的单纯形(三角形或四面体的概念到任意尺寸)来开始一组优化实验。Optimsimplex包通过使用 Spendley 方法提供了一种简单而有用的方法来实现这一点:
library('optimsimplex') #Paquete necesario
Ultra <- optimsimplex(method ='spendley',
x0=c(Vhno3=3,Vh2o2=1,Msample=300,Tsonic=15))
结果Ultra
是一个包含空间维度 (n) 的 optimsimplex 类对象,并且每个 (n+1) 顶点的 (n) 坐标。可以使用 len 选项指定单纯形的维度(长度):
len:单纯形的维度。如果长度是一个值,则在所有方向上都使用该唯一长度。如果 length 是具有 n 个值的向量,则每个长度都与相应的方向一起使用。仅在方法设置为 'axes' 或 'spendley' 时使用</p>
但是这个结果出现了一个我无法理解的错误:
Ultra <- optimsimplex(method ='spendley',
x0=c(Vhno3=3,Vh2o2=1,Msample=300,Tsonic=15),
len=c(pVhno3=0.5,pVh2o2=0.25,pMsample=50,pTsonic=5))
错误:optimsimplex:len 向量应为行矩阵,但当前形状为 1 x 4
那么,1 x 4 不是 {optimsimplex} 预期的行矩阵吗?这可能对应于包中的某种错误吗?提前致谢。