add plot4.R before readjustment