summary.pp<-function(pp1,pp2=NULL,method="MAP",cutoff=0.5)
  {
    tau0<-cutoff
    if(method=="MAP")
      {
        ind1<-pp1>cutoff
        if(length(pp2)>0)
          {
            ind2<-pp2>cutoff
          }
      }
    else if(method=="FDR")
      {
        tau0<-0.01
        ind1<-pp1>tau0
        
        while(sum(ind1)>0 && (sum(ind1*(1-pp1))/sum(ind1)>cutoff))
          {
            tau0<-tau0+0.01
            ind1<-pp1>tau0
          }
        
        if(length(pp2)>0)
          {
            tau0<-0.01
            ind2<-pp2>tau0
            while(sum(ind2)>0 && sum(ind2*(1-pp2))/sum(ind2)>cutoff)
              {
                tau0<-tau0+0.01
                ind2<-pp2>tau0
              }
          }
      }
    
    list(ind1=ind1,tau0=tau0)
  }