Observations of cosmic ray arrival directions made with the Pierre AugerObservatory have previously provided evidence of anisotropy at the 99% CL usingthe correlation of ultra high energy cosmic rays (UHECRs) with objects drawnfrom the Veron-Cetty Veron catalog. In this paper we report on the use of threecatalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3ptmethods, each giving a different measure of self-clustering in arrivaldirections, were tested on mock cosmic ray data sets to study the impacts ofsample size and magnetic smearing on their results, accounting for both angularand energy resolutions. If the sources of UHECRs follow the same large scalestructure as ordinary galaxies in the local Universe and if UHECRs aredeflected no more than a few degrees, a study of mock maps suggests that thesethree methods can efficiently respond to the resulting anisotropy with aP-value = 1.0% or smaller with data sets as few as 100 events. Using data takenfrom January 1, 2004 to July 31, 2010 we examined the 20, 30,..., 110 highestenergy events with a corresponding minimum energy threshold of about 49.3 EeV.The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. Inview of the multiple (correlated) scans performed on the data set, thesecatalog-independent methods do not yield strong evidence of anisotropy in thehighest energy cosmic rays.