First analysis script, work01.R. Update to readData.R -spk