r - 查找产生相关值的单个缺失值

Question

我正在尝试估算一个变量的缺失值，以便它产生给定的相关值

library(MASS)
mat <- mvrnorm(49, mu = c(0,5), Sigma = matrix(c(1,0.05,.05,1), ncol = 2), empirical = TRUE)

cor50row <-  function(x,y, rho){

  y_lnt <- length(y)

  x[length(x) +1] <- mean(x)

  val <- seq(-1000,0, .01)

  for(indx in val){
    y[y_lnt + 1] <- indx
    if(rho - cor(x,y) < 1e-6){
      break
    }
  }
  return(cbind(x,y))
}
a <- cor50row(x = mat[,1], y= mat[,2], rho = .06)

所以我们的想法是找到 y 变量的缺失值，该值将相关性增加 0.01

标签： rstatisticscorrelationimputation