以下是我最近兩天都解決不了的事，妻子覺得我在自找麻煩。 Macbook 在本周一突然壞 Harddisk 。至前天終於買了新的 harddisk 重新安裝。安裝時有麻煩，但對比以下的，只屬小事，噩夢在這之後才開始。 裝好 harddisk ，方發現 optical drive 也出了問題。測試幾次之後，才可以認出 Macbook 跟機的 Installation DVD 。趁這個珍貴機會，就將 Mac OS X Tiger 安裝進去。從此之後， optical drive 全死了。就算我試...
感謝讀者 Paladinoflab 提供了刪除 X1, X2, X3 的方法。跟據他提供的方法，我再想出更好的解決方案。 由於例子中的 table1 的第一行資料 (即 table1[1,]) 才是想用來當輸出 table 的第一行，故此要做以下的東西： 1. 將 table1[1,] 的資料變成 column name 2. 刪去 table[1,] 另外，其於 t() 的限制，要用到 data.frame(t()) 這樣樣衰的組合。 (( t 原意是用於 matrix ，但亦可用於 d...
根據非正式統計，本博講 R 或 Stat 的話 ((另外三個題目是小說、讀書及棋類遊戲)) ，讀者人數會銳減。就算本人講三腳貓的電腦編程，或者政治塗鴉，關注人數會較多。 但博客本來的精神，是忠於自己的偉大。對於一個下半生都應該會投入統計事業騙口苦飯吃的人來說，沒有比在博客分享統計電算之樂更正常的事。 ((其實是因為我這個人是沒有生活的，像個活死人，除了講這些東西之外，我的生命空虛得像個荒廢的防空洞。我的博客是沒有精彩如重遇當年暗戀的女孩或者飲飲食食之類的題才。)) 統計的其中一個重要部份，...
最近看了兩篇東西，第一篇是 Nature 期刊探討 PhD 問題的專題，第二篇是明報報道，研究生出路窄 增薪遠遜本科生。 最近有個覺悟，或者我此生都沒有機會讀 PhD ，現在這個假郎中的生涯，要麼是永久持續，要麼是隨時中斷投身地產霸權。由我來討論 PhD 問題，就有如高考物理揸槍講核輻射一樣，是沒有「資格」的。如果你仍想聽下去，記得我是「沒格」的評論員，本篇為無格的評論，故此內容全不可信。 Nature 專題指出各地大學濫發 PhD ，由其是新興國家如中國及印度，發達國家如美國及日本同樣濫...
愈沒有實力的人，在 A vs B 的討論中，會愈大聲，立場愈激進。因為他們以為參與這些討論，代表他的實力很高。 最常見的，是在網上抄錄一大段以前討論觀點，自己再說一次。 例如在 vi vs emacs 的討論中，很多人會找出 emacs 像個 OS 卻不像編輯器之類的觀點來講。這個觀點的確是對，但卻是 Cliché 中的 Cliché 。現代電腦用 Virtual Machine 同時行幾個 OS 已經不是新聞，難道行個 emacs 系統需求會高過行一個 VM ？還有那些 "EMACS ...
UbuntuHelpGuy 問， Ubuntu 用戶對 Unity 的意見。 我的意見是： Unity sucks! No doubt about it. ((最近甚至覺得 Mac OS X 麻麻地好用。因為 Mac OS X UI 阻住做野。)) 我在工作間的主電腦是用 Ubuntu ，最近已轉去 Xubuntu ，因為我不想轉成 11.04 時要用 Unity 。而家中 Macbook 安裝時已用 Xubuntu 。 Xubuntu 11.04 beta 的 XFCE 4.8 tu...
如果你能夠坐著一小時看完整個 lecture ，恭喜你，你的減肥計劃已經成功了一半。
研究一下那段 code ，是很簡單的結構，但卻寫得超長的。 這個 code 我想有六七年歷史，最初代版本是 SPSS Syntax ，後來根據 logic 在 R 重寫。這 code 不停重覆的結構是 # 今時不同住日，現在會就會想試下去改，當係教自己怎樣寫 program 也好。當年寫這些 code 應該用了不少次 copy and paste ，現在看應該要 generalize 那個 pattern 。 結構上的另一個問題，是數據和邏輯混在一起，要想個方法將兩者分拆。功能上也有...
今早出門口上班時，不小心用手掙撞到妻子的心口。我想她應該不太痛，除了連忙說對不起，腦中突然想到過期政治笑話一個：「要不要去瑪麗醫院驗傷？」但是這個新聞已經過了氣一個月，也不講出來了。 三月中旬發生了太多的事，三月頭特首報稱被襲事件，像是數年前發生似的。 每一單新聞都有「此日期前最佳」，期限有長有短。有些新聞過了氣，記者就不太會去跟進。但是，新聞正在有效日期時，就會有多人關注。近年的傳媒愛好 yellow journalism 以及甚麼 agenda setting ，政府也樂於大搞 doc...
查實用了 R 好耐，已經寫落不少 Code 。最近有點想將這些 Code 整理一下，將常用的東西全部儲在一個叫 csx.R ((別問我 csx 是甚麼，是亂改的。或者可以叫做 CSX: Statistical eXtension 吧。)) 的檔案內。下次要用時，可以用 source() 將所有常用的 function 存入。 上篇文才說過，識 R 只代表識 R 的規則，不代表寫得出好 Code 。我正是只識 R 規則的人。以下是我從 csx 找出來的一段 code 。 # = 16 &...
之前上完了 MIT 6.00 課，可謂獲益良多。教這課程的其中一位 Professor 叫 Eric Grimson ，最近榮升 MIT 校長。那課堂用的教學語言是 Python ，但是 Grimson 不停說自己是 old time Lisp programmer 。 從 MIT 6.00 結業，下一課叫 MIT 6.001 - Structure and Interpretation of Computer Programs (SICP)。 由於 MIT 的 Video lecture...
我的 research 生涯，做得最多的研究類型是 Retrospective case series (RCS, 也稱為 Chart review ) ， Prospective 的反而不多。 RCS 有其限制，但是 Clinicians 們卻以為 RCS 是平靚正的選擇，以為 reviewer 見到其 sample size 夠大就會刊登。 多年來，已經知道怎樣才能將 RCS 做好寫好，而且送到「較佳」期刊刊登的法門。講真，期刊連 Randomized Controlled Trial...
這不是用 R 畫的，是用 Pylab 。數據來自 98 至 09 年 HA Statistical report 。 圖中每格是 500 人。 Source code: from pylab import * medical = [3581,3804,3979,4222,4460,4617.5,4871.8,4858.8,4898.1,4935.35,5057,5160.5] nursing = [19614,20435,19880,19746,19710,19567.5,1...
新聞不停轉載 5 克 Plutonium-239 足以殺所有人。 這是怎樣計出來的呢？我也很想知道，雖然我只識得三年級的數學。 Pu-239 的 LD50 是 10 ng/kg 。一個 ng (納克, nanogram) ，即是 10-9 克。 LD50 是怎樣計出來？這個數值是指如果一群人都攝入這個劑量的 Pu-239 ，有 50% 人會死。 全球有 70 億人。但是怎樣算全球人的平均體重？我想，將全球人的平均體重設定在 80 kg ，是 overestimate 了。但由於今次是測上限...
曾在網友間討論某人是否好煩的男人。慢慢地覺得，我也是很煩的男人。 煩的人都有一個特點，就是指指點點，高高在上似的。友人說，這些人生活苦悶，很多意見，但老婆不聽他講耶蘇。 邊緣回望，我發現我也有這個特點。請讓我爆一下。 例如我不停指正廣東話讀音。某人打電話來，我次次都要和他說，我叫阿 /hɔŋ˥/ ，唔係阿 /hɔn˥/ ，已達三四年。去到一個位已經覺得唔好意思再話佢，覺得自己好撚婆媽，反思自己係唔係應該收聲。 "I DON'T CARE WHAT YOUR FUCKING NAME IS!...
Previous — 36 of 210 — Next