问题描述
我需要确定拟合数据的 GMM 的 99% 概率轮廓。按照这个示例，我希望能够指定要绘制哪些等高线，以及它们的 x、y。
mu1 = [1 2]; Sigma1 = [2 0; 0 0.5];
mu2 = [-3 -5]; Sigma2 = [1 0;0 1];
X = [mvnrnd(mu1,Sigma1,1000); mvnrnd(mu2,Sigma2,1000)];
GMModel = fitgmdist(X,2);
figure
y = [zeros(1000,1);ones(1000,1)];
h = gscatter(X(:,1),X(:,2),y);
hold on
gmPDF = @(x,y) arrayfun(@(x0,y0) pdf(GMModel,[x0 y0]),x,y);
g = gca;
fcontour(gmPDF,[g.XLim g.YLim])
title('{\bf Scatter Plot and Fitted Gaussian Mixture Contours}')
legend(h,'Model 0','Model1')
hold off
因此，在下图中，我希望能够用黑色虚线“k”绘制 99%。知道如何做到这一点吗？
解决方法
指定Enter your WRDS username [root]:.........
Enter your password:........
WRDS recommends setting up a .pgpass file.
You can find more info here:
https://www.postgresql.org/docs/9.5/static/libpq-pgpass.html.
There was an error with your password.
---------------------------------------------------------------------------
OperationalError Traceback (most recent call last)
/usr/local/lib/python3.7/dist-packages/sqlalchemy/engine/base.py in _wrap_pool_connect(self,fn,connection)
3210 try:
-> 3211 return fn()
3212 except dialect.dbapi.Error as e:
83 frames
OperationalError: fe_sendauth: no password supplied
The above exception was the direct cause of the following exception:
OperationalError Traceback (most recent call last)
OperationalError: (psycopg2.OperationalError) fe_sendauth: no password supplied
(Background on this error at: http://sqlalche.me/e/14/e3q8)
During handling of the above exception,another exception occurred:
OperationalError Traceback (most recent call last)
OperationalError: could not translate host name "12GO@wrds-pgdata.wharton.upenn.edu" to address: Name or service not known
The above exception was the direct cause of the following exception:
OperationalError Traceback (most recent call last)
/usr/local/lib/python3.7/dist-packages/psycopg2/__init__.py in connect(dsn,connection_factory,cursor_factory,**kwargs)
120
121 dsn = _ext.make_dsn(dsn,**kwargs)
--> 122 conn = _connect(dsn,connection_factory=connection_factory,**kwasync)
123 if cursor_factory is not None:
124 conn.cursor_factory = cursor_factory
OperationalError: (psycopg2.OperationalError) could not translate host name "12GO@wrds-pgdata.wharton.upenn.edu" to address: Name or service not known
(Background on this error at: http://sqlalche.me/e/14/e3q8)
的LevelList属性可以显示并获取给定轮廓线的坐标，然后读取轮廓句柄的ContourMatrix属性：
fcontour