Jun 21 2014

R 期末上機考上傳狀況

目前為止,尚未收到以下同學的上機考程式

400651195
400657036
499650306

Apr 20 2014

統計電腦應用程式練習題

資料檔: babies
http://netstat.stat.tku.edu.tw/show_samples.php#c_babies
資料檔: cancer
http://netstat.stat.tku.edu.tw/show_samples.php#c_cancer
0. 請下載 babies.csv, 並在程式中讀入, 儲存成 babies 資料框架變數
1. 請用 babies 資料檔做一個一因子(One-Way)設計分析
2. 請用兩種不同方法將 bwt 變數依照 Q1, Q2, Q3 為切割點,轉碼成 B1, B2, B3, B4 分類變數
3. 請將 cancer.csv 資料檔讀入 R, 並分別使用 stack 與 unstack 做資料轉換
(1) stack 處理後儲存成 cancer2
(2) 用 unstack 把 cancer2 轉成 cancer3

Jan 6 2014

週二 R 上機考教室

B204, B206  ( 6:10 pm 開始)

以下學號同學在 B206, 其餘在 B204

400656129

400657044

400657135

400657366

401657084

401657092

401657118

401657274

401657316

401658033

402658024

498650539

498650869

498651701

499490125

499490430

499650017

499650025

499650082

499650173

499650231

499650264

499650496

499650561

499650645

499650660

499650736

499650777

499650827

499650868

499650884

499650942

499651098

499651130

499651361

499651429

499651502

499651643

499651726

499651791

499651817

499656063

499657186


Jan 6 2014

統計計算補充講義

http://steve-chen.net/R/stcmp.zip (需密碼)

http://steve-chen.net/R/package.zip (需密碼)


Dec 9 2013

R 期中上機考成績

400...031    莊   35
400...049    許   15
400...064    柯   20
400...072    梁   70
400...122    王   100
400...130    莊   45
400...148    周   15
400...189    楊   35
400...213    楊   55
400...262    郭   20
400...270    王   85
400...296    江   25
400...304    陳   20
400...312    王   0
400...320    曹   50
400...353    賴   25
400...361    章   50
400...379    劉   70
400...387    周   45
400...437    葉   45
400...445    王   85
400...486    邱   20
400...569    翁   15
400...619    張   25
400...643    童   100
400...692    黃   40
400...700    黃   60
400...726    劉   25
400...734    徐   35
400...767    陳   10
400...791    曾   50
400...809    趙   60
400...833    曾   70
400...882    林   100
400...940    陳   55
400...957    姜   25
400...013    劉   20
400...054    謝   45
400...062    許   25
400...088    蘇   25
400...120    魏   45
400...146    張   40
400...153    陳   15
400...179    王   30
400...195    蔡   35
400...203    王   35
400...211    吳   40
400...252    洪   10
400...278    楊   60
400...286    李   30
400...344    鄭   35
400...369    王   30
400...427    王   55
400...484    劉   15
400...534    陳   25
400...542    陳   50
400...559    李   25
400...567    周   35
400...591    黃   25
400...625    楊   35
400...682    張   15
400...724    李   35
400...740    李   75
400...757    李   55
400...773    黃   15
400...020    黎   20
400...103    黃   40
400...129    黃   25
400...044    程   40
400...135    唐   15
400...366    鄭   50
401...084    蘇   55
401...092    陳   100
401...118    彭   15
401...274    吳   75
401...316    賴   20
401...033    韓   15
402...024    李   25
498...539    林   25
498...869    康   15
498...701    李   0
499...026    劉   0
499...125    蕭   45
499...315    王   0
499...398    江   0
499...406    吳   10
499...430    劉   40
499...463    吳   0
499...017    王   30
499...025    吳   100
499...074    蔡   25
499...082    曾   60
499...173    林   60
499...223    陳   20
499...231    許   40
499...264    林   75
499...496    黃   50
499...553    何   0
499...561    盧   20
499...645    朱   40
499...660    翁   30
499...736    陳   20
499...777    陳   10
499...827    吳   75
499...868    林   15
499...884    呂   75
499...942    許   25
499...098    張   80
499...130    陳   20
499...361    溫   35
499...429    黃   35
499...502    趙   40
499...643    許   10
499...726    賴   85
499...791    陳   70
499...817    張   20
499...063    李   25
499...097    黃   10
499...147    羅   10
499...186    程   15

