2010-01-01から1年間の記事一覧

回廊亭

人に勧められて回廊亭殺人事件を読んだ.推理小説とかほとんど読んだことがなかったから新鮮だった.っていうかこんな展開かよ!みたいなね.まぁ,登場人物全部が怪しいからあり得ると言えばあり得るけど.まぁ面白かったさ.回廊亭殺人事件 (光文社文庫)作…

温泉に行った

ずっと働きっぱなしだったので休みを取って温泉に行って来た.(といっても昼間は働いていましたけど) で,行ったところは箱根湯本の弥次喜多の湯というところ. http://www.kokoga-e.com/fureai/yajikita/ 評価は今まで行った温泉の中で最低ランクかな. …

うー

忙殺されていて猫の手も借りたい勢い.なんつうか,図を描いたりする時間さえない.やばいなぁこれは.

matplotlib.pylab で時系列データを扱う

pylab で時系列データを plot するにはどうしたらいいか小一時間悩んだ.結論は単に datetime オブジェクトを放り投げればいいだけだった. datetime.datetime(*strptime(row[0], "%Y-%m-%d %H:%M:%S")[0:6])) row[0] には datetime.datetime.now() で得た時…

シリアルポート増設カード

シリアルポート増設カードが欲しい.しかし,貸し出しをしてくれた Interface 社のカードは 64bit 非対応であった. そこで他社を調査することにした. 品名 Linux64bit 専用ケーブル 参考 AREA SD-PCI9845-4S ? 必要なし ラトックREX-PCI64 ○ Dsub44Dsub9 x…

Poisson 分布のエラーと Gaussian 分布のエラーの違い

http://d.hatena.ne.jp/youtsumi/20100109/1263050919

Poisson分布に従う少ない観測量に対する誤差統計量の近似式

Poisson 分布は観測値が大きければ中心極限定理により Guassian 分布に近づくわけですが,得てして天文業界では観測値が大きくない場合が多いです.そんなわけで Poisson 分布関数をそれなりにまじめに計算しないと本当はいけない場合が多いです.ざっと の…

忙殺されて殺気立っていた。

今日は本当に忙殺されていて殺気立っていた。 昨日、おじいさんの看病に広島に帰った母親が足を折って帰京。病院に行くために午前休。 飯を食いながら14時まで集中して昨日残った宿題をこなし、二、三件の議論をこなし、そして別件の仕事をあげて17時からア…

1/7

今日は会議日。まぁ、正月がんばったおかげで基本的に全部の会議は楽々こなせた。 帰ったら親が骨折してた。もう若くないのだから気をつけてよね。

積分

当たり前だけど import scipy.integrate print scipy.integrate.quad(lambda x: x**2, 0,1) で積分できちゃうのすごいな. print scipy.integrate.quad(lambda x: x**-0.5, 0,1) とかも

Rを使う

R

ベクトルの書き方 > x <- c(1,2,3) プリントの仕方 > x

宇宙年齢を求める

(知らない人には不親切な)理屈 宇宙論の教科書(Cosmological Physics, J. A. Peacock とか)に って書いてあるのでであることを使えば, とかって書き換えられますから.これをつまり,まで積分するとでる. hubble parameter を使うと はやりの宇宙論パラ…

家に帰ったらやること

/work03/CFHQS/image/noiseを見て昔書いたプログラムを解析する.そして,動径分布を調べる.

ARTDATA

EXAMPLES 1. Create a galaxy cluster with a power law distribution of field galaxies and stars as background/foreground. ar> gallist galaxies.dat 100 spatial=hubble lum=schecter egal=.8 ar> gallist galaxies.dat 500 ar> starlist galaxies.dat…

Kolmogorov-Smirnov test with python

KS 検定を勘違いしていなければ、これで一応 XX % の確率で分布が違うとか、一致する。と言える。KS 検定についてはもう少し勉強する。しかし、数が少なくてもいいとか便利すぎ。 #!/usr/bin/env python from scipy import stats dmax=0.37 n=7 f = lambda x…

Two dimensional Kolmogorov-Smirnov test

J. A. Peacock (1983) http://adsabs.harvard.edu/abs/1983MNRAS.202..615P

買い物

ちょっと用があってアメリカでそこそこ値段のするものを買ったんです。で、円高だったのもあったし、日本で買うのに比べどれくらい安かったのかぁ?と思って日本で調べてみたら定価にして1万円近く開きがあることがわかって、すごい得した気分になりました。…

やがて哀しき外国語

やがて哀しき外国語 (講談社文庫)作者: 村上春樹出版社/メーカー: 講談社発売日: 1997/02/14メディア: 文庫購入: 5人 クリック: 29回この商品を含むブログ (115件) を見る最初ははなしが散漫だなぁと思ってたいくつだったけど、だんだん自分が感じていること…

来た!

いっぱいあけおめ。新しい時代来たよ。電話も来てうれしかった。