0

一些数据：

a <- function(a) {3*a+12*a^2}
dom1 <- seq(-1,4,0.1)
vec1 <- a(dom1)[1:10]

c <- function(c) {-5*c^2+2*c^3}
dom2 <- seq(-0.1,0.2,0.01)
vec2 <- c(dom2)[1:10]

d <- function(d) {2*d^2+5*d^3+12*d^4}
dom3 <- seq(0.1,0.5,0.01)
vec3 <- d(dom3)[1:10]

w <- function(w) {7*w-3*w^2}
dom4 <- seq(0.5,2.5,0.05)
vec4 <- w(dom4)[1:10]

现在假设我们适合lm model更大的数据集lm(y~a+c+d+w)并且 lm 参数是c(-0.2,0.2,0.1,0.6)

fun.mean <- function(a,c,d,w) {-0.2*a+0.2*c+0.1*d+0.6*w}

我试过了，但没有按预期工作：

这是循环从相关函数（&domains）生成的向量

for(a in vec1) {
 for(c in vec2) {
  for(d in vec3) {
   for(w in vec4) {
    sol <- fun.mean(a,c,d,w)
    if (sol %% 1 > 0.40 & sol < 0.50) print(c(a,c,d,w))
  }}}}

所以我正在寻找的是找到 c(a,c,d,w) 的组合，它等于 0.5 或理想情况下等于 interval 0.4-0.5。

c(a,c,d,w)因此，将输出与函数相乘fun.mean不会得到所需的值（间隔 0.4-0.5）。我做错了什么？给定“目标”值，是否有更好的方法来找到 c(a,c,d,w) 值？什么是替代方案，for loop因为它非常慢。

4

1 回答 1

1

您使用 %%1 进行过滤，但在检查答案时不使用它。以下代码返回 9.00 -0.052000 0.02620000 2.7500，即 -0.15778。它符合您的标准（-0.15778 %% 1 = 0.85 > 0.4 和 -0.15 < 0.5）。

我猜你的标准是错误的，你可能应该这样做

sol %% 1 > 0.4 and sol %% 1 < 0.5

或者

sol > 0.4 and sol < 0.5.

在第一种情况下，也许您应该在 fun.mean 计算中添加 %% 1 。

于 2013-10-07T10:07:39.650 回答