statistics - python matplotlib - contour plot - confidence intervals -

-

i'm trying plot contours (doable) on grid of data using matplotlib.pyplot.contour, contours placed @ 1, 2 , 3 sigma away peak value. there neat way apart brute force? thanks!

python version is

python 2.7.2 |epd 7.2-2 (64-bit)| (default, sep 7 2011, 16:31:15) [gcc 4.0.1 (apple inc. build 5493)] on darwin

you can specify list of z-values contours drawn. have collect correct z-values distribution. here example '1, 2, , 3 sigma away peak value':

enter image description here

code:

import numpy np import matplotlib.cm cm import matplotlib.mlab mlab import matplotlib.pyplot plt  #set 2d gaussian: delta = 0.025 x = np.arange(-3.0, 3.0, delta) y = np.arange(-3.0, 3.0, delta) x, y = np.meshgrid(x, y) sigma = 1.0 z = mlab.bivariate_normal(x, y, sigma, sigma, 0.0, 0.0) #get z values contours 1, 2, , 3 sigma away peak: z1 = mlab.bivariate_normal(0, 1 * sigma, sigma, sigma, 0.0, 0.0) z2 = mlab.bivariate_normal(0, 2 * sigma, sigma, sigma, 0.0, 0.0) z3 = mlab.bivariate_normal(0, 3 * sigma, sigma, sigma, 0.0, 0.0)  plt.figure() #plot gaussian: im = plt.imshow(z, interpolation='bilinear', origin='lower',                  extent=(-50,50,-50,50),cmap=cm.gray) #plot contours @ whatever z values want: cs = plt.contour(z, [z1, z2, z3], origin='lower', extent=(-50,50,-50,50),colors='red') plt.savefig('fig.png') plt.show()

Comments

Post a Comment

Popular posts from this blog

django - How can I change user group without delete record -

-
on website user have 1 group. , user can change group. it's made user.groups.clear() and user.groups.add(new_group) but it's not efficient, because there 2 sql query: delete, insert. how can change group update query? user , group related each other using manytomanyfield . means intersection table exists relating both entities, , if don't specify model map (using through attribute) django creates 1 you. looking @ sources django.contrib.auth.models see that's case. fortunatly, can access intermediary model using through attribute of manager (in case, user.groups.through ). can use regular model. example: >>> alice = user.objects.create_user('alice', 'alice@example.com', 'alicepw') >>> employee = group.objects.create(name='employee') >>> manager = group.objects.create(name='manager') >>> alice.groups.add(employee) >>> alice.groups.all() [<group: employee...
Read more

java - Need to add SOAP security token -

-
i need set custom soap header attribute on jax-ws generated webservice client. case webservice calls must go through proxy server requiring specific token (recieved web request header) present in soap request header. e.g.: 1 carserviceservice service = null; 2 service = new carserviceservice(new url(url), new qname(qname); 3 carserviceendpoint port = service.getcarserviceport(); it seems in line 3 wsdl retrieved , call fails due missing security token. 1 point direction on how done? a detailed example has been mentioned here: creating , deploying jax-ws web service on tomcat 6 this article shows how create , use security token .
Read more

java - EclipseLink JPA Object is not a known entity type -

-
i have strange problem eclipselink , object want persist. have 1 object (keypointlistimpl) stores object keypointimpl in list. persisting keypointimpl objects works great if try persist keypointlistimpl object java.lang.illegalargumentexception says object keypointimpl isn't known entity type. here keypointimpl code: @entity @table(name="keypoints") public class keypointimpl implements keypoint { @id @generatedvalue(strategy = generationtype.identity) private int id; @enumerated(enumtype.string) private detectortype keypointtype; private float x; private float y; private float size; private float angle; private float response; private int octave; private int classid; ... } here keypointlistimpl code: @entity @table(name="keypointlists") public class keypointlistimpl implements keypointlist { @id @generatedvalue(strategy = generationtype.identity) private int id; @onetoone(cascade={ca...
Read more