我现在有省份变量 有行业变量 也有个体的收入变量
比如

id   province  industry  income    meanincome
1     11             1            135
2     11             2            234
3     13             1             333


我现在想 求出每一个个体在他所在的省份所处的行业的收入均值
比如对id1来说，求出省份11 行业1的收入均值 然后写入meanincome位置

那么应该如何去编程呢？