求帮忙解答一段回归分析程序

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

全部程序如下，目标是看一个遗传位点SNP与结局变量group（1=有病，0=没病）之间的关系.
已知exposure是暴露因素，SNP与exposure强相关。

reg exposure SNP
predict yhat, xb
logit group yhat  , robust

rename yhat y2hold
rename exposure yhat
predict double res, residual
rename yhat exposure                  
rename y2hold yhat  
replace res = res^2  
summarize res

scalar realmse = r(mean)*r(N)/e(df_r)                              
matrix bmatrix = e(b)
matrix Vmatrix = e(V)
matrix Vmatrix = e(V) * realmse / e(rmse)^2
ereturn post bmatrix Vmatrix, noclear
ereturn display

运行程序之后，出错在下面一步：
. ereturn post bmatrix Vmatrix, noclear
estimates post: matrix has missing values
r(504);

求高手帮忙解答

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：回归分析 分析程序 求帮忙 Estimates summarize 回归分析 程序

如能解答，小弟愿意奉送论坛币，不太会操作，到时候直接送给解答的人还是怎么操作啊

错误信息是矩阵里有缺失值。
你得发数据来，不然没法帮你试。

夏目贵志 发表于 2015-10-16 22:55
错误信息是矩阵里有缺失值。
你得发数据来，不然没法帮你试。

好的，能否留个邮箱，我发你邮箱如何？谢谢啦

夏目贵志 发表于 2015-10-16 22:55
错误信息是矩阵里有缺失值。
你得发数据来，不然没法帮你试。

datafor_renda.xlsx (35.56 KB)

2015-10-17 09:41:51 上传

我分析的数据

不知道是否上传成功，还请牛人帮忙搞搞

麦迪他哥哥 发表于 2015-10-17 09:42
不知道是否上传成功，还请牛人帮忙搞搞

你的程序的问题是logit模型回归之后不会设置e(df_r)值，即
scalar realmse = r(mean)*r(N)/e(df_r)      
会生成缺失值，导致之后Vmatrix缺失。

夏目贵志 发表于 2015-10-17 09:55
你的程序的问题是logit模型回归之后不会设置e(df_r)值，即
scalar realmse = r(mean)*r(N)/e(df_r)      ...

请教大牛如何解决呢？小弟平时用sas多点，也是为了解决一个分析才用stata，纯属新人
求大牛帮忙指点下。谢谢啦

麦迪他哥哥 发表于 2015-10-17 17:26
请教大牛如何解决呢？小弟平时用sas多点，也是为了解决一个分析才用stata，纯属新人
求大牛帮忙指点下。 ...

你得知道写程序的人他想通过这个程序干什么。然后根据这个目的自己算出相当于e(df_r)的值。

夏目贵志 发表于 2015-10-17 22:43
你得知道写程序的人他想通过这个程序干什么。然后根据这个目的自己算出相当于e(df_r)的值。

嗯，是这样：目标是看一个遗传位点SNP与结局变量group（1=有病，0=没病）之间的关系.
已知exposure是暴露因素，SNP与exposure强相关。

本来是做2sls回归的，但是因变量是二分类变量，因此用ivregress 做不了，所以单独自己来算。先算第一阶段回归的预测值，再用预测值来进行第二阶段的logistic回归。

但是这样得到的标准误差有偏差，也就是估计的sd，因此有人建议用上面的程序进行调整这个sd。

请问这样清楚了吗？


请牛人帮忙解答下。

麦迪他哥哥 发表于 2015-10-18 08:47
嗯，是这样：目标是看一个遗传位点SNP与结局变量group（1=有病，0=没病）之间的关系.
已知exposure是暴 ...

你这个用ivprobit不就能直接搞定了吗？