r程序設(shè)計(jì)ppt課件-在線瀏覽
2025-06-22 18:22本頁(yè)面
【正文】 n個(gè)人的生日可取值數(shù)為 365n ， 而 n個(gè)人彼此不同的可能數(shù)為 365中取 n個(gè)的排列數(shù) ， 彼此不同的概率為 P3n65 。我們可以盡量用向量運(yùn)算來(lái)實(shí) 現(xiàn)，速度要快得多： ? x numeric(365) for(n in 1:365) { x[n] =1 prod((365:(365n+1))/365) } 這段程序只用了 1秒。（ choose(n, k) 、 factorial(k) ） ?另外要注意使用 for(i in 1:n)格式的計(jì)數(shù)循環(huán)時(shí)要避免一個(gè)常見(jiàn)錯(cuò)誤 ， 即當(dāng) n為零或負(fù)數(shù)時(shí) 1:n是一個(gè)從大到小的循環(huán) ， 而我們經(jīng)常需要的是當(dāng) n為零或負(fù)數(shù)時(shí)就不進(jìn)入循環(huán) 。 ?while循環(huán)是在開(kāi)始處判斷循環(huán)條件的當(dāng)型循環(huán)，如： ? While (baeps) { c (a+b)/2。 還可以使用 ?repeat 表達(dá)式 循環(huán)，在循環(huán)體內(nèi)用 break跳出。 如 i=1。repeat{x=x+rnorm(1)。i=i+1}。 利用二分法求方程的根： fzero(f,1,2,1e6) fzerofunction(f, a, b, eps=1e5){ if (f(a)*f(b)0) list (fail=“finding root is fail!) else { repeat { if (abs(ba)eps) break x(a+b)/2 if (f(a)*f(b)0) bx else ax } list(root=(a+b)/2, fun=f(x)) } } ffunction(x) {x^3x1} fzero(f,1,2,1e6) 實(shí)際上： uniroot(f,c(1,2)) $root [1] $ [1] $iter [1] 7 $ [1] 控制語(yǔ)句的幾個(gè)簡(jiǎn)單示例 ? x=NULL。x=NULL。i=i+1} ? x=rnorm(1)。x=rnorm(1)。if(x3)break。print(c(i,x)) 怎么調(diào)出 Packages來(lái)使用 ? Packages:libraries 敲 library(),就知道有什么 libraries了 ,缺省 library是 mass,就敲 library(MASS) 每個(gè) library都有許多數(shù)據(jù) 在每個(gè) library敲 data(),就知道有什么數(shù)據(jù)了 比如敲 data(Titanic),就調(diào)出數(shù)據(jù) Titanic來(lái)了 . 注意 :R語(yǔ)言對(duì)大小寫(xiě)敏感 . R程序設(shè)計(jì) 對(duì)于復(fù)雜一些的計(jì)算問(wèn)題應(yīng)該編寫(xiě)成函數(shù)。 工作空間管理 ?前面我們已經(jīng)提到 ， R在運(yùn)行時(shí)保持一個(gè)變量搜索路徑表 ， 要讀取某變量時(shí)依次在此路徑表中查找 ，返回找到的第一個(gè)；給變量賦值時(shí)在搜索路徑的第一個(gè)位置賦值 。 ?用 ls()函數(shù)可以查看當(dāng)前工作空間保存的變量和函數(shù)，用 rm()函數(shù)可以剔除不想要的對(duì)象。定義函數(shù)可以在命令行進(jìn)行，例如 ? hello = function(){ cat(Hello, are you a good student?\n) 。} ? hello （查看函數(shù)具體內(nèi)容） function(){ cat(Hello, are you a good student?\n) 。} ? hello() （運(yùn)行函數(shù)） Hello, are you a good student? ?函數(shù)體為一個(gè)復(fù)合表達(dá)式，各表達(dá)式的之間用換行或分號(hào)分開(kāi)。 ?在命令行輸入函數(shù)程序很不方便修改，所以我們一般是打開(kāi)一個(gè)其他的編輯程序（如 Windows 的記事本），輸入以上函數(shù)定義，保存文件，比如保存到了 E:\R\，我們就可以用 ? source(“E:/R/) ?運(yùn)行文件中的程序。 參數(shù)（自變量） ?函數(shù)可以帶虛參數(shù)（形式自變量）。即調(diào)用時(shí)實(shí)參與虛參可以按次序結(jié)合，也可以直接指定虛參名結(jié)合。 ?函數(shù)在調(diào)用時(shí)可以不給出所有的實(shí)參， 這需要在定義時(shí)為虛參指定缺省值 。 作用域 ?函數(shù)的虛參完全是 按值傳遞 的，改變虛參的值不能改變對(duì)應(yīng)實(shí)參的值。x} f(x) x ? 函數(shù)體內(nèi)的變量也是局部的，對(duì)函數(shù)體內(nèi)的變量賦值當(dāng)函數(shù)結(jié)束運(yùn)行后變量值就刪除了，不影響原來(lái)同名變量的值。 x= 20。但是也要注意，函數(shù)中的顯示函數(shù)調(diào)用時(shí)局部變量 x還沒(méi)有賦值，顯示的是全局變量 x 的值。 程序調(diào)試 ?SPLUS和 R目前還不象其它主流程序設(shè)計(jì)語(yǔ)言那樣具有單步跟蹤 、 設(shè)置斷點(diǎn) 、 觀察表達(dá)式等強(qiáng)勁的調(diào)試功能 。 ?對(duì)任何程序語(yǔ)言 ， 最基本的調(diào)試手段當(dāng)然是在需要的地方顯示變量的值 。 例如 ， 我們?yōu)榱苏{(diào)試前面定義的 larger()函數(shù) ， 可以顯示兩個(gè)自變量的值及中間變量的值： ? larger function(x, y){ cat(39。, x, 39。) 。y = 39。\n39。 =(yx)。 = 39。\n39。x[] = y[] 。例如： ? larger= function(x, y){ (yx)。x[] =y[]。} ? 我們運(yùn)行此程序： larger(c(1,3,5),c(2, 4, 9)) ? 退出 R的 browser()菜單可用 c（在 SPlus中用return()）。 ? R提供了一個(gè) debug()函數(shù)， debug(f)可以打開(kāi)對(duì)函數(shù) f()的調(diào)試，執(zhí)行到函數(shù) f時(shí)自動(dòng)進(jìn)入單步執(zhí)行的 browser()菜單。 ? 試試 larger(c(1,5),c(2, 4, 9)) ? Warning message:長(zhǎng)的目標(biāo)對(duì)象長(zhǎng)度不是短的目標(biāo)對(duì)象長(zhǎng)度的整倍數(shù) Called from: larger(c(1, 5), c(2, 4, 9)) Browse[1] y [1] 2 4 9 Browse[1] x [1] 1 5 Browse[1] yx Warning in y x : longer object length is not a multiple of shorter object length [1] TRUE FALSE TRUE Browse[1] c Error: subscript (3) out of bounds, should be at most 2 程序設(shè)計(jì)舉例 ?設(shè)計(jì) R程序是很容易的 ， 在初學(xué)時(shí)我們只要使用我們從一般程序設(shè)計(jì)中學(xué)來(lái)的知識(shí)并充分利用 R中現(xiàn)成的各種算法及繪圖函數(shù)就可以了 。 ?用 R語(yǔ)言開(kāi)發(fā)算法 ， 最重要的一點(diǎn)是要記住R是一個(gè)向量語(yǔ)言 ， 計(jì)算應(yīng)該盡量通過(guò)向量 、矩陣運(yùn)算來(lái)進(jìn)行 ，
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1