######################################## #TELECOM Saint-Etienne - année 2009/2010 #Cours n°13 #L. Carraro - O. Roustant ######################################## #graphique du test t avec les niveaux x<-seq(-4,4,by=0.05) y<-dt(x,df=17) plot(x,y,type="l",xlab="",ylab="",main="test t de niveau 5%",cex.main=2) seuil=qt(0.975,df=17) xright<-x[x>seuil] yright<-y[x>seuil] polygon(c(xright,rev(xright)), c(yright,rev(rep(0,length(yright)))), col = 'blue', border = "black" ) xleft<-x[x<(-seuil)] yleft<-y[x<(-seuil)] polygon(c(xleft,rev(xleft)), c(yleft,rev(rep(0,length(yleft)))), col = 'blue', border = "black" ) #graphique du test t avec la p-valeur #exemple avec t_obs=1.5 t_obs<-1.5 x<-seq(-4,4,by=0.05) y<-dt(x,df=17) plot(x,y,type="l",xlab="",ylab="",main=c("test t","p-valeur pour t_obs=1.5")) seuil=qt(0.975,df=17) xright<-x[x>t_obs] yright<-y[x>t_obs] polygon(c(xright,rev(xright)), c(yright,rev(rep(0,length(yright)))), col = 'blue', border = "black" ) xleft<-x[x<(-t_obs)] yleft<-y[x<(-t_obs)] polygon(c(xleft,rev(xleft)), c(yleft,rev(rep(0,length(yleft)))), col = 'blue', border = "black" ) p.valeur<-2*(1-pt(t_obs,17))