Nov 11 2013

週三(11/13)上機考 B206 教室名單

註:沒有被列入以下名單的同學,考試教室在 B204

座號 學生證號碼

2 400650049
8 400650189
12 400650296
17 400650361
18 400650379
23 400650569
26 400650692
27 400650700
30 400650767
32 400650809
34 400650882
35 400650940
36 400650957
38 400651054
39 400651062
42 400651146
43 400651153
44 400651179
45 400651195
51 400651344
52 400651369
56 400651542
58 400651567
63 400651740
72 401657084
73 401657092
75 401657274
76 401657316
78 402658024
81 498651701
83 499490125
84 499490315
88 499490463
89 499650017
90 499650025
92 499650082
97 499650496
99 499650561
100 499650645
101 499650660
102 499650736
103 499650777
104 499650827
107 499650942
108 499651098
109 499651130
112 499651502


Jun 26 2013

2013/06 時間數列期末報告成績

以下為時間數列期末報告各組成績:

第 1 組: 80
第 2 組: 85
第 3 組: 74
第 4 組: 80
第 5 組: 78
第 6 組: 84
第 7 組: 73
第 8 組: 77
第 9 組: 78
第 10 組: 82
第 11 組: 82
第 12 組: 80
第 13 組: 82
第 14 組: 78
第 15 組: 82
第 16 組: 76
第 17 組: 82
第 18 組: 83
第 19 組: 83
第 20 組: 78
第 21 組: 40
第 22 組: 72
第 23 組: 90
第 24 組: 87
第 25 組: 80
第 26 組: 75
第 27 組: 82


Jun 14 2013

時間數列期末報告:各組最佳模型比較

以下模型均採用台北市流浪狗資料,期數 123 期(至 2013/03)

註:模型之好壞還得視殘差分析結果、以及預測能力而定

     Models                         AIC         BIC
ARIMA( 0 , 1 , 1 )x( 1 , 1 , 1 )  1478.682*   1489.931*
ARIMA( 0 , 1 , 1 )x( 0 , 1 , 1 )  1484.456    1492.892*
ARIMA( 0 , 1 , 1 )x( 2 , 1 , 1 )  1479.357*   1493.418
ARIMA( 1 , 1 , 0 )x( 0 , 1 , 1 )  1491.690    1500.127
ARIMA( 1 , 1 , 0 )x( 3 , 1 , 1 )  1483.222    1500.095
ARIMA( 2 , 1 , 1 )x( 0 , 1 , 1 )  1479.357    1493.418
ARIMA( 2 , 1 , 1 )x( 1 , 1 , 0 )  1488.721    1502.782
ARIMA( 2 , 1 , 2 )x( 2 , 1 , 2 )  1485.458    1510.768
ARIMA( 3 , 1 , 0 )x( 2 , 1 , 1 )  1478.870    1498.556
ARIMA( 3 , 1 , 3 )x( 2 , 1 , 2 )  1479.871    1510.805
ARIMA( 4 , 1 , 2 )x( 2 , 1 , 2 )  1481.293    1512.227
ARIMA( 4 , 1 , 2 )x( 2 , 1 , 1 )  1479.315    1507.437
ARIMA( 5 , 1 , 2 )x( 1 , 1 , 1 )  1479.315    1507.437

介入模式:
ARIMA( 0 , 1 , 1 )x( 1 , 1 , 1 )* 1472.350**  1486.409**

此模型是 ARIMA(0,1,1)x(1,1,1) + 97年1月開始的 Step function 介入模式
介入模式的R 程式碼:

X1 = rep(0,length(Z))  # 製造出 123 個 0
X1[61:length(Z)] = 1   # 第 61 期(97/01)開始到最後改成 1

