getwd()
setwd("/home/loic/R")
guichet=read.table("guichet.txt");guichet
load("guichet.rda")
typeof(guichet)
is.data.frame(guichet)
q=guichet$queue
q=guichet[[1]];q
attach(guichet);queue
summary(queue)
var(queue)
sd(queue)
sqrt(var(queue))
x11()
sum(table(queue))
boxplot(queue)
163/2000
(163+418)/2000
(163+418+500)/2000
queue
hist(queue)
table(queue)
typeof(table(queue))


tqueue=table(queue);tqueue # tableau des effectifs

#par(mfrow=c(2,2))
plot(tqueue) # diagramme en bâtons, effectifs
plot(tqueue/(sum(tqueue))) # diagramme en bâtons, fréquences

plot(cumsum(tqueue)/sum(tqueue),type="s") # courbe fréquences cumulées, à la main
plot(ecdf(queue)) # courbe fréquences cumulées, direct
x11()
plot(tqueue/(sum(tqueue))) # diagramme en bâtons, fréquences
lines(0:9,dpois(0:9,mean(queue),log=FALSE),col="red")
help(plot)
help(dpois)

# À retenir : load, read.table, var, sd, boxplot, summary, plot, ecdf, lines, dpois, help


# Exercice 3 ####

tenfants=matrix(c(1:8,10,510,248,132,55,29,15,2,8,1),ncol=2);tenfants
tenfants=c(510,248,132,55,29,15,2,8,0,1);tenfants
tenfants[3]
enfants=rep(1:10,tenfants);enfants
dev.off()
x11()
plot(tenfants,type="h")
plot(table(enfants)/)
which.max(tenfants)
prop.table(table(enfants))
table(enfants)/sum(table(enfants)) #tableau des fréquences
cumsum(table(enfants))/sum(table(enfants)) #tableau des fréquences cumulées
plot(ecdf(enfants))
summary(enfants)
boxplot(enfants)
mean(enfants)
var(enfants)
sd(enfants)
dev.off()
dev.new()
plot(table(enfants)/sum(table(enfants)))
lines(1:10,dgeom(1:10,1/(1+mean(enfants))),col="red")

X=rnorm(150,0,1);X 

# diagramme en batons (fréquences) + lines( dgeom ) (bon paramètre, p=1/(1+moyenne))