IDNote

aStat 實證醫學統計小工具介紹

2010-11-28T23:10:00.001+08:00

(部份內容已先行發表於 mobile01.com)

前言

因為工作/研究上的需求，常常得算些二維列表的統計量，比如說: 卡方檢定，費雪精密檢定等。雖然已經有太多的軟體可以計算這些東西，但是很多資料我都是整理成乾淨的 excel 列表再使用樞紐分析表的功能去檢視變數之間的相關性，網路上工具一堆，可是上網有時也不是頂方便的。後來更因為要計算敏感度，特異性的信賴區間，以及 positive likelihood ratio (PLR) 等等，發現常用的 SPSS 並不提供比例或是 PLR 的信賴區間估計功能，按計算機或是用 excel 表又有點麻煩，因而萌生自己寫個小工具的念頭。

aStat 統計小工具 (A Statistical Calculator for Android) 就是這樣的狀況下的產物。

特點

1. 二維列連表 (2x2 contingency table) 的相關統計量計算，包括:

統計意義檢驗 (Test of significance): 卡方檢定 (Chi-square test), 連續性校正的卡方檢定 (Yate's corrected Chi-square test), 費雪精密檢定 (Fisher's exact test)
相關性強度檢驗 (Test of association): 勝算比 (odds ratio), 相對危險性 (risk ratio or relative risk), 危險差 (risk difference), 包含 95% 信賴區間的數值。
檢驗的正確性評價 (Evaluation of diagnostic power): 計算敏感度 (sensitivity), 特異性 (specificity), 陽性預測值 (PPV: positive predictive value), 陰性預測值 (negative predictive value), 陽性相似比 (positive likelihood ratio), 陰性相似比 (negative likelihood ratio), 亦包含 95% 信賴區間的數值。

2. Person-time 統計推論。

比較不同的組別 (暴露 vs. 非暴露) 的 incidence rate 是不是有差?
計算 incidence rate ratio 以及 incidence rate difference

3. NNT/NNH (number needed to treat/harm)

