用SigmaPlot做一个多靶单击模型拟合剂量存活曲线,大家来讨论下怎么做拟合。
Sigmoidal里的Chapman 4 parameters的方程是f=y0+a*(1-exp(-b*x))^c,怎样将方程修改成多靶单击模型S=1-(1-e-D/Do)N
数据如
0
1
1
2
0.305
0.62
4
0.13
0.42
6
0.042
0.286
8
0.0012
0.138
10
0.0000001
0.022
12
0.00000001
0.00008
克隆接种率PE(plating efficiency)=(对照组克隆数/接种细胞数)×100%,存活分数SF(surviving fraction)=克隆数/(接种细胞数×PE)。
上面分别是0,2,4,6,8,10,12GY照射的SF值,要绘制对数图.
和算出简单的多靶单击模型S=1-(1-e-D/Do)N拟合剂量存活曲线,主要获得N,D0,Dq和SF2等多个放射生物学参数。S为剂量D对应的数值.
主要意义如下:目前应用较多的是fertil等建立的克隆形成率分析法,是用多个剂量的多个克隆存活分数进行数学模型拟合。可通过简单的多靶单击模型S=1-(1-e-D/Do)N拟合剂量存活曲线,主要获得N,D0,Dq和SF2等多个放射生物学参数。这些衡量放射敏感性的参数能较真实全面的反映细胞放射生物学特性,已被广泛应用于细胞放射敏感性分析。该种拟合获得的N、Dq、Do和SF2值,可从不同角度反应细胞放射敏感性:N值是曲线指数部分外延至Y轴的截距,称外推值,反映细胞对放射引起损伤的修复能力,N增大表示细胞修复能力增强,杀死细胞所需剂量增大,放射抗拒性增强;Dq为细胞受损所需准阈剂量,Dq增大,细胞存活曲线肩区增宽,放射抗拒性增强;Do为曲线指数区下降63%所需的剂量,Do越大,放射抗拒性越强;2Gy时存活分数SF2是代表细胞放射敏感性的重要指标,SF2越大,放射抗拒性越强。 式中S为受到剂量照射的存活分数;存活曲线的直线部分斜率的倒数为Do值,称为细胞的平均致死剂量,为使存活率下降63%所需剂量。若将直线部分外推与纵坐标相交点的数值称为外推值n,代表细胞内靶的个数或所需击中靶的次数。由纵坐标1.0处作一条与横坐标的平行线,与外推线的交点在横坐标上投影点的数值即为Dq,称准域剂量,代表细胞累积亚致死性损伤的能力,为克服肩区的剂量。
下载地址:
http://www.softsale.cn/demo/michael/多靶单击模型.rar