real fra1(26),err1(26),fra2(26),err2(26),fra3(26), +err3(26),fra4(26),err4(26),agn_zspec(26) real sum1,sum2,sum3,sum4,mean1,mean2,mean3,mean4 c open(1,file='../data/final_agn_gal.tbl',status='old') open(1,file='../data1/final_agn_gal.tbl',status='old') sum1=0.0 sum2=0.0 sum3=0.0 sum4=0.0 do 10 i=1,26 read(1,101)agn_zspec(i),fra1(i),fra2(i),fra3(i), +fra4(i),err1(i),err2(i),err3(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. sum1=sum1+err1(i) sum2=sum1+err2(i) sum3=sum1+err3(i) sum4=sum1+err4(i) 10 continue mean1=sum1/26.0 mean2=sum2/26.0 mean3=sum3/26.0 mean4=sum4/26.0 c101 format(92x,f5.3,100x,8f5.2) 101 format(92x,f5.3,102x,8f5.2) c call pgbegin(0,'f_z_group1.eps/vcps',1,1) call pgbegin(0,'?',1,1) call pgpaper(0.0,0.258) call pgscf(2) call pgsls(1) call pgslw(4.0) 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,-10.,30.) call pgbox('bcnst',0.0,0,'bcnst',0.0,0) call pglabel('\fiz\d\frspec\u','\fif\d\frgroup\u(%)','') call pgsch(1.3) call pgptxt(-0.13,6.5,90.,0.,'AGN') call pgsch(2.0) call pgtext(0.2,-8.,'\fir\d\frs\u=0.094;') call pgtext(0.5,-8.,'slope=0.012') call pgtext(0.05,26.,'sample 1 (\fiP\d\frmem\u>0.1)') c call pgsch(1.7) call pgtext(0.05,-8.,'(a)') call fit(agn_zspec,fra1,26,err1,0,a,b,siga,sigb,chi2,q) print*, 'a=',a,'b=',b,'siga=',siga,'sigb=',sigb call pearsn(agn_zspec,fra1,26,r,prob,z00) print*,'r=',r,'prob=',prob do i=1,26 call pgpt1(agn_zspec(i),fra1(i),23) c call pgerry(1,agn_zspec(i),fra1(i)+err1(i),fra1(i)-err1(i),1.2) end do call pgerry(1,0.1,22.+0.5*mean1,22.-0.5*mean1,1.5) c call pgsch(1.7) call pgtext(0.13,22.,'mean error bar') 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,-10.,30.) call pgbox('bcnst',0.0,0,'bcnst',0.0,0) call pglabel('\fiz\d\frspec\u','\fif\d\frgroup\u(%)','') call pgsch(1.3) call pgptxt(-0.13,6.5,90.,0.,'AGN') call pgsch(2.0) call pgtext(0.2,-8.,'\fir\d\frs\u=0.098;') call pgtext(0.5,-8.,'slope=0.022') call pgtext(0.05,26.,'sample 2 (\fiP\d\frmem\u>0.3)') c call pgsch(1.7) call pgtext(0.05,-8.,'(b)') call fit(agn_zspec,fra2,26,err2,0,a2,b2, + siga2,sigb2,chi22,q2) print*, 'a2=',a2,'b2=',b2,'siga2=',siga2,'sigb2=' + ,sigb2,'chi22=',chi22 call pearsn(agn_zspec,fra2,26,r2,prob2,z002) print*,'r2=',r2,'prob2=',prob2 do i=1,26 call pgpt1(agn_zspec(i),fra2(i),23) c call pgerry(1,agn_zspec(i),fra2(i)+err2(i),fra2(i)-err2(i),1.2) end do call pgerry(1,0.1,22.+0.5*mean2,22.-0.5*mean2,1.5) c call pgsch(1.7) call pgtext(0.13,22.,'mean error bar') 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,-10.,30.) call pgbox('bcnst',0.0,0,'bcnst',0.0,0) call pglabel('\fiz\d\frspec\u','\fif\d\frgroup\u(%)','') call pgsch(1.3) call pgptxt(-0.13,6.5,90.,0.,'AGN') call pgsch(2.0) call pgtext(0.2,-8.,'\fir\d\frs\u=0.163;') call pgtext(0.5,-8.,'slope=0.050') call pgtext(0.05,26.,'sample 3 (\fiP\d\frmem\u>0.5)') c call pgsch(1.7) call pgtext(0.05,-8.,'(c)') call fit(agn_zspec,fra3,26,err3,0,a3,b3, + siga3,sigb3,chi23,q3) print*, 'a3=',a3,'b3=',b3,'siga3=',siga3,'sigb3=' + ,sigb3,'chi23=',chi23 call pearsn(agn_zspec,fra3,26,r3,prob3,z003) print*,'r3=',r3,'prob3=',prob3 do i=1,26 call pgpt1(agn_zspec(i),fra3(i),23) c call pgerry(1,agn_zspec(i),fra3(i)+err3(i),fra3(i)-err3(i),1.2) end do call pgerry(1,0.15,22.+0.5*mean3,22.-0.5*mean3,1.5) c call pgsch(1.7) call pgtext(0.18,22.,'mean error bar') 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) end