Apan’s Notes

Create a NA matrix in R

2011-12-10T01:13:00.001-08:00

matrix(data=NA, nrow=3, ncol=5) can create/initiate a 3 × 5 NA matrix on the fly. Of course, matrix(data=0, nrow=3, ncol=5) may be useful in case.

Statistical computing by using R

2011-10-30T16:56:00.000-07:00

最近把東海大學生命科學系大學部的生物統計學實驗講義中的 SAS code 全改寫成 R code。內容是如何使用 R 計算常見的統計檢定，包括例題及 R code。該文件以 Creative Commons Attribution-ShareAlike 3.0 Unported License 授權。

View Statistical computing by using R on Scribd

A function of Pearson/Deviance goodness-of-fit for (generalized) linear regression

2011-08-04T05:57:00.000-07:00

After an analysis of (generalized) linear regression, a Pearson/Deviance goodness-of-fit test is useful to test if the model fitted reasonably. I wrote a simple function in R to do it.

Source Code

# Usage example:
#   model.name <- glm (...)
#   lm.fit.test(model.name)

lm.fit.test <- function (my.model) {
  df <- my.model$df.residual;
  chisq.pearson <- sum(resid(my.model , type="pearson")^2);
  chisq.deviance <- sum(resid(my.model , type="deviance")^2);
  p.pearson <- pchisq( chisq.pearson , df , lower.tail=F );
  p.deviance <- pchisq( chisq.deviance , df , lower.tail=F );
  ratio.pearson <- chisq.pearson / df;
  ratio.deviance <- chisq.deviance / df;
  cat(
    "Pearson chisq = " , chisq.pearson ,
    ", df = " , df ,
    ", p = " , p.pearson ,
    ", chisq / df = " , chisq.pearson/df ,
    ".\n" ,
    "Deviance chisq = " , chisq.deviance ,
    ", df = " , df ,
    ", p = " , p.deviance ,
    ", chisq / df = " , chisq.deviance/df ,
    ".\n"
  );
}

Usage and example

After defining a glm() model as a custom name, copy the source code into R interpreter and then input lm.fit.test(your.model.name) where your.model.name is the custom glm model name. See the fallowing example

y <- c(rpois(20,1) , rpois(20,2));
x <- gl(2 , 20);
my.model <- glm(y~x , family=poisson);
summary(my.model);
lm.fit.test(my.model);

Games-Howell post-hoc in R

2011-07-07T08:56:00.001-07:00

Mr. Shigenobu AOKI displayed a quick R function for both of Tukey and Games-Howell post-hoc comparisons after an ANOVA. The Games-Howell post-hoc is famous for non-equal variances between treatments. Read more on the webpage about Games-Howell post-hoc in R.

邏輯斯迴歸、勝算與勝算比的關係 (logistic regression, odds and odds ratio)

2011-06-22T23:21:00.000-07:00

公式一覽

\[ \log ( \text{odd} ) = \text{logit} (p) = \log \left( \frac{\Pr (y=1) }{1 - \Pr (y=1)} \right) = \beta_0 + \sum_{i=1}^{n} \beta_i x_i \] \[ \text{odd} = \frac{\Pr (y=1) }{1 - \Pr (y=1)} = \exp \left( \beta_0 + \sum_{i=1}^{n} \beta_i x_i \right) \] \[ \Pr (y=1) = \frac{\text{odd}}{1+\text{odd}} = \frac{1}{1+\text{odd}^{-1}} = \frac{1}{1+\exp (-( \beta_0 + \sum_{i=1}^{n}\beta_i x_i)) } \] \[ \log (\text{odds ratio}) = \log \left( \frac{\text{odd}_1}{\text{odd}_2} \right) = \log(\text{odd}_1) - \log(\text{odd}_2) \]

寫在前面

如何解讀邏輯斯迴歸 (logistic regression) 的迴歸式? 因為我常搞錯邏輯斯迴歸結果與勝算或勝算比的關係, 特書此筆記. 本文最主要的內容是唸過 “FAQ: How do I interpret odds ratios in logistic regression?” 一文後的筆記。

名詞翻譯

logistic regression: 邏級式迴歸;
odds: 勝算;
odds ratio: 勝算比.

請注意 “勝算” 與 “勝算比” 是兩回事.

僅包含常數項的迴歸式

在模型 \[ \text{logit}(p) = \log(\frac{1}{1-p}) = \beta_0,\] 其中 $p$ 為依變數 $y = 1$ 之機率, 而 $y$ 為布林資料 (0 或 1) 以表示肯定或否定. 此時, $\beta_0$ 即為 \[ \log(\frac{\Pr (y = 1)}{\Pr (y = 0)}).\] 說成白話文, $y = 1$ 的機率除以 $y = 0$ 的機率為 $\exp(\beta_0)$.

我以下列資料舉例. 抽樣的 10 名學生中, 有 6 人有手機而 4 人沒有. 計算 “有手機” 對 “無手機” 的勝算為 $\frac{6/10}{4/10} = 1.5$. 在迴歸式中可得 $\text{logit}(p) = \log( \frac{1}{1-p}) = 0.4055$, 可驗證 $\exp(0.4055) = 1.500052$ (因為 $y = 1$ 為有手機而 $y = 0$ 為無手機). 以下為驗證之 R code:

> isPhone<-c(rep(1,6),rep(0,4))
> m<-glm(isPhone~1,family=binomial);summary(m)

Call:
glm(formula = isPhone ~ 1, family = binomial)

Deviance Residuals: 
   Min      1Q  Median      3Q     Max  
-1.354  -1.354   1.011   1.011   1.011  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)
(Intercept)   0.4055     0.6455   0.628     0.53

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 13.46  on 9  degrees of freedom
Residual deviance: 13.46  on 9  degrees of freedom
AIC: 15.46

Number of Fisher Scoring iterations: 4

> exp(0.4055)
[1] 1.500052

包含一個二元類別自變數的迴歸式

在模型 \[ \text{logit}(p) = \beta_0 + \beta_1 x_1, \] 其中 $p$ 為依變數 $y = 1$ 之機率而 $y$ 為布林資料 (0 或 1) 以表示肯定或否定, $x_1$ 為二元類別資料並轉換為布林資料 (0 或 1). 此時, $\beta_0$ 即為 $x = 0$ 時 $\log ( \frac{\Pr (y = 1)}{\Pr (y = 0)})$, 而 $\beta_1$ 為 $x_1 = 1$ 對 $x_1 = 0$ 在 $y$ 是否為 1 的勝算比. 說成白話文, $x_1 = 0$ 時 $y = 1$ 的機率除以 $y = 0$ 的機率為 $\exp( \beta_0)$; $x_1 = 1$ 對 $x_1 = 0$ 的勝算比為 $\beta_1$.

我以下列資料舉例. 抽樣的 10 名學生中, 有 6 名男生 4 名女生, 且每人分別擁有手機或沒有, 資料如下:

有無手機	男生 (x₁ = 1)	女生 (x₁ = 0)
有 (y = 1)	5	1
無 (y = 0)	1	3
小計	6	4

“男生有手機” 對 “男生無手機” 的勝算為 $\frac{5/6}{1/6} = 5$; “女生有手機” 對 “女生無手機” 的勝算為 $\frac{1/4}{3/4} = 0.3333$; 二者的勝算比為 $5 / 0.3333 = 15$.

在男生 = 1, 女生 = 0, 有手機 = 1, 無手機 = 0 的設定下, 迴歸式中可得 $\text{logit}(p) = -1.099 + 2.708(\text{性別})$, 可驗證 $\exp (-1.099) = 0.3332041$ 為 “女生有手機” 對 “女生無手機” 的勝算 (因為女生的編碼為 $y = 0$); $\exp (2.708) = 14.99925$ 為性別與有無手機的勝算比, 可解釋作 “男生有手機之勝算為女生的 15 倍”. 以下為驗證之 R code:

> isMale<-c(1,1,1,1,1,1,0,0,0,0)
> isPhone<-as.factor(c(1,1,1,1,1,0,0,0,0,1))
> m<-glm(isPhone~isMale,family=binomial);summary(m)

Call:
glm(formula = isMale ~ isPhone, family = binomial)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.8930  -0.7585   0.6039   0.6039   1.6651  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)   -1.099      1.155  -0.951   0.3414  
isMale         2.708      1.592   1.701   0.0889 .
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 13.4602  on 9  degrees of freedom
Residual deviance:  9.9054  on 8  degrees of freedom
AIC: 13.905

Number of Fisher Scoring iterations: 4

> exp(-1.099)
[1] 0.3332041
> exp(2.708)
[1] 14.99925

包含一連續變數作為自變數的迴歸式

在模型 \[\text{logit}(p) = \beta_0 + \beta_0 x_1,\] 其中 $p$ 為依變數 $y = 1$ 之機率而 $y$ 為布林資料 (0 或 1) 以表示肯定或否定, $x_1$ 為連續型變數. 此時, $\beta_0$ 即為 $x = 0$ 時 $\log (y = 1 \text{對} y = 0 \text{的勝算})$, 而 $\beta_1$ 為 $x_1$ 對勝算的貢獻量. 說成白話文, $x_1=0$ 時 $y = 1$ 的機率除以 $y = 0$ 的機率為 $\exp(\beta_0)$; 若 $x_1 \neq 0$ 時 $y = 1$ 的機率除以 $y = 0$ 的機率為 $\exp (\beta_0 + \beta_1 x_1)$, 而 $\beta_1$ 為正數或負數則影響了勝算的正向或負向影響。

