python - PyMC2 and PyMC3 give different results...? -


i'm trying simple pymc2 model working in pymc3. i've gotten model run models give different map estimates variables. here pymc2 model:

import pymc theta = pymc.normal('theta', 0, .88)  x1 = pymc.bernoulli('x2', p=pymc.lambda('a', lambda theta=theta:1./(1+np.exp(-(theta-(-0.75))))), value=[1],observed=true) x2 = pymc.bernoulli('x3', p=pymc.lambda('b', lambda theta=theta:1./(1+np.exp(-(theta-0)))), value=[1],observed=true)  model = pymc.model([theta, x1, x2]) mcmc = pymc.mcmc(model) mcmc.sample(iter=25000, burn=5000) trace = (mcmc.trace('theta')[:]) print "\nthe map value theta is", trace.sum()/len(trace)

that seems work expected. had sorts of trouble figuring out how use equivalent of pymc.lambda object in pymc3. came across deterministic object. following code:

import pymc3  pymc3.model() model:      theta = pymc3.normal('theta', 0, 0.88)     x1 = pymc3.bernoulli('x1', p=pymc3.deterministic('b', 1./(1+np.exp(-(theta-(-0.75))))), observed=[1])     x2 = pymc3.bernoulli('x2', p=pymc3.deterministic('c', 1./(1+np.exp(-(theta-(0))))), observed=[1])      start=pymc3.find_map()     step=pymc3.nuts(state=start)     trace = pymc3.sample(20000, step, njobs=1, progressbar=true)  pymc3.traceplot(trace)

the problem i'm having map estimate theta using pymc2 ~0.68 (correct), while estimate pymc3 gives ~0.26 (incorrect). suspect has way i'm defining deterministic function. pymc3 won't let me use lambda function, have write expression in-line. when try use lambda theta=theta:... error:

astensorerror: ('cannot convert <function <lambda> @ 0x157323e60> tensortype', <type 'function'>)

something theano?? suggestions appreciated!

it works when use theano tensor instead of numpy function in deterministic.

import pymc3 import theano.tensor tt  pymc3.model() model:      theta = pymc3.normal('theta', 0, 0.88)     x1 = pymc3.bernoulli('x1', p=pymc3.deterministic('b', 1./(1+tt.exp(-(theta-(-0.75))))), observed=[1])     x2 = pymc3.bernoulli('x2', p=pymc3.deterministic('c', 1./(1+tt.exp(-(theta-(0))))), observed=[1])      start=pymc3.find_map()     step=pymc3.nuts(state=start)     trace = pymc3.sample(20000, step, njobs=1, progressbar=true)  print "\nthe map value theta is", np.median(trace['theta'])  pymc3.traceplot(trace);

here's output:

enter image description here


-

Comments

Post a Comment

Popular posts from this blog

c# - Binding a comma separated list to a List<int> in asp.net web api -

the default binder in web api expecst http://url.com/webapi/report/?pageids=3243&pageids=2365 to bind public ihttpactionresult report(list<int> pageids){ // exciting webapi code} i wish bind http://url.com/webapi/report/?pageids=3243,2365 as running out of space in url get . i have created class public class commaseparatedmodelbinder : system.web.http.modelbinding.imodelbinder { public bool bindmodel(httpactioncontext actioncontext, modelbindingcontext bindingcontext) { //binding in here } } and registered in webapiconfig.cs var provider = new simplemodelbinderprovider( typeof(list<int>), new commaseparatedmodelbinder()); config.services.insert(typeof(modelbinderprovider), 0, provider); i have altered method signature use model binder so public ihttpactionresult report( [modelbinder] list<int> pageids){ // exciting webapi code} however break point in binder not being ...
Read more

how to prompt save As Box in Excel Interlop c# MVC 4 -

i want export excel sheet using excel interloop... i have datatable convert excel sheet... now there way save file using save dialogue box in mvc 4 here code public string datatabletoexcel(datatable dt, string htmlheading, string key) { try { microsoft.office.interop.excel.application excel; microsoft.office.interop.excel.workbook excelworkbook; microsoft.office.interop.excel.worksheet excelsheet; microsoft.office.interop.excel.range excelcellrange; excel = new microsoft.office.interop.excel.application(); excel.visible = false; excel.displayalerts = false; excelworkbook = excel.workbooks.add(type.missing); excelsheet = (microsoft.office.interop.excel.worksheet)excelworkbook.activesheet; excelsheet.name = "sheet1"; dataset ...
Read more

xslt 1.0 - How to access or retrieve mets content of an item from another item? -

my use case have item has link in item. example, item 123456789/152 has metadata field dc.relation.hasversion=123456789/717 . within item view of 123456789/152 , how can retrieve values of metadata of 123456789/717 ? example, want retrieve dc.language.iso value of 123456789/717 123456789/152 . i looked @ related items feature in dspace don't know how metadata displayed in related items list pulled items in list. i using dspace version 5.3 mirage 2 theme . edit based on schweerelos answer, below actual code. using custom metadata field dc.relation.languageversion . basically, want link other version instead of displaying value, display dc.language.iso of other version in item-view.xsl . i've incorporated answer of schweerelos question in code just displaying value of dc.relation.languageversion . sample values of dc.relation.languageversion 10665.1/9843 ; 10665.1/9844 ie not complete uri handle. thanks in advance! actual code <xsl:template name=...
Read more