real fra1(10),err1(10),fra2(10),err2(10),fra3(10),
     +err3(10),fra4(10),err4(10),agn_zspec(10)

c        open(1,file='../data/agn_ratio_bin.tbl',status='old')
        open(1,file='../data1/agn_ratio_bin2.tbl',status='old')

        do 10 i=1,10
          read(1,101)agn_zspec(i),fra1(i),err1(i),fra2(i),err2(i),
     +fra3(i),err3(i),fra4(i),err4(i)
          fra1(i)=fra1(i)*100.
          fra2(i)=fra2(i)*100.
          fra3(i)=fra3(i)*100.
          fra4(i)=fra4(i)*100.
          err1(i)=err1(i)*100.
          err2(i)=err2(i)*100.
          err3(i)=err3(i)*100.
          err4(i)=err4(i)*100.
10     continue

        if(i.eq.1)print*,agn_zspec(i),fra1(i),err1(i),fra2(i),
     +err2(i),fra3(i),err3(i),fra4(i),err4(i)

101     format(f4.2,8e11.4)
   
************** begin to plot *******************

c      call pgbegin(0,'f_z_bin.eps/vcps',1,1)
      call pgbegin(0,'?',1,1)
      call pgpaper(0.0,0.258)
      call pgscf(2)
      call pgsls(1)
      call pgslw(4.)
      call pgsch(2.0) 
      call pgsci(1)

************ plot the p>0.1**********

        call pgvport(0.04,0.32,0.19,0.95)
        call pgwindow(0.0,1.1,-1.,3.)
        call pgbox('bcnst',0.0,0,'bcnst',0.0,0)
        call pglabel('\fiz','\fif\d\frAGN\u(%)','')
        call pgtext(0.05,2.4,'sample 1 (\fiP\d\frmem\u>0.1)')
        call pgtext(0.05,-0.8,'(a)') 
        call pgtext(0.25,-0.8,'\fir\d\frs\u=0.394;')
        call pgtext(0.55,-0.8,'slope=0.00356')

	call fit(agn_zspec,fra1,10,err,0,a,b,siga,sigb,chi2,q)
	print*, 'a=',a,'b=',b,'siga=',siga,'sigb=',sigb

	call pearsn(agn_zspec,fra1,10,r,prob,z00)
	print*,'r=',r,'prob=',prob

	
      do i=1,10
        call pgpt1(agn_zspec(i),fra1(i),23)
	call pgerry(1,agn_zspec(i),fra1(i)+err1(i),fra1(i)-err1(i),1.2)
	end do
            
	call pgmove(0.0,a+b*0.0)
	call pgdraw(1.2,a+b*1.2)

************ plot the p>0.3 **********

        call pgvport(0.37,0.65,0.19,0.95)
        call pgwindow(0.0,1.1,-1.,3.)
        call pgbox('bcnst',0.0,0,'bcnst',0.0,0)
        call pglabel('\fiz','\fif\d\frAGN\u(%)','')
        call pgtext(0.05,2.4,'sample 2 (\fiP\d\frmem\u>0.3)')
        call pgtext(0.05,-0.8,'(b)')
        call pgtext(0.25,-0.8,'r\ds\u=0.387;')
        call pgtext(0.55,-0.8,'slope=0.00416')

	call fit(agn_zspec,fra2,10,err2,0,a2,b2,
     +siga2,sigb2,chi22,q2)
	print*, 'a2=',a2,'b2=',b2,'siga2=',siga2,'sigb2=',
     +sigb2,'chi22=',chi22

	call pearsn(agn_zspec,fra2,10,r2,prob2,z002)
	print*,'r2=',r2,'prob2=',prob2

	do j=1,10
       call pgpt1(agn_zspec(j),fra2(j),23) 
        call pgerry(1,agn_zspec(j),fra2(j)+err2(j),fra2(j)-err2(j),1.2)
	end do

	call pgmove(0.0,a2+b2*0.0)
	call pgdraw(1.2,a2+b2*1.2)

************ plot the p>0.5 **********

        call pgvport(0.7,0.98,0.19,0.95)
        call pgwindow(0.0,1.1,-1.,3.)
        call pgbox('bcnst',0.0,0,'bcnst',0.0,0)
        call pglabel('\fiz','\fif\d\frAGN\u(%)','')
        call pgtext(0.05,2.4,'sample 3 (\fiP\d\frmem\u>0.5)')
        call pgtext(0.05,-0.8,'(c)')
        call pgtext(0.25,-0.8,'\fir\d\frs\u=0.359;')
        call pgtext(0.55,-0.8,'slope=0.00490')

        call fit(agn_zspec,fra3,10,err3,0,a3,b3,
     +siga3,sigb3,chi23,q3)
        print*, 'a3=',a3,'b3=',b3,'siga3=',siga3,'sigb3=',
     +sigb3,'chi23=',chi23

        call pearsn(agn_zspec,fra3,10,r3,prob3,z003)
        print*,'r3=',r3,'prob3=',prob3

        do j=1,10
        call pgpt1(agn_zspec(j),fra3(j),23) 
        call pgerry(1,agn_zspec(j),fra3(j)+err3(j),fra3(j)-err3(j),1.2)
        end do

	call pgmove(0.0,a3+b3*0.0)
	call pgdraw(1.2,a3+b3*1.2)

************ plot the p>0.7 **********


*************************************
        call pgend
        
        close(1);close(2);close(3);close(4)
        end