fit2 = arima(Z,order=c(0,1,1),seasonal=list(order=c(1,1,1),12),xreg=X1)
fit2
k = 1+2+1+1 # MA(1)+SAR(1)+SMA(1)+w0+sigma
bic = -2*fit2$loglik+k*log(length(Z))
bic

May 14 2013

R 期中上機考成績

http://steve-chen.net/R/2013mids.txt


Dec 26 2011

SAS 期末上機考範例

1 . 請下載 iris data , 並讀入 SAS, data set 名稱為 iris, 變數名稱依序為 SL, SW, PL, PW, species

2. 請使用 SAS 方法與 PROC SQL 方法將整組 data set 依照 species 切割成 3 份新的 data sets,名稱各為 setosa, versicolor, virginica

3. 請使用 SAS 方法與 PROC SQL 方法將 2. 的 3 組 data set 合併成新的 data set,名稱為 iris2

4. 請自行創造一筆新資料,各以 SAS 方法及 PROC SQL 方法增加到 iris2 資料集中


Oct 24 2011

SAS 考古題

 請從以下網頁下載資料檔

http://steve-chen.net/sas/iris.txt

資料檔內的變數及位置如下:

x1 : 第 1 格~ 第 3 格

x2:第 4 格 ~ 第 6 格

x3:第 7 格 ~ 第 9 格

剩下兩個變數為 x4 跟 species

(1) 請使用「至少」兩種輸入法來讀入資料,並僅從第 4 列讀到第 110 列

(2) 請依照 x2 然後 x3 的順序從小到大排序,並使用 PUT 指令,將所有變數內容及 FIRST.x2 , LAST.x2 , FIRST.x3 , LAST.x3 , 及 INFILE 尾端所指定的 COLUMN 變數內容寫到 LOG 視窗


May 26 2011

進學班 R 軟體期末上機考考古題

1. 有 10 個資料檔,名稱分別為 d001.txt, d002.txt, … d010.txt, 資料檔裡面都有相同的兩個變數 IQ 與 scores。請用迴圈方式讀入 R 程式,每個資料檔畫出一個 IQ vs. scores vs.  的 X-Y 散佈圖

2. 請用 R 程式作出 t 分配(自由度 30, 31, …, 48,49, 50) 的查表表格(請參考一般統計學後面 t 分配表格的格式與配置)

3. 請將 2 寫成一個 function,function 中必須有內部的 function,function 的參數也必須有預設值

4. 請畫出 beta(1.0,0.5) 分配的機率密度圖 (可用 plot 函數或 Chapter 12 中的 curve 函數來畫)


Nov 9 2010

「R軟體:應用統計方法」修訂版

我寫的「R軟體:應用統計方法」由於第一版存貨已經賣完,目前修訂版已經出版了,有興趣的朋友們,可以在各網路書店買到,例如 博客來書局網頁

這次修訂版僅是微幅修正一些打字排版錯誤之處,實際內容跟初版相差極小,但東華書局重新幫我設計了封面、並重新排版。其實我還是比較喜歡第一版的封面,因為是我自己設計的 :)

之前購買第一版的朋友們,麻煩下載以下的堪誤表 ( 此檔副檔名為 .docx,若下載後變成 .zip 檔,請直接更改為 .docx ):

http://steve-chen.net/wp-content/uploads/2010/11/correction.docx

很感謝大家對這本書的支持,也對書局那邊往往無法及時供貨在此道歉。我會盡量請東華書局那邊出書能夠更為順暢。


Jul 1 2010

期末收到的幾封 emails

每個學期末的期末考完之後,我都會收到一堆關於成績問題的 emails。這些類似的 emails 每年每個學期都會出現,我想藉這篇文章來談談這些問題。

Continue reading


Jun 30 2010

時間數列期末報告成績

時間數列:期末報告成績網址如下:

http://steve-chen.net/TS/tsreport.html

註:
1. 以上網址僅列出期末報告成績,並非學期成績
2. 本文一週後自動刪除