Gibbs point estimation

May 23, 2010

From Lehmann and Casella (1998):
\[ \frac{1}{M}\sum_{i=1}^ME(\Lambda|x, b_i) = \frac{1}{M}\sum_{i=1}^M\frac{b_i}{1+b_i}(a+x)\],
the posterior estimate:
\[\hat\pi(\lambda|x) = \frac{1}{M}\sum_{i=1}^M\pi(\lambda|x, b_i) = \frac{\lambda^{a+x-1}}{M\Gamma(a+x)}\sum_{i=1}^M\left(\frac{1+b_i}{b_i}\right)\exp\left(-\lambda\frac{1+b_i}{b_i}\right)\]

Advertisements

下面的Python碼:

import timeit
from scipy import *
from scipy.linalg import *
from scipy.stats.distributions import *

def svdTime():
    return svd(norm.rvs(size=[1e3, 1e3]))

t = timeit.Timer('svdTime()', 'from __main__ import svdTime')
print(t.timeit(number=1))

會得到:

gong-yi@21:17 ~ $ python test.py 
8.96158790588

這還是在Intel core i5上得到的結果咧!
R可以在2秒左右, 而Octave更是在1秒左右,SciPy的開發團隊要注意一下了

這台T410是最讓我心酸的一台小黑

那種大國氣勢已經不在了,全身上下都是cost down的傷㾗

過氣了嗎?

(原文同時po於mobile01)

個人感想: 我上週才拿到我的 T410:

規格: intel core 5i
memory: ddr3 2g
HD: 250gb 5400rpm (WD, 非常安靜)
LCD: LED, 1440×900, roll-cage 已取消
顯卡: intel HD (在Linux 看 1080p HD flash, 例如 big bunny buck 不會lag)
centrino 1000 N wireless, bluetooth, finerprint reader, smartcard reader, SD reader,
3xusb (1 powered usb) 1xdisplayport, 1xVGA port, RJ-45 RJ-11. 9-cell 電池
內建陣列mic, webcam. 七列標準鍵盤, DVD rw, combo headphone/mic 插孔

跟 T61
(T61 規格:
intel T7500 (core duo)
memory: ddr2 3g
HD: 100gb 7200rpm (segate)
LCD: 傳統燈管背光式, 1440×900
顯卡:Nvidia Quadro MV 140, 在Linux 看 1080p HD flash, 例如 big bunny buck 會lag
Intel wireless, bluetooth, finerprint reader, pccard slot (寛窄皆有),
3xusb 1xVGA port, RJ-45 RJ-11. 9-cell 電池
 內建mic, webcam. 七列標準鍵盤,DVD rw, headphone 孔,mic 孔)

比較:

優點:

  1. LCD 較亮, 變薄,螢幕對稱了 
  2. 續航力大輻增強 (2.7 hour -> 5.5 hour, 在Linux下, 以編寫程式為主)
  3. 比較輕,附的AC adaptor 也變輕變小,全機整體略薄於T61
  4. 鍵盤手感較佳,造型亦較有質感, ESC 和 DEL 的設計對我而言真的是非常有用
    ( 常用 vi 的人就知道我在說什麼, 常要處理spreadsheet的人也會常用到DEL )
    P.S.: T61的三個銀色按鈕不知道Yamato Lab是不是看了老夫子才放那三個醜的要命的圈圈上去
  5. WD 硬碟很安靜
  6. 指紋開機
  7. 裝第二條memory 變方便了
  8. 喇叭音量變大(但音質沒有說好很多)
  9. 燈號變少更加簡潔

缺點 :

  1. LCD A蓋和B蓋很軟,用手加壓LCD會漾出水波紋,這在有roll-cage的 T61上是不會出現的
  2. 不支援 bay battery
  3. 主機體變厚(不含LCD的部份)
  4. 鍵盤打鍵音變大,可能會吵到已經在睡的室友或家人
  5. CPU 風扇的高頻音非常明顯,在夜中工作時有點惱人(在office聽不到) 
  6. Linux不支援這個 fingerprint reader的晶片, 寫信給Lenovo沒人回
  7. 機底的 memory slot 塑膠蓋子是全機最脆弱的部份
  8. combo 耳機要另外買新的 (或至少得買個轉接頭), 價格比較高
  9. 整體的”塑膠感變重”,原來的”真小黑”的感覺走味了,有降格成R系列的感覺
  10. Lenovo 字樣變大了

    

感想:

  1. Lenovo再這樣走cost down的路下去,我的下一台小黑(如果我還是買Thinkpad) 絕不會是T系列
    我相信原T系列的用戶會更加重視耐撞耐操的持性(較之於輕薄的持性)。 
    取消T410的roll cage是個讓我非常不滿的改變,這不該是T系列該作的讓步。
    我寧可多付一些錢,但我要耐操耐撞的T系列
  2. 設計得再好,用料不佳,品管不嚴也是惘然
  3. Lenovo的客服人員有很大的比例是印度人,有的腔調重到我聽不懂他在講殺小。

P.S. T61在台灣買的, T410在米國買的,兩台價差約 NTD$11000

以上。

今天把 GotoBLAS2 (1.13) 裝在我的新筆電上,遇上了點麻煩。

在編譯 kernel/x86_64/gemm_ncopy4.S 這個檔案的時候, 發現一直編不過, 一直
遇到這個錯誤: RPREFETCHSIZE: undefined

實在令人很火大, 之前在Core-Duo平台上編都沒問題, 為什麼今天就踩到雷?

後來發現是 Makefile bug, 如果沒有在 Makefile.rule
中指定 TARGET = NEHALEM, make 會真的不知道 RPREFETCHSIZE 要設成什麼。

原因很簡單,棍,他們懂個屁政治。

要制定政策不是看看政論節目放放嘴砲就作得出來的,如果你連公開的法條都懶得看,然後就自詡高級選民(這點台北人很擅長)
別人都是理盲濫情的人渣,那麼, garbage in, garbage out.

穿衣品味這回事,有眼睛就可以看, 不見得需要腦思考,同時還可以練練刻薄的毒舌, 這應該很合台北人的味口。

反正沒腦的台北人喜歡背誦垃咭老師的鬼扯也懶得思考,就讓他們繼續穿著破布袋去笑人沒品味吧,台北人應得的。

If $a, b, p, q > 0$ and $\frac{1}{p} + \frac{1}{q} = 1$ then
\[ ab \leq \frac{a^p}{p} + \frac{b^q}{q}\]
why? Because: let $X$ be a random variable with support $\{\{p\log a\}, \{q\log b\}\}$, and since $f(\cdot) = e^\cdot$ is convex, then
\[ ab = \exp(\frac{p\log a}{p} + \frac{q\log b}{q}) = f(E(X)) \leq E(f(X)) = \frac{\exp(p\log a)}{p} + \frac{\exp(q\log b)}{q} = \frac{a^p}{p} + \frac{b^q}{q}\]