Basic Example with Synthetic Data¶
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 |
|---|---|---|---|---|
|
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 |
|
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 |
|
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.