Comments (2)
It's cause the argument is to_pobs
import numpy as np
import pandas as pd
import scipy.stats as st
from copulae import ClaytonCopula
np.random.seed(100)
v1 = np.sort(st.norm.rvs(loc=0, scale=1.5, size=10000))
v2 = st.norm.rvs(loc=0, scale=1, size=10000)
v3 = np.sort(st.norm.rvs(loc=0, scale=1.5, size=10000))
v4 = st.norm.rvs(loc=0, scale=1, size=10000)
total_sims = pd.DataFrame([0.1 * (v1 + v2), 0.1 * (v3 + v4)], index=['V1', 'V2']).T
gumbel_ca_unif = pd.DataFrame((st.gumbel_l.cdf(total_sims['V1'], *st.gumbel_l.fit(total_sims['V1']))),
columns=['gumbel_V1'])
gumbel_es_unif = pd.DataFrame((st.gumbel_l.cdf(total_sims['V2'], *st.gumbel_l.fit(total_sims['V2']))),
columns=['gumbel_V2'])
gumbel_unif = np.asarray(pd.concat([gumbel_ca_unif, gumbel_es_unif], axis=1))
for to_pobs in [True, False]:
cop = ClaytonCopula(2)
cop.fit(data=gumbel_unif, to_pobs=to_pobs)
print(f"""
to_pobs = {to_pobs}
{cop.__class__.__name__} theta: {cop.params}
{cop.__class__.__name__} Kendall : {cop.tau}")
{cop.__class__.__name__} LogLikihood : {cop.log_lik(gumbel_unif, to_pobs=to_pobs)}")
""".strip())
from copulae.
Ahh, dumb mistake. The results now are consistent with SAS. Thank you!
from copulae.
Related Issues (20)
- Copula Fit Replication with SAS using Uniform Marginals HOT 2
- Error after update HOT 2
- the summary( ) got converging error HOT 2
- ModuleNotFoundError: No module named 'copulae.special._specfunc' HOT 2
- Python 3.6 not supported? HOT 4
- How to Inverse Normalisation?[Quantile Function] HOT 2
- Fixing Correlation Matrix if copulae.fit() HOT 2
- test_gmc failures HOT 1
- Feature Request - Plot a copula
- Gaussian Mixture Copula Model gives wrong sampling
- Frank and Gaussian Copula give Error
- PEP 517 builds HOT 1
- Ver 0.7.7 - Kendall Tau not supported for Student Copula HOT 1
- Fit Summary HOT 1
- Support for discrete marginals
- Thank you! HOT 2
- An error occurred while importing copulae HOT 2
- GitHub project documentation link is broken HOT 1
- Reference for gaussian mixture copula
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from copulae.