如何使用PROC GENMOD (或者其他PROC）来分析我的数据？

我想看看，消费者的购买行为与其教育年限，收入，性别，民族之间的关系，现有如下四个示意性质的方程：
方程1：Y苹果（0，1） - 教育年限，收入，性别，民族
方程2：Y樱桃（0，1） - 教育年限，收入，性别，民族
方程3：Y梨类（0，1） -教育年限，收入，性别，民族
方程4：Y桃类（0，1） - 教育年限，收入，性别，民族
假如，部分数据如下：

Fruit	Y-是否购买	X1-Education	X2-Income	X3-Gender	X4-Race
Apple	0	2	4	F	MEXICO
Apple	0	3	5	M	CHINA
Apple	1	4	6	F	MEXICO
Apple	1	7	2	M	CHINA
Cherry	0
Cherry	1	9	6	F	CANADA
Cherry	0	2	5	M	MEXICO
Cherry	1	9	4	F	JAPAN
Pear	0	2	4	M	CHINA
Pear	0	7	6	M	MEXICO
Pear	1	3	9	F	JAPAN
Pear	0	1	5	F	MEXICO
Peach	1				MEXICO
Peach	0	8	1	M	ENGLAND
Peach	0	4	2	F	MEXICO
Peach	1	2	4	M	ENGLAND

现在，假定每个方程均为PROBIT模型，需要4个方程联合估计，并且还有方程之间参数限制：四个方程的收入参数之和为1，四个方程的截距之和为0.
请问，如何使用PROC GENMOD进行建模，并且伴随参数限制。
或者提供其他PROC 建模代码，或者提供例子，均可，非常感谢。

问题补充：
如果四个方程按照PROBIT模型联合估计比较困难，逐个估计各个方程亦可；如果估计过程加入参数限制比较困难，忽略参数限制亦可；如果感觉PROC GENMOD比较复杂，推荐其它PROC 亦可。总之，为了达到我的目标，什么要求都可妥协，想要选择一个能够达到目标，并且最优的方案。

如果四种水果是互斥的话，可以用
proc Logistic；
model a=education gender income race / link=mlogic;
run;
这是多元logistic 回归，保证四种概率之和为1，但截距没注意过，楼主可以试一试
自变量要先变成dummy variable