Bonjour,
je vais faire une simulation de jeux de donnée de 25 individu de binormal corréler de 1.
le problème que R, il simule de donnée corréler de 0.2 à 0.3, et moi je vais un data corréler à 1 (ou bien 0.9).
Merci en avance de me aidée si vous avez des idées pour résoudre cette problème(même avec d'autre logiciel).
Le code R utilisée est:
mean1 <- 0.5
mean2 <- 0.5
var1 <- 0.5
var2 <- 0.5
covariance = correlation*var1*var2
sigma <- matrix(data=c(var1,covariance,covariance,var2),ncol=2)
valeur_propre <- eigen(sigma)
if (valeur_propre$values[1] > 0 & valeur_propre$values[2] > 0){
x <- rmvnorm(n=nb_studies, mean=c(mean1,mean2), sigma=sigma)
}