200币|悬赏 求问 简单logistic问题（内附stata数据）宗教有否提升提高了二胎生育意愿

悬赏不够可继续追加。每条都奖论坛币！越快解决，奖励越高。

待解问题：宗教有否提升提高了二胎生育意，即考察“希望生育子女数”是否受宗教信仰的影响（附件为数据和问卷）
因变量： 希望生育子女数（数据a37a）
自变量： 有否宗教信仰（a511-a520有一个为“是”） 。
范围：   限定为年龄小于42岁（即数据a3a >= 1971），并且现只有一个孩子（数据a681加a682刚好等于1），排除所有伊斯兰教信仰者，即a514为“是”的均排除

请按照序号，写出stata的命令行，logistic 和 probit 都需要 ，并且显示双尾检验统计显著度P值

1 基准模型 什么都不控制
2 控制变量 性别（a2） 年龄（a3a） 婚姻 （a69 以已婚（含同居）为参照，分别考察“未婚”和“离异或丧偶）  民族（a4 以“非汉族为”为参照）
3 上一条的所有控制变量，另控制 政治面貌 （a10 以“不是党员不是团员”为对照）
户口性质（a18 以“农业户口”为对照） 教育年限 （a7a） 家庭收入 （a62） 个人收入 （a8a）
4 上一条的所有控制变量，另控制 现只有一个男孩（数据a681等于1 并且a682等于0）


附件链接如下
CGSS2013原始数据.rar
https://bbs.pinggu.org/a-2341698.html

#1 基准模型 什么都不控制#
gene cons=1
ologit a37a cons
oprob a37a cons
#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 a2 a3a single divorce han
oprob a37a a2 a3a single divorce han
#上一条的所有控制变量，另控制 政治面貌 （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 a2 a3a single divorce han party regis a7a a62 a8a
oprob a37a a2 a3a single divorce han  party regis a7a a62 a8a
#4 上一条的所有控制变量，另控制 现只有一个男孩（数据a681等于1 并且a682等于0）#
gene boy=0
replace boy=1 if a681==1&a682==0
replace boy=. if a681==.|a682==.
ologit a37a a2 a3a single divorce han party regis a7a a62 a8a boy
oprob a37a a2 a3a single divorce han  party regis a7a a62 a8a boy
我已经运行了~可以跑通，有什么不懂的继续问哈~

马上给附件哈~正在帮您写code~

马上给附件哈~正在帮您写code~

molicoolgirl 发表于 2017-10-17 08:56
马上给附件哈~正在帮您写code~

就是这个，论坛的链接， 不是病毒  https://bbs.pinggu.org/a-2341698.html
这个也行，https://bbs.pinggu.org/thread-6027009-1-1.html

你的因变量希望生育子女数是有序变量，所以需要用有序（ordered）logit和probit模型，在stata里面用ologit和oprob模型解决，你可以在stata里help一下ologit和oprob看一下具体的用法

molicoolgirl 发表于 2017-10-17 09:19
你的因变量希望生育子女数是有序变量，所以需要用有序（ordered）logit和probit模型，在stata里面用ologit和 ...

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

61603161 发表于 2017-10-17 09:22
马上试试代码！你简直是救星。需要的话马上就给币。大救星

不客气，有问题随时联系

molicoolgirl 发表于 2017-10-17 09:23
不客气，有问题随时联系

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