真値 0 を持ち，分散 1 で観測される x という変数を 5 回観測したとする．その5回の観測量から，元の分布関数の特長を調べたいので bootstrap する． import numpy import pylab n = 5 trial = 1000 a=numpy.random.normal(0,1,n) b=numpy.array([a[map(lam…

lambda で演算する方 import numpy import time f = lambda x: x**2 def test0(n): """空""" t1=time.time() t2=time.time() print t2-t1 def test1(n): """内包表現""" t1=time.time() a=[ f(x) for x in range(n) ] t2=time.time() print t2-t1 def test2(…

ちょっとしたテスト import numpy import time def test0(n): """空""" t1=time.time() t2=time.time() print t2-t1 def test1(n): """内包表現""" t1=time.time() a=[ x for x in range(n) ] t2=time.time() print t2-t1 def test2(n): """append""" t1=tim…

最近考えた TSV (tab separated values) 形式で書かれたデータファイルを読んで辞書のリストに登録する． import re tsvfile = "filename" f=open(tsv,"r") s=re.compile("[\t\s]") objects=map( lambda x: { "id" : int(x[0]), # ここら辺はそれぞれのデー…

たぶんどれも一緒だけど http://www.cv.nrao.edu/fits/documents/wcs/wcs.all.ps http://www.aoc.nrao.edu/~egreisen/AnAI.pdf 細かく描いてある．FITS の手引きもこれを参考にしたのではないだろうか．

python を使った天文用 plotter として大変便利な astLib ですが，scalebar の長さが違うという重要なバグがありましたのでメモしておきます． なにが間違っているかと言うと，scalebar を描く時の終点を決めるところで，ra 方向にも角度をそのまま足してし…

いつも忘れてしまうのでメモしておく matplotlib を使って図を作ってるときに，tex 形式のテキストをいれたいときに行う環境変数の設定． from matplotlib import rc rc('text', usetex=True) これで使いたいところで r"$\log\left( h^{-1} M/ M_{\odot}\rig…

これが分かりやすいかは別として，テキストファイルととかから読み込んだデータを配列に格納する． p=re.compile("[\t\n]") data = numpy.array( \ map( lambda x: map(lambda y: float(y), p.split(x)), \ [ \ filter( lambda x: x is not "", x) \ for x i…