我可以在Excel中轻松完成此操作,但我的数据集已经太大了。在excel中,我会使用求解器。
Column A,B = random numbers
Column C = random number (which I want to maximize the correlation to)
Column D = A*x+B*y where x,y are coefficients resulted from solver
在一个单独的单元格中,我会有相关的(C,D)
在求解器中,我会将correl(C,D)的目标设置为max,方法是更改变量x,y并设置某些约束(例如x,y都必须是正数)。
我怎样才能在R中这样做?谢谢你的帮助。
答案 0 :(得分:9)
在R中,您创建一个function
,其输出是您想要最大化或最小化的值。基础R中包含的一个优化器称为optim()
:
set.seed(1)
A <- runif(100)
B <- runif(100)
C <- runif(100)
# these are your x and y to optimize
pars <- c(x=1,y=1)
OptPars <- function(pars,A,B,C){
D <- A*pars[1]+B*pars[2]
-cor(C,D)
}
# optim is one of many R-ish ways to do Excel's solver
# it minimizes by default (though you can tell it not to)
# and that's why I told it to take -cor()
optim(pars,OptPars,A=A,B=B,C=C)
如果您希望x
和y
有约束,请在您要优化的函数中包含约束,例如abs(x)
代替x
。