数据数据总共包含5年，解释变量为happiness，如果想区分高收入阶层与低收入阶层以及城镇居民和农村居民两组分别进行有序probit回归，还包括年龄、性别等控制变量，应该怎么做？求大神指点！！！