Suppose we have a joint distribution p(x_1,x_2), and we know x_1,x_2,p. Both are discrete, (x_1,x_2) is scatter, its contour could be drawn, marginal as well. I would like to show the area of 95% quantile (a scale of 95% data will be contained) of the joint distribution, how can I do that?
How to calculate a 95 credible region for a 2D joint distribution?
507 Views Asked by TxWang At
2
There are 2 best solutions below
0
fakufaku
On
As the other points out, there are infinitely many solutions to this problem. A practical one is to find the approximate center of the point cloud and extend a circle from there until it contains approximately 95% of the data. Then, find the convex hull of the selected points and compute its area.
Of course, this will only work if the data is sort of concentrated in a single area. This won't work if there are several clusters.
Related Questions in PYTHON
- How to store a date/time in sqlite (or something similar to a date)
- Instagrapi recently showing HTTPError and UnknownError
- How to Retrieve Data from an MySQL Database and Display it in a GUI?
- How to create a regular expression to partition a string that terminates in either ": 45" or ",", without the ": "
- Python Geopandas unable to convert latitude longitude to points
- Influence of Unused FFN on Model Accuracy in PyTorch
- Seeking Python Libraries for Removing Extraneous Characters and Spaces in Text
- Writes to child subprocess.Popen.stdin don't work from within process group?
- Conda has two different python binarys (python and python3) with the same version for a single environment. Why?
- Problem with add new attribute in table with BOTO3 on python
- Can't install packages in python conda environment
- Setting diagonal of a matrix to zero
- List of numbers converted to list of strings to iterate over it. But receiving TypeError messages
- Basic Python Question: Shortening If Statements
- Python and regex, can't understand why some words are left out of the match
Related Questions in NUMPY
- Why numpy.vectorize calls vectorized function more times than elements in the vector?
- Producing filtered random samples which can be replicated using the same seed
- Numpy array methods are faster than numpy functions?
- When I create a series of spectrograms from a long audio file, the colour intesities vary noticably
- How do I fix a NumPy ValueError for an inhomogeneous array shape?
- How should I troubleshoot "RuntimeWarning: invalid value encountered in arccos" in NumPy?
- Unravel by multi-index/group
- Calculating IRR Using Numpy
- Integrating with an array of upper limits without sacrificing time efficiency
- Why doesn't this code work? - Backpropagation algorithm
- How to remove integers from a mixed numpy array containing sub-arrays and integers?
- How to transfer object dataframe in sklearn.ensemble methods
- Rust cannot borrow as mutable
- Why does the following code detect this matrix as a non-singular matrix?
- How to detect the exact boundary of a Sudoku using OpenCV when there are multiple external boundaries?
Related Questions in SCIPY
- Integrating with an array of upper limits without sacrificing time efficiency
- Hourly and annual optimization problem over matrix
- Checking Event in solve_ivp
- Is polar decomposition commutative for diagonal matrices?
- loading variables from a python script in matlab
- cov2corr() for scipy sparse matrices
- Solving double integration in python using scipy.integrate
- Saving a scipy.sparse matrix directly as a regular txt file
- How can I apply scipy.interpolate.RBFInterpolator on an image / ndarray?
- Truncated normal distribution doesn't agree with untruncated normal distribution?
- Incorrect world-to-camera matrix with euler angles (yaw, pitch, roll)
- Scipy solve_ivp extremely slow/freezing
- Interpolating 3D volumetric data with scipy's RegularGridInterpolator
- Hilbert envelope in scipy matches original signal
- How to fill the area for segmentation using fill_holes function?
Related Questions in CONFIDENCE-INTERVAL
- How to bootstrap indirect relationships within a multilevel sem model?
- Bootstrapped confidence intervals for non-linear mixed effects model
- ECharts confidence bands in legend as one item
- Confidence intervals on predictions using predictor variables for a non-linear mixed effects model (nlme)
- How Can I Get Bootstrapped BCa Confidence Intervals for Spearman's Correlation Matrix in R?
- Confidence Interval for Some (But Not All) Interaction Terms in a Polynomial Interacted with a Binary Treatment
- Bootstrap for a variable after group by (getting confidence intervals for each group)
- How to change the confidence interval for marginal effects at levels of a moderator variable in R
- Add confidence intervals to a plot
- Error band using plt.fill_between is giving two areas instead of one
- Calculating confidence intervals using the expss package
- Confidence Interval for a Chi-Square & T-Student in R
- Plotting confidence or prediction intervals on point plots ggplot2
- Negative confidence interval for geom_smooth() in ggplot
- Problem: pooling confidence intervals when calculating SIRs?
Related Questions in CREDIBLE-INTERVAL
- reorder community order using multiSpeciesTP, siberdensityplot
- Analytic Highest Density Interval in Python (preferably for Beta distributions)
- How to find and plot different credible intervals?
- Calculate 95 % credible interval for multiple columns in data frame R
- Plotting both 50% quantile and credible interval from package 'mcp' output in R
- Extract credible intervals for robust correlations in R
- determine equal-tail credible interval
- Extracting posterior modes and credible intervals from glmmTMB output
- Gaussian Process Confidence vs Credible Intervals
- Extract posterior estimate and credible intervals for random effect for lme4 model in R
- Highest Density Interval (HDI) for Posterior Distribution Pystan
- Coverage probability estimation for Bayesian Credible interval (Normal Distribution)
- How to calculate a 95 credible region for a 2D joint distribution?
- Error propagation in a Bayesian analysis of a Markov chain
- Reduce credible intervals in Causal Impact model
Trending Questions
- UIImageView Frame Doesn't Reflect Constraints
- Is it possible to use adb commands to click on a view by finding its ID?
- How to create a new web character symbol recognizable by html/javascript?
- Why isn't my CSS3 animation smooth in Google Chrome (but very smooth on other browsers)?
- Heap Gives Page Fault
- Connect ffmpeg to Visual Studio 2008
- Both Object- and ValueAnimator jumps when Duration is set above API LvL 24
- How to avoid default initialization of objects in std::vector?
- second argument of the command line arguments in a format other than char** argv or char* argv[]
- How to improve efficiency of algorithm which generates next lexicographic permutation?
- Navigating to the another actvity app getting crash in android
- How to read the particular message format in android and store in sqlite database?
- Resetting inventory status after order is cancelled
- Efficiently compute powers of X in SSE/AVX
- Insert into an external database using ajax and php : POST 500 (Internal Server Error)
Popular # Hahtags
Popular Questions
- How do I undo the most recent local commits in Git?
- How can I remove a specific item from an array in JavaScript?
- How do I delete a Git branch locally and remotely?
- Find all files containing a specific text (string) on Linux?
- How do I revert a Git repository to a previous commit?
- How do I create an HTML button that acts like a link?
- How do I check out a remote Git branch?
- How do I force "git pull" to overwrite local files?
- How do I list all files of a directory?
- How to check whether a string contains a substring in JavaScript?
- How do I redirect to another webpage?
- How can I iterate over rows in a Pandas DataFrame?
- How do I convert a String to an int in Java?
- Does Python have a string 'contains' substring method?
- How do I check if a string contains a specific word?
If you are interested in finding a pair x_1, x_2 of real numbers such that P(X_1<=x_1, X_2<=x_2) = 0.95 and your distribution is continuous then there will be infinitely many of these pairs. You might be better of just fixing one of them and then finding the other