start on MAR data_prep.r script borrowing from NP script