#運(yùn)算符:的返回值是一個(gè)向量峻村,也就是一個(gè)一位數(shù)組
1:6
#屏幕上出現(xiàn)的6個(gè)數(shù)字并沒(méi)有保存在計(jì)算機(jī)內(nèi)存里
#想要把他們保存下來(lái)需要建立R對(duì)象
#對(duì)象object
die <- 1:6
#如果想查看已經(jīng)命名了哪些R對(duì)象读存,可以使用ls命令
ls()
#Rstudio的環(huán)境面板也可以看到
#擲骰子
#骰子所有面的值隨機(jī)選擇一個(gè)
#需要建立一個(gè)函數(shù)(function)
#函數(shù)
#R自帶函數(shù)round可以實(shí)現(xiàn)數(shù)字的四舍五入
#factorial函數(shù)可以實(shí)現(xiàn)階乘操作
round(3.13145)
factorial(3)
#傳遞到函數(shù)中的數(shù)據(jù)稱為該函數(shù)的參數(shù)(argument)。
#參數(shù)可以是原始數(shù)據(jù)、R對(duì)象、另一個(gè)R函數(shù)的返回結(jié)果
round(3.5)
round(mean(1:6))
round(mean(die))
#R使用函數(shù)嵌套時(shí),從最內(nèi)層的運(yùn)算開(kāi)始解析褐着,知道最外層的運(yùn)算為止
#R首先找到對(duì)象die,然后算mean(1:6)托呕;最后四舍五入
#sample(x,size)從向量x中抽取size個(gè)元素并返回
sample(x = 1:4, size = 2)
#擲骰子
sample(x = die, size = 1)
#一個(gè)R函數(shù)可以設(shè)置多個(gè)參數(shù)含蓉,只要用逗號(hào)將它們隔開(kāi)
#對(duì)數(shù)據(jù)對(duì)象賦值給該函數(shù)的哪個(gè)參數(shù)
#將這個(gè)數(shù)據(jù)對(duì)象的名稱與參數(shù)名用等號(hào)連起來(lái)
#R用戶大都不會(huì)指定R函數(shù)的第一個(gè)參數(shù)名稱频敛。
sample(die, size = 1)
#args函數(shù)可以查看某函數(shù)的所有參數(shù)名
args(round)
#在調(diào)用一個(gè)包含多個(gè)參數(shù)的函數(shù)時(shí),從第二個(gè)參數(shù)或者第三個(gè)參數(shù)開(kāi)始馅扣,應(yīng)該寫(xiě)出每個(gè)參數(shù)的名稱
#可放回抽樣
#sample函數(shù)在抽樣時(shí)默認(rèn)使用了不可放回抽樣
sample(die, size = 2)
#額外設(shè)定參數(shù)replace = TRUE
#模擬一對(duì)骰子的正確方法
sample(die, size = 2, replace = TRUE)
#可放回抽樣是創(chuàng)建獨(dú)立隨機(jī)樣本的一種簡(jiǎn)單方法
#如果想要知道兩個(gè)骰子的總點(diǎn)數(shù)
#將模擬結(jié)果交給sum函數(shù)
#若有一個(gè)對(duì)象在每次調(diào)用時(shí)都重新擲一次骰子斟赚,需要編寫(xiě)一個(gè)自定義的R函數(shù)
#編寫(xiě)自定義函數(shù)
die <- 1:6
dice <- sample(die, size = 2, replace = TRUE)
sum(dice)
#我們要寫(xiě)一個(gè)函數(shù)roll(), 每次調(diào)用函數(shù)roll(),R都會(huì)返回兩個(gè)骰子的總和
#函數(shù)構(gòu)建器
#任何一個(gè)函數(shù)都包含三個(gè)部分:函數(shù)名差油、程序主體拗军、參數(shù)集合
my_function <- function(){}
# function的作用是將打括號(hào)中的所有代碼構(gòu)建一個(gè)函數(shù)
# 例如,若要將擲骰子的模擬構(gòu)成為一個(gè)函數(shù)
roll <- function(){
? die <- 1:6
? dice <- sample(die, size = 2, replace = TRUE)
? sum(dice)?
}
# 請(qǐng)注意大括號(hào)內(nèi)部的每一行代碼都做了縮緊操作蓄喇,R會(huì)忽略所有空格和換行符
#左大擴(kuò)號(hào)后发侵,在每一行結(jié)束時(shí)按回車鍵就可以調(diào)入下一行的編輯
#使用方法很簡(jiǎn)單,只需要鍵入這個(gè)對(duì)象的名稱妆偏,跟上一對(duì)括號(hào)即可
roll()
#如果不帶括號(hào)刃鳄,R會(huì)展示這個(gè)函數(shù)的代碼
#大括號(hào)內(nèi)的內(nèi)容叫函數(shù)的主體
#R會(huì)順序執(zhí)行該函數(shù)的主體中所有代碼,并且返回最后一行代碼的運(yùn)行結(jié)果
#如果最后一行代碼不返回任何值钱骂,整個(gè)函數(shù)也不會(huì)返回任何值
#因此在編寫(xiě)函數(shù)時(shí)應(yīng)該確保最后一行代碼具有明確的返回值
# 參數(shù)
# 如果函數(shù)中沒(méi)有定義對(duì)象bones叔锐,在函數(shù)調(diào)用過(guò)程中就會(huì)粗線錯(cuò)誤
# 但是如果將這個(gè)對(duì)象bones作為這個(gè)函數(shù)的參數(shù)名,就可以將bones作為參數(shù)傳遞給sample函數(shù)
roll2 <- function(bones) {
? dice <- sample(die, size = 2, replace = TRUE)
? sum(dice)?
}
#我們?cè)谑褂眯潞瘮?shù)roll2時(shí)见秽,在其后的括號(hào)總鍵入bones作為該函數(shù)的參數(shù)
#這樣我們可以擲不同類型的骰子
#如果不給roll2提供bones的值愉烙,那么roll2函數(shù)會(huì)報(bào)錯(cuò)
#避免這種情況,可以給一個(gè)初始默認(rèn)值給bones
#只需要在定義roll2函數(shù)時(shí)解取,將bones與相應(yīng)的默認(rèn)值用等號(hào)連接即可
roll2 <- function(bones = 1:6) {
? dice <- sample(die, size = 2, replace = TRUE)
? sum(dice)?
}
#可以在自定義函數(shù)中設(shè)置任意數(shù)量的參數(shù)步责,只要在function后的括號(hào)中列出名稱,用逗號(hào)隔開(kāi)
#讓這對(duì)骰子擲出大點(diǎn)數(shù)的概率稍高于擲出小點(diǎn)數(shù)的概率
#對(duì)骰子進(jìn)行加權(quán)
#首先要確保這對(duì)骰子原本是均勻的
#重復(fù)(repetition)和可視化(visualization)
library(ggplot2)
#使用c函數(shù)創(chuàng)建數(shù)值向量禀苦,每個(gè)數(shù)字用逗號(hào)隔開(kāi)
#散點(diǎn)圖:描述兩個(gè)變量之間的關(guān)系
#直方圖:可視化單一變量的分布情況
#檢查骰子有沒(méi)有被做過(guò)手腳
# replicate()可以快速重復(fù)運(yùn)行一段R命令
#replicate(重復(fù)運(yùn)行的次數(shù)蔓肯,想要重復(fù)運(yùn)行的R代碼)
replicate(3, 1 + 1)
replicate(10, roll())
#模擬投擲10000次骰子作圖
rolls <- replicate(10000, roll())
qplot(rolls,binwidth = 1)
#將兩個(gè)骰子出現(xiàn)點(diǎn)數(shù)為6的概率提高
#其實(shí)只需要在sample函數(shù)中添加一個(gè)新參數(shù)就可以實(shí)現(xiàn)加權(quán)
roll <- function(){
? die <- 1:6
? dice <- sample(die, size = 2, replace = TRUE,
? ? ? ? ? ? ? ? prob = c(1/8, 1/8, 1/8, 1/8, 1/8, 3/8))
? sum(dice)?
}