R 程式設計/示例會話
< R 程式設計
本頁是 R 程式語言的簡介。它展示瞭如何使用 R 執行非常簡單的任務。首先，您需要安裝 R（請參閱設定頁面）。如果您使用 Windows 或 Mac OS，最簡單的解決方案是使用 R 圖形使用者介面（單擊其圖示）。如果您使用 Linux，請開啟終端並在命令提示符下鍵入 R。
通常，當您開啟 R 時，您會在控制檯中看到類似於以下的訊息
R version 3.5.1 (2018-07-02) -- "Feather Spray"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
Natural language support but running in an English locale
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
[Workspace loaded from ~/.RData]
>
您可以在尖括號 > 後鍵入您的程式碼。
R 可以用作一個簡單的計算器，我們可以執行任何簡單的計算。
> # Sample Session
> # This is a comment
>
> 2 # print a number
[1] 2
> 2+3 # perform a simple calculation
[1] 5
> log(2) # natural log
[1] 0.6931472
我們還可以使用賦值運算子 <- 儲存數字或字串物件。
> x <- 2 # store an object
> x # print this object
[1] 2
> (x <- 3) # store and print an object
[1] 3
>
> x <- "Hello" # store a string object
> x
[1] "Hello"
我們還可以儲存向量.
> Height <- c(168, 177, 177, 177, 178, 172, 165, 171, 178, 170) #store a vector
> Height # print the vector
[1] 168 177 177 177 178 172 165 171 178 170
>
> Height[2] # Print the second component
[1] 177
> Height[2:5] # Print the second, the 3rd, the 4th and 5th component
[1] 177 177 177 178
>
> (obs <- 1:10) # Define a vector as a sequence (1 to 10)
[1] 1 2 3 4 5 6 7 8 9 10
>
> Weight <- c(88, 72, 85, 52, 71, 69, 61, 61, 51, 75)
>
> BMI <- Weight/((Height/100)^2) # Performs a simple calculation using vectors
> BMI
[1] 31.17914 22.98190 27.13141 16.59804 22.40879 23.32342 22.40588 20.86112
[9] 16.09645 25.95156
我們還可以用以下方法描述向量length(), mean()和var().
> length(Height)
[1] 10
> mean(Height) # Compute the sample mean
[1] 173.3
> var(Height)
[1] 22.23333
我們還可以定義一個矩陣.
> M <- cbind(obs,Height,Weight,BMI) # Create a matrix
> typeof(M) # Give the type of the matrix
[1] "double"
> class(M) # Give the class of an object
[1] "matrix"
> is.matrix(M) # Check if M is a matrix
[1] TRUE
> is.vector(M) # M is not a vector
[1] FALSE
> dim(M) # Dimensions of a matrix
[1] 10 4
我們可以使用以下方法繪製資料plot().
> plot(Height,Weight,ylab="Weight",xlab="Height",main="Corpulence")
我們可以定義一個資料框.
> mydat <- data.frame(M) # Creates a dataframe
> names(mydat) # Give the names of each variable
[1] "obs" "Height" "Weight" "BMI"
> str(mydat) # give the structure of your data
'data.frame': 10 obs. of 4 variables:
$ obs : num 1 2 3 4 5 6 7 8 9 10
$ Height: num 168 177 177 177 178 172 165 171 178 170
$ Weight: num 88 72 85 52 71 69 61 61 51 75
$ BMI : num 31.2 23 27.1 16.6 22.4 ...
>
> View(mydat) # Look at your data
>
> summary(mydat) # Descriptive Statistics
obs Height Weight BMI
Min. : 1.00 Min. :165.0 Min. :51.00 Min. :16.10
1st Qu.: 3.25 1st Qu.:170.2 1st Qu.:61.00 1st Qu.:21.25
Median : 5.50 Median :174.5 Median :70.00 Median :22.70
Mean : 5.50 Mean :173.3 Mean :68.50 Mean :22.89
3rd Qu.: 7.75 3rd Qu.:177.0 3rd Qu.:74.25 3rd Qu.:25.29
Max. :10.00 Max. :178.0 Max. :88.00 Max. :31.18
>
您可以儲存 R 會話（記憶體中的所有物件）並載入會話。
> save.image(file="~/Documents/Logiciels/R/test.rda")
> load("~/Documents/Logiciels/R/test.rda")
我們可以定義一個工作目錄。注意 Windows 使用者：R 在目錄中使用斜槓（“/”）而不是反斜槓（“\”）。
> setwd("~/Desktop") # Sets working directory (character string enclosed in "...")
> getwd() # Returns current working directory
[1] "/Users/username/Desktop"
> dir() * Lists the content of the working directory
R 中有一些特殊字元
- NA : 不可用（即缺失值）
- NaN : 非數字（例如 0/0）
- Inf: 無限大
- -Inf : 負無窮大。
例如，0 除以 0 會得到一個NaN但 1 除以 0 會得到
> 0/0
[1] NaN
> 1/0
[1] Inf
我們可以使用以下方法退出 Rq(). 這no引數指定不儲存 R 會話。
q("no")