Question

我有两个数组。例如，第一个是A，第二个是B：

A <-c（2,5,6,10,11） B <-c（13,2,6,8,12）

我的常数teta = 1。我做了一个数组：

D = B-teta

以及A和D之间的相关系数

koeff <-cor（A，D）

面临的挑战是通过更改常数teta来使koeff达到最大值。据我了解，我可以使用nloptr库。我应该如何使用它？

最诚挚的问候！

P.S。我在excel中使用求解器一般减小的梯度进行了此工作。我在R中找不到此功能。

Answer 1

说明log(B+ϑ)不会更改顺序。

> # ranks for different theta
> theta <- -1; rank(log(B+theta))
[1] 5 1 2 3 4
> theta <- 3; rank(log(B+theta))
[1] 5 1 2 3 4
> theta <- 10; rank(log(B+theta))
[1] 5 1 2 3 4
> # theta=-2 is a border case: we get -infinity
> theta=-2; log(B+theta)
[1] 2.397895     -Inf 1.386294 1.791759 2.302585
> rank(log(B+theta))
[1] 5 1 2 3 4
>

因此，所有Spearman相关性都应该相同：

> A <- c(2,5,6,10,11)
> cor(A,B,method="spearman")
[1] 0
> theta <- -1; cor(A,log(B+theta),method="spearman")
[1] 0
> theta <- 3; cor(A,log(B+theta),method="spearman")
[1] 0
> theta <- 10; cor(A,log(B+theta),method="spearman")
[1] 0
> theta <- -2; cor(A,log(B+theta),method="spearman")
[1] 0
>

如何使用nloptr库最大化两个数组之间的Spearman相关性

1 个答案: