Индекс вне границ на языке R - PullRequest
0 голосов
/ 25 февраля 2019
genBi <- function(rho, mu1, mu2, s1, s2){
  library(MASS)
  mu <- c(mu1, mu2) #mean
  sigma <- matrix(c(s1^2, s1*s2*rho, s1*s2*rho, s2^2),
                  2) #covariance matrix
  bvn1 <- mvrnorm(4000, mu = mu, Sigma = sigma )
  colnames(bvn1) <- c("x","y")
  return(bvn1)
}

#get samples from Z with sample size 20,
getSlice <- function(data){
  Z <- seq(0, 1, length.out = 200) #initialize Z
  for(i in 0:199){
    temp <- data[i * 20 + 1 : (i + 1) * 20, ]
    R <- cor(temp[,1], temp[,2])
    Z[i] <- 0.5 * log((1 + R)/(1 - R))
  }
  return(Z)
}
data <- genBi(0.6, 1, 1, 2, 2)
Z <- getSlice(data)

Возврат с ошибкой Ошибка в temp [, 2]: индекс за пределами.Пожалуйста, помогите определить проблему!

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...