我以下列資料舉例. 抽樣的 10 名學生中, 每人有不同每日零用錢及是否擁有手機, 資料如下:

有無手機 (y, 1 為有, 0 為無)	每日零用錢 (x₁)
1	100
1	120
1	90
1	110
1	110
0	80
0	100
0	90
0	70
0	110

上述資料可得回歸式 $\text{logit}(p) = -9.71278 + 0.09840(\text{每日零用錢})$. 由迴歸式可觀察出, 每日零用錢每增加 1 單位則對有手機之勝算增加 $\text{exp}(0.09840) = 1.103404$, 亦表示有無手機在多拿 1 單位零用錢與沒有者的勝算比是 $\exp(\beta_1 \times 1) = 1.103404$，多拿 2 單位則為 $\exp(\beta_1 \times 2)$, 以此類推.

我們可以估計每日零用錢為 70 時有手機對無手機的勝算: $\exp (-9.71278 + 0.09840 \times 70) = 0.05932171$, 可解釋為每日零用錢為 70 的學生中有手機的機率為沒有手機的 0.05932171 倍. 要估算每日零用錢為 75 的勝算亦如法炮製: $\exp( -9.71278 + 0.09840 \times 75) = 0.09702564$. 因此, 二者的勝算比為 $0.09702564 / 0.05932171 = 1.635584$, 也就是 $\exp (0.09840 \times 5)$.

由於$\text{勝算} = p / (1 - p)$, 故可由勝算回推 $p$, 公式為 $p = \text{勝算} / (1 + \text{勝算})$. 前面有談到零用錢為 70 時有手機對無手機的勝算為 0.05932171, 可推得零用錢為 70 時有手機的機率為 $0.05932171 / (1 + 0.05932171) = 0.05599971$. 同理, 每日零用錢為 75 時有手機的機率為 $0.09702564 / (1 + 0.09702564) = 0.08844428$.

> isPhone<-c(1,1,1,1,1,0,0,0,0,0)
> money<-c(100, 120, 90, 110, 110, 80, 100, 90, 70, 110)
> m<-glm(isPhone~money,family=binomial);summary(m)
Call:
glm(formula = isPhone ~ money, family = binomial)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-1.67078  -0.76661   0.07113   0.75438   1.55604  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)
(Intercept) -9.71278    6.59477  -1.473    0.141
money        0.09840    0.06574   1.497    0.134

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 13.863  on 9  degrees of freedom
Residual deviance: 10.481  on 8  degrees of freedom
AIC: 14.481

Number of Fisher Scoring iterations: 4

結論

邏輯斯迴歸的迴歸式 $\log ( \text{odd} ) = \text{logit} (p) = \log \left( \frac{\Pr (y=1) }{1 - \Pr (y=1)} \right) = \beta_0 + \sum_{i=1}^{n} \beta_i x_i$ 所加總得到的數值是 log(勝算), 故 exp(logit(p)) = 勝算;
將二個勝算相除可以得到勝算比;
可由勝算回推機率, 公式為機率 = 勝算 / (1 + 勝算);
因為指數與對數在運算上的特性, log(勝算比) 常成為二個 logit(p) 相加減之結果.
若不解, 一切回歸到 “log(勝算) = log (p/(1 - p)) = 迴歸式” 及 “勝算比 = 勝算₁ / 勝算₂” 的基本概念, 再往下推導.

Least-significant difference (LSD) in R

2011-06-20T02:48:00.000-07:00

How to apply LSD (as a kind of multiple comparisons after ANOVA) in R?

A dataset for typical one-way ANOVA is given:

y <- rnorm(40)
x <- gl(4,10)
m <- aov(y~x)

Applying LSD by using function pairwise.t.test:

pairwise.t.test(y,x,p.adj="none")

Applying LSD by using function LSD.test from package agricolae:

library(agricolae)
LSD.test(m,"x",group=T) # with grouping
LSD.test(m,"x",group=F) # without grouping

Specifying a reference level in R: function relevel

2011-06-15T04:12:00.000-07:00

How to specify a specific level as a reference level in GNU R while applying (generalized) linear regression?

A dataframe dt is given:

y<-rpois(40,1)
x1<-gl(4,10)
x2<-gl(4,1,40)
dt <- data.frame(cbind(y,x1,x2))

Applying poisson regression:

m1<-glm(y~factor(x1)+factor(x2),data=dt)
summary(m1)

The output shows x1=1 and x2=1 are both reference level defined by R automatically. To switch the reference level from x1=1 to x1=2, the function relevel could be used:

m2<-glm(y~relevel(factor(dt$x1),"2")+factor(x2),data=dt)
summary(m2)

阿盤的 Beamer 中英字型配方

2011-05-20T01:23:00.000-07:00

本文介紹我自己的 Beamer 中英文字型配置。

詳文請閱阿盤的 Beamer 中英字型配方。另歡迎至討論串伺教。

學習 HTML, PHP 及 Database

2010-06-25T15:38:00.000-07:00

在 telnet://ptt.cc 的 Web_Design 板有板友剛入門寫網頁. 他有為自己訂立了學習目標, 包括學習架站, 寫討論區等等, 但又不知道如何入門, 也不知道學會了這些技術有什麼幫助, 於是發表標題為 “[問題] 新手上路，卻不知路在何處...” 的文章. 有些板友給予他回應, 大多認為他想學的東西太基礎, 早已經有許多現成的技術可以取代, 例如 “手寫討論區? 何不用 XOOPS?” 等等. 我看了還蠻多感觸的, 於是寫了篇文章回應, 內容如下.

我願意分享一些自己的心得給原 PO.

我不是本科生, 也完全沒有相關背景, 頂多小時候學過 BASIC 語言罷了. 上大學之後, 為了當時流行每個人有自己的個人網頁, 就開始入門. 一開始, 當然也只是用 WYSIWYG 的軟體來製作網頁. 不知不覺, 就對 HTML 有了興趣. 畢竟這才是做網頁的基礎.

HTML 不只是做網頁罷了, 其內涵是讓你使用正確的 tag 來展現文稿結構. 當時看了一些書, 主要是教導大家屏棄舊式瀏覽器支援的 HTML 語法 (如 <font>), 而採用 HTML 4 + CSS. 學這些沒為了什麼, 就是讓自己的網頁在各種瀏覽器都可以正確解析罷了. (於是加入了 “凡是網頁設計師必痛恨 IE 6” 的族群.)

我就利用自己學會的 “標準 HTML + CSS” 開始寫自己的日記. 做到這裡, 就很有成就感了: 每天寫寫內容, 畫畫 CSS.

好日子過沒多久, 發現我的日記越來越多. 由於我是以 “一篇日記就是一個 HTML file” 的方式進行. 要讓每篇日記都共享相同的網頁結構變得非常困難. 看來, 是該把我的日記做得像 blog 才行. 於是我就開始學習 PHP.

當然, PHP 必竟是一種程式語言, 學習上會吃力一些. 於是我就從最簡單的 PHP 教本開始唸. 只唸我有需要的部份. 有時候簡單的教本還沒唸完, 就又買了進階級的教本來唸. 就這些唸了又練, 練了又唸, 花了不少時間, 終於把網頁以 PHP 改寫完.

沒多久, 我才認知, 原來我還需要學會 database 才能幫助我的日記網站有更好的維護方式. 那時我連 database 是幹嘛的都不懂, 只覺得大概就像個 excel 檔吧! 好在 PHP 的教本都是連帶著 MySQL 一起教, 也就學了 database 了.

從學 HTML 到 PHP + MySQL, 也花課餘時間二三年了呢. 雖然程式不到專業級, 但給自己的成就感是很豐富的了. 有了一些些實做的經驗, javascript, python, perl, 正規表示式也都學了點皮毛. 我幫老師寫過實驗室網站, 也利用這些知識幫助我目前的工作. 但我明白我是不可能用這些基本功養家的.

也許原 PO 會覺得, 花了這些多時間學會這些不能拿來養家的技術, 何苦? 我倒完全不這麼覺得啦.

我學會了 HTML, 也就學會了 “文稿結構” 這件事. 這對寫任何文章都有幫助, 也使我學習 LATEX 非常快.

為了學 CSS, 就探究了平面設計與繪圖軟體操作. 這在許多工作也都用得著.

學會了 PHP, 就等於你是個 “會寫程式的人” (即使只是半調子.) 會寫點小程式對做很多事都很有幫助, 例如別人要用手按一星期的工作, 而你只要半天就完成. 此外, 只要學會第一種程式語言, 第二種就容易許多了.

學了 database, 我連計帳本都做正規化了! 有時候老闆交代幫他整理資料或實驗數據, 懂 database 概念的人真的得心應手許多.

好吧, 也許我是在催眠自己. 但我沒後悔過花這麼多時間學過這些技術. 它們雖然沒直接幫助我賺什麼錢, 但附加價值是無窮的, 也大大地豐富我的精神生活. 不過, 假如現在我再重新寫一個中型網站專案, 我沒空也懶惰就是了!

MacTeX + CJKutf8 + DVIPDFMx 的新增中文 TTF 字形

2008-07-27T08:25:00.000-07:00

