Solutions des TPs
633
• Fichier Poids naissance.xls
3.1- url <- "http://www.biostatisticien.eu/
springeR/Poids_naissance.xls"
poid.nais <- read.xls(url)
PTL1 <- poid.nais$PTL
PTL1[poid.nais$PTL>=2] <- 2
poid.nais <- cbind(poid.nais,PTL1)
3.2- FVT1 <- poid.nais$FVT
FVT1[poid.nais$FVT>=2] <- 2
poid.nais <- cbind(poid.nais,FVT1)
3.3- poid.nais[order(poid.nais$BWT),]
3.4- poid.nais[poid.nais$RACE<=2 & poid.nais$SMOKE==1,]
B- Gestion des valeurs manquantes
Installez le package gdata. Et puisque ce package fonctionne avec PERL,
installez aussi PERL (http://www.biostatisticien.eu/springeR/Rtools.
exe). Maintenant, on lit les donn´ ees :
url <- "http://www.biostatisticien.eu/springeR/Infarct.xls"
read.xls(url)
On constate que les valeurs manquantes sont cod´ ees avec des ".". On utilise
donc l’instruction suivante :
infarct <- read.xls(url,na.strings=".")
3.1- indNA <- which(apply(is.na(infarct),FUN=any,MARGIN=1)==TRUE)
3.2- sup1 <- function(x) {sum(x)>1}
indNA <- which(apply(is.na(infarct),FUN=sup1,MARGIN=1)==TRUE)
infarct$NUMERO[indNA]
3.3- colnames(infarct)[which(apply(is.na(infarct),FUN=any,
MARGIN=2) == TRUE)]
3.4- a) infarct[as.logical(apply(!is.na(infarct),1,prod)),]
b)infarct[!apply(apply(infarct,1,is.na),2,any),]
c) infarct[apply(!apply(infarct,1,is.na),2,all),]
d)infarct[complete.cases(infarct),]
e) na.omit(infarct)