程序如何修改，或如何实现更好！

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

  比如rankext 3 6 结果有问题。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：如何实现 Rank ext RAN 程序 如何

这个程序的目的，可否简介一下？

sungmoo 发表于 2012-12-27 08:39
这个程序的目的，可否简介一下？

样本含量小的时候 计算配对符号秩检验的确切概率。

*上面的程序可以大幅减化（首数是样本量，末数是可能的秩和，求给定首数时末数对应的概率）
clear
cap pr drop rankext
pr rankext
args n e
if wordcount("`*'")!=2|int(`e')!=`e'|int(`n')!=`n'|`e'<0|`n'<0|`e'>`n'*(`n'+1)/2 {
n di as error "Two and only two nonnegative integers, of which"
n di as error "the last is smaller than the accumlated sum by the first, are allowed"
exit
}
qui{
clear
set ob `=2^`n''
forv i=1/`n'{
g v`i'=`i'*mod(int((_n-1)/2^(`n'-`i')),2)
}
egen exp=rowtotal(v*)
collapse (count)t=v1,by(exp)
egen p=pc(t),prop
su p if exp<=`e'
n di
n di in yellow "exact P-value for the ranksum of `e' given the sample size of `n' = " r(sum)
clear
}
end

*上面的程序可以大幅减化：

clear
cap pr drop rankext
pr rankext
args e n
if wordcount("`*'")!=2|int(`e')!=`e'|int(`n')!=`n'|`e'<0|`n'<0|`e'>`n' {
n di in red "Two and only two nonnegative integers, of  which the 2nd is larger, are allowed"
exit
}
clear
qui{
set ob `=2^`n''
forv i=1/`n'{
g v`i'=0
replace v`i'=`i' if mod(int((_n-1)/2^(`n'-`i')),2)
}
egen exp=rowtotal(v*)
collapse (count)t=v1,by(exp)
egen p=pc(t),prop
keep if exp==`e'
n di
n di in yellow "exact P Value for `e'/`n' = " p
}
clear
end
谢谢！

sungmoo 发表于 2012-12-28 23:41
*上面的程序可以大幅减化：

clear

运行 rankext 9 10.5 出错

运行 rankext 9 10.5 出错

上面的程序已经规定只许取非负整数。

sungmoo 发表于 2012-12-29 00:00
上面的程序已经规定只许取非负整数。

这样不完整啊！

这样不完整啊！

秩和可以是非整数？

sungmoo 发表于 2012-12-29 01:04
秩和可以是非整数？

可以，相同等级或取值取平均秩次，实际秩和中有可能出现0.5。