給定實驗組與對照組的 events rate (i.e., bad outcomes) 計算 absolute risk reduction/increase
提供 NNT/NNH 信賴區間 (Newcombe's method)

4. 單一比例的信賴區間估計

使用三個方法計算比例的信賴區間: exact method, Wilson's score method, Agresti-Coull method

5. 自由設定您所要的信賴區間: 90%, 95%, 99%，能符合大多數人的需求
6. 設定顯示非零小數位數: 2-8 位數。
7. 除了提供一般教科書讓所使用的信賴區間估計法，另外還提供諸如 Agresti-Coull, Wilson, Newcombe 等人所提出的信賴區間方法，特別適用於小樣本數的狀況 (如: N < 40-50)。
8. 免費 (自由使用)

對象
研究人員，公衛/統計相關人員，實證醫學領域相關人士

當前版本
0.6

系統需求
Android 1.6 以上

近期計畫
1. Chi-square test for trend。之後應該暫時不會再更新。
中文介面的部份考量到詞彙的統一性，以及還在增加功能的過程，暫時還沒有做出來，還請見諒。

下載
您可以在 android market 搜尋關鍵字 aStat 就可以了，或是掃描以下的 QR code。

開發網頁
http://aStat.twbbs.org

展示圖片

結語

aStat並不同於一般 android market 上的統計工具，其目的並非要取代 PC 上的統計軟體，而是期待能提供一個快速，方便而準確的評估數據的方法。雖然是小眾市場，但是秉持著好東西要分享的精神，希望對您能有所幫忙! 若您覺得好用，請在 market 上給個好評吧! 此外，若您對區間估計的 algorithms 有興趣，程式碼亦可提供給您。目前所使用的機率分佈函式庫為 JSci，如果有人找到更小的歡迎再告訴我，謝謝!

致謝
1. Google。提供了這麼一個開放的平台，讓想像與時間是唯一的限制。
2. gasolin, 就是因為您的教材: 深入淺出 android ，讓非資訊相關產業的我有辦法完成這個小品!
3. 我們的研究團隊跟老闆。

可轉債自動報價查詢 II

2010-08-02T00:25:00.005+08:00

可轉債, CB, 自動化查詢, 股票, 金融

改進了之前的 script ，全部改成 python 來寫。可由標準輸入檔案讀入想查詢的台灣股市上市上櫃公司代號，完畢後輸出成 html table 並發信至指定的電郵信箱。

usage:

python thisscript.py < stocklist.txt

這支 script 基本上做了下面的事情:

1. 依代號查網頁 (tw.stock.yahoo.com)

2. 解析網頁資料

3. 輸出成 html table

4. 寄信 with attachment

#! /usr/bin/env python
# -*- coding: big5 -*-
import urllib2, re, sys, csv
import smtplib
from datetime import datetime
from email.mime.multipart import MIMEMultipart
from email.mime.text import MIMEText
from email import encoders
from email.utils import COMMASPACE

html = '<html><table border="1" style="border-collapse:collapse;" borderColor=black>\n'
html = html + '<tr><th>代碼</th><th>HTTP</th><th>名稱</th><th>成交價</th><th>漲跌</th><th>買進</th><th>賣出</th>'
html = html + '<th>張數</th><th>昨收</th><th>開盤</th><th>最高</th><th>最低</th></tr>\n'
data = []
for line in sys.stdin:
F1, F2, F3, F4, F5 = '', '', '', '', ''
F1 = re.match('\d+', line).group(0)
url = 'http://tw.stock.yahoo.com/q/q?s=' + F1
myrequest = urllib2.Request(url)
myrequest.add_header('User-Agent', 'Mozilla 5.0')
try:
content = urllib2.urlopen(myrequest)
F2 = 'OK'
except:
F2 = 'Failed'
c = content.read()
m = re.search(r'href="/q/bc\?s=\d+">(\d+.*?)</a>', c)
F3 = m.group(1) if m else '－'
m = re.search(r'nowrap><b>(\d+.?\d+|\d+|－)</b></td>', c)
F4 = m.group(1) if m else '－'
m = re.search(r'nowrap><font color=#......>(▽\d+.\d+|△\d+.\d+|\d+.\d+|－)', c)
F5 = m.group(1) if m else '－'
F6 = re.findall(r'wrap>(\d+.?\d+|\d+|－)</td>', c)
data.append([F1, F2, F3, F4, F5]+F6[1:])
for line in sorted(data, key=lambda x: float(x[3]) if (x[3] <> '－') else x[3]):
result = ''
for x in line:
result = result + x + ','
result = re.sub(',' , '</td><td>', result)
result = '<tr><th>'+result+'</tr>\n'
result = re.sub(r'<th>(\d+)</td>', r'<th>\1</th>', result)
html = html + result

html = html + '</table></html>\n'

tt = datetime.now().timetuple()

me = "sender@email.address"
you = ['recipient1@email.address', 'recipient2@email.address']

# Create message container - the correct MIME type is multipart/alternative.
msg = MIMEMultipart('alternative')
msg['Subject'] = "CB Report " + str(tt[0]) + "-" + str(tt[1]) + "-"  + str(tt[2]) + " "  + str(tt[3]) + ":"  + str(tt[4]) + ":"  + str(tt[5])
msg['From'] = me
msg['To'] = COMMASPACE.join(you)

# Create the body of the message (a plain-text and an HTML version).
# text = """ <html><b>test</b></html> """
# html = """ this is html data """
# Record the MIME types of both parts - text/plain and text/html.
# part1 = MIMEText(text, 'plain')
part2 = MIMEText(html, 'html', _charset='Big5') # _charset is important to avoid bad characters

# Attach parts into message container.
# According to RFC 2046, the last part of a multipart message, in this case
# the HTML message, is best and preferred.
#msg.attach(part1)
msg.attach(part2)

# Send the message via local SMTP server.
s = smtplib.SMTP('localhost')
# sendmail function takes 3 arguments: sender's address, recipient's address
# and message to send - here it is sent as one string.
s.sendmail(me, you, msg.as_string())
s.quit()

可轉債自動報價查詢

2010-07-21T14:53:00.003+08:00

關鍵字: perl, python, 可轉債 (convertible bond)

上市公司及上櫃公司的可轉債約有200多卷，可由公開資訊觀測站查詢相關的訊息。以可轉債的特性而言，比較有用的是跌到票面價格以下的債卷，但必須要常常追蹤價格以及市場消息，對於非專業人士實在力有未逮。如果要一筆一筆輸入 Yahoo 的投資組合來 monitor，也是耗時又費力。因此在一台 linux 上面寫了個 perl 的查詢器，再放到 crontab 來定時自動查詢，方便性增加不少。在這邊做個紀錄。

下面是一開始用的 perl script，可以用標準輸入讀入可轉債代號文字檔案 (一個代號一列)，查詢後輸出到標準輸出 (csv tab 分隔文字格式)。以下的方法也可以用來查詢一般的股票。

#!/usr/bin/perl -w
use LWP::UserAgent;
use HTTP::Request;
use HTTP::Response;
$ua = new LWP::UserAgent; # 產生 UserAgent 物件
my $year = (localtime)[5] + 1900; # year
my $month = (localtime)[4] + 1; # month
my $day = (localtime)[3]; # day
my $date = $year."-".$month."-".$day;
print "代碼"."\t"."HTTP"."\t"."名稱"."\t"."成交價"."\n";
foreach (<>) {
chomp($_);
$request = new HTTP::Request('GET', 'http://tw.stock.yahoo.com/q/q?s='.$_); # 產生 Request 物件
$response = $ua->request($request); # 開始抓取網頁，並將結果傳會 $response
if ($response->is_success) { # 若抓取網頁成功，則印出 HTML 原始碼
   my $html = $response->content;
   if ($html =~ m{href="/q/bc\?s=\d+">(\d+.*?|－)}) {
   print $_."\tOK\t".$1."\t";}
   else {print $_."\tOK\t"."抓不到\t";}
   if ($html =~ m{nowrap>(\d+.?\d+)}) {
   print $1."\n";}
   else {print "抓不到\n";}
   } else { # 若抓取網頁不成功，則印出錯誤訊息
   #print $response->error_as_HTML;
   print $_."\t"."網頁讀取有問題"."\n";
}
}

配合上 crontab 執行以下的 perl script 可以固定每天定時查詢後發信至指定的信箱，包括日期跟時間

#!/usr/bin/perl -w
my $year = (localtime)[5] + 1900; # year
my $month = (localtime)[4] + 1; # month
my $day = (localtime)[3]; # day
my $hour = (localtime)[2]; # hour
my $minute = (localtime)[1]; # minute
my $second = (localtime)[0]; # second
my $cbdate = $year."-".$month."-".$day." ".$hour.":".$minute.":".$second;
my $output = "/tmp/"."cblist-".$year."-".$month."-".$day."-".$hour."-".$minute."-".$second.".txt";

`/path-to-script/queryscript.pl < /convertible_bond_list.txt > $output`;
sleep(25);
open IN, ("< $output");

$to='your@email.address';
$from= 'your@email.address';
$subject='CB Report '.$cbdate;

open(MAIL, "|/usr/sbin/sendmail -t");

## Mail Header
print MAIL "To: $to\n";
print MAIL "From: $from\n";
print MAIL "Subject: $subject\n\n";
## Mail Body
foreach () {
print MAIL $_;
}

close(MAIL);

查完後發現還是不夠完美，我希望可以針對成交價由小至大做排序，這樣只要看前面的幾個可轉債就好了，所以把查詢的 script 用 python 來改寫，輸出 csv 逗號分隔檔至標準輸出裝置，另外增加了 yahoo 上可查詢到的其他欄位。

#! /usr/bin/env python
# -*- coding: big5 -*-
import urllib2, re, sys, csv
print '代碼,HTTP,名稱,成交價,漲跌,買進,賣出,張數,昨收,開盤,最高,最低'
data = []
for line in sys.stdin:
   F1, F2, F3, F4, F5 = '', '', '', '', ''
   F1 = re.match('\d+', line).group(0)
   url = 'http://tw.stock.yahoo.com/q/q?s=' + F1
   myrequest = urllib2.Request(url)
   myrequest.add_header('User-Agent', 'Mozilla 5.0')
   try:
   content = urllib2.urlopen(myrequest)
   F2 = 'OK'
   except:
   F2 = 'Failed'
   c = content.read()
   m = re.search(r'href="/q/bc\?s=\d+">(\d+.*?)', c)
   F3 = m.group(1) if m else '－'
   m = re.search(r'nowrap>(\d+.?\d+|\d+|－)', c)
   F4 = m.group(1) if m else '－'
   m = re.search(r'nowrap>(▽\d+.\d+|△\d+.\d+|\d+.\d+|－)', c)
   F5 = m.group(1) if m else '－'
   F6 = re.findall(r'wrap>(\d+.?\d+|\d+|－)', c)
   data.append([F1, F2, F3, F4, F5]+F6[1:])
for line in sorted(data, key=lambda x: float(x[3]) if (x[3] <> '－') else x[3]):
   result = ''
   for x in line:
   result = result + x + ','
   print result[:-1]

這個 script 有幾個問題
1. 在比大小的時候，python 內建的排序遇到浮點數會被當成字串，因為一開始我們查詢回來的資料就是字串，因此比大小的時候要 type casting。
2. 有用到 key，所以必須要用 python2.4 以上。

單一比例的信賴區間估計 (Confidence Interval Estimation for Proportions)

2010-05-27T01:12:00.012+08:00

整理一下幾種單一比例的信賴區間估計方法。

令

X 代表成功次數

N 代表試驗次數

p = X/N 代表比例

z (alpha/2) = 1.960 at alpha level = 0.05 (i.e., 95% confidence interval)

sqrt = 開平方根

I. Wald method (Wald interval)

教科書裏面最常提到的為 Wald method，又稱為 Wald Interval。

下限 (LL) = p + z * sqrt (p(1-p)/N)

上限 (UL) = p - z * sqrt (p(1-p)/N)

Wald Interval 雖便於計算，但還是有使用上的限制。由於是透過常態分佈近似而來的區間，在樣本數不足，或是 p 接近 0 或 1 時，計算得出的並不是真正的 95% 信賴區間，有 overshooting (大於一) 或 degeneracy (小於零) 的問題。甚至有人證明即使樣本數很大，也不能保證其區間機率 (coverage probability) 近似 95%。但是好算好教好懂是它的最大優點，大部分教科書也是以 Wald Interval 為例子。

II. Exact method (Clopper-Pearson Interval)

由 Clopper 與 Pearson 提出的方法，計算方法較為複雜，可由查表或數值方法或跟 F distribution 的關係作計算。現在電腦很方便，excel 運算功能強大，利用 excel 的 F 分佈反函數即可算出特定單一比例的上下限。以 excel 的內建函數 FINV 則 95% 上下限可用下面的式子表示:

LL = IF(X=0,0,X/(X+(1+N-X)*FINV(0.025,2*(1+N-X),2*X)))
UL = IF(X=N,1,(X+1)*FINV(0.025,2*(X+1),2*(N-X))/(N-X+(X+1)*FINV(0.025,2*(X+1),2*(N-X))))

對於一般的醫學研究，可以用這個方法快速求出某個比例的 95% 信賴區間，例如二維列聯表的 sensitivity 或 specificity，這也是 Graphpad 統計軟體計算所使用的方法。缺點是可能 over conservative。

III. Wilson's score method (Wilson interval)

LL = (2*n*p + z*z - z*sqrt(z*z+4*n*p*(1-p)))/(2*(n+z*z))
UL = (2*n*p + z*z + z*sqrt(z*z+4*n*p*(1-p)))/(2*(n+z*z))

這個 interval 被近來許多的統計學家用模擬的方法證實，coverage probability 較為接近 95%，即使在樣本數少的狀況下。另外，計算上也還不算太繁複，相較於 exact method 所估計出來的區間比較窄一些。

IV. Modified Wald method (Agresti-Coull Interval)

在 Newcombe 比較七種方法和推崇 Wilson's score method 的差不多時間 (1998)，UF 的一個統計學教授 Alan Agresti 發表文章，原本的 Wald method 經由分子加二，分母加四 (作者稱為 pseudo-observations) 可以大大地改進原本的 Wald method 的 coverage probability，使之接近 95% 名目機率。作者還進一步闡述，經過變換，其實加上 pseudo-observations 使得原本的 Wald interval 可以作為 Wilson's score interval 的近似。也就是說，大部分的狀況，都可以用這個改變後的 p' 數值代入原本的 Wald method 得到接近 95% coverage probability 的信賴區間。這個方法有可能產生 overshooting，也就是區間大於一或小於零的狀況，這時候要把他裁掉。

令 p'= (X+2)/(N+4)

下限 (LL) = p' + z * sqrt (p'(1-p')/N)

上限 (UL) = p' - z * sqrt (p'(1-p')/N)

其實加上的 pseudo-observations:
2 是 z*z/2 的近似值
4 是 z*z 的近似值
因此也可以拿來套用在不同的 alpha level。

在 Agresti 與 Caffo 合著的文章 (2000) 中提到，即使沒有資料也可以計算信賴區間，因為分子為二分母為四可以帶入運算，這是 Agresti 提出的方法一個有趣的地方。

CONCLUSIONS

以應用上來說，若要好算又嚴謹，可以直接使用 Agresti-Coull interval。若要名字好聽，那就用所謂的 exact confidence interval，這是多數統計軟體提供的，也是所謂的金標準。雖然它已經被證實不比起 Agresti-Coull interval 跟 Wilson's interval 精確，而且它的信賴區間過度保守 (conservative ~ 比較寬)，不過一般的接受度應該是很高的。

個人是傾向 Newcombe 提倡的 Wilson's interval 或 Agresti-Coull interval，除了用 excel 方便計算，區間在小樣本的時候看起來窄一點之外，更重要的是它們都可以應用在建立兩個比率的差異的信賴區間。最後，忘了在哪一篇文章中曾提到，樣本數若偏小 (N<40)，Wilson's interval 可能是最好的選擇。

REFERENCES:

Newcombe RG. Two-sided confidence intervals for the single proportion: comparison of seven methods. Stat Med. 1998;17(8):857-872.
Clopper-Pearson Interval. Available at: http://en.wikipedia.org/wiki/Binomial_proportion_confidence_interval#Clopper-Pearson_interval [Accessed May 26, 2010]
Agresti A, Coull BA. Approximate Is Better than "Exact" for Interval Estimation of Binomial Proportions. The American Statistician. 1998;52(2):119-126.
Agresti A, Caffo B. Simple and Effective Confidence Intervals for Proportions and Differences of Proportions Result from Adding Two Successes and Two Failures. The American Statistician. 2000;54(4):280-288.

A 2x2 contigency table calculator

2010-05-24T19:59:00.000+08:00

http://sites.google.com/a/idman.twbbs.org/linkoucgmhid/files/2x2.zip

Updated. 2010/05/24

1. Confidence interval estimation for two proportions (unpaired samples) using Newcombe's procedure (1998), which was suitable for both equal n or unequal n.

2. Confidence interval estimation for a single proportion with three methods: Exact method (reverse F distribution), modified Wald method (Agresti, 1998), and Wilson's score method.

Cox 比例風險回歸 (Cox Proportional Hazards Model) 的應用

2010-01-17T23:26:00.001+08:00

The prototype of hazard function：h(t,X) = h0(t) exp (b1X1 + b2X2 + b3X3 + ... + bkXk)

h0(t): baseline hazard function
t: time to event
X: covariates
b: coefficients
exp: exponential function

例子：加護病房病人產生黴菌尿後，有些接受治療有些不治療。經過一段時間後，他們的預後有沒有差別？換句話說，他們的死亡跟有沒有用藥到底有沒有關聯？

此例並非隨機雙盲試驗，我們無法只比較 Kaplan-Meier curves 來得到結論。用不用藥物受到病人的本身身體狀況，以及疾病嚴重度的影響。有沒有用藥又進一步受到時間因素的干擾：有些人來不及用藥，有些剛開始沒有用藥，後來因為狀況不好又把藥物加上去。在這裡我們可以用 Cox proportional hazards model 來考慮這些干擾因素的影響，以評估藥物與病人預後（死亡）的相關性。

Cox porportianl hazards model 的特性 [1-3]：
1. 不需事先知道 ho(t) 這個 baseline hazard function，也不要求服從特定的機率分佈。因此它有 non-parametric model 的特性。

2. 解釋變數或共變數（covariates）的部份有 parametric model 的特性。所以有些人說 Cox 回歸屬於 semi-parametric models。

3. 解釋變數的風險（hazards）是一個定值，不隨著時間變化。亦即，符合比例風險假定（proportinal hazards assumption）。

4. 解釋變數可以是連續變項，類別變項或是有次序的變項。

應用 Cox proportional hazards model 來做回歸我們必須注意[1]：

1. 研究樣本數要足夠，一般需要自變數個數的十倍以上。若研究得出的結果違背常理，則必須考慮到模型設計上面的缺陷，不一定是科學上的新發現。

2. 要符合比例風險假定 (proportional hazards assumption)。任意兩個個體在同一個變數的風險比 (risk ratio) 是一個定值，不隨著時間變化。不符合比例風險假定則可分層或是作 Cox regression with time-depedent covariates。以前面黴菌尿的例子，治療會受到時間因素的影響，因此必須作時間協變的 Cox regression。
3. covariates 若是連續變數，可考慮轉成類別變數再分析。比如說血球數目，回歸的時候單位是用顆來計算，所得到的結果意義就代表，血球數每增加或減少一個所帶來的風險變化。
4. covariates 若是沒有次序性的，例如居住地等等，則需加入 dull variables，以免造成解釋上的困難。
5. 變數的篩選需考慮：臨床上有意義的重要的變數要放入模型，無關的不擺進去，或是先作單變數分析再篩選等等。

判斷是否符合風險比例假定有幾個方法 [1-3]：
1. 圖形法（graphic methods）。

直接比較 survival curves 是最直觀的方法之一。就在 SPSS 裡面直接作 Cox regression 把懷疑受時間影響的變數放進 Strata，covariates 不放，然後畫圖。如果圖形有明顯交錯就不符合風險比例假定。
Log-minus-log plots
Partial residual plots

2. Goodness-of-fit test。查不到怎麼實做。
3. Computational methods, like time-interaction test. 一個符合比例風險假定的 covariates，和時間交互作用後的變項，將其加入 Cox 回歸模型，其回歸係數應該很接近零（和 0 沒有顯著差異）[3]。若 covariate 對時間的交互作用項在 Cox 回歸模型中有統計意義，covariate 就不符合比例風險假定 [1]。

REFERENCES:
1. 生存分析。醫學統計學，第二十三章。人民衛生出版社 2005 年第一版。
2. Proportional hazards model. Available at: http://en.wikipedia.org/wiki/Proportional_hazards_models [Accessed January 16, 2010]
3. Cox regression. Available at: http://faculty.chass.ncsu.edu/garson/PA765/cox.htm [Accessed January 17, 2010]

費雪精密檢定（Fisher's Exact Test）

2010-01-01T20:41:00.000+08:00

好文章：
1. http://www.math.nsysu.edu.tw/outstanding/use/report/h.pdf
2. http://en.wikipedia.org/wiki/Fisher%27s_exact_test
計算工具：

1. http://www.langsrud.com/fisher.htm
2. http://cgmh.idman.twbbs.org/Home/tools

布瓦松分佈與區間估計 (Poisson Distribution and Estimation of Confidence Interval)

2009-09-30T13:04:00.046+08:00

Poisson Distribution 可視為二項式分佈 (Binomial Distribution) 的極限情形。亦即，當試驗次數 (n) 趨近無窮大，p 很小，但一段時間裏面 np 為一不大不小的定值時。從品質管理的角度，如：每日不良產品的數目，不良事件的發生數目等等，每週醫療照顧相關感染的發生次數等等，皆適用於布瓦松分佈。[1]

依從 Poisson 分布的群體，最佳的點估計即為平均值 u。當我們只做了一個單次的觀察，可利用卡方分佈 (Chi square Distribution) 與 Poisson 分布的關係，很快的求出對應的 95% 信賴區間估計 (exact confidence interval)。[2-4]

令 x 為單次的觀察結果 (發生次數)，InvChiSqare (v,y) 函數代表自由度 v 的 Chi Sqare 分布下面積為 y 的機率分布對應的 Chi Sqare 值。則依從 Poisson 分布的 95 % 信賴區間上下限 (exact confidence interval; UL; LL) 與卡方分布有如下的關係：

LL = 1/2 * InvChiSqare (2x, 0.025)
UL = 1/2 * InvChiSqare (2x + 2, 0.975)

當觀察次數不止一次時 (observations = N)，所有的觀察次數合 (n) 依從 Poisson 分布。其 95% 信賴區間上下限可如此表示：[5]

LL = 1/2 * InvChiSqare (2n, 0.025) / N
UL = 1/2 * InvChiSqare (2n + 2, 0.975) / N

References:

Poisson 分布. Available at: http://episte.math.ntu.edu.tw/articles/sm/sm_16_07_1/index.html [Accessed Sep 30, 2009]
Confidence intervals for the mean of a Poisson distribution. Available at: http://www.math.mcmaster.ca/peter/s743/poissonalpha.html [Accessed Sep 30, 2009]
Poisson confidence interval. Available at: http://www.statsdirect.com/help/parametric_methods/pest.htm [Accessed Sep 30, 2009]
An exact method for calculating a confidence interval of a poisson parameter. Mulder et al. American Journal of Epidemiology 117 (3): 337. Available at: http://aje.oxfordjournals.org/cgi/pdf_extract/117/3/377 [Accessed Sept 30, 2009]
Estimating the mean of a Poisson population. Chapter 6.3. Available at Google Books [Accessed Sep 30, 2009]

在 Excel 內增加色票 (Add new colors to the default excel palette)

2009-07-03T13:19:00.002+08:00

Refer to http://www.bloggpro.com/add-new-colors-to-the-default-excel-palette/

選擇顏色的工具: http://colorschemedesigner.com/

在 ubuntu 下播放 HiChannel 頻道

2009-02-15T16:13:00.002+08:00

有人寫了個 script http://fourdollars.blogspot.com/2007/11/hichannel.html

大陸有人改成 gui 在此

近30年美股指數

2009-02-01T16:06:00.000+08:00

Yahoo::Finance::Interactive chart

舊電腦何處去?

2009-02-01T15:20:00.001+08:00

想把舊電腦送給有需要的人...
http://www.3c-dr.com.tw/About-Us.html
http://blog.hfu.edu.tw/blog/hfufreehw

Obama's Inaugural Speech

2009-01-21T23:52:00.001+08:00

Read the full article from CNN.com

Vanguard GNMA

2009-01-18T16:43:00.001+08:00

Google::Finance::VFIIX

GDIPLUS.DLL dependency of AdvStringGrid

2009-01-05T14:29:00.002+08:00

This is an excerpt from the website of TMS Component suite
---
Q22: Users of older operating systems have an error message on application startup related to a missing gdiplus.dll

A22: Either redistribute the Microsoft GDIPLUS.DLL (explained in README.TXT) or remove the gdiplus.dll dependency by commenting the line {$DEFINE TMSGDIPLUS} in TMSDEFS.INC

A Free PDF Splitter and Merger

2008-12-29T12:07:00.001+08:00

http://sourceforge.net/projects/pdfsam/
Written in Java

APACHE II Score and Predicted Mortality Rate

2008-12-06T15:03:00.004+08:00

The following result was based on original APACHE II (Acute Physiology and Chronic Health Evaluation II) score, published by Knaus et al. in 1985.
----
score = predicted mortality reate
0 = 0.0288323804597415
1 = 0.0332141827839264
2 = 0.0382356916987282
3 = 0.0439818435946882
4 = 0.0505461568401609
5 = 0.0580307274067857
6 = 0.0665458910150112
7 = 0.0762094426987653
8 = 0.0871452904871902
9 = 0.0994814113171375
10 = 0.113346978732479
11 = 0.12886854890866
12 = 0.146165230272512
13 = 0.165342828461315
14 = 0.186487056810361
15 = 0.209656033245229
16 = 0.234872441095319
17 = 0.262115898694886
18 = 0.291316235325043
19 = 0.322348474912558
20 = 0.355030345948782
21 = 0.389123033352462
22 = 0.424335648915215
23 = 0.460333532118856
24 = 0.496750045768218
25 = 0.533201073292596
26 = 0.569301043624338
27 = 0.604679084771601
28 = 0.638993887971074
29 = 0.671946054939952
30 = 0.70328705911182
31 = 0.732824402144069
32 = 0.760423002108014
33 = 0.78600323342325
34 = 0.809536301764071
35 = 0.831037763630525
36 = 0.85056000246354
37 = 0.868184382009093
38 = 0.884013650945879
39 = 0.898165005907147
40 = 0.910764060077794
41 = 0.921939828549029
42 = 0.93182073729233
43 = 0.940531590674707
44 = 0.948191389177323
45 = 0.954911868587507
46 = 0.960796628155722
47 = 0.96594072236451
48 = 0.970430604427108
49 = 0.974344325979221
50 = 0.977751914291392
51 = 0.980715864241253
52 = 0.983291696442564
53 = 0.985528544995485
54 = 0.987469748266567
55 = 0.989153424071641
56 = 0.990613016854848
57 = 0.991877809199794
58 = 0.992973393538575
59 = 0.993922102483462
60 = 0.994743398004825
61 = 0.995454220894235
62 = 0.996069302727867
63 = 0.99660144299736
64 = 0.9970617542933
65 = 0.997459878479777
66 = 0.997804176739593
67 = 0.998101896237682
68 = 0.998359315973963
69 = 0.998581874196749
70 = 0.99877427953855
71 = 0.998940607827617

醫事人員積分系統查詢

2008-11-13T13:55:00.001+08:00

http://chs-www.doh.gov.tw/csnpt/index.html
帳號: 身分證號碼
密碼: 生日末四碼加身分證號碼

王永慶給兒女的一封信 (2004-05-20)

2008-11-08T14:58:00.001+08:00

子女們
財富雖然是每個人都喜歡的事物，但它並非與生俱來，同時也不是任何人可以隨身帶走。人經由各自努力程度之不同，在其一生當中固然可能累積或多或少之財富，然而當生命終結，辭別人世之時，這些財富將再全數歸還社會，無人可以例外。因此如果我們透視財富的本質，它終究只是上天託付作妥善管理和支配之用，沒有人可以真正擁有。面對財富問題，我希望你們每一個人都能正確予以認知，並且在這樣的認知基礎上營造充實的人生。
我本出身於貧困家庭，歷經努力耕耘，能夠有所成就。在一生奮鬥過程中，我日益堅定的相信，人生最大的意義和價值所在，乃是藉由一己力量的發揮，能夠對於社會作出實質貢獻，為人群創造更為美好的發展前景，同時唯有建立這樣的觀念和人生目標，才能在漫長一生當中持續不斷自我期許勉勵，永不懈怠，並且憑以締造若干貢獻與成就，而不虛此生。
基於這樣的深刻體會，因此我希望所有子女也都能夠充分理解生命的真義所在，並且出自內心的認同和支持，將我個人財富留給社會大眾，使之繼續發揮促進社會進步，增進人群福祉之功能，並使一生創辦之企業能達到永續之經營，長遠造福員工與社會。與此同時，我也殷切期盼所有子女，在創業與日常生活中，不忘以服務奉獻社會、造福人群為宗旨，而非只以私利做為追求目標，如此才能建立廣闊和宏偉的見識及胸襟，充分發揮智慧力量，而不負於生命之意義。

本文由經濟日報提供。經濟日報為紀念王永慶先生辭世，出版「王永慶紀念特刊」，今天隨報附贈。

採用 UCT 演算法的圍棋程式 Mogo

2008-09-28T09:32:00.001+08:00

http://www.lri.fr/~gelly/MoGo_Download.htm

Unpredictable Fortune

2008-09-18T12:38:00.003+08:00

以下為 2008/07/23 舊聞

【鉅亨網記者陳慧琳‧台北】台股昨(22)日小跌，融資波段減幅儘管昨日已首度超越台股跌幅，但是融資還是持續殺出16億元，餘額來到2741億元今(2008)年低水位，並創近17個月的新低，法人表示，今日台股開盤大漲，除了因美股走揚，最主要就是反映籌碼換手成功，底部訊號確定。

融資餘額持續每日每日的遞減，相較於 5月22日最高峰的3686億元，波段洗出高達945億元，幅度 25.64% 。至於大盤指數則從 5月20日當天最高點的9309點下修，因前天暴漲 270點，融資減幅就已首度超越大盤跌幅，再計算至昨天收盤的7065點，波段跌幅為24.1%。

融資餘額已降到2741億元，此水位是從去年 3月以來新低的水位，代表散戶動向的融資餘額減肥成功，市場籌碼逐漸穩定，有利股市反彈。

若分析近 1個月融資洗出最多個股，法人表示，融資減肥主要集中在跌深的面板、低價電子及部分金融股，奇美電 3009(TW) 、友達 2409(TW) 累計20個交易日融資分別減58240、39757張；其餘錸德 2349(TW) 、歌林 1606(TW) 、華邦電 2344(TW) 、新纖 1409(TW) 及彰銀 2801(TW) 、國泰金 2882(TW) 。

至於向來是散戶最愛的IC設計族群，前波跌勢兇猛，融資猛殺，不少個股股價腰斬、甚至還有打到 3折，不過法人表示，近期IC設計族群率先止跌，包括聯家軍的聯陽 3014(TW) 、原相 3227(TW) 、揚智 3041(TW) 等，隨著籌碼沉澱，股價都已經在6月底、7月初領先大盤落底反彈，就已是浮現台股底部訊號。

法人表示，整體來看台股 7月見低點、8月整理，9 月台股伴隨電子業擺脫低迷，旺季來臨，預期可望回升。電子族群中，NB相關族群業績成長佳，逢低仍可分批承接； TFT面板相關族群在報價尚未止穩之前，宜維持觀望。

台塑砸150億護盤分析師：短多

2008-09-14T09:57:00.000+08:00

【經濟日報╱記者黃俊苔／台北報導】

2008.09.13 01:40 am

台塑集團投入資金護盤，宣示穩定投資人信心意味十足。不過，分析師認為，護盤僅能代表股價來到長期相對低點，不代表股價立刻有表現機會，台塑集團動作可能帶動短期反彈效果，股價仍需回歸基本面看待。

台塑集團動用150億元資金護盤，為民國93年之後第二度進行護盤。美嘉資產管理總經理杜富國表示，這次護盤雖對股價沒有強烈鼓舞作用，但以台塑集團過去所建立的誠信形象，還是會受到市場肯定，也代表台塑集團股的投資買點浮現。

聯邦投信副總經理黎方國表示，企業經營階層理是最了解公司營運與產業狀況的人，買回自家的股票，代表肯定公司的經營面，具有宣示與振奮人心的效果，就理性投資人而言，確實可視為買點。

不過，四年前台塑集團買回自家公司股票，當時台灣的經濟成長率為6.07%，而今年則是備受質疑的4.3%，且台塑集團93年買進股票後，石化業景氣進入旺季，整體經濟環境優於現在，這次台塑集團再進場護盤，實質效益仍待觀察。

國際大通投顧分析師連乾文指出，去年亞洲乙烯新增產能僅有486萬噸，而今年第四季與明年上半年皆預期產能大增，導致乙烯價格快速滑落，自然壓縮塑化類股的獲利空間。從產業面來看，台塑集團確實面臨較艱難的困境。

連乾文說，短線大盤如果出現止跌，台塑四寶有反彈機會，不過，就中期表現而言，基本面確實不佳，股價難有太好表現。

阿特拉斯的軟體倉庫

2008-08-08T21:07:00.001+08:00

有很多自製好用的影音相關軟體

譬如: flv 轉檔工具

Refer to http://www.atlas-zone.com/software/index.html

Delphifans

2008-07-28T16:42:00.001+08:00

http://www.delphifans.com/

Proud of You

2008-07-13T23:42:00.003+08:00

Keywords: song, English
原唱: Fiona Fong (馮翠樺)
作詞/作曲: 陳光榮

References

http://shenez.blogspot.com/2006/11/proud-of-you-fiona-fung.html 有正確的歌詞
http://www.wretch.cc/blog/tangmf/2488900
http://blog.xuite.net/ronnyqa00/test/7705138 原唱 mp3

Premenopausal osteoporosis

2008-07-09T09:20:00.002+08:00

keywords: premenopause, woman, osteoporosis, synposium
refer to http://www.postgradmed.com/issues/2006/06/008/

Google Apps + twbbs.org

2008-07-06T13:36:00.000+08:00

refer to http://blog.pixnet.net/s25655/post/19608497

裝潢

2008-07-04T23:26:00.000+08:00

參考一下
全能住宅改造王--東王漢宮 @ 室內設計家居‧人文雅致‧風采權釋 PIXNET 痞客邦

寶來 ETF 相關資訊

2008-07-02T23:46:00.002+08:00

宜蘭礁溪之旅

2008-07-02T12:38:00.000+08:00

2008/07/12-13

住宿

宜蘭礁溪蘭陽溫泉大飯店

交通 - 客運

火車

新玩具 - JVC GZ-MG275

2008-07-01T21:18:00.003+08:00

NT$ 21,000 at Nova, Taipei
Standard equip. with a second Li-ion battery, a photo bag, one UV len, and one cleaning set.

Refer to Mobile01.com 的開箱文

People said there:
....
這台會高光溢出...(????)

如果錄影時，太亮或太暗，可使用

1.背光補償功能
2.單點曝光功能
3.調整快門速度
4.調整光圈比例

最大光圈F1.2

1/3.9CCD感光元件、10倍光學變焦。採用JVC獨家研發晶片，影像清晰又銳利，在黯暗的環境下也能拍攝出清晰的畫面喔。最大光圈F1.2，比一般F1.8明亮約2倍，室內拍攝最佳幫手

2008/07/03 補充
發現 JVC 拍出來的檔案格式為 .mod, google 後發現這個格式以 16:9 拍攝後在 xp 下直接播放竟會變成 4:3 格式! 必須再轉檔, 或使用 sdcopy 工具處理.
參考

數位照片處理工具 - 光影魔術手

2008-07-01T02:01:00.001+08:00

這個大陸人有夠牛!
原站: http://www.neoimaging.cn/tw/index1.htm
功能
正片效果
數碼補光
晚霞渲染
褪色舊相
黑白效果
正片負沖
負片效果
夜景效果
白平衡一指鍵
CCD死點修復
界面：對比模式和處理模式
色階和曲線
人像褪黃
水印
花樣邊框
撕邊效果
數碼減光
變形校正
証件照排版
制作組合圖
高ISO去噪
柔光鏡模擬
包圍曝光三合一
嚴重白平衡錯誤校正
紅飽和衰減
人像美容
IE魔朮圖
去紅眼
影樓風格人像
素描

Bookstart - 台北市立圖書館::嬰幼兒父母學習講座

2008-06-30T20:15:00.002+08:00

報名

資料來源: 台北市立圖書館
一、活動目標
（一）為提倡閱讀向下紮根理念，讓嬰幼兒能夠從出生即接觸閱讀的世界，以提升兒童對閱讀的興趣及喜好，並增進親子互動。（二）讓新生兒父母及社會大眾瞭解Bookstart的意義，並學習如何引導新生兒藉由閱讀認識新的事物，並提升其學習能力。
二、活動內容
（一）嬰幼兒父母學習講座：開啟寶寶的閱讀之門
1.講座內容以「閱讀禮袋」中的圖畫書、導讀手冊、推薦書單為課程內容。經由說明、投影片介紹，讓家長了解Bookstart閱讀起步走的精神、價值及意義；並以閱讀禮袋中的推薦書單為例，介紹嬰幼兒閱讀行為發展、親子共讀策略。 2.對象：以設籍臺北市且家有6~18個月的嬰幼兒家長為對象。3.報名時間：自即日起受理報名，至各梯次報名額滿為止。 4.報名方式：採現場、傳真或線上方式報名。（1）現場報名總館：每週一至週五，上午9：00至12：00，下午2：00至5：00，持戶口名簿正本或三個月內謄本親至推廣課（臺北市建國南路2段125號）填寫報名表並辦專屬借閱證。分館：於開館時間內，持戶口名簿正本或三個月內謄本親至民生分館、萬華分館、景新分館、東湖分館、士林分館、永建分館及大直分館綜合服務臺填寫報名表並辦證。（2）傳真報名填妥報名表後與戶口名簿一併傳真至推廣課或辦理之分館。（3）線上報名於本館網站「市圖首頁」中，進入「線上申辦」網頁，並於「活動線上報名」網頁申請，並將戶口名簿傳真至推廣課或辦理之分館。（二）贈送「閱讀禮袋」1. 凡申辦「Bookstart嬰幼兒借閱證」，且參加「嬰幼兒父母學習講座」之嬰幼兒父母，可於講座活動結束後領取「閱讀禮袋」1份（閱讀禮袋內有2本專家推薦優質圖畫書、1本導讀手冊及推薦書單）。2. 每位寶寶以領取1份閱讀禮袋為限，年滿6個月且未超過18個月之年齡限制以領取閱讀禮袋當天為計算基準。
三、講座梯次、日期及辦理地點一覽表（PDF檔 DOC檔)
四、注意事項
（一）參加者請事先申辦嬰幼兒借閱證或當天攜帶戶口名簿正本於活動現場申辦，以便領取「閱讀禮袋」1份。（二）如有相關問題可逕洽本館推廣課02-27552823#2117；民生分館02-27138083；東湖分館02-26323378；萬華分館02-23391056；景新分館02-29331244；永建分館02-22367448；士林分館02-28361994；大直分館02-25336535。
下載報名表（PDF檔 DOC檔)
線上報名

TAdvStrigGrid's dependency of gdiplus

2008-06-30T16:43:00.000+08:00

Modify TMSDEFS.INC

Comment {$DEFINE TMSGDIPLUS} to disable GDI+ use in TAdvStrigGrid
----
//{$DEFINE TMSGDIPLUS} // comment this line to disable GDI+ use in TAdvStringGrid

更改圖檔名稱工具 - Tools for modifying filenames with the datetime of the images

2008-06-22T17:11:00.007+08:00

經過一段時間的尋找與試用, 有以下的幾個選擇:

Image Renaming

Xnview

Siren

Image Renaming

優點: 可以簡單的選擇多轉換的日期格式

缺點: 要安裝 .Net 2.0 Framework, 轉換檔名的速度有點慢

XnView

優點: 功能強大, 連 mp3 都可以轉, 非常多樣的檔名版型組合

缺點: 要花一點時間找一下轉檔名的功能表單

Siren

優點: 支源音樂檔圖檔, 專門的改檔名免費工具
缺點:

台東 - 海草健康輕食

2008-06-19T15:34:00.001+08:00

關鍵字: 餐飲, 健康, 台東, 評價

特色

口味相較於台東一般的餐廳較為清淡
創意料理
手工麵條 (海草)
廚房, 內裝乾淨整齊
研磨咖啡免費
備嬰兒座椅

部落格

eth0 problems

2008-06-13T16:43:00.003+08:00

Keywords: fixed IP, eth0, network config, runlevel

After B2D linux (B2D pureKGB 20080513) was installed, eth0 was not up automatically during booting up. I do not know exactly the reason causing this. However, the problem could be solved by setting the networking service manually.

bash#cd /etc/rc5.d/
bash#ln -s ../init.d/networking S16networking

Install PIL (Python Image Library) at Bluehost.com

2008-06-11T23:15:00.003+08:00

I tried to install PIL1.1.6 at Bluehost and encoutered the following error messages:


bash#python setup.py install --prefix=$HOME
......
......
_imagingtk.c:20:16: tk.h: No such file or directory
_imagingtk.c:23: error: syntax error before '*' token
_imagingtk.c:31: error: syntax error before "Tcl_Interp"
_imagingtk.c:31: warning: no semicolon at end of struct or union
_imagingtk.c:32: warning: data definition has no type or storage class
_imagingtk.c: In function `_tkinit':
_imagingtk.c:37: error: `Tcl_Interp' undeclared (first use in this function)
_imagingtk.c:37: error: (Each undeclared identifier is reported only once
_imagingtk.c:37: error: for each function it appears in.)
_imagingtk.c:37: error: `interp' undeclared (first use in this function)
_imagingtk.c:45: error: syntax error before ')' token
_imagingtk.c:50: error: `app' undeclared (first use in this function)
_imagingtk.c:50: error: syntax error before ')' token
_imagingtk.c: At top level:
_imagingtk.c:55: warning: parameter names (without types) in function declaration
_imagingtk.c:55: error: conflicting types for 'TkImaging_Init'
_imagingtk.c:23: error: previous declaration of 'TkImaging_Init' was here
_imagingtk.c:55: error: conflicting types for 'TkImaging_Init'
_imagingtk.c:23: error: previous declaration of 'TkImaging_Init' was here
_imagingtk.c:55: warning: data definition has no type or storage class
_imagingtk.c:57: error: syntax error before '&' token
error: command 'gcc' failed with exit status 1

The problem can be corrected by not importing tkinter. Modify the program block of the file "setup.py"


try:
    import _tkinter
except ImportError:
    _tkinter = None

to the following


_tkinter = None

With this, PIL would not be compiled with tkinter support.

TPMath - A Scientific Library Written in Pascal

2008-03-30T22:27:00.001+08:00

http://www.unilim.fr/pages_perso/jean.debord/tpmath/tpmath.htm

IT is a library of scientific programs written in Pascal. Available in several versions according to the compiler, it is proposed as an alternative to the famous Numerical Recipes, for which the Pascal version is no longer developed, or to the Borland Numerical Methods Toolbox which is no longer available.

An Opensource Equation Editor

2008-03-13T20:07:00.000+08:00

http://mathcast.sourceforge.net/home.html

白崇亮「勇於真實」

2008-03-13T00:04:00.003+08:00

http://blog.roodo.com/josephpai

A Collection of Delphi Code Libraries - Fundamentals

2008-03-11T19:04:00.001+08:00

http://fundementals.sourceforge.net/index.html

The collection consists of:

Utility functions:
String operations, Dynamic array operations, System functions and Date & Time operations.
Unicode functions:
Unicode codecs and Unicode character and string functions.
Streams:
Common stream implementations that includes parsing support and binary packing.
Datastructures:
Commonly used data structures, including efficient array and dictionary implementations.
Sockets:
UDP and TCP server and client classes. Support for asynchronous or threaded modes.
Mathematics:
Units included for vectors, matrices, rational numbers and complex numbers.

A Cross-platform Packer For Executables（跨平台的執行檔壓縮工具）

2008-03-06T13:41:00.007+08:00

在 Lazarus 網站的 FAQ 中提到，可使用 UPX 這個工具，加上移除 GNU debugger 的 symbols（使用 strip）使得 freepascal 編譯出的 windows 執行檔變小。

經測試壓縮效率比 PECompact 差一些，不過差別在於：

UPX 為 GPL。
UPX 跨平台。
UPX 支援許多不同平台的執行檔格式。（win32, linux, linux kernel, mac, even Playstation 1 !)

