Statistiques

tout sur les statistiques
 
AccueilAccueil  FAQFAQ  RechercherRechercher  S'enregistrerS'enregistrer  MembresMembres  GroupesGroupes  Connexion  

Partagez | 
 

 Genérer des covariables fortement corrélées sur R

Voir le sujet précédent Voir le sujet suivant Aller en bas 
AuteurMessage
smile24



Nombre de messages : 1
Date d'inscription : 23/04/2011

MessageSujet: Genérer des covariables fortement corrélées sur R   Sam 23 Avr - 0:17

Bonjour,

J'effectue une simulation de données pour une analyse de survie avec le modèle Cox. Mon étude porte sur un cas spécial où des covariables sont fortement corrélées. Lorsque je genere X iid N(0,1) elles ne le sont pas puisque normales, comment pourrais-je générer des variables pour les forcer à être corrélées?
J'ai inclus le code R effectué ci-dessous.
Merci !

n = 10000
p=1000
beta= c(2,-1,3, rep(0,997))
lambdaT = .002
lambdaC = .004
X<-matrix(rnorm(n*p), nrow=n, ncol=p) # Matrice de covariables Xi iid N(0,1)


T = rweibull(n, shape=1, scale=lambdaT*exp(X %*% beta))
C = rweibull(n, shape=1, scale=lambdaC) #Temps censuré
time = pmin(T,C) #Temps observé
evenement = time==T # 1 si l'évènement est observé

table(evenement)
evenement
FALSE TRUE
4083 5917


library(survival)
fit1 <-coxph(Surv(time, evenenement)~ X, method="breslow")
Revenir en haut Aller en bas
Voir le profil de l'utilisateur
 
Genérer des covariables fortement corrélées sur R
Voir le sujet précédent Voir le sujet suivant Revenir en haut 
Page 1 sur 1
 Sujets similaires
-
» Cedric Villani : Le respect dû par une société à ses enseignants est très fortement corrélé aux résultats
» Corrélation entre 2 séries de données, quel test?
» Corrélation et significativité
» Comparaison de corrélation sur 2 population indépendantes
» test de pearson et coefficient de corrélation [RESOLU]

Permission de ce forum:Vous ne pouvez pas répondre aux sujets dans ce forum
Statistiques :: Statistiques-
Sauter vers: