统计程式语言.ppt

載入套件【方法二】 套件查詢 Currently, the CRAN package repository features 5816 available packages. Package 查詢：/web/packages/ Package “leaps”： /web/packages/leaps/index.html 套件查詢 變數命名及語法 例子：定義一個新的變數叫temp，並且指定其值為24 紅字是輸入，藍字是輸出 大小寫代表不同物件 等號可用箭頭(-)代替 變數名稱必須為字母或是點(.)開頭 變數名稱由字母、數字、點(.)或下底線( _ )組成，不能含有特殊符號，如：＄＆！。 正確的命名：x5; x_5; x.5; .x5; x.y.z.5 錯誤的命名：x-5; 5x; x$5 c, q, t, C, D, F, I, T; diff, df, pt, pi, …不可以是物件名稱 語法 程式以分號（;）或是換行做分隔。 註解：R不執行#後之內容，故常將程式註解置於#之後。 Recall : 上下鍵。 基本運算： 可以直接輸入數值計算。 若物件屬性與維度相同，可以進行四則運算。 若不能計算，會回傳NaN (not a number)，如：log(-1)。 */43 R as a calculator 可以直接輸入數值計算 若物件屬性與維度相同，可以進行四則運算 常用的數學函數，如log, exp, sin,..可直接計算 若不能計算，會回傳NaN (not a number) * 查詢、移除已定義物件 ls() rm(“A”) rm(list=ls()) * Environment commands — ls() rm() 紅色字代表輸入的指令、藍色字代表輸出的結果。 R的註解符號為「#」，當程式遇到「#」，後面的文字將不會執行 查詢已定義物件：尚無物件 查詢已存在的package：共九個package 定義新物件(定義新變數為Temp01，並指定其值為9)：大小寫代表不同變數 查詢已定義物件：兩個物件temp01 Temp01 移除已定義物件：移除temp01 Practice 請用R計算 請用R計算 請比較log(10); log2(10); log10(10) * Hint: 可用x!=gamma(x+1)或是x!=prod(1:x) 輸入向量 * c(…) 可以是數值 可以是字串 還可以把每一個元素取名字 Practice Weight: 50, 62, 47, 48, 72, 81, 45, 58, 65, 55 Height: 1.58, 1.77, 1.43, 1.52, 1.84, 1.72, 1.45, 1.62, 1.70, 1.55 請計算每個人的BMI及平均BMI * 將物件存檔 save.image(C:\\Users\\user1\\Desktop\\統計軟體\\temp.RData) * Input a sequence 用“:”表示從1到10 等距的數列可以用seq (sequence) 重複生成相同的數/數列，用rep (replication) * Practice: 如何生出 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 這個數列? Q:如何生出 10 9 8 7 6 5 4 3 2 1 呢? 輸入矩陣 * 直接輸入矩陣的元素，再指定其維度 利用matrix這個function。 1:12為矩陣內元素 nrow=number of row Note: ncol=number of column 矩陣運算 維度相同才能加減 x*y: x，y相同位置的元素直接相乘 x%*%t(y): 需維度符合才可相乘/除 t(y): transpose of y * 其他相關函數 cbind/rbind solve(A) rownames/colnames is.matrix/as.matrix, is.numeric/as.numeric, … colSums/rowSums, colMeans/rowMeans eigen(x) * Operations Practice 輸入矩陣 計算 計算 * Operations 定義變數 Function Example Output 數字運算 1+1 2 序列 -2:2 -2 -1 0 1 2 c(1:5) 1 2 3 4 5 rep(1, 5) 1 1 1 1 1 rep(1:2, 3)