以下轉貼自原始網站 http://upx.sourceforge.net/

UPX is a versatile executable packer with the following features:

excellent compression ratio: typically compresses better than WinZip/zip/gzip, use UPX to decrease the size of your distribution !
very fast decompression: ~10 MB/sec on an ancient Pentium 133, ~200 MB/sec on an Athlon XP 2000+.
no memory overhead for your compressed executables because of in-place decompression.
safe: you can list, test and unpack your executables. Also, a checksum of both the compressed and uncompressed file is maintained internally.
universal: UPX can pack a number of executable formats.
portable: UPX is written in portable endian-neutral C++
extendable: because of the class layout it's very easy to add new executable formats or new compression algorithms
free: UPX is distributed with full source code under the GNU General Public License, and may be used freely even with commercial programs.

You probably understand now why we call UPX the "Ultimate Packer for eXecutables".

UPX is based on experience with our previous packers DJP and lzop and uses the NRV compression library.

Pascal Tutorial

2008-03-06T01:26:00.001+08:00

http://www.taoyue.com/tutorials/pascal/
Pascal 的入門教學。

隨身碟病毒 - kavo kava tavo tava etc.

2008-02-14T22:45:00.006+08:00

最近家中電腦感染了隨身碟病毒，起先沒有發現（Norton 2005）。notebook 中安裝的是 Avira Antivir Personal Edition 才偵測出隨身碟已受到感染。參考了幾篇文章之後動手清除。其中比較有價值的文章分別是：

