Fitting Gaussian DAG Models with Hidden Variables in R