***************************************************************;
* Wine 2_2                                                    *;
* "eclatement" de données correspondant à un traitement       *;
***************************************************************;
* Pierre Desmet 21/11/2010 ;

goptions devmap=winansi keymap=winansi ;    /* prise en compte des caractères accentués */

Title1 		H=2 j=L "Master Marketing" J=R "Universite Paris-Dauphine";
Title2 		H=1 j=L "Tests et Experimentation" ;
Title3		H=2 j=c "Cas R&W (2) " ;
Footnote1 	h=1 j=L "(c)2010, Pierre DESMET" ;
Data in ;
	input trait $ nb  ;
	do i=1 to nb ;
		input temps_restau @@;
		output ;
	end ;
Label temps_restau="Temps passé dans le restaurant" ;

cards ;
white 19
  80    98  115    89   103
  91   119   31   109    95
  71   105   66   141    79
 113    69  120    87    
red 23
  93   87   106   76   121   
 143   81    97   74   107
 112   67   107   72   116
  99  104    91  132    78
 107   101   92
;
run ;
proc print data=in ;
run ;
proc format ;
	value $ traitfmt "red"="Vin rouge"  "white"="Vin blanc" ; /* signaler un format alpha avec "$" */
run;
proc tabulate data=in ;
	class trait ;
	var temps_restau ;
	table (trait all),						/* lignes */
			N temps_restau*(mean std); 		/* colonnes */
	format trait $traitfmt. ;				/* signaler un format alpha avec "$" */
run ;

	ods graphics on ;						/* sorties graphiques automatiques */
proc univariate data=in normal;
	class trait ;
	var temps_restau ;
	histogram temps_restau / normal ;
run ;
proc ttest data=in ;
	class trait ;
	var temps_restau ;
run ;

proc glm data=in ;
	class trait ;
	model temps_restau= trait / SS3 solution;
	means trait / hovtest=levene ;
	contrast "effet type de vin" trait -1 1 ;
run ;
ods graphics off ;
***************************************************************;
* test de l'équivalence TOST ;
***************************************************************;
title3  "Test TOST de l'équivalence" ;
proc ttest data=in tost (-20, 20) ;
	class trait;
	var temps_restau ;
run ;