說明 MacTeX + CJKutf8 + DVIPDFMx 的環境下，新增中文 TrueType 字形的詳細步驟。

本文可能已經過時，請斟酌參考。

DVIPDFMx 的優點

在利用 CJK package 進行中文文稿由 DVI 編譯為 PDF 的時候，若採用 DVIPDFMx 就可以直接利用 TrueType 字型檔而不用使用 Type 1 字型檔。除此之外，DVIPDFMx 也是藉由 Ghostscript 的管路來產生 PDF，所以 Ghostscript 的參數也是可以在編譯時一併使用，功能仍不輸給較傳統的 latex + dvips + ps2pdf 方法。另外，利用自行撰寫的 shell script 也可以使 latex + dvipdfmx 的過程在 TeXShop 中方便運行。

運行環境及所需軟體

OS: MAC OS X 4 Tiger
TeX system: MacTeX-2007 Package
原始中文字型: TrueType 之 TTF 或 TTC 字型檔
（選用）純文字編輯器: Smultron
（選用）造字改字應用程式: FontForge

開工

以下我們以 cwming.ttf 為例。

利用 `ttf2tfm` 產生 TFM 檔

一個中文字型檔必定會產生大量的 TFM 檔。我們先把 cwming.ttf 放在桌面上 cwming 資料夾中，再到 Terminal 中輸入

cd ~desktop/cwming
ttf2tfm cwming.ttf cwming@Unicode@

以到達工作路徑並以生成 TFM 檔。經過等待，你就會在該資料夾中得到一大堆 TFM 檔。其中，cwming@Unicode@ 表示生成的 TFM 是以 cwming 做檔名頭，而以 Unicode 編碼運行。

若對於 ttf2tfm 想了解更多，請在 Terminal 中輸入 man ttf2tfm 或 ttf2tfm --help 以取得其說明文件。

如果你的原中文字型檔是 TTC，請了解一下 ttf2tfm -f 的用法。

如果你要動用到 FontForge 先修改字型，那也可以在 FontForge 中在輸出新 TTF 檔時一併輸出 TFM 檔。

把 TFM 搬到正確的位置

在 /usr/local/texlive/2007/texmf-dist/fonts/tfm/ 中先建一個 cwming 資料夾，再把方才所有的 TFM 檔搬到這裡，這都需要管理者權限。在 Terminal 中輸入

cd /usr/local/texlive/2007/texmf-dist/fonts/tfm
sudo mkdir cwming
sudo mv ~/desktop/cwming/*.tfm /usr/local/texlive/2007/texmf-dist/fonts/tfm/cwming/

即可。

把 TTF/TTC 搬到正確的位置

在 /usr/local/texlive/2007/texmf-dist/fonts/truetype/ 中先建一個 cwming 資料夾，再把方才的 TTF 或 TFM 檔搬到這裡，這都需要管理者權限。在 Terminal 中輸入

cd /usr/local/texlive/2007/texmf-dist/fonts/truetype
sudo mkdir cwming
sudo mv ~/desktop/cwming/cwming.ttf /usr/local/texlive/2007/texmf-dist/fonts/truetype/cwming/

即可。

修改 `cid-x.map`

利用純文字編輯器開啟 /usr/local/texlive/2007/texmf/fonts/map/dvipdfm/dvipdfmx/cid-x.map，對在新的一列中加入
cwming@Unicode@ unicode :0:cwming.ttf
即可。若需說明，請仔細看看 cid-x.map。

製作 FD 描

利用純文字編輯器製作 c70cwming.fd，其中 c70 是因為你用 Unicode。c70cwming.fd 的內容就定義了 LaTeX 指令如何控制這個新字型，例如

\ProvidesFile{c70cwming.fd}
\DeclareFontFamily{C70}{cwming}{\hyphenchar \font\m@ne}
\DeclareFontShape{C70}{cwming}{m}{n}{<-> CJK * cwming}{}
\DeclareFontShape{C70}{cwheib}{b}{n}{<-> CJKb * cwheib}{}
\DeclareFontShape{C70}{cwheib}{bx}{n}{<-> CJKb * cwheib}{}
\DeclareFontShape{C70}{cwming}{m}{it}{<-> CJK * cwkai}{}
\DeclareFontShape{C70}{cwming}{bx}{it}{<-> CJKb * cwkai}{}

就是定義了以 cwming 為 CJK 字型時的 Normal、Bold、Italic 與 Bold Italic 是什麼字型了。當然，其它的字型你也要一一裝上去才有的用。

之後把 c70cwming.fd 放到 /usr/local/texlive/2007/texmf-dist/tex/latex/CJK/UTF8/ 中即可。

更新 TeX 索引

在 Terminal 中輸入 sudo texhash 。以更新 TeX 索引。這樣就完成了。現在已經可以使用 latex filename;dvipdfmx filename 了。

TEX 檔範例

\documentclass[10pt,a4paper]{article}
\usepackage[T1]{fontenc}
\usepackage{CJKutf8}
\begin{document}
\begin{CJK}{UTF8}{cwming}
\CJKtilde
測試看看。
\clearpage
\end{CJK}
\end{document}

與 TeXShop 的配合

請見陳念波老師所寫的「如何在 Mac OS X 上裝設 LaTeX 以及 CJK 套件 (TeX Live 2007)」中的「在 TeXShop 裡方便的使用 latex + dvipdfmx 工作流程」。

與 Ghostscript 的配合

請見 /usr/local/texlive/2007/texmf-config/dvipdfm/config/dvipdfmx.cfg 中的 D "gs..." 及 Ghostscript 官網上 ps2pdf 的說明文件。簡單來說，就是把 ps2pdf 中的參數寫在 D "gs..." 裡。

CJK 全型標點符號「對齊」行尾

2008-02-23T05:57:00.000-08:00

在 LATEX 中對中文全型標點符號重新給予新命令，並達到「避頭點」、「貼近文字基線」及「可對齊行尾」之目的。

先前的做法

過去在網路上有些人為了使 CJK 中的中文全型標點符號更美觀，以自定義新命令來取代標點符號，如下例：

\def\pe{%
  \nolinebreak[4]% A
  \hspace{-0.25em}% B
  \raise -0.75em\hbox{。}% C
  \hspace{0.4em}% D
}

其中 A 列表示句點前不可斷行；B 列使文字往左移；C 列使文字下降後輸出句點；而 D 列代表文字後之距離。之後再將文中所有之句點以 \pe 來取代，即可得到貼近基線及左側並可避頭點之中文全型句點，非常美觀。

很抱歉我並不知道這個方法是誰先提出來的，在此向原創者至上感謝之意。

先前的做法有些問題

不過這個方法似乎有些問題，而且仍不及吳老師等人之 cwTEX 的效果。以下列出幾點。

由於使用了 \nolinebreak[4] 達到避頭點，但其效力仍貫穿了整個 \pe 命令，故使句號無法在行末出現。
中文的全型句點基乎都沒有所謂的 kerning 設定，故即使某行之行未為句點，且 LATEX 也貼心地做到左右對齊，但句點的右邊仍留下了原本的小小空間。

修正後的辦法

所以我想了一些方法，可以修正這二項問題，一樣也是以中文全型句點為例：

\def\pe{%
  {%
    \nolinebreak[4]% A
    \hspace*{-0.25em}% B
    \raise-0.75ex\hbox{。}% C
    \hspace*{-0.25em}% D
  }%
  \allowbreak% E
  \hspace{0.5em plus 0.1em}% F
  \ignorespaces% G
}%

其中與以往不同之處在下列表說明。

我特別為句點之後加了「必須執行」的負水平伸縮距離（B 與 D 列），使句點的左右兩邊空白都被迫減少。這二個伸縮距離就是該中文句點「字」本身左右所占據的空白水平長度。如此一來，B、C、D 三列可以使你得到一個「沒有左右空白」的中文句點。「必須執行」的用意在於就算行末是個句點，也必須執行 \hspace*。若僅使用 \hspace ，在行末的句點輸出後，LATEX 不理會 \hspace。
在輸出句號（C 列）及去除句號本身左右空白（B 與 D列）之前，為了避頭點，一樣不可以斷行（A 列），但這四列已被群組，表示不可斷行之功能僅限於這四列。這樣就可以達成句點後斷行。
在不可斷行之群組功能失效後，即容許斷行（但並非要求斷行，E 列）。這只是鼓勵 LATEX 可以在此斷行。
F 列的「不必須執行」水平拉伸是指「非行末」之句點右方之空白空間。基於 LATEX 的精神，句子結束又不斷行的情況下，水平空白應該比較大，所以在 F 列我給予了較大的空白。
至於 G 列，只是個習慣，必不必要我就不清楚了。

新辦法的模型

所以我把我的想法寫成一個模型，來表示類似句點這類必須避頭點之自定命令（如逗號及問號等，但左括號或左引號則不適用）：

\def\yourcmd{%
  {%
    \nolinebreak[4]%
    \hspace*{length1}%
    \raise length2 \hbox{×}%
    \hspace*{length3}%
  }%
  \allowbreak%
  \hspace{length4}%
  \ignorespaces%
}%

其中，

yourcmd 為取代後的 LATEX 命令；
length1 為該標點符號「文字本身」左方空白長度，應為負值，不過若你覺得文字與之後的標點符號太近了，可以自行修正一下；
length2 為該標點符號往下降的長度，應為負值；
× 為原本的中文全型標點符號；
length3 為該標點符號「文字本身」右方空白長度，應為負值；
length4 為該標點符號若並非出現於頁尾的情況下，後之後文字的距離。

尚未達成的目標

這項方法的缺點，是真的挺煩人的。要先抓到不同字體中每個標點符號的數值，還要記得在文中以新命令取代原本的符號。此外，吳老師等人也提到，像『「謝謝你。」』中有二個以上標點符號前後連續出現時，我的方法也還是會有失誤。

此外，我也正在想「避尾點」的標點符號的模式。

我練習 LATEX 並不久，文中若有任何錯誤，望大家多多指教。

替 MacTeX 新增 CJK 之 Type 1 字型

2007-12-27T09:30:00.000-08:00

本文可能已經過時，請斟酌參考。

本文說明如何在 MacTeX-2007 下為 CJK 新增中文 Type 1 字型

前言

首先感謝李果正先生的資料，其中對於新手的指點非常充足。

以 LaTeX 系統排中文，最有名的方式就是 CJK 了。不過，在 MacTeX-2007 中，只提供了 bsmi（為一種明體）作為正體中文之 CJK 字體，確實不太夠用。而網路上大部份說明如何新增 CJK 字型的中文資料，針對 MacTeX 者少。

我用 xelatex 也好一陣子了。不過 xelatex 最大的缺點就是，中英文字型會用同一種字型，除非每次遇到中英文時分別使用更換字型的命令（也許真的有更好的辦法，但小弟並不知道）。

以下我說明替 MacTeX-2007 新增 Arphic bkai00mp-20071106.tar.bz2 為 CJK 字型的步驟。若要從 TTF 字型轉成 TYPE 1 字型開始，我也不會，請參考上述李果正先生的資料。

工作環境

我的條件是在 OS 10.4 下灌好 MacTeX-2007，並未其它任何其它 TeX 或ghostscript 系統。

進行步驟

下載 bkai00mp-20071106.tar.bz2 並在 OS X 中直接解壓縮。
解壓縮後，將相關的資料夾或檔案丟到正確位置。過程中一定會要求密碼驗證。不過太擔心丟錯地方會怎樣，其實新字型資料夾的結構和要丟去的地方有相似之處。如果你在 Finder 中找不到怎麼到這些目的地，請試著用 apple + shift + g。
- .../fonts/afm/arphic/bkaiu
  這個資料夾丟到
  /usr/local/texlive/2007/texmf-dist/fonts/afm/arphic/
  裡。
- .../fonts/map/arphic/dvips/arphic/bkaiu.map
  這個檔丟到
  /usr/local/texlive/2007/texmf-dist/fonts/map/dvips/arphic/
  裡。
- .../fonts/tfm/arphic/
  內的
  bkaimp 與 bkaiu
  二個資料夾丟到
  /usr/local/texlive/2007/texmf-dist/fonts/tfm/arphic/
  裡。
- .../fonts/type1/arphic/bkaiu
  這個資料夾丟到
  /usr/local/texlive/2007/texmf-dist/fonts/type1/arphic/
  裡。
- .../fonts/vf/arphic/bkaimp
  這個資料夾丟到
  /usr/local/texlive/2007/texmf-dist/fonts/vf/arphic/
  裡。
確定
/usr/local/texlive/2007/texmf-dist/tex/latex/CJK/UTF8/
中有沒有
c70bkai.fd 及 c70bkai.fdx
這二個檔。如果沒有，到 CTAN 裡來抓。
在 terminal 中，執行：
sudo texhash
並輸入管理者密碼。此動作可以更新
/usr/local/texlive/2007/**/
下數個 Is-R 檔。
在 terminal 中，執行：
sudo updmap-sys --enable Map bkaiu.map
之後可以在
/usr/local/texlive/2007/texmf-config/web2c/updmap.cfg
這個檔案中找到
Map bkaiu.map
這行敘述。
（此步驟不知道是否必要）在
/usr/local/texlive/2007/texmf-dist/dvips/arphic/config.arphic
中加入
p +bkaiu.map
一行。
到此步驟為止，已經可以用 latex + dvipdfmx 成功編譯，不過檔案很大，且 latex + dvips + ps2pdf 在 dvips 時仍會找不到新字型。我搞了很久，終於在
/usr/local/texlive/2007/texmf-config/dvips/config/config.ps
（記得用純文字編輯器打開）中，加入
p +bkaiu.map
一行即可。之後，latex + dvips + ps2pdf 可行，且 PDF 檔案也小。（但我不知道為什麼 MacTeX-2007 內建的 bsmiu.map 與 gbsnu.map 就不用。）

