R Programming - Note

Week 1

http://www.reibang.com/p/1be2821f3929

Week 2

1 控制結(jié)構(gòu)語句(Control Structures)

1.1 If-else

If語句有兩種形式菱鸥，其中else語句不是必須的

if (<condition>) {
          ## do something
} else {
          ## do something else
}


if (<condition>) {
          ## do something
} else if (<condition>) {
          ## do something different
} else {
          ## do something different
}

可以直接給變量用if-else語句賦值

y <- if(x > 3) {
     10
} else {
     0
}

1.2 For

基本范式

for(i in 1: 10) {
    print(i)
}

seq_along()返回和對象等長的數(shù)列

for(i in 1:10) {
    print(x[i])
}
for(i in seq_along(x)) {
    print(x[i])
}
for(letter in x ) {
    print(letter)
}
for(i in seq_len(nrow(x))) {
    for(j in seq_len(ncol(x))) {
        print(x[i, j])
    }
}

嵌套for循環(huán)（最好不要超過三層）

x <- matrix(1:6, 2, 3)
for (i in seq_len(nrow(x)) {
    for (j in seq_len(ncol(x))) {
        print(x[i, j])
    }
}

1.3 While

基本范式

count <- 0
while(count < 10) {
      print(count)
      count = count + 1
}

rbinom()正態(tài)分布隨機(jī)數(shù)的生成函數(shù)：rnorm（n,mean=0,sd=1）其中n表示生成的隨機(jī)數(shù)數(shù)量，mean是正態(tài)分布的均值非剃，默認(rèn)為0始衅，sd是正態(tài)分布的標(biāo)準(zhǔn)差，默認(rèn)時(shí)為1

z <- 5
while (z >=3 && z <=10) {
    print(z)
    coin <- rbinom(1, 1, 0.5)
    if (coin == 1) { ##random walk
        z <- z+1
    } else {
        z <- z-1
    }
}

1.4 Repeat

退出一個(gè)重復(fù)循環(huán)唯一的方式是調(diào)用break

1.5 Next

next用來跳過一些迭代

for (i in 1:100) {
    if (i <=20) {
        ##skip the first 20 iterations
        next
    }
    ## do something here
}

1.6 Return

退出循環(huán)并返回一個(gè)值

2 編寫R函數(shù)(Functions)

所有R函數(shù)會自動返回最后一個(gè)表達(dá)式的值
參數(shù)定義時(shí)可以設(shè)定缺省值

above <- function(x, n = 10) { ##默認(rèn)參數(shù)等于10
    use <- x >n
    x[use]
}
columnmean <- function(y, removeNA = TRUE) {
    nc <- ncol(y) ## 數(shù)據(jù)類型integer
    means <- numeric(nc) ## 轉(zhuǎn)化為numeric
    for (i in 1:nc) {
        means[i] <- mean(y[,i], na.rm = removeNA)
    }
    means
}

2.1 參數(shù)(Arguments)

2.1.1 基本方法

形式參數(shù)是包含在函數(shù)定義里的參數(shù)，formals()查看或設(shè)置函數(shù)的形式參數(shù)
args()查看參數(shù)名稱和參數(shù)缺省值
可通過命名參數(shù)來快速調(diào)用參數(shù)
R函數(shù)檢查參數(shù)匹配的順序：名稱完全匹配、名稱部分匹配侥祭、位置匹配

2.1.2 "..."參數(shù)

用來表明一些可以傳遞給另一個(gè)函數(shù)的參數(shù)叁执，不用重復(fù)創(chuàng)建參數(shù)（另一用法是用在泛型函數(shù)中）

myplot <- function(x, y, type = "1", ...) {
       plot(x, y, type = type, ...)
}

無法確定參數(shù)數(shù)量的時(shí)候也用"..."

> args(paste)
function (..., sep = " ", collapse = NULL) 
> args(cat)
function (..., file = "", sep = " ", fill = FALSE, labels = NULL, 
    append = FALSE)

任何出現(xiàn)在"..."之后的參數(shù)需要明確地給出名稱，而且不能進(jìn)行部分匹配

2.1.2 惰性求值(Lazy Evaluation)

即僅在使用這個(gè)參數(shù)時(shí)求值

> f <- function(a, b) {
    a^2
}
> f(2)
[1] 4
> f <- function(a, b) {
    print(a)
    print(b)
}
> f(45)
[1] 45
Error in print(b) : argument "b" is missing, with no default

3 作用域規(guī)則(Scoping rules)

搜索順序全局環(huán)境->包

make.power <- function(n) {
    pow <- function(x) {
        x^n
    }
    pow
}
> cube <- make.power(3)
> square <- make.power(2)
> cube(3)
[1] 27
> square(3)
[1] 9

> ls(environment(cube))
[1] "n" "pow"
> get("n", environment(cube))
[1] 3

y <- 10
f <- function(x) {
    y <- 2
    y^2 + g(x)
}
g <- function(x) {
    x*y
}
> f(3)
[1] 34

作用域規(guī)則在統(tǒng)計(jì)學(xué)中的優(yōu)化也很有用（暫略矮冬，具體看課程視頻）

4 時(shí)間與日期

R里面存儲日期的數(shù)據(jù)類型是Date

> x <- as.Date("1970-01-01")
> x
[1] "1970-01-01"
> class(x)
[1] "Date"
> unclass(x)
[1] 0
> unclass(as.Date("1970-01-02"))
[1] 1

R里面存儲時(shí)間的數(shù)據(jù)類型是POSIXct和POSIXIt

> x <- Sys.time()
> x
[1] "2017-11-18 15:10:32 CST"
> p <- as.POSIXlt(x)
> names(unclass(p))
 [1] "sec"    "min"    "hour"  
 [4] "mday"   "mon"    "year"  
 [7] "wday"   "yday"   "isdst" 
[10] "zone"   "gmtoff"
> p$sec
[1] 32.1765

對時(shí)間和日期起作用的函數(shù)：weekdays()返回一周中的某一天谈宛、months()返回月份、quarters()返回季度（"Q1", "Q2", "Q3", "Q4"）
strptime()字符串轉(zhuǎn)換為日期與時(shí)間格式

最后編輯于：2019.09.26 16:32:01

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者