Subsampling the artifical PNe data

Our goal here is to use the same number of PNe as Feldmeier/Okamura observed in each field and see if we can still detect the correlation between PNe and light. We need to bring our artifical catalogs down to the size of the observations.

Simply put, we randomly sampled out from the catalog the number of PNe found by Feldmeier/Okamura.

Code --> subs.sm  -  subs  -  loads catalogs of artifical PNe from PNe8alpha program
                                             sets desired sizes of catalogs [ 75 for 3.fits, 11 for 4.fits, 77 for 7.fits, and 34 for Sub.fits ]
                                             generates vector of random numbers the same length as our artifical data
                                             sorts our data in order of the random vector and takes the first however many values, 75, 11, 141, or 34
                                             concatenates the PNe data from each field into full catalog as it goes along
                                             writes PNeFAKE.dat with all RA and DEC


Examples of Selected Catalogs of PNe (compare to complete catalogs) overlaid on the fields:
[these were generated with conv.sm (program to get PNeFAKE.dat into PNeFAKE.reg that ds9 can read) and ds9.sm (quickie program to use command line ds9 to open images and overlay regions and save them, quickly]

3mm.fits with PNeFAKE overlay



4mm.fits with PNeFAKE overlay


7mm.fits with PNeFAKE overlay


Submm.fits with PNeFAKE overlay



Coremm.fits with PNeFAKE overlay



LPCmm.fits with PNeFAKE overlay



FCJmm.fits with PNeFAKE overlay




Back to main Virgo page