STATA 中如何求得 LD50

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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(.) 是正态分布的累积概率函数。

谢谢！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Stata tata Probit trial Event Stata

. 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.

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

谢谢版大

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

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

再次谢谢版大……

请问LD50和ED50可以用同一种算法求得吗