Path: csiph.com!usenet.pasdenom.info!news.redatomik.org!newsfeed.xs4all.nl!newsfeed1a.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail
To: python-list@python.org
From: Terry Reedy <tjreedy@udel.edu>
Subject: Re: Question about 'x' in pymc.invlogit(a+b*x)
Date: Sat, 07 Mar 2015 13:13:59 -0500
References: <c17c3949-32dc-4f02-850d-7534abbf9897@googlegroups.com>
Mime-Version: 1.0
Content-Type: text/plain; charset=utf-8; format=flowed
Content-Transfer-Encoding: 7bit
User-Agent: Mozilla/5.0 (Windows NT 6.1; WOW64; rv:31.0) Gecko/20100101 Thunderbird/31.5.0
In-Reply-To: <c17c3949-32dc-4f02-850d-7534abbf9897@googlegroups.com>
Precedence: list
Newsgroups: comp.lang.python
Message-ID: <mailman.157.1425752060.21433.python-list@python.org>
Lines: 43
NNTP-Posting-Host: 2001:888:2000:d::a6
Xref: csiph.com comp.lang.python:87119

On 3/7/2015 11:44 AM, fl wrote:
> Hi,
>
> I once learnt Python for a few weeks. Now, I try to using a Python package
> pymc. It has the following example code:
>
>
>
>
> import pymc
> import numpy as np
> n = 5*np.ones(4,dtype=int)
> x = np.array([-.86,-.3,-.05,.73])

x is defined here as a module ('global') name

> 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)

x is used here.  Names in functions can be resolved in the local, 
(nonlocal, when present, but not here), global, and builtin namespaces, 
in that order.  They are resolved when the function is called.

> d = pymc.Binomial('d', n=n, p=theta, value=np.array([0.,1.,3.,5.]),\
> observed=True)

> I don't understand the 'x' in pymc.invlogit(a+b*x). Could you help me on it?

See comments, and contemplate this simplified example

 >>> def f(a): return a+b

 >>> b = 3
 >>> f(4)
7

-- 
Terry Jan Reedy