行政院資安論壇中的一個討論串。其中有一篇說明 kavo virus 的詳細清除步驟。甚至提到病毒會自我更新病毒碼！
台灣論壇中的一篇文章。包括其他的隨身碟病毒的各種名稱：tava/kava/avpa/amva/mmva/mnsa。但解法不如 1. 的完整。
您不可不知的KAVO病毒與解法。內有提到使用 HiJackThis 工具，來掃瞄並修復系統。此項工具可從趨勢網站下載。

系統是 XP Home，步驟與資安論壇中的文章類似但略有不同，紀錄如下：

關閉系統還原。
關閉網路。不一定要進入安全模式。
修復登錄檔，把顯示隱藏檔案的功能打開。執行 regedit 把 HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\Windows\CurrentVersion\Explorer\Advanced\Folder\Hidden\SHOWALL 值改成 1。
移除 HKEY_CURRENT_USER\Software\Microsoft\Windows\CurrentVersion\Explorer\MountPoints2
先執行 cmd.exe 後用工作管理員把 explorer.exe 這個 process 殺掉，之後在 cmd 視窗下再執行 explorer。因為其中有 kavox.dll 會被 explorer 調用。
在 cmd 視窗下執行 explorer 以開啟 explorer，如：explorer c:。將資料夾選項中的隱藏檔案功能取消，如此可看到病毒寫入的隱藏檔案，將之清除。如：

