bioinformatique

Bonnjour,

quelqu'un pourrait m'expliquer comment on fait pour calculer les probas a l'exercice 4 du TP 3-4 de bioinfo ???

Merci d'avance !

> summary(data) affiche un résumé des données que tu as lu à partir de ton fichier
Taille Chirurgie Succès
Grande:343 Ouverte :350 Non:138
Petite:357 Ultrasons:350 Oui:562
> t<-table(data$Succès,data$Chirurgie,data$Taille) Ici tu crées un objet pour y fixer ta table
> t
, , = Grande


Ouverte Ultrasons
Non 71 25
Oui 192 55

, , = Petite


Ouverte Ultrasons
Non 6 36
Oui 81 234
> t[2,1,1] ici tu demandes à R de te donner la valeur de la deuxième ligne, le première colonne, la première table
[1] 192
> t[2,1,1]/sum(t[,1,1]) Ici tu calcule la probabilité d'avoir cette valeur (la case divisée par la somme des valeurs de la première colonne de la première table)
[1] 0.730038
> p<-vector (length=6) Tu crées un vecteur de longueur 6 ( pour 6 données)
> p[1]<-t[2,1,1]/sum(t[,2,1]) Tu introduits chaque proba dans le vecteur
> p[2]<-t[2,2,1]/sum(t[,1,1])
> p[3]<-t[2,1,2]/sum(t[,1,2])
> p[4]<-t[2,2,2]/sum(t[,2,2])
> p[5]<-sum(t[2,1,])/sum(t[,1,])
> p[6]<-sum(t[2,2,])/sum(t[,2,])
> p
[1] 0.7300380 0.6875000 0.9310345 0.8666667 0.7800000 0.8257143
> sum(t[,,])
[1] 700

Désolée mon code n'est pas tout à fait bien affiché mais c'est du copie-colle

Merci

De rien biquette