這樣就完成了。我以
\usepackage{CJKutf8}......\begin{CJK}{UTF8}{bkai}......\end{CJK}
再加上 latex + dvips + ps2pdf 可以完成很不錯且體積小的 PDF 檔，且與 hyperref 這個 package 在 PDF 中也配合得很好。

後記

這些過程我試很多次了，所以前前後後也許有遺漏之處，盼前輩能夠指正。希望這篇說明對其它 MAC 愛好者能有所小小幫助。

若要加入其它中文字型給 MacTeX-2007 內的 CJK 用，我想大同小異吧。反正改錯了再改回原來的樣子就好了，勇敢試下去吧！

以 MATLAB 轉換 OOIBase32 原始檔

2007-09-22T01:21:00.000-07:00

本文分享一個 MATLAB 之 m 函數原始檔，可將 OOIBase32 來源之原始文字檔轉換成二個經過三次平滑之一維陣列，並生成一個純文字（.dat）之結果檔案。

OOIBase32 所產生的原始檔，通常長這個樣子。

OOIBase32 Version 2.0.6.5 Data File
++++++++++++++++++++++++++++++++++++
Date: 01-24-2007, 17:39:35
User: Valued Ocean Optics Customer
Spectrometer Serial Number: USB4C02907
Spectrometer Channel: Master
Integration Time (msec): 35
Spectra Averaged: 5
Boxcar Smoothing: 0
Correct for Electrical Dark: Disabled
Time Normalized: Disabled
Dual-beam Reference: Disabled
Reference Channel: Master
Temperature: Not acquired
Spectrometer Type: S4000
ADC Type: USB4000
Number of Pixels in File: 3648
Graph Title: 
>>>>>>Begin Spectral Data<>>>>
178.59 0.000
178.81 0.000
179.02 88.000
179.24 -2350.000
179.46 -14.286
179.67 34.043
179.89 60.000
180.10 -3600.000
180.32 100.000
180.54 1266.667
.
.
.
887.29 110.120
887.46 110.528
887.63 111.647
887.80 109.825
887.97 111.218
888.13 110.559
888.30 109.965
>>>>>>End Spectral Data<>>>>

以下的函式可以將上述原始檔進行轉換（假如已經更改過 OOIBase32 原始檔，那我不保証這個函式能正確運作）。請複製以下原始碼，並儲存為 trans_data.m。

function [trans_x,trans_y] = trans_data(source_file,extension,begin,inter,endd)
% Load the source_file.extension,
% then smooth the data and generate a 2*4001 array,
% and creat a source_file.dat.
% out put: trans_x and trans_y: the transformed and smoothe data

%load the source file
source_file_ext = [source_file '.' extension];
fid = fopen(source_file_ext, 'r');

% sweap several beginning lines
for n = 1:30
 s1=fgetl(fid);
end

%scan wavelength
mydata=fscanf(fid,'%g',[2 inf]);
fclose(fid);
xx=mydata(1,:);
yy=mydata(2,:);
trans_x=begin:inter:endd;
trans_y=interp1(xx,yy,trans_x,'cubit');

%write source_file.dat
fid=fopen([source_file '.dat'],'w');
fprintf(fid,'%9.9f %9.9f\n',[trans_x;trans_y]);
fclose(fid);

trans_data 函式共有 5 個引數，分別為來源檔之檔名、來源檔之副檔名、起始值、平滑間距及終止值。舉例來說，我取得一個由 OOIBase32 得到的來源檔，名叫 color.Master.transmission ，其中包括許多 x 與 y 之對應值。我希望將該檔內的資料由 x=300 取至 x=700，且平滑間距為 0.1，轉換成新的 x 叫 new_x 且新的 y 叫 new_y。方法如下：將 trans_data.m 與 color.Master.transmission 置於同一個資料夾，再從 MATLAB 中到達該資料夾，並輸入：

[new_x,new_y]=trans_data('color','Master.transmission',300,0.1,700)

即可，而且也會在該資料夾中生成一個 color.dat 檔。由於其中 color 與 Master.transmission 為字串，故需以單引號括起來。若想觀看結果，可在 MATLAB 中輸入：

plot (new_x,new_y)

以觀看新產生的結果。

以上例為例，trans_data 函式的輸出有二個變數（new_x與new_y）及一個輸出文字檔（color.dat）。你可以根據你的需要，利用這二種不同類型的資料再引入其它函式中操作。

請注意。由於 trans_data 函式會將 OOIBase32 原始檔前 30 列去除（個人的經驗中前 30 列包括不需要的文字說明及無意義資料），若有需要，請修改將函式中 for n = 1:30 的 30 修改成其它數值。

如何產生 PDF 文件？

2007-09-01T14:49:00.001-07:00

本文說明 PDF 文件之重要性，以及如何產生 PDF 文件。

