Another useful resource for learning the causal relationships between variables is the Tetrad program. Tetrad contains implementations of many major causal structure learning algorithms, including PC, PC-Stable, GES, FCI, ICA-LiNGAM, and others.
The Tetrad website has many resources explaining how to use Tetrad, although it is still not a particularly intuitive interface. Richard Scheines has recorded a tutorial on Tetrad (and causal structure learning in general), which might be more helpful.
The choice of search algorithm depends on the assumptions you are willing to make about your data. For example, if you think there are confounding variables that you have not observed, which might produce an association between the variables you have observed, then you should use the FCI algorithm (which allows for the presence of unobserved confounders). If you think your data is (at least approximately) multivariate Gaussian, then you can apply the default conditional independence test in FCI (the Fisher-Z test). However, if you think the data are not jointly Gaussian, you should probably use an independence test that allows for non-Gaussianity, such as the Conditional Correlation test. Can you tell us a little more about your data?