根目錄下面的 autorun.inf, ntdelect.com (NOT ntdetect.com!), 數字/亂碼.bat/com。這些檔案不一定都存在。
C:\windows\system32\kavo.exe
C:\windows\system32\kavo1.dll
C:\windows\help 下也要檢查。（？）

使用 cmd 視窗進入使用者的家目錄下的 Local Settings 之後執行 "rd /s /q Temp" 清除暫存目錄下的檔案。
每個磁碟根目錄都用 cmd 視窗下執行 explorer 磁碟機代號清除病毒檔案。
最後用 HiJackThis scan 之後勾選 kava/tava/等等將之修復。
重新開機。

VirusTotal - 免費線上病毒和惡意軟體掃瞄

2008-02-04T13:03:00.000+08:00

網站 http://www.virustotal.com/

dvd 轉 divx

2008-02-03T16:53:00.000+08:00

xmpeg

Binary Logistirc Regression (SPSS)

2008-01-27T15:05:00.001+08:00

基本算式: Ln (Odd) = Ln (Y'/1-Y') = a + bX
Y' represents the predicted probability of the event which is coded with 1 (one).

For dichomotous variables, the calculated odds ration (OR) of a logistic model is

(the odd in favor of the event coded with 1) / (the odd in favor the event coded with 0)

