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