2017-06-15 101 views
0

假设在下面的方程(1)d = .99。另外,sd = 1.2。理想的是,5+(.1*5) <= m1 <= 5+(.5*5)5 <= m2 <= 5+(.3*5)用2个未知数求解一个方程,得到R中的可能答案的具体范围

公式(1):d = (m1-m2)/sd

问题

有一定的许多m1m2可能的答案。但R中,我怎么能获得为我上​​面指定m1m2下降的m1的范围内且m2可能的答案(这就是为什么我用“理想”上面)?

+0

如果答案没有确定,因为您似乎暗示那么有无数的可能的解决方案。绘制定义平等我会有帮助。不幸的是,你的不平等限制对任何变量都没有依赖性。 –

+0

@ 42-,谢谢你的洞察力。如果可以的话,你可能会提供一个简短的回应,说明如何实施你的方法吗? – rnorouzian

+0

只有当你对不平等的限制取决于一个变量时,它才对我有意义。目前看起来像m1和m2变量位于矩形框中。相当微不足道。 'sd'和'mean'之间的联系完全不清楚。 –

回答

1

解决您的方程m1产量是m1 = m2 + d*sd,所以:。

m1 = m2 + 1.188 

你的不平等

5.5 <= m1 <= 7.5 
5.0 <= m2 <= 6.5 

如果我们m2 + 1.188在第一个不等式更换m1和简化,我们得到了两个不等式:

4.312 <= m2 <= 6.312 
    5.0 <= m2 <= 6.5 

要让他们两个真的,我们需要

max(4.312,5.0) <= m2 <= min(6.312,6.5) 

so:

5.0 <= m2 <= 6.312 

在R你可以做例如

> m2 <- seq(5.0,6.312,length.out = 10) 
> m1 <- m2 + 1.188 
> cbind(m1,m2) 
      m1  m2 
[1,] 6.188000 5.000000 
[2,] 6.333778 5.145778 
[3,] 6.479556 5.291556 
[4,] 6.625333 5.437333 
[5,] 6.771111 5.583111 
[6,] 6.916889 5.728889 
[7,] 7.062667 5.874667 
[8,] 7.208444 6.020444 
[9,] 7.354222 6.166222 
[10,] 7.500000 6.312000 
2

我会尝试一些数据可视化,以查看是否在评论点可以明确:

> d = .99; sd = 1.2; png(); plot(x=seq(-10,10), -10:10) 
> abline(h= c(5+(.1*5), 5+(.5*5))) 
> abline(v=c(5+(.1*5) , 5+(.5*5))) 

enter image description here

(我不知道这是不是一个答案..随意downvote或西红柿扔我不会在意)

相关问题