Fitting multivariate α-stable distributions to data is still not feasible in higher dimensions since the (non-parametric) spectral measure of the characteristic function is extremely difficult to estimate in dimensions higher than 2. This was shown by  and . α-stable sub-Gaussian distributions are a particular (parametric) subclass of the multivariate α-stable distributions. We present and extend a method based on  to estimate the dispersion matrix of an α-stable sub-Gaussian distribution and estimate the tail index α of the dis-tribution. In particular, we develop an estimator for the off-diagonal entries of the dispersion matrix that has statistical properties superior to the normal off-diagonal estimator based on the covariation. Furthermore, this approach allows estimation of the dispersion matrix of any normal variance mixture distribution up to a scale parameter. We demonstrate the behaviour of these estimators by fitting an α-stable sub-Gaussian distribution to the DAX30 components. Finally, we conduct a stable principal component analysis and calculate the coefficient of tail dependence of the prinipal components.