24小时热门版块排行榜    

查看: 10460  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

propamidine

铁杆木虫 (著名写手)


[求助] 关于使用SPSS进行EC50的计算问题 已有1人参与

本人之前进行杀菌剂EC50的计算,是根据师兄留下来的Excel表,照猫画虎的填数据计算,从未仔细思考过其计算原理或公式来源。
近日看了篇文献"J.-L. Li, X.-Y. Liu, J.-T. Xie, et al. A Comparison of Different Estimation Methods for Fungicide EC50 and EC95 Values. Journal of Phytopathology, 2015. 163(4):239-244."
这文章对常用的生测数据统计软件(如SPSS、DPS、GraphPad、Excel等)进行了比较,指出进行EC50计算和EC95时,Excel很不准确,应使用SPSS、DPS和GraphPad。
文献里对生测数据的统计也以SPSS为多,所以本人想用SPSS进行数据分析。
但搜不到计算杀菌剂EC50的中文文献,而多为杀虫剂LD50的文献,如:"齐军山, 辛志梅, 李林, 等. 应用SPSS软件进行农药试验数据的统计分析. 山东农业科学, 2008. 7:100-104."
按理说EC50和LD50的算法相似,大概知道怎么在SPSS输入数据进行计算了。
之后又查了不少资料,却越发糊涂。
(1)好像EC50的估算方法有很多种,见:"赵劲松, 于书霞. 利用概率单位法估算EC50的研究. 生态毒理学报, 2010. 420-425."
(2)Probit、ANOVA的关系和区别也搞的我很糊涂?
(3)外文文献的很多结果,就是毒力方程、EC50、95%置信限和相关系数。但我们老师说不要相关系数,要卡方值。卡方检验是否应该用适合性检验,以比较测定值与理论值是否符合,即毒力方程是否可信?

如果需要得到毒力方程、EC50、95%置信限、相关系数和卡方值,在SPSS进行数据处理是否分为不同部分进行?
哪位虫友可热心制作一份操作指南?定重金酬谢。

附上几组校正抑制率数据:
药剂浓度(mg/L)
320
160
80
40
20
10
5

药剂A校正抑制率
81.08
75.94
67.16
56.03
48.11
14.91
5.71

药剂B校正抑制率
73.37
70.38
66.09
62.67
56.24
27.12
8.06

[ Last edited by propamidine on 2015-6-16 at 21:56 ]
回复此楼

» 收录本帖的淘帖专辑推荐

农药制剂加工 感悟人生 【经验-心得】 真菌

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

propamidine

铁杆木虫 (著名写手)


引用回帖:
2楼: Originally posted by zhenguoyan at 2015-06-17 08:19:15
既然齐军山等发表的文章,那恭喜你已经找到了方法。
首先,SPSS计算害虫的LC50需要的变量有:药剂浓度、总虫数、死虫数
那么对应到菌上,抑菌试验的话,可以总的菌当做100%,然后抑菌率当做死亡的部分。所以得到以 ...

你好!
按照文献方法,我下面如图的操作是否正确?
关于使用SPSS进行EC50的计算问题
0001.jpg


关于使用SPSS进行EC50的计算问题-1
0002.jpg


关于使用SPSS进行EC50的计算问题-2
0003.jpg


关于使用SPSS进行EC50的计算问题-3
0004.jpg


关于使用SPSS进行EC50的计算问题-4
0005.jpg

3楼2015-06-17 16:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

zhenguoyan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
propamidine: 金币+30, ★★★很有帮助, 感谢帮助,请麻烦再看一下我数据统计的结果是否正确。另外,相关系数在哪里找的到? 2015-06-17 16:17:30
既然齐军山等发表的文章,那恭喜你已经找到了方法。
首先,SPSS计算害虫的LC50需要的变量有:药剂浓度、总虫数、死虫数
那么对应到菌上,抑菌试验的话,可以总的菌当做100%,然后抑菌率当做死亡的部分。所以得到以下对应:
药剂浓度
总虫数-------100
死虫数-------抑菌率
比如药剂A,浓度输在第一列,第二列总虫数为100,第三列死虫数为相应的抑菌率。就可以计算出来了。
另外,毒力方程、EC50、95%置信限、相关系数和卡方值,是一次性出来的。
如果是棵小草,即使在最好的地方,你也长不成大树。果真如此,不如历经风雨,把自己培养成名贵花卉。
2楼2015-06-17 08:19:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

propamidine

铁杆木虫 (著名写手)


上面的数据处理以药剂A为例。

毒力方程:y=-2.072+1.283x
EC50=41.247 mg/L  95%置信限24.150~69.915
卡方值X2=20.501  当自由度为5时,查X2表,p=0.05时,X2=11.07        20.501>11.07,表示差异不显著,所得方程y=-2.072+1.283x是符合实际的。

以上对结果的基本表述没有问题吧?

另外,相关系数在probit结果里找不到啊?
4楼2015-06-17 16:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenguoyan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
propamidine: 金币+100, ★★★很有帮助, 非常感谢,顿感清晰了不少。冒出个新问题,在SPSS中双击图,添加总计拟合线,生成的“R2线性”就是相关系数吧?小数点后怎么保留4位? 2015-06-18 09:34:06
引用回帖:
4楼: Originally posted by propamidine at 2015-06-17 16:15:11
上面的数据处理以药剂A为例。

毒力方程:y=-2.072+1.283x
EC50=41.247 mg/L  95%置信限24.150~69.915
卡方值X2=20.501  当自由度为5时,查X2表,p=0.05时,X2=11.07        20.501>11.07,表示差异不显著,所 ...

方法对了,但是卡方检验是拟合度不好,sig.(p)要大于0.05,卡方值要小于临界值,拟合优度才好。
卡方检验表那,sig的值0.001后面不是标了b么,你看下面的注释,肯定是说方程存在异质性的。要是拟合度好,它会提示不存在异质性。
总之,就是卡方值越小,P越大,方程的拟合度越好,P越接近于1,方程越符合实际。
鉴于此,你看对数据进行一定的取舍。
从正态分布的角度,概率在16%-84%直接,误差最小。
从图中可知,可以把较小的数去掉后再拟合。
相关系数可以在SPSS中双击图,添加直线,可以得到。
如果是棵小草,即使在最好的地方,你也长不成大树。果真如此,不如历经风雨,把自己培养成名贵花卉。
5楼2015-06-18 08:27:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见