PDF 文件與其優點

所謂的 PDF 文件，基礎上來說，就是給「印表機」看的文件。只要有相同的 PDF 文件，在印表機看來，都是完全相同的列印內容。也因此，PDF 可以不限平台、不限作業系統、不限原生軟體之下，都有相同的列印格式。

PDF 有許多好處。最大的好處就是可以讓所有的人都毫無偏差地看見同一份文件。過去我們習慣使用 MS office word 所產生的 .doc 檔來分享文件，但在不同的版本的 office 所產生的 .doc，在另一個版本的 office 所開啟之後，格式都會跑掉。而且，誰規定每個人都必須購買 MS office 呢？沒有購買 MS office 的人就沒有權利接受別人生產的文件嗎？在 PDF 的世界，因為 PDF 是給印表機看的檔案，所以不會有這樣的問題。另外，通常 PDF 是不可以再更改的。而且，PDF 也可嵌入字型，以確保文字正確顯示。這些都是 PDF 足以做為分享分件的最大優點。

如何產生 PDF 文件

一般在 Windows 系統上，若使用 MS office 的話，並不能直接生產 PDF 文件。但是解決辦法有許多種：

購買 Adobe 公司的 Acrobat Professional。很貴。
使用其它的印表機外掛程式。我引述「戴嘉賢 (Nelson Tai)」先生所整理的免費軟体。請大家選用。該網頁中所介紹的「CutePDF Writer」、「PrimoPDF」等免費軟體可能是你所需要的。當然也可以自行 google 這些軟體再自行下載。
如果你很熟電腦，阿盤最推薦的辦法是，不要使用 MS office 而改用開放原始碼的 OpenOffice。OpenOffice 可以直接生成 PDF。

如果是使用 Unix-like 作業系統，既然你已經玩 Unix-like，那產生 PDF 應該是小 Case 囉！至於 Mac 使用者，內建有產生 PDF 的程式。

XHTML 鑲嵌 QuickTime 媒體

2007-08-31T08:24:00.000-07:00

本文討論了三種在 XHTML 中鑲嵌 QuickTime 媒體的有效語法。

引言

在 XHTML 網頁中鑲嵌影音媒體一直是個備受討論的議題。最主要是由於 <embed> 不是 W3C 標準，加上 IE6 又引入了 ActiveX，使得大家有點不知所措。因為小弟在網頁中常常會鑲嵌 QuickTime 媒體，所以對這個問題做了一些 review，希望對網頁創作者有些幫助。

以我所知，有效的方法有三種：

<object><embed></embed></object> 巢狀結構語法
以 Javascript 插入 <object><embed></embed></object>
<object><object></object></object> 巢狀結構

實際做法

以下將分別討論這三種方法。

1. `<object><embed></embed></object>` 巢狀結構語法

此方法是最常見的方式，但不能通過 W3C 檢測。請見以下的例子。

<object
  classid="clsid:02BF25D5-8C17-4B23-BC80-D3488ABDDC6B"
  width="320"
  height="256"
  codebase="http://www.apple.com/qtactivex/qtplugin.cab#version=6,0,2,0"
  align="middle">
  <param name="src" value="sample.mov" />
  <param name="autoplay" value="true" />
  <embed src="sample.mov" width="320" height="256"
    pluginspage=http://www.apple.com/quicktime/download/
    align="middle" autoplay="true" bgcolor="black" >
  </embed>
</object>

2. 以 Javascript 插入 `<object><embed></embed></object>` 的巢狀結構語法

此方法是 Apple 官方說法。只是將第一個方法改用 Javascript 寫出來，以逃避 W3C 的檢測。但在本質上，是完全相同的。請參考 Preparing Websites with Active Content for Upcoming Changes to Internet Explorer for Windows，其中介紹了二種方法。一種為直接撰寫 Javascript 以取代 <object><embed></embed></object> 巢狀結構語法，另一種則是以 Apple 官方提供的 Javascript Function 來完成。個人比較喜好後者，因為該 Function 很容易使用且語法簡單。無論哪種方式，都可以避開 <embed> 的存在而通過 W3C 檢測。但是一旦用戶端關閉 Javascript，則完全沒有作用。

3. `<object><object></object></object>` 巢狀結構

此方法由 realise.com 提出，請參考 XHTML Quicktime Object 以了解作法。本方法不使用 <embed> 故可通過 W3C 檢測，且在各種新版瀏覽器都有不錯的表現。這項方法有一特色是必須配合 CSS 以使 IE 有正常表現。但是依自己的經驗，若在 <object> 中需要 target 等標籤，仍不能通過 W3C 檢測，因為某些 <object> 內的標籤在 W3C 規範中也不被接受。

結論

如果您不在乎 W3C 標準，那其實什麼方法都不錯。如果您在乎 W3C 標準，請參考第二項及第三項解決方案。雖然第二項解決方案還蠻完美的，但若用戶端關閉 Javascript，則完全失效。所以，真的沒有一項完全 100 分的解決方案就是了！囧。

利用 PowerPoint 製作學術海報

2007-08-13T08:01:00.000-07:00

本文說明如何利用 PowerPoint 製作學術海報，以及大圖輸出時之注意事項。

寫在前面

對於第一次做學術海報的同學來說，千頭萬緒可能不知道怎麼開始才好。大家熟悉的 MS office 的 PowerPoint 可以幫助同學完成看起來還蠻正點的海報。受老闆之託付，阿盤在卓逸民老師的教學網頁中撰寫本文，之後再移植至 Google Page。在接下來的說明中，阿盤將一步步帶領你使用 PowerPoint 來製作海報（以 office PowerPoint 2003 為例）。

如果你不使用 PowerPoint，也可以使用其它的 office 軟體（例如免費的 openoffice）中所提供的投影片軟體或繪圖軟體來製作海報，或是直接利用專業的向量圖製作軟體（ illustrator、coreldraw 等），但難度較高。

阿盤在這個網頁中的一些海報基礎知識，大多是經由東海大學生命科學系卓逸民老師在「生物學研究導論」中的上課內容。在此感謝。

開工了！

第一步：先前工作

先在紙上完成草稿。
在製作數位海報之前，請同學已經在紙上完成海報的草稿（包括文案、版型、引用圖片等），並大致決定海報的配色（例如白底黑字或黑底白字）。請記得，草稿有許多重點要注意，例如「1/8 面積的標題區」、「有意義的大網」、「使用條列式段落」等等。
收集解析度夠大的圖片。
如果放在海報上的點陣圖片解析度不夠，輸出成大海報時可能會出現「馬賽克」的真失現象。但是解析度要多大才行？先讓我們了解何謂「dpi」。所謂的dpi就是指「每英吋有幾個點」。dpi值越高，表示單位長度內的點數越多，也就有越好的印刷品質。通常而言，300 dpi 以上的印刷品可以達到肉眼無法分辨的品質。也就是說，如果你有一張 1200×900 相素的照片，在 300 dpi 的標準下，最大可以印成 4 吋寬 3 吋高（1 吋為 2.54 公分）的印刷品。超過這個輸出尺寸會造成照片上有格子的瑕疵。

第二步：PowerPoint 基本設定

設定海報尺寸。
在台灣學術海報習慣採用 90 公分寬 120 公分高的尺寸。在 PowerPoint 中打開「檔案」→「版面設定」以設定海報尺寸。假如版面設定無法達到這麼大的版面大小，其實採用與輸出大小相同的比例即可，但要注意之後 PowerPoint 中尺規的數字大小要經過換算才是真正的輸出大小。
多多利用尺規。
「右鍵」→「尺規」可以拉出尺規，讓我們可以大致想像海報中每項物件的實際大小。尺規中的對齊標記也很實用，讓每個物件可以對齊。
設定背景顏色。
「右鍵」→「背景」可以達成。值得注意的是，背景顏色不要太豔，以免模糊主題焦點。

第三步：製作海報中的文字

拉出一個文字方塊，將大小做適當調整，並放置在合適的位置。
在文字方塊中輸入文字，並設定字型大小與顏色，甚至文字的特效。字型建議使用「無襯線字體」，例如英文的 Arial 系列、Helvetica 系列（ Windows XP 中並無內建中文的無襯線字體，可用標楷體代替）。使用無襯線字體的好處是非常容易被閱讀。
經驗談：可以利用 PowerPoint 內鍵的「項目符號及篇號」來建構你的條列式文字段落，而盡量不要自行輸入項目符號或篇號，否則之後排版會很辛苦。

第四步：插入圖片

「插入」→「圖片」→「從檔案」，再選擇檔案後插入圖檔，再將圖案大小做適當拉伸。最大可以拉多大才不會有馬賽克？記得參考「第一步：先前工作」中的「收集解析度夠大的圖片」。
在圖片上連按二下，可以設定圖片的許多參數，例如大小、排列、邊框等等。
在「繪圖工具列」中的「對齊或等距分布」可以讓多張圖片完美地排列整齊，而「圖片工具列」中可以快速地設定「亮度」、「對比」及「裁剪」，是優化圖片的最佳工具。

第五步：插入 Excel 圖

做法很簡單，只要在 excel 中的圖案複製，再到 PowerPoint 中貼上即可。因為 excel 圖為向量圖，不用考慮 300 dpi 的限制，想把圖拉多大就拉多大。請注意該 Excel 圖是否已經符合科學論文之格式，例如標示座標軸名稱及單位等等。

