my Poisson simulation

2010-03-19 12:21:02

# -*- coding:utf-8 -*-
"""
Created on 2010-3-18

@author: Cross
"""
from scipy.stats import poisson

p1 = poisson(1.725)
p2 = poisson(0.935)
margin_HAD = 1.12
margin_HILO = 1.08

psize = 100000


a1 = p1.rvs(psize)
a2 = p2.rvs(psize)

h=0
d=0
a=0
HI = 0
LO = 0

for i in range(psize):
    if a1[i]>a2[i]:
        h = h + 1
    elif a1[i]        a = a + 1
    else:
        d = d + 1

    if a1[i] + a2[i]<2 .5:=".5:" br="br">        LO = LO + 1
    else:
        HI = HI + 1


print "======================="
print "HDA"
print 1/margin_HAD/(h/float(psize))
print 1/margin_HAD/(d/float(psize))
print 1/margin_HAD/(a/float(psize))
print "======================="
print "HILO 2.5"
print 1/margin_HILO/(LO/float(psize))
print 1/margin_HILO/(HI/float(psize))

possion function example (from scipy)

>>> from scipy.stats import poisson, lognorm >>> myShape = 5;myMu=10 >>> ln = lognorm(myShape) >>> p = poisson(myMu) >>> ln.rvs((10,)) generate 10 RVs from ln array([ 2.09164812e+00, 3.29062874e-01, 1.22453941e-03, 3.80101527e+02, 7.67464...... http://bit.ly/9MkbZC

导出 豆瓣 日记 到其他Blog和twitter

导出 日记 到其他blog： 豆瓣日记 -> twitterfeed -> ping.fm -> blogger/wordpress 导出 广播 到twitter： 豆瓣广播 -> feedburner -> twitter (->facebook) 豆瓣广播 到 新浪微博： 豆瓣广播 -> google reader (add star) stared -> 新浪微博（关联博客） http://bit.ly/dpOErL

[Blog post] python reads excel

python reads excel crossz | January 11, 2010 at 1:12 am | Tags: python tutorial howto | Categories: life | URL: http://wp.me/pthU9-hH

Use the excellent xlrd package, which works on any platform. That means you can read Excel files from Python in Linux! Example usage:

Open the workbook

import xlrd wb = xlrd.open_workbook('myworkbook.xls')

Check the sheet names

wb.sheet_names()

Get the first sheet either by index or by name

sh = wb.sheet_by_index(0) sh = wb.sheet_by_name(u'Sheet1')

Iterate through rows, returning each as a list that you can index:

for rownum in range(sh.nrows):     print sh.row_values(rownum)

If you just want the first column:

first_column = sh.col_values(0)

Index individual cells:

cell_A1 = sh.cell(0,0).value cell_C4 = sh.cell(rowx=3,colx=2).value

(Note Python indices start at zero but Excel starts at one)

Add a comment to this post

--

WordPress.com | Thanks for flying with WordPress! Manage your subscriptions | Publish text, photos, music, and videos by email using our Post by Email feature. | Signup with WordPress.com