筆記：

Continuous variables in a logistic regression model - the importance of the "scale" - description 對應的某單位變動的勝算比，使用的單位不同，勝算比也不同。如：早產兒的體重使用 1 gram vs. 1 kilogram，odds ratio 有不同的大小。

參考資料：

生物谷的教學
Multivariate Analysis with SPSS - 一個教學資訊豐富的網頁，包括如何使用 spss 做 binary logistic regression analysis。

DBGrid 加強

2008-01-10T01:25:00.000+08:00

從 delphi.ktop.com.tw 看到的 SMDBGrid component freeware。
可補原 TDBGrid 之不足。

Firefox 的 Menuitem 問題。

2008-01-09T17:07:00.001+08:00

不知為何在瀏覽器的視窗下面出現了一排紅色的 menuitem 錯誤。Google 之後發現是 JAVA 寫入語系的時候「寫錯了」。

修改 C:\Program Files\Mozilla Firefox\extensions\{CAFEEFAC-0016-0000-0001-ABCDEFFEDCBA}\chrome.manifest 這個檔案，先關閉 firefox，再把 zh_TW 改成 zh-TW 之後重新啟動 firefox 即可解決。

投資金律 (The Four Pillars of Investing)

2007-12-30T21:18:00.001+08:00

I. The First Pillar - Investing Theory