在 PowerPoint 中的 excel 圖可以直接按二下以進行編輯，包括資料、顏色、樣式、字型等。但是調整這些參數需要大量的經驗才能得心應手（尤其是黑底白字的海報上，難度更大），同學只能靠多多練習。

第六步：插入表格

表格的製作方式其實與 word 完全相同，但要達到正式科學格式也是需要經驗才能順手。也要注意表格是否符合科學論文之格式，例如沒有縱線等等。

第七步：最後排版

完成所有內容之後，最終排版才是最難之處。要怎麼排才能排得均勻、整齊、好看？以下阿盤分享幾個技巧。

多多利用「繪圖工具列」中的「對齊或等距分布」將物件排列到整齊且均勻。
雙按圖片可以加入邊框，以增加圖片的明顯及立體。
小心地使用分隔線，以區分不同的大網或主題，但用太多只會讓觀眾覺得混亂。
設定文字方塊的底色，可以讓重點更搶眼。這招常常用在大標題及次標題。
若是對外比賽，可以把校徽放在海報上，大家都喜歡這麼做。
如果excel圖或是表格怎麼調都調不好，那就先休息一下吧。這是需要耐心的。
在整張大海報的四周也可加入邊框（利用方塊圖），單純地增加美感和視覺集中度。
如果你喜歡用「半透明」或「漸層」，那就快去學 illustrator 或 coreldraw 吧！PowerPoint 在這方面是廢物（你印出來就會明白了），請小心使用。
如果你懂得將照片去背只留下主體，那不彷做一張去背照放在海報中，因為這會讓大家都想去看你的海報

完成作品並輸出

將海報先輸出為 PDF 檔

把作品送到印刷廠，最怕的就是送回來的海報和自己設計的長不一樣！PDF 檔在這時候可以幫你的忙。PDF 檔案其實是給「印表機」看的檔案。換句話說，理論上 PDF 檔的樣子就是印出來的樣子。PDF 也會將字型鑲嵌在內，所以不用怕是不是用了很特別的字型造成別台電腦不認識你的字型。

要輸出 PDF 檔，不一定非得購買 Adobe Acrobet Professional，其實有許多免費軟體可以達成：

CutePDF Writer
PrimoPDF
OpenOffice內建有輸出 PDF 檔的功能
Mac OS X中即內建 PDF 輸出功能

到印刷廠的注意事項

記得一定要仔細對廠商說明輸出大小。在給檔案之前一定要在廠商的電腦上預覽過你的檔案，仔細確定所有細節後才能給廠商檔案。

有一項非常容易發生的問題，就是輸出的顏色和你在螢幕上看到的不太一樣。這是很正常的，真的。這世界上沒有一台螢幕的顏色是準的。

有一點很重要，就是輸出 PDF 的設定一定要很仔細檢查。因為 PDF 生成的過程可能會壓縮點陣圖（如照片）的 dpi。假如最終輸出的 dpi 小於 300 將使點陣圖失真。這一點一定要好好檢查，最好的辦法就是不要壓縮點陣圖。

以試算表計算簡單的統計檢定

2007-08-13T07:54:00.000-07:00

本文分享阿盤製作的 Excel 與 Open Document 檔案，可進行數種基本之統計檢定。

為什麼要做這件事？

在 2005 年的春天，阿盤參與了東海大學生命科學系的低年級課後輔導工作。那時候我教了學弟妹一些統計。為了讓他們可以更喜歡「玩」統計，所以阿盤以 Excel 上設計了七個統計檢定，目的只是希望他們可以因為玩統計而喜歡統計。當時這些檔案也分享在 PTT 站的 Share 版，接受網友下載。但由於這個 Excel 檔放置於學生空間，隨著阿盤大學畢業就消失了。後來陸續又有 PTT 的網友希望取得該 Excel 檔，為了一勞久逸，阿盤將該 Excel 檔放置於 Google Page 中，以永久性的保存。

為了正義與公平起見，另外製作了 Open Document 格式之試算表，可於 OpenOffice 2.0 以上使用。

作者聲明

本網頁及附屬Excel檔案僅供教學使用，且對於統計結果僅為參考，請不要做為學術研究及商業營利使用。若造成損失請使用者自行負責，與本人無關。

下載檔案

歡迎下載 Excel 統計檢定檔案或下載 Open Document 統計檢定檔案（可用 OpenOffice 2.0 以上開啟）。

統計項目

卡方比例檢定
Chi-square ratio goodness-of-fit test
單樣本T-test檢定
One-sample T-test
雙樣本T-test成偶檢定
Two-sample paried T-test
變異數相等之雙樣本T-test非成偶檢定
Two-sample unparied T-test with equal variance
單因子變異數分析
One-way ANOVA
簡單線性相關
Simple linear correlation
簡單線性回歸
Simple linear regression

撰寫環境

該 Excel 檔以 Microsoft Office 2004 Excel on Mac 的環境撰寫，並儲存成「Excel 97-2004 & 5.0/95 Wordbook」以供最大相容性。而 Open Document 檔則以 NeoOffice 2.2.2 之下製作，相信 OpenOffice 也可以開啟並運作（已於 OpenOffice 2.0 on Mac 測試）。這七項不同的檢定儲存在不同的工作頁，請在 Excel 或 OpenOffice 左下方列表切換。

該 Excel 檔中所運用的函式都是內建於 Excel 中的基本函式，不包括任何套件，以供最大相容性。

使用說明

怎麼用這個檔案？: 在每一個檔案中，灰色儲存格供使用者輸入數值資料。在輸入之後，Excel 會自動立即計算結果，所有重要的結果會以紅色、綠色及藍分表示。另外在下方會自動作圖供使用者參考。該圖的數值為組內平均數或頻率，誤差為組內 SD。在該 Excel 檔中已經簡單表示如何使用，相信只要粗淺明白統計檢定的朋友都可以明白如何使用。
重覆數太多了，灰色格子不夠用！: 每項統計檢定檔案中，僅供有限的重覆或組別數量。若有需要增加，請自行修改該 Excel 檔。
灰色格子沒填滿，可以嗎？: 在本 Excel 檔中，已經儘量以使用者方便為前題的方式製作完成。所以只要有填數字在灰色格子裡，Excel 程式應該都可以正常運作，不一定非得填滿或從最上面的灰格子開始填。如果違反了統計公式運作而造成錯誤（例如使分母為零，或是配對的數值沒有對齊），不會造成可怕的後果（如 Excel 當掉），只會使最終結果算不出來或不正確而已。
計算結果不太對勁……: 請注意，這個 Excel 或 Open Document 統計檢定檔案中，並不包括任何對於該統計檢定之前題。若資料與統計檢定之前題不符，則該檢定結果自然也不會正確，但 Excel 中並不會告知。除此之外，所有檢定的 alpha 值皆設定為 0.05。若需改變請自行修改語法。當然，該檔案也可能有錯誤，請原諒阿盤。

常用的免費 OS X 軟體

2007-08-13T05:38:00.000-07:00

本文可能已經過時，請斟酌參考。

本文介紹阿盤在 Mac OS X 上所使用各種功能之免費軟體。

瀏覽器、BBS

FireFox — 功能強大且安全的瀏覽器: Firefox 無論在各種作業系統之下都是非常叫好的瀏覽器。開放原始碼，而且具有許多套件可以使用，大大增加瀏覽網頁的效率。Firefox 更新地非常頻煩，感覺蠻安全的。目前 OS X 之下的 Firefox 常常會有「中文掉字」或「中文粗細體混雜」的問題，一直困擾著許多使用者，但最近聽說有越來越正常的趨勢。和 Safari 相比，Firefox 雖然功能強大，但速度並比不上 Safari。在 2.0 版本之後，Firefox 內建的功能可以更容易地將 Safari 的 Bookmarks 匯入到 Firefox 中，非常方便。
Camino — Firefox 的兄弟: Camino 使用的核心和 Firefox 相同，但是完全替除 Firefox 的外掛擴充功能，另外 Camino 的外觀也更相似於 OS X 的應用程式，執行速度也比 Firefox 快。簡單來說，Camino 算是 Firefox 的陽春版。假如用不上 Firefox 外掛擴充功能，Camino 真的非常適合你。在 Safari 無法開啟的網頁，我也必定先以 Camino 試試。
SafariStand — 讓Safari進化吧！: SafariStand是Safari的外掛，提供Safari許多強力的功能，其中最吸引人的是「縮圖標籤側邊欄位」。它讓Safari的側邊產生一個新欄位，其中排列了所有標籤的網頁之縮圖，方便常常使用大量標籤頁的Safari使用者。另外我非常喜歡SafariStand可以將整個網頁以「不分頁」的方式轉成PDF，而且只要按一個右鍵就行了。此外SfariStand提供使用者極容易地下載.flv（例如Youtube）的影片檔。SafariStand並不是一個應用程式，只是Safari的外掛，所以不是直接拉到應用程式中來使用。請在下載之後請仔細閱讀其說明文件，才能正確地讓他出現在Safari之中。
Vienna — 功能完善的RSS/ATOM瀏覽器: 雖然Safari的RSS/ATOM瀏覽功能非常好用，但你也可以試試功能完善的Vienna。Vienna也可以下載免費的CSS樣式，使你觀看RSS/ATOM更順心。
AlienBBS — 麻雀雖小、五臟俱全: 在BBS軟體之中，AlienBBS可稱是最快、最簡單、最實用、也最美觀的佼佼者了！而且作者也是台灣人呢！雖然有人認為AlienBBS並不能做太多的設定，例如自定字體、自動login等功能，但其實非常實用，可以說上簡單確完善。AlienBBS也有內建的防呆功能，且快速鍵也很有邏輯，所以容易記憶與使用。AlienBBS也有標籤列，也支援文中URL自動擷取。目前美中不足的是並不支援雙位元文字辨視，所以刪除中文時要辛苦一點。不過由藍永倫所製作的MulderPatch（一個AlienBBS的外掛）已經為AlienBBS解決了雙位元文字的問題，也增加了許多功能，使得AlienBBS算是非常強大且實用了。

