Я делаю некоторый анализ образцов MCMC и использую пакет Python GetDist для создания своих контурных графиков.Однако контурные графики - это только часть всего анализа, и я хотел бы показать некоторые другие графики вместе с контурным графиком на том же рисунке.
Я использую matplotlib для генерации всех остальных графиков,так что мой вопрос: есть ли способ иметь график GetDist в подспоте matplotlib, чтобы у меня была фигура matplotlib с несколькими графиками и график GetDist в ней?
Я использую GridSpec для разделенияфигурировать в подсюжетах (а также разбивать подсюжеты в подсюжетах).Я попытался установить конкретный субплот в качестве текущей оси перед созданием графика треугольника, а также попытался взглянуть на исходный код GetDist, чтобы найти способ передать требуемый субплот в качестве аргумента в GetDist, но безуспешно.
Прямо сейчас мой код выглядит примерно так:
import matplotlib.pyplot as plt
import getdist.plots
from matplotlib import gridspec
import numpy as np
import numpy.random
N = 4
# Generate mock random data
chain = np.array([numpy.random.normal(loc=(i+0.5), scale=(i+0.5), size=100000) for i in xrange(N)])
Names = [str(unichr(i+97)) for i in xrange(N)]
Labels = [str(unichr(i+97)) for i in xrange(N)]
Sample = getdist.MCSamples(samples=chain.T, names=Names, labels=Labels)
#Set up plot layout
fig = plt.figure(figsize=(10.5,12.5))
gs = gridspec.GridSpec(3, 2, width_ratios=[2,3], height_ratios=[10,0.5,4], wspace=0.2, hspace=.05)
Lplot = gridspec.GridSpecFromSubplotSpec(N, 1, subplot_spec=gs[0,0], hspace=0.)
Rplot = gridspec.GridSpecFromSubplotSpec(2,1,subplot_spec=gs[0,1], height_ratios=[4,1.265])
Dplot = plt.subplot(gs[2,0:2])
axL = [plt.subplot(Lplot[i]) for i in xrange(N)]
axR = plt.subplot(Rplot[1])
GD = plt.subplot(Rplot[0])
for i in axL+[axR]+[GD]+[Dplot]: i.set_xticks([]); i.set_yticks([])
plt.sca(GD)
# Generate triangle plot
g = getdist.plots.getSubplotPlotter()
g.triangle_plot(Sample, filled=True)
plt.savefig("outfile.pdf", bbox_inches='tight')
Я бы хотел, чтобы мой контурный график находился в подпункте "GD".Любая помощь?