PRO UNI_STATS
device,decomposed=0
loadct,39

restore,'unitemp.sav'

data = read_ascii('results.csv',template=temp)

;SET_PLOT,'ps'

;DEVICE, /ENCAPSUL, BITS_PER_PIXEL=8, /LANDSCAPE,/COLOR,FILENAME='output.eps'

window,0,xsize=1024,ysize=800

!p.multi = [0,2,2,1,1]

plot,data.ucas,data.first,psym=1,color=0,background=255,yrange=[0,50],title='1st Class Degrees',xtitle='UCAS Points',ytitle='Percentage'
fit = LINFIT(data.ucas,data.first)
x = indgen(600)
fity = fit[0] + x*fit[1]
oplot,x,fity,psym=0,color=250
;print,fit[0],fit[1]
print,correlate(data.ucas,data.first)
plot,data.ucas,data.two_one,psym=1,color=0,background=255,yrange=[0,50],title='2.1 Degrees',xtitle='UCAS Points',ytitle='Percentage'
fit = LINFIT(data.ucas,data.two_one)
x = indgen(600)
fity = fit[0] + x*fit[1]
oplot,x,fity,psym=0,color=250
;print,fit[0],fit[1]
print,correlate(data.ucas,data.two_one)
plot,data.ucas,data.two_two,psym=1,color=0,background=255,yrange=[0,50],title='2.2 Degrees',xtitle='UCAS Points',ytitle='Percentage'
fit = LINFIT(data.ucas,data.two_two)
x = indgen(600)
fity = fit[0] + x*fit[1]
oplot,x,fity,psym=0,color=250
;print,fit[0],fit[1]
print,correlate(data.ucas,data.two_two)
plot,data.ucas,data.pass,psym=1,color=0,background=255,yrange=[0,50],title='Pass Degrees',xtitle='UCAS Points',ytitle='Percentage'
fit = LINFIT(data.ucas,data.pass)
x = indgen(600)
fity = fit[0] + x*fit[1]
oplot,x,fity,psym=0,color=250
;print,fit[0],fit[1]
print,correlate(data.ucas,data.pass)

meanf = mean(data.first)
sdf = stddev(data.first)
meanuc = mean(data.ucas)
sduc = stddev(data.ucas)

print,meanf,sdf,meanuc,sduc

;device,/close


END