楼主: 笑翠鸟
2086 1

[编程问题求助] stata编程显示invalid ·2' [推广有奖]

  • 4关注
  • 0粉丝

已卖:62份资源

大专生

85%

还不是VIP/贵宾

-

威望
0
论坛币
948 个
通用积分
0.0051
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
381 点
帖子
16
精华
0
在线时间
93 小时
注册时间
2013-5-16
最后登录
2020-3-31

楼主
笑翠鸟 学生认证  发表于 2017-2-24 13:32:57 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
初学stata,编写程序,运行显示 invalid ·2',求纠错
capture program drop nemenyi
program define nemenyi
version 12.0
local g =`1'
if (`g'== 3 | `g'== 4){
if `g' ==3{
local n1 = `2'
local n2 = `3'
local n3 = `4'
local Rs1 = `5'
local Rs2 = `6'
local Rs3 = `7'
local H =`8'
local Hc = `9'

local N = `n1'+`n2'+`n3'
local r1 = `Rs1' / `n1'
local r2 = `Rs2'/ `n2'
local r3 = `Rs3' / `n3'
local C = `H'/ `Hc'

local x12 =( `r1' -`r2' ) 2/( (`N' *(`N' +1) / 12) * ( 1/ `n1' +1/ `n2')*`C')
local x13 =( `r1' -`r3' ) 2/( (`N' *(`N' +1) / 12) * ( 1/ `n1' +1/ `n3')*`C')
local x23 =( `r2' -`r3' ) 2/( (`N'*(`N' +1) / 12) * ( 1/ `n2 ' +1/ `n3')*`C')
local p12 = chiprob( 2, `x12')
local p13 = chiprob( 2, `x13')
local p23 = chiprob( 2, `x23')

di newline(1)
di " Nemenyi test of 3 Groups"
di "--------------------------- "
di "Group chi2 P value "
di "--------------------------- "
di " 1 and 2" in ye %9. 4f `x12' in ye %9. 4f `p12'
di " 1 and 3" in ye %9. 4f `x13' in ye %9. 4f `p13'
di " 2 and 3" in ye %9. 4f `x23' in ye %9. 4f `p23'
di?--------------------------- "
}
if `g'==4{
local n1 = `2'
local n2 = `3'
local n3 = `4'
local n4 = `5'
local Rs1 = `6'
local Rs2 = `7'
local Rs3 = `8'
local Rs4 = `9'
local H = `10'
local Hc = `11'

local N = `n1'+ `n2' + `n3' + `n4'
local r1 = `Rs1'/ `n1'
local r2 = `Rs2'/ `n2'
local r3 = `Rs3'/ `n3'
local r4 = `Rs4' / `n4'
local C = `H' / `Hc'

local x12 =( `r1'-`r2') 2/( (`N'* (`N'+1) / 12) *( 1/ `n1'+1/ `n2') *`C')
local x13 =( `r1'-`r3') 2/( (`N'* (`N'+1) / 12) *( 1/ `n1'+1/ `n3') *`C')
local x14 =( `r1'-`r4') 2/( (`N'* (`N'+1) / 12) *( 1/ `n1'+1/ `n4') *`C')
local x23 =( `r2'-`r3') 2/( (`N'* (`N'+1) / 12) *( 1/ `n2'+1/ `n3') *`C')
local x24 =( `r2'-`r4') 2/( (`N'* (`N'+1) / 12) *( 1/ `n2'+1/ `n4') *`C')
local x34 =( `r3'-`r4') 2/( (`N'* (`N'+1) / 12) *( 1/ `n3'+1/ `n4') *`C')
local p12 = chiprob( 3, `x12')
local p13 = chiprob( 3, `x13')
local p14 = chiprob( 3, `x14')
local p23 = chiprob( 3, `x23')
local p24 = chiprob( 3, `x24')
local p34 = chiprob( 3, `x34')

di new line(1)
di "Nemenyi test of 4 Groups"
di "--------------------------- "
di "Group chi2 P value"
di "--------------------------- "
di " 1 and 2" in ye %9. 4f `x12' in ye %9. 4f `p12'
di " 1 and 3" in ye %9. 4f `x13' in ye %9. 4f `p13'
di " 1 and 4" in ye %9. 4f `x14' in ye %9. 4f `p14'
di " 2 and 3" in ye %9. 4f `x23' in ye %9. 4f `p23'
di " 2 and 4" in ye %9. 4f `x24' in ye %9. 4f `p24'
di " 3 and 4" in ye %9. 4f `x34' in ye %9. 4f `p34'
di "--------------------------- "
}
}
else{
di new line(1)
di "Error ! The number of groups is not equal to 3 or 4."
}
end
二维码

扫码加我 拉你入群

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

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

关键词:stata编程 invalid Stata Valid tata 编写程序 invalid local

沙发
夏目贵志 发表于 2017-2-24 23:18:12
把数据和你执行的时候用的命令也一起发上来

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群
GMT+8, 2026-1-5 05:01