Question

我有一个（对我而言）非常复杂的问题。我有两个载体：

vectora <- c(111, 245, 379, 516, 671)
vectorb <- c(38, 54, 62, 67, 108)

此外我有两个变量

x = 80
y = 0.8

第三个向量基于变量x和y，方式如下：

vectorc <- vectora^y/(1+(vectora^y-1)/x)

目标是通过改变x和y来最小化vectorb和vectorc的偏差。偏差由以下函数定义：

deviation <- (abs(vectorb[1]-vectorc[1])) + (abs(vectorb[2]-vectorc[2])) + (abs(vectorb[3]-vectorc[3])) + (abs(vectorb[4]-vectorc[4])) + (abs(vectorb[5]-vectorc[5]))

我怎样才能在R？中做到这一点？

Answer 1

您可以使用optim程序！

以下是它的工作方式：

vectora <- c(111, 245, 379, 516, 671)
vectorb <- c(38, 54, 62, 67, 108)
fn <- function(v) {
    x = v[1]
    y = v[2]
    vectorc <- vectora^y/(1+(vectora^y-1)/x);
    return <- sum(abs(vectorb - vectorc))
}
optim(c(80, 0.8), fn)

输出是：

$par
[1] 91.4452617  0.8840952

$value
[1] 37.2487

$counts
function gradient 
     151       NA 

$convergence
[1] 0

$message
NULL

R中的Excel求解器 - 最小化函数

1 个答案: