马上试试代码！你简直是救星。需要的话马上就给币。大救星

不客气，有问题随时联系

logistic a37a a511 a512 a513 a514 a515 a516 a517  if a3a >= 1971
这是我的很幼稚的代码。R平方是0.0062  为什么你代码的第一部分的R平方是0 呢
方便加我微信吗 13828422627

重贴一下代码：
#1 基准模型 什么都不控制#
gene religion=0
replace religion=1 if a511==1|a512==1|a513==1|a514==1|a515==1|a516==1|a517==1|a518==1|a519==1|a520==1
gene children=a681+a682
ologit a37a religion if a3a>= 1971&children==1&a514!=1
oprob a37a religion if a3a>= 1971&children==1&a514!=1
#2控制变量 性别（a2） 年龄（a3a） 婚姻 （a69 以已婚（含同居）为参照，
分别考察“未婚”和“离异或丧偶）  民族（a4 以“非汉族为”为参照） #
gene single=0
replace single=1 if a69==1
replace single=. if a69==.
gene divorce=0
replace divorce=1 if a69==6|a69==7
replace divorce=. if a69==.
gene han=0
replace han=1 if a4==1
replace han=. if a4==.
ologit a37a religion a2 a3a single divorce han if a3a>= 1971&children==1&a514!=1
oprob a37a religion a2 a3a single divorce han if a3a>= 1971&children==1&a514!=1
#3 上一条的所有控制变量，另控制 政治面貌 （a10 以“不是党员不是团员”为对照）
户口性质（a18 以“农业户口”为对照） 教育年限 （a7a） 家庭收入 （a62） 个人收入
（a8a）#
gene party=0
replace party=1 if a10==1|a10==2|a10==3
replace party=. if a10==.
gene regis=1
replace regis=0 if a18==1
replace regis=. if a18==.
ologit a37a religion a2 a3a single divorce han party regis a7a a62 a8a if a3a>= 1971&children==1&a514!=1
oprob a37a religion a2 a3a single divorce han  party regis a7a a62 a8a if a3a>= 1971&children==1&a514!=1
#4 上一条的所有控制变量，另控制 现只有一个男孩（数据a681等于1 并且a682等于0）#
gene boy=0
replace boy=1 if a681==1&a682==0
replace boy=. if a681==.|a682==.
ologit a37a religion a2 a3a single divorce han party regis a7a a62 a8a boy if a3a>= 1971&children==1&a514!=1
oprob a37a religion a2 a3a single divorce han  party regis a7a a62 a8a boy if a3a>= 1971&children==1&a514!=1