Basic Example with Synthetic Data

Python Version License

This basic example contains synthetic data created for the purpose of testing the tool during its development


Files in the data Folder

  • synthetic_data_ETP.txt: Synthetic data.

  • age_model_ETP.txt: Synthetic age-depth model that is adapted to the data.

  • La2010d_ecc3L.txt: Astronomical solution data (eccentricity component).


Data Used

File

Description

Reference

Transformations / Pre-Processing

Version

synthetic_data_ETP.txt

Final synthetic dataset produced by combining astronomical signals with AR(1) correlated noise and white noise; represents ~1 Myr sequence with 1001 data points.

Marin & Robert (2014) for AR(1) reference

Added AR(1) noise with ρ = 0.5; noise variance adjusted so astronomical component is ~50% of total variance.

1.0

age_model_ETP.txt

Implicit age model computed from a predefined sedimentation rate function centered on 1 Myr/m, used to transfer the synthetic ETP signal from time to depth domain.

Synthetic model described in the paper

Computed using time–depth transfer function ( t(D) ); sampling every 0.1 cm for 1 m.

1.0

La2010d_ecc3L.txt

Synthetic astronomical signal created by summing eccentricity, obliquity, and climatic precession components from the Laskar et al. (2004) solution, each normalized to zero mean and unit variance.

Laskar et al. (2004)

Normalization of each orbital component to zero mean and unit variance before summation.

1.0


Execute the Example

Note

  • The folder examples has to be in the same directory as the AstroGeoFit_tool.py file. If not, this will not work.

To run the example, open a terminal in the AstroGeoFit folder. Once this is open execute the following

python AstroGeoFit_tool.py --basic_example

This command will execute the significance test, the genetic algorithm fitting, the MCMC, and the MCMC weights calculation.

If you only want to run specific components (e.g., the genetic algorithm fitting and the MCMC), you can use:

python AstroGeoFit_tool.py --basic_example -fit -mcmc

The full list of execution options can be found in AstroGeoFit Script Execution.


Visualization of the Results

To visualize the results, you can use the two Jupyter notebooks provided. To obtain the results of the basic_example data, just modify the variable configuration_file_path in the following way:

configuration_file_path = "basic_example"

And then the visualization notebook will be ready to be executed.

For a detailed guide on using the notebooks and interpreting the figures, please refer to Visualization of the Results Guide.