公卫人

 找回密码
 立即注册

QQ登录

只需一步,快速开始

不劳无获:如何获取钢镚? 因为论坛,所以相逢。 捐赠百科答题至尊

公卫考研:一起风雨兼程 因为梦想,所以努力。 真题答案政治英语

职称考试:诸君逢考必过 因为热爱,所以执着。 模拟考场技能执医中级

查看: 3715|回复: 8

[分享] LD50计算程序

[复制链接]
admin 发表于 2007-5-9 19:59:52 | 显示全部楼层 |阅读模式

注册后推荐绑定QQ,之后方才可以使用下方的“用QQ帐号登录”。

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
有关IC50,LC50,LC50,ED50的讨论IC50, z; M$ w3 x# K/ C! r
(From Wikipedia, the free encyclopedia)" T8 q/ s$ a1 `7 G0 w: d
IC50, the half maximal inhibitory concentration, represents the concentration of an inhibitor that is required for 50% inhibition of things like an enzyme, a cell, a cell receptor or a microorganism.5 ]  E6 {& ?" P$ C: h

$ N* X% {) r3 J) h, nIt's commonly used as a measure of drug effectiveness. Another measure of drug efficacy is EC50 which represents the concentration of a compound that is required to obtain 50% of the maximum effect. (Ex. EC100 would be the minimum concentration of a compound to obtain 100% of its effect)." q/ t. X/ z  F* T7 f( j

