• make sure to use par(no.readonly = TRUE) TH.data