with R analysis