Para generar muestras con distribución t multivariada, uso esta función:

def multivariatet(mu,Sigma,N,M):
    '''
    Output:
    Produce M samples of d-dimensional multivariate t distribution
    Input:
    mu = mean (d dimensional numpy array or scalar)
    Sigma = scale matrix (dxd numpy array)
    N = degrees of freedom
    M = # of samples to produce
    '''
    d = len(Sigma)
    g = np.tile(np.random.gamma(N/2.,2./N,M),(d,1)).T
    Z = np.random.multivariate_normal(np.zeros(d),Sigma,M)
    return mu + Z/np.sqrt(g)

pero lo que estoy buscando ahora es eldistribución t multivariada de estudiantes por sí mismo para poder calcular la densidad de los elementos dondedimension > 1.

Eso será algo comostats.t.pdf(x, df, loc, scale) del paqueteescoria pero en espacio multidimensional.