No project description provided
Project description
VineCopulas
VineCopulas
is a Python package that is able to:
- Fit both bivariate and vine copulas
- Simulate from both bivariate and vine copulas
- Allow for both discrete as well as continuous input data
- Draw conditional samples for any variables of interest with the use of bivariate copulas and different vine structures
Installation
pip install vinecopulas
Getting Started
Get started by testing checking out the package functionality using the Abalone example data.
import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
from vinecopulas.marginals import *
from vinecopulas.bivariate import *
datapath = 'https://raw.githubusercontent.com/VU-IVM/vinecopula/develop/doc/sample_data.csv'
df = pd.read_csv(datapath)
df.head()
Transform the data to pseudo data and fit a survival gumbel copula between two variables. Use the fitted copula to generate random samples.
x = np.array(df)[:,:-1]
u = pseudodata(x) # computes the pseudodata
cop = 4 # copula 4 is the gumbel copula with 180 degree rotation
par = fit(cop, u[:,:2]) # fit the variables in the first 2 columns
n = len(u) # number of samples to generate
ur = random(cop, par, n) # generate random samples
# plot
plt.scatter(u[:,0],u[:,1], label = 'Data')
plt.scatter(ur[:,0], ur[:,1], alpha = 0.5, label = 'Random')
plt.xlabel('$u_1$')
plt.ylabel('$u_2$')
plt.legend()
Fit a vine copula between multiple variables in the data, considering all possible copulas available in the package.
cops = list(range(1,16)) # fit vine copula according to these copulas
M, P, C = fit_vinecop(u, cops, vine = 'R') # fit R-vine
plotvine(M,variables = list(df.columns[:-1]), plottitle = 'R-Vine') # plot structure
For more examples, please have a look at the documentation.
Contribution Guidelines
Please look at our Contributing Guidelines if you are interested in contributing to this package.
Asking Questions and Reporting Issues
If you encounter any bugs or issues while using VineCopulas
, please report them by opening an issue in the GitHub repository. Be sure to provide detailed information about the problem, such as steps to reproduce it, including operating system and Python version.
If you have any suggestions for improvements, or questions, please also raise an issue.
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for VineCopulas-0.2.5-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 1b4e63fe70b6e89a163da4f279f21f401832b031aee0e00f2998d6e87216bf1f |
|
MD5 | 96d520c66c2cb4075685e43d0f059658 |
|
BLAKE2b-256 | 7816d86ab694f1241c028dd68388ce7d852cd26e07ea77451e626e73e28c9a46 |