! re-analysis of incoming Y1 student data
GRAV[1]=`<70%'
GRAV[2]=`70-74%'
GRAV[3]=`75-79%'
GRAV[4]=`80-84%'
GRAV[5]=`85-89%'
GRAV[6]=`90-99%'
YEAR=[2005:2012]

BPHY[2005,[1:6]]={0;0;0;0;0;0}
CAST[2005,[1:6]]={0;0;0;0;0;0}
PHYS[2005,[1:6]]={0;0;1;0;0;0}
 FMS[2005,[1:6]]={3;42;61;60;56;23}

BPHY[2006,[1:6]]={0;1;0;0;1;0}
CAST[2006,[1:6]]={0;0;0;0;0;0}
PHYS[2006,[1:6]]={0;0;2;1;0;1}
 FMS[2006,[1:6]]={0;43;58;43;50;31}

BPHY[2007,[1:6]]={0;0;1;0;0;0}
CAST[2007,[1:6]]={0;0;0;0;0;0}
PHYS[2007,[1:6]]={0;5;2;2;0;0}
 FMS[2007,[1:6]]={2;72;89;70;47;27}

BPHY[2008,[1:6]]={0;1;0;0;0;0}
CAST[2008,[1:6]]={0;0;1;1;1;0}
PHYS[2008,[1:6]]={0;1;0;1;0;0}
 FMS[2008,[1:6]]={0;37;69;54;39;21}

BPHY[2009,[1:6]]={0;0;0;0;0;0}
CAST[2009,[1:6]]={0;0;0;0;0;0}
PHYS[2009,[1:6]]={0;3;3;0;3;1}
 FMS[2009,[1:6]]={1;65;73;54;55;23}

BPHY[2010,[1:6]]={0;0;0;0;0;0}
CAST[2010,[1:6]]={0;1;0;0;0;0}
PHYS[2010,[1:6]]={0;1;0;0;2;1}
 FMS[2010,[1:6]]={0;59;81;73;49;22}

BPHY[2011,[1:6]]={0;0;0;0;1;0}
CAST[2011,[1:6]]={0;0;0;0;0;0}
PHYS[2011,[1:6]]={0;0;2;2;0;1}
 FMS[2011,[1:6]]={0;70;106;81;56;37}

BPHY[2012,[1:6]]={0;0;0;0;1;2}
CAST[2012,[1:6]]={0;0;1;0;0;0}
PHYS[2012,[1:6]]={0;2;3;5;0;3}
 FMS[2012,[1:6]]={2;127;126;84;55;39}

PTTL[[2005:2012],[1:6]]=BPHY+CAST+PHYS

clear 
defaults

!window 5
!surface\hist PTTL 45 45
!window 7
!surface\hist FMS 45 45
orientation portrait

window 3
label\x `Average High School grade'
set font triumf.2
set %yuaxis 97
set %xlaxis 8
scales 65 97 4 0 6 3
scales 65 97 4 0 35 7
set xvmin 70
set xvmax 90
set nsyinc 2
set histyp 2
legend on
legend frame 10 70 25 90
legend title `Physics'
legend title `FMS %'
legend frame off
legend auto off

x={67.5;72.5;77.5;82.5;87.5;92.55}

do i=[1:len(YEAR)]
  set colour i
  set pchar 2+i 0.08
  stat FMS[2004+i,*] norm\sum
!  graph\noax\hist rchar(YEAR[i]) x-2+i*0.45,PTTL[2004+i,*]
  graph\noax\hist rchar(YEAR[i]) x-2+i*0.45,FMS[2004+i,*]/norm*100
enddo
set colour 1
graph\ax PTTL[2005,*]

window 4
terminal
scales 65 97 4 0 150 3
set xvmin 70
set xvmax 90
set nsyinc 5
legend on
legend frame 10 70 25 90
legend frame off
legend title `FMS'
legend title `FMS count'
legend auto off

do i=[1:len(YEAR)]
  set colour i
  set pchar 2+i 0.08
  graph\noax\hist rchar(YEAR[i]) x-2+i*0.45,FMS[2004+i,*]
enddo
set colour 1
graph\ax FMS[2005,*]