通訊、網路

Adium — 通訊全都包！: Adium相容於MSN Messenger、Yahoo Messenger、.Mac、AOL等等通訊媒介，可以說是全都包了！Adium最大的好處就是使用者設定非常強大，無論是帳號、聯絡人清單、對話視窗、甚至於Dock中的圖示都可以自定。除此之外，官方網站也提供了超級多的資源（總稱為Adium Xtras）讓使用者下載。當然，Adium也支援多帳號同時登入，無論是不是相同的通訊系統。Adium也支援Growl提示（本文也有介紹），讓每個訊息都不會漏掉！真的是「全都包了！」
Preteus — 簡單好用的文字通訊軟體: 與Adium相比，Preteus的功能比較陽春，但是更容易上手。個人的經驗也認為Preteus比Admin穩定許多啦。事實上，我很少用Adium，反而是以Preteus來使用MSN。
Skype — 音訊視訊都可以: Skype在PC及MAC上的功能幾乎沒有差異，目前也已經推出視訊功能。使用上也非常簡單啦，所以阿盤有時候還蠻喜歡到PTT的Wanted版跟團到Skype上唱歌呢！
Cyberduck — 功能齊全支援多的FTP軟體: FTP軟體中，不用錢當中就屬Cyberduck最有名氣。Cyberduck支援FTP、FTP-SSL與SFTP協定，當然也有書籤之類的功能。我喜歡Cyberduck的原因是他支援Growl（本文中也有介紹）顯示、也支援某些文字編輯程式（如Smultron，本文中也有介紹）的線上存取，讓寫網頁的阿盤非常方便。但Cyberduck最大的缺點就是不夠不穩定，常常莫名斷線或是自動終結。
URLcollector — 功能齊全的續傳軟體: URLcollector是功能非常齊全的續傳軟體。它也支援從網頁擷取可下載的元件，愛抓圖的朋友一定會喜歡它！目前的缺點是太久沒有更新了，而且只有日文介面可以看。一開始可能不易上手，但看看漢字和用快速鍵猜猜看也可以明白怎麼用的。
Google Earth — 全世界的影像都在電腦裡: Google Earth，這麼強大且出名的軟體，應該不用多說什麼了。

音訊、視訊、多媒體

Sing that iTune! — 同時展示歌詞、封面與控制iTunes: Sing that iTune!是個自動上網搜尋歌詞的Widget，可以搜英文、日文及中文歌。搜尋到的歌詞也可以自動匯入iTunes資料庫，也可以自重儲存在一個資料夾中成為你自己的「歌詞資料檔案庫」。它在iTunes執行時會自動從iTunes資料庫中讀取專輯封面，若無法取得則也可以從Amazon自動下載，非常方便。它也能簡單地控制iTunes的動作，算是功能齊全的iTunes相關Widget。
iTunes Artwork — iTunes的CD盒: iTunes Artwork是一個可以展示CD封面的Widget。在畫面中，CD封面會被放在一個CD盒子裡，看起來就像真的一樣。若iTunes資料庫中沒有該歌曲的CD封面，iTunes Artwork也能自動從Amazon下載。這個Widget的功能並沒有上述的Sing that Itune!多，但阿盤還是非常喜歡它，最主要是因為阿盤喜歡收集CD，而它讓阿盤在Dashboard裡看到像真實CD般的畫面。
Hinet Radio — 讓MAC也能收聽台灣的數位廣播: Hinet所提供的線上數位廣播非常豐富，只要有網路想聽什麼就聽什麼，但播出的界面必須依靠Windows的IE與Windows media player。但Hinet Radio打破了這個規定。這個簡單的Widget可以收聽Hinet所提供的線上數位廣播，讓MAC族非常受用！印像中作者也是台灣人呢！非常感謝他的貢獻。
GuitarChords — 吉它和弦資料庫: 這是個專門設計給彈吉它朋友的Widget，功能與設定都很齊全。你可以在GuitarChords上找到所有調性的和弦組合，並看清楚有哪些不同的按法。除此之外，也可以設定每條琴弦的空弦音、左撇子或右撇子，最後還能聽聽看midi合成的和弦音，算是一個很方便的吉它和弦工具。
Perfect Pitch — 弦樂器的調音程式: Perfect Pitch是一個偵測MIC傳入聲音頻率的應用程式。簡單來說，用來調音就是了！不過當阿盤想調吉它音高而使用iBook的內建MIC輸入聲音時，十分容易受到其它噪音影響而造成干擾，所以並沒有真的比一般的電子調音器好用。但手邊沒有調音器或其它樂器時，還是能大概調出正確音高。
WMV Player — 以QuickTime開啟WMV與WMA: WMV Player可以算是QuickTime的一個外掛，目前被Microsoft拿來MAC上開啟WMV與WMA格式的解決方案。灌了WMV Player之後，就連FInder都可以預覽WMV或WMA了。安裝WMV Player之後，它的設定會出現在系統偏好設定的其它那類群裡，這一點比較特別。
Perian — 讓QuickTime開啟更多種影音檔: 這個QuickTime的外掛可以讓QuickTime開啟FLV、3ivX、DivX、Flash Screen Video、MS-MPEG4、Sorenson H.263、Truemotion VP6、Xvid等等的影音檔案。有了它，從Youtube下載的影片就可以直接由QuickTime觀看了。大力推薦！
MPEG Streamclip — 多媒體檔案的轉檔高手: MPEG Streamclip可以將許多不同格式的視訊檔做簡單的剪輯，再輸出成其它的壓縮檔案格式。也就是說，如果沒有QuickTime Pro的話，利用MPEG Streamclip也可以製作H.264的.mov檔案、將視訊檔的聲音擷取成音訊檔、把VCD轉成mpeg1或.mov檔……MPEG Streamclip也可以處理音訊檔案。總而言之，對於多媒體檔案的擷取與轉檔之功能強大！
ffmpegX — 超強地影音轉檔工具: ffmpegX是基於鼎鼎大名的ffmpeg的影音轉檔工具，能輸入/輸出的檔案種類非常多。不過在使用之前，必須先下載一些函式庫（軟體本身會有所提示）。不過，ffmpegX並不是免費軟體，但是可以無限時間/次數地使用，值得推薦給大家。
MacTheRipper — 原封不動地備份DVD: MacTheRipper用來把DVD原封不動地拷進別的磁碟機，並不參與轉檔。換句話說，它幫你直接把DVD拷備到別的地方去罷了。MacTheRipper也可以針對DVD的章節來拷備資料，但通常電影DVD不會只看某個章節吧。如果你有非常珍愛的DVD，又捨不得拿出來看，那你一定需要MacTheRipper。
Handbreak — 將DVD轉成高壓縮的多媒體檔案: Handbrake可以將DVD轉成別的壓縮格式，包括H.264。使用介面可以選定DVD章節、字幕、聲音、影片尺寸等等。如果想收集大量DVD影片，可以利用HandBrake將DVD壓成H.264再儲存於硬碟之中。通常DVD原尺寸大小而言，壓至1800 kbps（相當於只有1/4的空間）都能看起來「沒有失真似的」。不過生成H.264需要極大量的運算，所以想使用HandBrake來壓整部電影對notebook會非常辛苦的。
Max — 強大的CD/音訊轉檔工具: 雖然iTunes已經非常完善了，但你仍然可以試試Max。Max可以把CD音訊轉成其它的音訊格式，例如MP3。Max也提供轉檔功能，包括MP3、Ogg、FLAC、AAC、Apple Lossless、Monkey's Audio、WavPack、Speex、AIFF、及WAVE，其中許多格式是不常在MAC上使用，但PC上卻常見的。我習慣利用Max將iTunes中的AAC轉成MP3，以讓不支援AAC的MP3隨身聽讀取音訊。雖然iTunes也可以完成這項工作，但在iTunes中轉成MP3的檔案也會被收錄在iTunes Library中，我還要花時間把這些多餘的MP3清理掉是件煩麻的事，所以才交給Max來做。Max編碼MP3的方式是基於LAME函式庫，所產生的MP3品質也算是優異了。
VLC — 讓MAC族看影片不再跛腳: 這個軟體可以應付九成以上的多媒體檔，連看DVD也成。除了接受的檔案豐富，對於使用者自定也非常豐富，例如DVD的字幕、DVD音訊、DVD去條紋等等。總而言之，預到QuickTime開不了的多媒體就試試VLC，再不行就試試MPlayer OS X（在此並沒有介紹），大多不成問題。

靜態影像

