全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4068 4
2008-06-26
LD50 是生物实验中常见的毒性计量标准,常用 probit 模型获得,
不过我看了 stata 中 probit 的帮助,还是弄不明白该如何计算。数据如下:

我知道在 SAS 中可以使用 proc probit 来求解,不知道 stata 中应该如何做。

dose event trial
0.0 3  15
1.1 4  15
1.3 4  15
2.0 3  15
2.2 5  15
2.8 4  15
3.7 5  15
3.9 9  15
4.4 8  15
4.8 11 15
5.9 12 15
6.8 13 15

其中因变量应该是 event/trial,自变量是 dose 的对数值。模型表达式为:

p = Pr(Y = 0) = C + (1 - C)F(x'\beta)

其中 F(.) 是正态分布的累积概率函数。

谢谢!


二维码

扫码加我 拉你入群

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

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

全部回复
2008-6-26 10:55:00


. gen y= event/trial

. gen x=ln( dose)
(1 missing value generated)

. probit y x

outcome does not vary; remember:
                                  0 = negative outcome,
        all other nonmissing values = positive outcome
r(2000);

Since no observation has a y-vlaue of zero, probit model cannot be built here.

二维码

扫码加我 拉你入群

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

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

2008-6-26 12:18:00

clear

input dose event trial
0.0 3  15
1.1 4  15
1.3 4  15
2.0 3  15
2.2 5  15
2.8 4  15
3.7 5  15
3.9 9  15
4.4 8  15
4.8 11 15
5.9 12 15
6.8 13 15
end

bprobit  event trial dose
gprobit  event trial dose

blogit  event trial dose
blogit  event trial dose,or

glogit  event trial dose
glogit  event trial dose,or

二维码

扫码加我 拉你入群

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

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

2008-6-26 14:14:00
谢谢版大

一开始没有注意到这些命令。

现在的问题是,这些命令在某些情况下,比如,在实验中,对照组没有致死率的情况下,是可以的。但如果对照组有死亡率时,似乎就有问题了,因为它不能给出我在第一个帖子中提到的 C。

再次谢谢版大……
二维码

扫码加我 拉你入群

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

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

2020-5-12 17:14:41
请问LD50和ED50可以用同一种算法求得吗
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群