Hi John,

Do you want to minimize correlation or |correlation|?

Here is a simple code, with example, that shows how to optimize the magnitude 
(that is corr^2) of correlation between y and x^p, where p is a positive real 
number.  You can modify it to suit your needs.

corxy.opt <- function(x, y, interval=c(0,4), maximum=TRUE){

obj <- function(p, x, y) cor(y, x^p)^2

ans <- optimize(f=obj, interval=interval, maximum=maximum, x=x, y=y)

ans
}

set.seed(321)

x <- runif(100)

y <- rnorm(100, mean=1-2*x, sd=1)

corxy.opt(x, y)

#  Why minimizing does not make sense
obj2 <- function(p, x, y) sapply(p, function(p) cor(y, x^p)^2)

plot(seq(0,5,length=1000), obj2(seq(0,5,length=1000),x,y), type="l")

Best,
Ravi

Ravi Varadhan, Ph.D.
Assistant Professor
The Center on Aging and Health
Division of Geriatric Medicine & Gerontology
Johns Hopkins University
rvarad...@jhmi.edu<mailto:rvarad...@jhmi.edu>
410-502-2619


        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to