Screenshot Plus — 全功能的螢幕擷取Widget: 這是一個功能強大的螢幕擷取Widget！可以擷取單一視窗、自選區域、定時全螢幕擷取、甚至連Dashboard裡的Widget長什麼樣都抓下來！擷取之後提供簡單的預覽，並且可選擇儲存位置或是匯入什麼應用程式以進行後製。非常非常強大！
Xee — 簡單又方便的圖片瀏覽程式: 簡單來說，Xee的功能就是早期的ACDsee，作用於擊點單一圖形檔之後，還可以上下瀏覽同資料夾之下的其它圖形檔。雖然iPhoto讓我們減少了這樣的機會，但有了Xee，找一些圖形資料還是很方便的。阿盤的習慣是，先在Finder之下設定所有的GIF、JPG及PNG都預設以Xee開啟，接下來就如果過去用ACDsee的經驗完全一樣了。

文書、開發

Smultron — 功能齊全的文字編輯程式: Smultron是標準的文字編輯程式，可以讓開發者有個順手的工作園地。Smultron支援的語系完善，對於各種語法的高亮度提示也很完備，也提供HTML預覽，最近也加入了英文拚音校正，最棒的是它可以與Cyberduck（本文也有介紹）這個FTP程式一同工作，讓我在線上存取檔案且直接編輯都沒有問題，就連阿盤現在也是用Smultron來寫這個網頁呢。我曾經寫信給作者，表達阿盤對他的感謝之意，他的回信只有四個字：Thank you very much。哈哈。
TranslateIt!（原Mdict） — 非常好用的字典平台: TranslateIt!是一個字典的平台，並不包括中文字典。所以要先去下載免費的中文字典資料庫，再利用TranslateIt!這個平台來使用字典資料庫。TranslateIt!可以同時接受許多不同字典資料庫的資料，所以可以同時有英英、英漢、漢英、英日……字典，也可以有牛津、朗道……字典，而且並不需作任何切換動作。中文相關字典資料庫可以到stardict.sourceforge.net來找找。
PDFLab — 合併PDF的實用工具: PDFLab的功能就是合併多個PDF檔案。雖然內建的Automoter也可以合併檔案，但利用PDFLab卻較直接也簡單。在合併PDF也可做一些選項，例如只合併該PDF中的某些頁面等等，小小工具大大有用。
MacTeX — MAC上的TEX完整套件: MacTeX是TEX在MAC上的一個套件，內容非常完整。2007版的MacTeX已經加入了CJK及XETEX，對中文文件的編排更是方便。需要使用TEX系統的朋友，可以不用再以i-installer一項一項的下載，而使用MacTex一次幫你安裝好。

系統、檔案、其它

OnyX — 強大的系統維護程式: OnyX是很知名的系統維護與設定的程式，裡頭包含許多調整OS X隱藏設定與系統維護的選項。雖然OS X每天都會自動維護自己，但往往是在半夜進行的，若關機就錯過了呢。有時候用OnyX整理整理系統，讓MAC的效能維持下去吧！值得注意的是，這類系統維護的程式都是針對特定的作業系統版本而設計的，例如Panther就有給Panther用的OnyX版本，Tiger就有給Tiger用的，千萬不要用錯了！此外，每當OS X有重大更新的時候，個人會停用這類的系統維護程式，等過幾天該程式推出改版，才去下載新版來用，以免造成系統被舊版的系統維護程式搞壞了。
Yasu — 介面簡單的系統維護程式: 與Onyx相似，Yasu也是個系統維護程式，但並不包括調整OS X隱藏設定的功能。使用的介面非常簡單，這是我最喜歡它的原因。有一些人認為Yasu跑過之後OS X真的很順，我也沒有實驗過到底順不順，這種順不順的問題就見人見志囉！一樣的，請下載適用於目前使用OS X版本的系統維護程式，原因和OnyX的介紹一樣。
Growl! — 隨時提醒有事發生: Growl是一個提示程式。它本身沒有介面，必須讓其它程式來呼叫它。舉個例子，只要先在Adium中開啟Growl提示功能，若有人用MSN給我訊息，則可以在螢幕上顯示一個神奇的提示文字區域。這就和Windows中使用MSN有事件發生時在螢幕右下角會有個提示訊息欄跳出來的情況非常相似。隨著支援Growl的程式越來越多，例如Cyberduck（一種FTP軟體，本文也有介紹）、MDict（字典平台軟體，本文也有介紹）等等，你將會發現Growl真的非常好用。另外，Growl產生的提示訊息外觀也可以自己選擇，都蠻美觀的，至少我朋友都說比MSN的好看多了。安裝Growl之後，它的設定會出發在系統偏好設定的其它那類群裡，這一點比較特別。
The Unarchiver — 比內建BOMArchiveHelper更好用的解壓縮程式: The Unarchiver可以解許多種類的壓縮檔類型，功能完善。我最喜歡The Unarchiver是因為它可以把解壓後的壓縮檔自動丟到Trash裡，蠻方便的。第一次使用前請先到The Unarchiver的Preferences裡做好所有的設定。
iStat nano — 小而美的系統狀態Widget: iStat nano是一個小小的Widget，但可以顯示許多系統現況，例如中央處理器、記憶體、磁碟機、溫度、電池、網路……。iStat nano有個哥哥叫iStat pro，功能與iStat nano一樣，不過內容較多，卻也比較佔Dashboard的空間。我比較愛小小的iStat nano啦。
Youpi Optimizer — 刪除不必要的語系檔案: 在OS X之下的應用程式中，其它都藏著不同語文的語系檔案。而「系統偏好設定→國際設定→語言」的順序就決定了該應用程式被啟動時所載入本身語系檔的順序。然而，有太多的語系檔也許使用者一輩子都用不著吧。Youpi Optimizer可以把應用程式中的語系檔找出來，並且讓使用者勾選想要刪除的語言。以阿盤為例，包留了英語、日語、正體中文及簡體中文，其它所有的語系檔都刪除，竟也生了半GB出來，很神奇吧！
UNO — 轉換程式的金屬與水樣底色: 在作業程式中，有「金屬」（例如Finder與Safari的工具列）及「水樣』（例如系統偏好設定的底圖）二種視窗版面，還有一種稱為「Unified（例如系統偏好設定上方部分）」的版面，受到許多人的喜愛，但偏偏許多應用程式並沒有選用Unified版面。UNO可以改變OS X中所有應用程式預設的版面，將金屬或水樣改成Unified。改過之後，OS X真的變美很多。
Stellarium — 看星星就像真的一樣: Stellarium算是一個教育軟體，就是模擬天空的星象，而且非常逼真。在Stellarium中可以設定經緯度與時間，讓你無論晴天雨天都可以看星星和認識星座。

Apan’s Notes

Create a NA matrix in R

Statistical computing by using R

A function of Pearson/Deviance goodness-of-fit for (generalized) linear regression

Source Code

Usage and example

Games-Howell post-hoc in R

邏輯斯迴歸、勝算與勝算比的關係 (logistic regression, odds and odds ratio)

公式一覽

寫在前面

名詞翻譯

僅包含常數項的迴歸式

包含一個二元類別自變數的迴歸式

包含一連續變數作為自變數的迴歸式

結論

Least-significant difference (LSD) in R

Specifying a reference level in R: function relevel

阿盤的 Beamer 中英字型配方

學習 HTML, PHP 及 Database

MacTeX + CJKutf8 + DVIPDFMx 的新增中文 TTF 字形

DVIPDFMx 的優點

運行環境及所需軟體

開工

利用 ttf2tfm 產生 TFM 檔

把 TFM 搬到正確的位置

把 TTF/TTC 搬到正確的位置

修改 cid-x.map

製作 FD 描

更新 TeX 索引

TEX 檔範例

與 TeXShop 的配合

與 Ghostscript 的配合

CJK 全型標點符號「對齊」行尾

先前的做法

先前的做法有些問題

修正後的辦法

新辦法的模型

尚未達成的目標

替 MacTeX 新增 CJK 之 Type 1 字型

前言

工作環境

進行步驟

後記

以 MATLAB 轉換 OOIBase32 原始檔

如何產生 PDF 文件？

PDF 文件與其優點

如何產生 PDF 文件

XHTML 鑲嵌 QuickTime 媒體

引言

實際做法

1. <object><embed></embed></object> 巢狀結構語法

2. 以 Javascript 插入 <object><embed></embed></object> 的巢狀結構語法

3. <object><object></object></object> 巢狀結構

結論

利用 PowerPoint 製作學術海報

寫在前面

開工了！

第一步：先前工作

第二步：PowerPoint 基本設定

第三步：製作海報中的文字

第四步：插入圖片

第五步：插入 Excel 圖

第六步：插入表格

第七步：最後排版

完成作品並輸出

將海報先輸出為 PDF 檔

到印刷廠的注意事項

以試算表計算簡單的統計檢定

為什麼要做這件事？

作者聲明

下載檔案

統計項目

撰寫環境

使用說明

常用的免費 OS X 軟體

瀏覽器、BBS

通訊、網路

音訊、視訊、多媒體

靜態影像

文書、開發

利用 `ttf2tfm` 產生 TFM 檔

修改 `cid-x.map`

1. `<object><embed></embed></object>` 巢狀結構語法

2. 以 Javascript 插入 `<object><embed></embed></object>` 的巢狀結構語法

3. `<object><object></object></object>` 巢狀結構