Added R analysis script