Linear non-Gaussian models with latent variables for casual discovery