Commit 9327f867 authored by pdaude's avatar pdaude

Exercie 1 module 3 Done

parent 4aaa1ea3
...@@ -28,8 +28,19 @@ Nous commençons donc par charger ces données: ...@@ -28,8 +28,19 @@ Nous commençons donc par charger ces données:
```{r} ```{r}
data = read.csv("shuttle.csv",header=T) data = read.csv("shuttle.csv",header=T)
data data
```
```
```{r}
data$m<-ifelse(data$Malfunction>1,1,0)
plot(data=data, m~Temperature)
plot(data=data, Malfunction/Count ~ Temperature, ylim=c(0,1))
logistic_reg = glm(data=data, Malfunction/Count ~ Temperature, weights=Count,
family=binomial(link='logit'))
summary(logistic_reg)
library(ggplot2)
ggplot(data,aes(x=Temperature,y=Malfunction/Count))+
geom_smooth(method="glm",method.args=list(family="binomial"))
```
Le jeu de données nous indique la date de l'essai, le nombre de joints Le jeu de données nous indique la date de l'essai, le nombre de joints
toriques mesurés (il y en a 6 sur le lançeur principal), la toriques mesurés (il y en a 6 sur le lançeur principal), la
température (en Farenheit) et la pression (en psi), et enfin le température (en Farenheit) et la pression (en psi), et enfin le
...@@ -44,8 +55,12 @@ moins un joint a été défectueux. ...@@ -44,8 +55,12 @@ moins un joint a été défectueux.
```{r} ```{r}
data = data[data$Malfunction>0,] data = data[data$Malfunction>0,]
data data
``` ```
Très bien, nous avons une variabilité de température importante mais Très bien, nous avons une variabilité de température importante mais
la pression est quasiment toujours égale à 200, ce qui devrait la pression est quasiment toujours égale à 200, ce qui devrait
simplifier l'analyse. simplifier l'analyse.
......
This diff is collapsed.
...@@ -43,10 +43,16 @@ Voici l'explication des colonnes donnée sur le [sur le site d'origine](https:// ...@@ -43,10 +43,16 @@ Voici l'explication des colonnes donnée sur le [sur le site d'origine](https://
La première ligne du fichier CSV est un commentaire, que nous ignorons en précisant `skip=1`. La première ligne du fichier CSV est un commentaire, que nous ignorons en précisant `skip=1`.
### Téléchargement ### Téléchargement
```{r} ```{r}
data = read.csv(data_url, skip=1) #data = read.csv(data_url, skip=1)
data_local='incidence-PAY-3.csv'
data = read.csv(data_local, skip=1)
``` ```
Regardons ce que nous avons obtenu: Regardons ce que nous avons obtenu:
```{r} ```{r}
head(data) head(data)
......
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment