Example: LISA data

from slipper.example_datasets.lisa_data import lisa_noise_periodogram, lisa_wd_strain
from slipper.fourier_methods import get_periodogram
from slipper.sample import LogPsplineSampler
import os
import numpy as np

np.random.seed(42)
OUTDIR = "out"
os.makedirs(OUTDIR, exist_ok=True)

---------------------------------------------------------------------------
ModuleNotFoundError                       Traceback (most recent call last)
Cell In[1], line 1
----> 1 from slipper.example_datasets.lisa_data import lisa_noise_periodogram, lisa_wd_strain
      2 from slipper.fourier_methods import get_periodogram
      3 from slipper.sample import LogPsplineSampler

ModuleNotFoundError: No module named 'slipper'

Noise

noise_mcmc = LogPsplineSampler.fit(
    data=lisa_noise_periodogram()[::5],
    sampler_kwargs=dict(Ntotal=200, burnin=100),
    spline_kwargs=dict(k=50, knot_locator_type='data_peak')
);
fig = noise_mcmc.plot_model_and_data(add_legend=True)
fig.savefig(f"{OUTDIR}/noise.png", dpi=100)

White dwarf background

import matplotlib.pyplot as plt
data = lisa_wd_strain()[0:10000]
pdgrm = get_periodogram(timeseries=data)[10:1000]
wdb_mcmc = LogPsplineSampler.fit(
    data=pdgrm,
    sampler_kwargs=dict(Ntotal=2000, burnin=1000),
    spline_kwargs=dict(k=30, knot_locator_type='data_peak')
);
fig = wdb_mcmc.plot_model_and_data(add_legend=True)
fig.savefig(f"{OUTDIR}/wdb.png", dpi=100)