estadisticas-stats + python

#Cree una función donde me retorne datos estadísticos,existe un modulo stats en python pero no tenia para la desviación estándar ponderada(pesada). Espero que les ayude y se diviertan


def stats(list_d,list_w):

    if type(list_d) and type(list_w) is list :

        lnd = len(list_d)
        lnw = len(list_w)

        if lnd == lnw and lnd > 0 :
          
            list_dxw = [list_d[i]*list_w[i] for i in xrange(len(list_d))]
            weimean = sum(list_dxw)/sum(list_w)
            list_rdw = [list_w[i]*((list_d[i] -weimean)**2) for i in xrange(lnd)]
            weistdev = ((sum(list_rdw))/(sum(list_w)))**0.5
            min_ = min(list_d)
            max_ = max(list_d)
            samnum = len(list_d)
            coefvar = weistdev/weimean
          
          
            return weimean, weistdev, min_, max_, samnum, coefvar