如何计算标准差?
本人是一枚生物学的学生,由于对生物信息学特别感兴趣,于是想自学生物信息学(新手莫怪)。了解到生物信息学要有编程基础,尤其是要会一门编程语言,例如:R语言、Python、Perl等,还要熟悉Linux系统,作为生信小白,听说Python挺简单的,于是就自学了Python,花了两天时间了解了Python的基础语法后,今天想做个练习题试试手(实践是检验真理的唯一标准),下面是练习题:(试题来源:《PYTHON生物信息学数据管理》)
今日练习题:如何计算标准差?
计算标准差有些复杂,因为需要 for 循环来计算每个值的平方差,必须先有预先计算出 的平均值。
然后将每个值减去平均值((value-average) * *2)。 所有平方差要加在一起,除以结果的数量,最后计算结果的平方根。平方差求和可以设置一个变量为 0.0,每得到一个平方差就加上去。
公式如下:
书中给出的计算的脚本是:
import math
data = [3.53,3.47,3.51,3.72,3.43]
average = sum(data)/len(data)
total = 0.0
for value in data:
total += (value - average)**2
stddev = math.sqrt(total/len(data))
print(stddev)
参