最近跑multinomial-logit 分析的时候遇到一些问题，麻烦各位高手解答。

下面是分析data
dependent variables
employstatusF：已婚妇女就业型态,家庭主妇＝0,full-time work worker＝1,part-time worker＝2

independent variables
childno：小孩数
educyrsF：妇女受教育年数
age20_29,age30_39,age40_49,age50_59：妇女年龄（dummy，分吸时用age20_29为基准）
educyrsM：丈夫受教育年数
wrhrsM：丈夫工作一周时间（时）
fulltimeM：丈夫是full-time worker（dummy）
parttimeM：丈夫是part-time worker（dummy）
selfemloyedM：丈夫职业是selfemployed （dummy）
retireM：丈夫退休（dummy，分析时以selfemloyedM和retireM为基准）
incomeM：丈夫年收（1为低于70万日圆，2为70万日圆到100万日圆，3为1百万日圆到130万日圆…以此类推）
freqmeal：丈夫做家事频度
howi：对于丈夫出外赚钱，妇女在家工作的赞成程度（程度1～7，7为完全赞成）
helphw：丈夫应该帮忙做家事的赞成程度（程度1～7，7为完全赞成）

目的：想知道小孩数，丈夫的收入，年纪，做家事频率对已婚妇女就业选择的影响
以职业是家庭主妇的妇女为基准

stata分析结果：
bysort employstatusF:sum childno educyrsF  age20_29 age30_39 age40_49 age50_59 educyrsM wrhrsM  fulltimeM parttimeM  selfemloyedM retireM incomeM freqmeal howi helphw

mlogit  employstatusF childno educyrsF  ageF educyrsM wrhrsM fulltimeM parttimeM  incomeM freqmeal howi helphw,baseoutcome(0)
分析结果  http://imageshack.us/photo/my-images/402/30012419.jpg/

想问的问题有三个
1. 分析出来好像有bais，但是想不出来哪里有问题，希望各位帮忙看一下哪里有问题
2. incomeM（丈夫的收入）可以直接像这样放进去mlogit分析吗？还是要另外用其他方法？
3.mlogit分析完后要做说明的话是不是还要用margin再跑一次？

我是用下面的指令去跑，不知道有没有错误，麻烦指正，谢谢。
margins , dydx(childno educyrsF   age30_39 age40_49 age50_59 educyrsM wrhrsM fulltimeM parttimeM  incomeM freqmeal howi helphw) predict(outcome(0))
Average marginal effects                          Number of obs   =        281
Model VCE    : OIM
Expression   : Pr(employstatusF==0), predict(outcome(0))
dy/dx w.r.t. : childno educyrsF age30_39 age40_49 age50_59 educyrsM wrhrsM fulltimeM
               parttimeM incomeM freqmeal howi helphw


margins , dydx(childno educyrsF age30_39 age40_49 age50_59 educyrsM wrhrsM fulltimeM parttimeM  incomeM freqmeal howi helphw) predict(outcome(1))
Average marginal effects                          Number of obs   =        281
Model VCE    : OIM
Expression   : Pr(employstatusF==1), predict(outcome(1))
dy/dx w.r.t. : childno educyrsF age30_39 age40_49 age50_59 educyrsM wrhrsM fulltimeM
               parttimeM incomeM freqmeal howi helphw

margins , dydx(childno educyrsF   age30_39 age40_49 age50_59 educyrsM wrhrsM fulltimeM parttimeM  incomeM freqmeal howi helphw) predict(outcome(2))
Average marginal effects                          Number of obs   =        281
Model VCE    : OIM
Expression   : Pr(employstatusF==2), predict(outcome(2))
dy/dx w.r.t. : childno educyrsF age30_39 age40_49 age50_59 educyrsM wrhrsM fulltimeM
               parttimeM incomeM freqmeal howi helphw

麻烦各位高手帮忙，谢谢。