Stata 显示ranksum的结果
我正试图修改@Nick Cox在我的网站上发布的代码,但我有一些问题 我已经设置了varlist和group变量。我还更改了col选项以适合我的varname。我想为每组rN 1/rN 2添加观察次数,并在结果列表上添加一些标题 我正试图研究display命令,但无法找到解决方案 我的代码如下:Stata 显示ranksum的结果,stata,stata-macros,Stata,Stata Macros,我正试图修改@Nick Cox在我的网站上发布的代码,但我有一些问题 我已经设置了varlist和group变量。我还更改了col选项以适合我的varname。我想为每组rN 1/rN 2添加观察次数,并在结果列表上添加一些标题 我正试图研究display命令,但无法找到解决方案 我的代码如下: foreach v of var BVCAlogMAR AvgSSI L1DensityWholeImage { quietly ranksum `v', by(G6PDcarente) po
foreach v of var BVCAlogMAR AvgSSI L1DensityWholeImage {
quietly ranksum `v', by(G6PDcarente) porder
scalar pval = 2*normprob(-abs(r(z)))
di "`v'{col 34}" %05.3f pval " " %6.4e pval " " %05.3f r(porder) ///
" " %05.3f r(N 1) " " %05.3f r(N 2)
}
我无法在结果列表中输入rN 1和rN 2的值。此外,我不知道如何显示具有以下标题的列标题:
p值,PValue2,Porder-Observ。第1组-观察。第2组
你能帮帮我吗 您错误地将rN_1和rN_2分别称为rN 1和rN 2 在您上一篇文章中提供的玩具示例中,我已更正了此错误,并插入了六行粗体字,实现了您的目标: sysuse auto, clear local i = 0 foreach v of var mpg price weight length displacement { local ++i quietly ranksum `v', by(foreign) porder scalar pval = 2*normprob(-abs(r(z))) if `i' == 1 { display %20s "P-Value", %5s "PValue2", %5s "Porder1", %5s "group 1", %5s "group 2" display "" } display "`v'{col 14}" %05.3f pval " " %6.4e pval " " %05.3f r(porder) /// " " %05.3f r(N_1) " " %05.3f r(N_2) }
有关使用display命令设置输出格式的详细信息,请在Stata的命令提示符下键入help format。请发布您的代码。特别是你尝试过的和不起作用的。谢谢!它是完美的。我尝试使用display命令,但对我来说太复杂了。
P-Value PValue2 Porder1 group 1 group 2
mpg 0.002 1.9e-03 0.271 52.000 22.000
price 0.298 3.0e-01 0.423 52.000 22.000
weight 0.000 3.8e-07 0.875 52.000 22.000
length 0.000 9.4e-07 0.862 52.000 22.000
displacement 0.000 1.1e-08 0.921 52.000 22.000