Risk and Reward

Growth stock vs. value stock

A Random Walk Down Wall Street

2007-12-25T20:51:00.001+08:00

2007-09-10

2007-09-10T08:44:00.000+08:00

音樂劇
little Jekyll and Hyde
Jekyll and Hyde: 改編音樂劇介紹
Chien-Ming Wang
Into the Woods 介紹

An Index to Numerical Methods in Pascal

2007-07-14T06:43:00.001+08:00

http://www-rab.larc.nasa.gov/nmp/nmpIndex.htm

Random Numbers With A Good Uniform distribution

2007-07-11T00:32:00.001+08:00

http://www.codepost.org/browse/snippets/105

Javascript Tutorials

2007-07-09T20:50:00.001+08:00

http://www.webdevelopersjournal.com/articles/jsintro1/js_begin1.html

qt opensource windows visual studio expression microsoft

2007-07-08T18:39:00.001+08:00

http://www.linuxext.com/bbs/simple/index.php?t12127.html

Using visual studio express edition with the microsoft platform SDK

visual studio express c++ download (zh-tw)

宮崎日式涮涮鍋

2007-07-02T22:57:00.001+08:00

宮崎日式涮涮鍋
壹蘋果美食王的介紹

Ubuntu With Acrobat Reader (Traditional Chinese)

