Analyses of primary remanence vector data from a large collection of lava flows: towards improved methodology in paleo-geomagnetism