X-Received: by 10.107.135.21 with SMTP id j21mr32053534iod.2.1450147388630; Mon, 14 Dec 2015 18:43:08 -0800 (PST)
X-Received: by 10.50.85.111 with SMTP id g15mr32318igz.2.1450147388589; Mon, 14 Dec 2015 18:43:08 -0800 (PST)
Path: csiph.com!au2pb.net!feeder.erje.net!2.us.feeder.erje.net!enother.net!enother.net!news.glorb.com!mv3no16300748igc.0!news-out.google.com!l1ni3579igd.0!nntp.google.com!mv3no12382167igc.0!postnews.google.com!glegroupsg2000goo.googlegroups.com!not-for-mail
Newsgroups: comp.lang.python
Date: Mon, 14 Dec 2015 18:43:08 -0800 (PST)
Complaints-To: groups-abuse@google.com
Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=50.100.117.144; posting-account=SZ_svQkAAACWRFG2bDA-zgq8ILyl4-vo
NNTP-Posting-Host: 50.100.117.144
User-Agent: G2/1.0
MIME-Version: 1.0
Message-ID: <ba4b4615-2739-49b3-9bb3-eafc527ddebc@googlegroups.com>
Subject: Question about figure plot
From: Robert <rxjwg98@gmail.com>
Injection-Date: Tue, 15 Dec 2015 02:43:08 +0000
Content-Type: text/plain; charset=ISO-8859-1
Xref: csiph.com comp.lang.python:100445

Hi,

When I run the following code, there is no figure shown in the end.


//////////
import pymc
import numpy as np

n = 5*np.ones(4,dtype=int)
x = np.array([-.86,-.3,-.05,.73])

alpha = pymc.Normal('alpha',mu=0,tau=.01)
beta = pymc.Normal('beta',mu=0,tau=.01)

@pymc.deterministic
def theta(a=alpha, b=beta):
    """theta = logit^{-1}(a+b)"""
    return pymc.invlogit(a+b*x)
    
d = pymc.Binomial('d', n=n, p=theta, value=np.array([0.,1.,3.,5.]),\
                    observed=True)
....
import pymc
import mymodel

S = pymc.MCMC(mymodel, db='pickle')
S.sample(iter=10000, burn=5000, thin=2)
pymc.Matplot.plot(S)



I find that the figures are shown after these two lines by myself:
*************
import matplotlib.pyplot as plt
plt.show()

I have searched around and have not found some explanation about it.
The plot function here is different from Matlab's. Is there better ways than
my last two lines? (I am not confident whether my last two lines is the 
only choice.

Thanks,