如何求解R中的一个不等式方程



我使用Tim Gerrodette的不等式方程(1987(来研究视觉调查数据的幂分析,但我很难找到用R.解决x的问题

我有以下方程式

r^2*n^3>12*cv^2(za/2+zB(^2

其中r=0.02,n=40,cv=0.67,za/2=1.645,我想求解zB(x,本质上(翻译过来就是:

0<zB^2+3.3zB-2.04其中-3.4<zB<3.4

我该如何在R中解决这个问题?我看过quadprog和limsolve,但我不确定我是否走在了正确的轨道上,因为许多人似乎正在使用这些函数来求解多个/联立方程。

我正在寻找一种在R中可以使用的方法,并在R、n、CV等的一系列值中使用sub。

如果问题是如何获得zB满足3个不等式的值集,我们可以使用蛮力将其确定为小数点后3位(或通过更改.001来确定其他(,从而使下面显示的区间满足不等式。

zB <- seq(-3.4, 3.4, .001)
range(zB[zB > -3.4 & zB < 3.4 & zB^2 + 3.3 * zB > 2.04])
## [1] 0.533 3.399

最新更新