Skip to main content

Posts

Showing posts from May, 2016

Simulation de maladie (Petit exemple avec r)

Dans une population la probabilité de loccurrence d'une maladie (Y) depend du sexe  (S) et de la presence d'une mutation (M). en plus on sait qe sexe et mutation sont des fcteurs independants, P(S=M) = 0.5 P(M=1) = 0.25, on note      M = 1 quand la maladie est presente sinon M=0                   Y = 1 si l'individu est malade Y=0 sinon                    S = 1 si ll'individu est un male S=0 sinon On suppose que la facon dont  probabilite de la maladie est relie aux autres variable peut s'ecrire Logit(P(Y=1/ S=m,M=1) ) = log(1.1) + log(1.5) * s + log(2) * m logit(p) = $\quad \log { \frac { p }{ 1-p }  } $ Nous allons coder sur r et obtenir un echantillon aleatoir de 1000 000  (Si. Mi. Yi )  Nous allons obtenir des graphiques et verifier les resultats theoriques avec ca CODE set.seed(100) n = 1000000 # le vecteur (1,0) ...