2007-07-02T16:28:00.000+08:00

1. Download this

2. tar xfvz this file

3. Run install script: sudo ./INSTALL

4. Conflicting with input method


sudo vim /usr/local/Adobe/Acrobat7.0/bin/acroread
GTK_IM_MODULE=xim

Delphi Game Programming With DelphiX

2007-06-30T16:18:00.001+08:00

http://www.savagesoftware.com.au/delphi/articles/article1.html
http://dn.codegear.com/article/26642

Install Oracle Client For Ubuntu/Debian

2007-06-28T17:38:00.000+08:00

1. Modify/etc/apt/sources.list with adding the following lines:
----
# oracle
deb http://oss.oracle.com/debian unstable main non-free
----
2. Update the repository with "sudo apt-get update"
3. Install the following software from the oracle repository:
-libaio
-oracle xe client
4. Update environmental variables in /etc/bash.bashrc << oracle_env.sh
5. Create tnsnames.ora and export TNS_ADMIN=/file/path/to/tnsnames.ora

Strip HTML Tags From Text

2007-06-25T13:06:00.001+08:00

use HTML::Strip
http://cpan.uwinnipeg.ca/htdocs/HTML-Strip/HTML/Strip.html

Convert Unicode String to XML Character Reference

2007-06-25T13:05:00.001+08:00

Python way

http://www.rhymingpanda.com/weblog/2005/10/20/21_02_33/index.html