If I understand correctly, rpart() will pick predictor at each node automatically. I am wondering if there is a way to force rpart() including a specific predictor. The reason I am asking is that I'd like to use rpart() to detect interaction terms for some variables. Thanks.