/* comfort experiment MJ 5,30 */
options nocenter ps=64 ls=80;
data a;
   infile 'comfort.dat';
   input temp chamb m1 m2 f1 f2;
data b; set a;
   score = m1; sex = 'm'; person = 1; output;
   score = m2; sex = 'm'; person = 2; output;
   score = f1; sex = 'f'; person = 1; output;
   score = f2; sex = 'f'; person = 2; output;
   drop m1--f2;
proc glm;
   class temp chamb sex;
   model score = temp chamb(temp) sex temp*sex / ss1;
   test h=temp e=chamb(temp);
   lsmeans temp / pdiff stderr e=chamb(temp);
   lsmeans sex temp*sex / pdiff stderr;