, {" N1 d9 q4 o+ J9 x* x4 o0 E8 h----------------------------------------------------------, w! E9 i$ {8 a' X+ }( }1 ?8 G) S
ED50: `8 l! l) S7 }, h6 i5 k
(From Wikipedia, the free encyclopedia)- d" g1 L1 ^- C3 I# g
ED 50 (European Datum 1950) is a geodetic datum which was defined after World War II for the international connection of geodetic networks.- `: L9 z7 ^( i7 Z* M: i, D9 q0 D

, t9 o8 y9 ?6 `: J: [/ W) MSome of the important battles of World War II were fought on the borders of Germany, the Netherlands, Belgium and France, and the mapping of these countries had incompatible latitude and longitude positioning. This led to the setting up of ED50 as a consistent mapping datum for much of Western Europe. It was, and still is, used in much of Western Europe apart from Great Britain, Ireland, Sweden and Switzerland, which have their own datums.
% l4 S$ y0 G; F8 a- H7 y7 o+ v9 j! M$ ^/ N* W; m& V0 ^7 h
It was based on the international Ellipsoid of 1909 resp. 1924 ("Hayford"-Ellipsoid) ( radius of the Earth's equator 6378,388 km, flattening 1:297) and widely used all over the world up to the 1980s, when GRS80 and WGS84 were established.
' P; I! o8 g* L$ ^3 v3 a
! K) c  O' q3 y6 `- i$ zMany national coordinate systems of Gauss-Krüger are defined by ED50 and oriented by means of Geodetic Astronomy. Up to now it has been used in data bases of gravity field, cadastre, small surveying networks in Europe and America, and by some developing countries with no modern baselines.2 Y! t6 E0 Z. B8 K' Q9 a6 c
, K0 u* b' ?$ C  g$ _
The geodetic datum of ED50 is centred at the Frauenkirche of Munich in Southern Germany, where the approximate centre of the Western Europe national networks was situated in the years of the cold war. ED50 was also part of the fundamentals of the NATO coordinates (Gauss Krüger and UTM) up to the 1980s.1 F0 ^% H! R- T- E- G: Q
7 {) o$ v9 t- o- r( A8 t, v0 b

  G" p1 t/ L# O: y: G& G1 c0 PDatum shift between ED50 and WGS84
2 u! D$ }( V6 u8 q0 iThe longitude and latitude lines on the two datums are the same in the Archangel region of north-west Russia. As one moves westwards across Europe, the longitude lines on ED50 gradually become further west than their WGS84 equivalents, and are around 100 metres west in Spain and Portugal. Moving southwards, the latitude lines on ED50 gradually become further south than the WGS84 lines, and are around 100 m south in the Mediterranean Sea. (NB. If the lines are further west, the longitude value of any given point becomes more easterly. Similarly, if the lines are south, the values become northerly.)
2 D! o2 ~0 Z7 C0 s( G3 [$ Z! M1 o9 i, H
The datum shift for the Universal Transverse Mercator grid is different. The eastings vary between 0 m in Eastern Europe and 100 m in the far west of the continent, roughly similar to the longitude shift, but the northings are about 200 m different across Europe, with the ED50 UTM northing lines south of the WGS84 UTM northing lines, as UTM northings are measured from the Equator, and the theoretical ED50 Equator is about 200 m south of the WGS84 one.
, e$ j0 p6 m7 w8 m-----------------------------------------! c# F& W$ T1 N
LD50
% j* w' C! C  n  c& h/ M7 nFrom Wikipedia, the free encyclopedia! M$ k+ I6 Y: R0 {/ b) Z
Jump to: navigation, search
- K& P: m7 e. o/ X- }( N$ O9 ^For other uses, see LD50 (disambiguation).
5 T1 |- I  u8 i2 }7 f' T1 N6 x' L" }: q* W
An LD50 test being administeredIn toxicology, the LD50 or colloquially semilethal dose of a particular substance is a measure of how much constitutes a lethal dose. The test was created by J.W. Trevan in 1927. In animal testing studies, and those of other organisms, the dose administered that kills half the test population is referred to as the LD50, for "Lethal Dose, 50%".
" d! s% s+ T5 Z2 ~' A
( j/ k/ [7 @, wAnimal-rights and animal-welfare groups have campaigned against LD50 testing on animals in particular, citing it as cruel and unnecessary, and in the case of some substances, causing the animals to die slow, painful deaths. Several countries, including the UK, have taken steps to ban the oral LD50, and the Organization for Economic Co-operation and Development (OECD) abolished the requirement for the oral test in 2001 (see Test Guideline 401, Trends in Pharmacological Sciences Vol 22, February 22, 2001).* U- b, ~$ q7 h3 C# U
----------------------------------------------------------. j3 d0 z! Y+ x! c( ^- \2 d" w- I( K. E
ASTM American Standard of Testing and Materialsls / |3 a; L' r* f4 k3 o/ Z
9 ~2 i' G/ a6 F2 P6 i
AFNOR French Association of Standardization 1 R* U/ c$ }* P' E8 {% |& e
1 k- j. I( I. \# J/ B2 E  U
bioindicator a bioindicator indicates biological or animal species which, because of their ecological characteristics, constitute the early index of biotic or abiotic modifications of the environment due to human activities (ex: lichen).
) o% {: ]5 t! Z5 N  H
. x5 R! A$ _5 x8 ?' c6 }& DBiomarquor a biomarquor is a parameter (bio)chimic whose measurement reflects the interaction between a biological system and an environmental agent. # p, m! A) S3 B& @1 o

" w) C% e  s, x( c4 T9 s' D- IDIN Deutsches Institut für Normung
! K) m# Z' k# J' ]
) w# r0 X- K- z9 C5 XECx Concentration having an effect on x% of the population compared to control
5 m2 E, l- Z! n/ s( Z2 }
, j0 {4 P7 y2 q! _1 qEC50 " Effect Concentration 50% " or concentration having 50% of effect compared to control
1 E' Y9 a  Z; L! A4 S
5 |$ w" R3 \" |2 y' }IC50 " Inhibition Concentration 50% " or concentration inhibiting of 50% a parameter like the growth or luminescence compared to control : ^3 J% u, a* Z) }4 F
$ n0 j! T. G& J( E
ISO International Standardization Organization
, \, o9 D& r9 Y' x
" t9 F6 W( o: \# O* X! x9 VLCx concentration having caused the death of X% of the population tested compared to control
! d# t4 ]% g* Q$ M
) _0 _& c$ x: r* ^0 L2 uLC50 concentration having caused the death of 50% of the population tested compared to control
+ R  l5 [$ B7 ]# C+ f8 S: k' l6 u, v- z! k! K6 ]6 S1 l- W
LD50 " Lethal Proportions 50% " or proportions having caused the death of 50% of the population tested.
; {4 u- k" Q- Q/ C/ [- A
. y0 I: {5 ]  I: E& J; e  f( {LOEC " Low Observed Effect Concentration " or concentration with a very weak observable effect. It is about the lowest concentration for which the effect is different from that of controls. It is the first concentration tested after the NOEC.
2 b1 z4 t: K7 X, w% B  t3 ^3 n) s7 x/ f9 x% D
NOEC " No Observed Effect Concentration " or concentration without observable effect. It is about the highest concentration tested for which the effect is significantly not different from controls
( `) j0 j7 R  b3 N- I' C. C
; J- U5 ]: A$ w9 B' _OECD Organisation for Economic Cooperation and Developement
; J: U2 Z1 E0 ~% g  s
+ l4 ?: S. Q5 A3 rToxicity Self capacity of a substance to cause harmful effects on living organisms # e# N: w' s; H

5 t" {* V3 B( _) y  ~: I
4 c: ^6 ~! Z! o) L; W1 y* S: U----------------------------------------------------
9 ?! P4 [* V( c7 H9 A* h' \: T7 Y其实IC50的计算包括很多方法,有正态概率转换法,有寇氏法,有序贯法,还有最为准确的加权直线法。不同的方法都对实验设计有不同要求的,要根据具体情况而使用。6 ~* @* m7 z4 b: O

( K# [) T( j$ g( x/ z- S+ H$ ]我习惯使用寇氏法,下面的帖子中有算法,excel可以比较方便的作出(其实也挺折腾):
$ s0 g% m( _9 M9 s( Q2 ohttp://www.dxy.cn/bbs/post/view? ... =1&age=30#96758
8 F4 I! |0 e# u9 V8 b- B+ C序贯法是用于急性实验的,我没用过,但如果掌握了原理,应该用excel也可以做。* i* u1 R6 k* _; d7 K4 ?, S/ ^
正态概率转换法可以用SPSS做,用regression中的probit程序即可,挺方便。2 b+ y4 I8 P/ E1 T. D9 C
加权直线法最准确,但数学公式极其复杂,我也没深入研究过,tboy兄最近好象在开发软件,是不是应该由tboy兄完成这个mission:impossible?2 [6 ]# r! A1 D7 K$ f
http://www.dxy.cn/bbs/post/view? ... amp;tpg=2&age=0 # n! ?. v# f2 J$ I: [

. S) _# s  l# A* j----------------------------------
8 {' r6 m& P8 z- a2 z但是对于您的实验来说,应当具体情况具体分析,因为不同方法的要求有所差异,对于我来说寇氏改良法比较适合,而IC50计算软件尽管方便,但是没有人提供该软件的确切理论依据,我不敢使用。还有,个人认为做直线回归求IC50不太好,因为实验结果并非总是如您所料,有时候很难做出一条满意的直线来。 ) C" K) P9 C5 \. _/ V
转载请注明:转载自《散人笔记》[http://www.eryi.org/blog/]
 楼主| admin 发表于 2007-5-9 20:01:01 | 显示全部楼层

自编LD50计算程序Excel版

本帖最后由 admin 于 2010-8-11 21:29 编辑
3 C  N# z' w& k  \
8 e) ~/ y) N0 x看到有不少战友没有合适的LD50计算软件,将以前编的一个小程序与大家分享。
5 g7 u( x/ X- {LD50、IC50、ED50等XX50半数量计算方法都相同,只是意义不同,本程序都可以算。算法为Bliss正规法,参考孙瑞元主编《定量药理学》一书相应部分,计算结果同NDST基本一致。
5 S2 N6 a& a5 q2 T: J本程序采用Excel的宏(VBA)编写,大量应用Execl提供的函数,故程序简短,可读性强,有爱好编程的朋友可以参考。# ^2 l' I- G7 B% T; S. Y6 Y
用法:打开表格,选择启用宏,任意位置输入相邻的三列数,鼠标选中,选择菜单“工具”-宏-宏(或按alt+F8),执行main即可。熟练的朋友可在工具栏上添加一个按钮,指向main,使用方便。$ @8 E) L9 R/ _4 R1 s  F4 X1 g
第一列:剂量(实际值,不用对数转换);第二列:每组总动物数;第三列:死亡/有效动物数。如:
+ _( N( z4 j6 j/ M: J+ p! V+ a0 o350  20  20
$ U) |5 s4 {( x, n7 h245  20  17
9 k1 P' Q) K/ p$ e; E. f) a! _172  20  16
  x' X5 ?2 g: {, c; }120  20  62 Z$ g* }$ T& t, X+ U+ l5 Z. G
84  20  2, p; F1 P  k( R# ]$ l9 s9 Y
59  20  0. e( e, e& f4 w( R

' `% Y" k; Z9 J( ], j5 GCITED FROM shex @DXY

LD.xls

104 KB, 下载次数: 30

回复

使用道具 举报

 楼主| admin 发表于 2007-5-9 20:03:33 | 显示全部楼层

LD50计算程序 SPSS软件的方法

关于LD50,SPSS软件的方法比较简单,具体方法介绍如下:4 y* j" N# o: G& k. g# z
1 SPSS——Analyze > regression>Probit…>Probit Analysis对话框:4 D: m, P9 D' r! E8 t9 X6 e
2 死亡数——Response Frequency框; 动物数——Totol Observed框; 剂量变量(死亡数)——Covariate框。
, b! S' S3 Y3 p" u9 K8 Z+ G& f3 Transform下拉菜单:选择Log base 10(以10为底的对数转换
: W# T( C8 `' s) M4 Model栏中:Probit 为概率单位模型
( M$ B4 o5 j( [6 r( @6 l- v3 p5 点击options…钮——Probit Analysis: Options对话框:
/ S7 G' S4 R7 I5 k1 {2 D选择 Calculate from data:依据现有数据所得的在缺少刺激条件下的响应率来估计自然响应率;3 d4 }: g( X( U$ C/ z
其他默认设置
9 |3 ^) k& d8 e  C1 R7 }这样点几下便可求出LD50。4 b. r6 @8 u( s. c. e
仅供参考,如有其他方法请相互交流。 0 f; W! F" a6 ~% I' D5 O: w
6 C" s4 A9 U; z  r# U8 K
CITED FROM 思痛@DXY
回复

使用道具 举报

 楼主| admin 发表于 2007-5-9 20:05:54 | 显示全部楼层

LD50计算程序 统计软件BLISS

本帖最后由 admin 于 2010-8-11 21:30 编辑
8 s0 s5 B. `; o+ e& ?
% ~( \, Z5 x' v% B. T统计软件BLISS,在DOS和win2000下均能使用。
; j' y: h% O# K1 G/ R可计算LD50和95%可信限,使用方便。
& F$ h2 q& d& @1 |9 `( V! H4 @( l! }
我的程序使用方法:
$ v' W; _. f4 S& _( d运行BLISS后,点1-input输入数据:% L" t6 t! ]- q3 U7 N. }
第一列为剂量,第二列为死亡数,第三列为动物总数。
2 \/ W$ K) y# q$ T7 f  N- s" ?9 i) V输入过程要仔细,不要弄错了哦。输入完毕,点F1退出,它会询问你是否继续编辑,如果数据无误,可点n,然后回车保存数据。
4 y" l# _, d7 Z$ y(注意序贯法适合要求各组动物数量必须一致,Bliss法则没有这样的要求)
6 X- \1 c0 |3 Z- D$ @点2-computing进行曲线拟和以及运算
9 i0 @" S4 w# G" L5 ?5 G  c9 i点3-dispaly 显示统计结果! t% K4 F0 d3 \3 V# R1 ^2 w
点4-graphics为显示拟和的曲线
3 \: W$ ~& {& J; ]  N) I2 ~$ \+ M点5-print为打印+ {# y  S. u2 V& o5 D" \
===========
* [  O! z9 W% h本软件的缺点在于不能将结果输入粘贴板,以便整理到word中。

BLISS.rar

126.06 KB, 下载次数: 19

回复

使用道具 举报

 楼主| admin 发表于 2007-5-9 20:10:19 | 显示全部楼层

Re:LD50计算程序 其他程序

本帖最后由 admin 于 2010-8-11 21:35 编辑 ' d9 b5 E1 {; C5 d

9 A1 x3 ^' b: P小软件汇总,有空自己也编写一个,赫赫
回复

使用道具 举报

 楼主| admin 发表于 2007-5-19 16:37:00 | 显示全部楼层

一种简单、精确计算EC50的方法

本帖最后由 admin 于 2010-8-11 21:34 编辑 : m: p" i, G8 a5 x5 C3 @" t
. ~& l6 Q* Z( K$ ]6 f
李晶晶  杨 璞  叶方立
. n4 X, j, X+ R- l(武汉科技大学医学院 武汉430062)6 J/ f% h* S9 W
摘 要:  提出了一种计算药物量- 效曲线参数半数效量(EC50) 的新方法, 可精确计算落在两已知浓度间的EC50, 不需绘图测
4 H, u4 W( y) H3 b3 @算或计算机曲线拟合。该方法比常用的绘图计算的方法更加快速、精确和简单, 能广泛应用于药效动力学研究。
. H; n2 j" ~4 R关键词:  EC50;  量- 效曲线

一种简单、精确计算EC50的方法.pdf

259.26 KB, 下载次数: 11

回复

使用道具 举报

 楼主| admin 发表于 2010-8-11 21:37:43 | 显示全部楼层
附件已修复:)
回复

使用道具 举报

cruise0303 发表于 2017-9-5 21:03:04 | 显示全部楼层
好东西 谢谢!
回复

使用道具 举报

harrybzy 发表于 2018-9-13 13:04:15 | 显示全部楼层
谢谢,感谢楼主奉献
回复

使用道具 举报

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

本版积分规则

提现|充值|关于|接种|公卫人 ( 沪ICP备06060850号-3 )

GMT+8, 2018-9-21 05:19 , Processed in 0.346948 second(s), 36 queries , Gzip On.

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表