版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
R軟件基本操作R軟件和軟件包安裝數(shù)學(xué)運(yùn)算向量運(yùn)算矩陣運(yùn)算列表、數(shù)據(jù)框概率函數(shù)軟件安裝(Windows)當(dāng)前版本3.2.2。至,點(diǎn)擊左邊download下鏈接CRAN
,任意選擇一個(gè)鏡像站點(diǎn),比如USA下http://r.B
。在DownloadandInstallR中選擇DownloadRforWindows
,點(diǎn)擊base,選擇文件DownloadR3.2.2forWindows
,點(diǎn)擊鼠標(biāo)右鍵另存為保存到指定目錄中。下載好后安裝R。軟件包安裝比如安裝“BRugs”,可以使用命令>install.packages("BRugs",repos="http://r.B",dep=TRUE)調(diào)用軟件包
>require(BRugs)#>library(BRugs)升級軟件包>update.packages(repos="http://r.B",dep=TRUE)刪除軟件包
>remove.packages("BRugs",lib=file.path("path","to","library"))命令的執(zhí)行主窗口REditor窗口(File-Newscript)require(Rcmdr)常用的小技巧查看工作目錄:>getwd()[1]"D:/ProgramFiles/R/R-2.11.1”設(shè)定工作目錄:>setwd("d:/ProgramFiles/R")>getwd()[1]"d:/ProgramFiles/R”取得幫助:>help.start()#打開幫助文檔>help.search('paneldata')>RSiteSearch('KalmanFilter')>?lm查看對象類型:>typeof()>class()單引號和雙引號的區(qū)別:>c("single'qutoe",'double"quote',"double\"quote",'single\'quote')[1]"single'qutoe""double\"quote""double\"quote""single'quote"更改選項(xiàng)設(shè)定>
options(prompt='R>',continue='R')運(yùn)算符號+-*/^<<=>>===!=%%余數(shù)%/%商%in%判斷>17%in%1:100[1]TRUE基本數(shù)學(xué)運(yùn)算加減乘除乘方運(yùn)算>5+6*7/12^(1/3)[1]23.34519指數(shù)對數(shù)運(yùn)算>exp(2)#e^2[1]7.389056>log(2,10)#lg(2)[1]0.30103>log(exp(1))#ln(e)[1]1基本數(shù)學(xué)運(yùn)算(續(xù))求導(dǎo)>D(expression(x^2),"x")2*x>D(expression(sin(x^2)),"x")cos(x^2)*(2*x)積分:利用函數(shù)integrate()>integrate(dnorm,-1.96,1.96)$value[1]0.9500042>integrate(dnorm,-Inf,Inf)$value[1]1#對于多重積分計(jì)算,利用adapt庫中函數(shù)adapt,不允許用無窮界限。f<-function(x){1/sqrt(2*pi)*exp(-x^2/2)}integrate(f,-1.96,1.96)$value復(fù)數(shù)運(yùn)算‘Re’-實(shí)部‘Im’-虛部‘Mod’-?!瓵rg’-角度‘Conj’-共軛
r=Mod(z)=√(x^2+y^2),φ=Arg(z),x=r*cos(φ)
y=r*sin(φ)
解方程(組)>polyroot(c(-3,2,1))[1]1+0i-3-0i>A<-matrix(c(2,3,1,1),ncol=2,byrow=T)>A[,1][,2][1,]23[2,]11>B<-c(1,-1)>solve(A,B)[1]-43x^2+2x-3=02x+3y=1x+y=-1集合運(yùn)算并>union(x,y)交>intersect(x,y)補(bǔ)>setdiff(x,y)>x<-1:20>y<-3:15>intersect(x,y)[1]3456789101112131415>union(x,y)[1]1234567891011121314151617181920>match(y,x)[1]3456789101112131415>z<-c(23,15,12)>match(y,z)[1]NANA
NA
NA
NA
NA
NA
NA
NA3NANA2>setdiff(x,y)#從x中消去y[1]121617181920排列組合choose(n,k)-組合數(shù)combn(n,k)-列出所有組合factorial(n)-階乘require(e1071)>permutations(3)[,1][,2][,3][1,]123[2,]213[3,]231[4,]132[5,]312[6,]321序列和向量
z<-seq(-1,10,length=100)#z<-seq(-1,10,len=100)z<-seq(10,-1,-1)#z<-10:-1>x<-rep(3,1:3)Errorinrep(3,1:3):invalid'times'argumentx<-rep(3:5,1:3)>x[1]344555x<-rep(c(1,10),c(4,5))w<-c(1,3,x,z);w[3]向量運(yùn)算
x<-rep(0,10);z<-1:3;x+z
[1]1231231231Warningmessage:longerobjectlengthisnotamultipleofshorterobjectlengthin:x+z
x*z[1]0000000000Warningmessage:longerobjectlengthisnotamultipleofshorterobjectlengthin:x*zrev(x)z<-
c("no
cat","has","nine","tails")z[1]=="nocat"[1]TRUE向量名和元素附加
x<-1:3;names(x)<-LETTERS[1:3]xABC123set.seed(10);append(x,runif(3),after=2)
ABC1.00002.00000.50750.30680.42693.0000向量賦值
z<-1:5z[7]<-8;z[1]12345NA8z<-NULLz[c(1,3,5)]<-1:3;z[1]1NA2NA3rnorm(10)[c(2,5)]z[-c(1,3)]#去掉第1、3元素.z[(length(z)-4):length(z)]#最后五個(gè)元素.向量的大小次序
set.seed(10);z<-sample(1:100,10);z
#比較sample(1:100,10,rep=T)
[1]5131426892226945740
order(z)#給出數(shù)據(jù)從小到大的位置,即最小的是第5個(gè)
[1]56721031948z[order(z)]
[1]9222631404251576894sort(z)
[1]9222631404251576894which(z==max(z))#給出下標(biāo)[1]8矩陣操作>(A<-matrix(1:12,3,4))[,1][,2][,3][,4][1,]14710[2,]25811[3,]36912>(A<-matrix(1:12,3,4,byrow=T))
[,1][,2][,3][,4][1,]1234[2,]5678[3,]9101112>dim(A)[1]34矩陣操作(續(xù))>nrow(A)[1]3>dim(A)[1][1]3>ncol(A)[1]4>dim(A)[2][1]4>A[3,][1]9101112>A[,3][1]3711>A[3,3][1]11矩陣操作(續(xù))>A[-1,]#去除A的第一行
[,1][,2][,3][,4][1,]5678[2,]9101112>A[,-2]#去除A的第二列[,1][,2][,3][1,]134[2,]578[3,]91112>(AA<-matrix(1:9,3,3))[,1][,2][,3][1,]147[2,]258[3,]369>AA[-2,-2]#刪除AA矩陣的第二行第二列
[,1][,2][1,]17[2,]39矩陣操作(續(xù))>A>3#判斷A中元素是否大于3
[,1][,2][,3][,4][1,]FALSEFALSEFALSETRUE[2,]TRUETRUETRUETRUE[3,]TRUETRUE
TRUE
TRUE>A==3
[,1][,2][,3][,4][1,]FALSEFALSETRUEFALSE[2,]FALSEFALSEFALSEFALSE[3,]FALSEFALSEFALSEFALSE>A[A>6][1]910711812>length(A[A>6])[1]6矩陣操作(續(xù))>(A<-matrix(0,3,3))[,1][,2][,3][1,]000[2,]000[3,]000>diag(A)<-1>diag(A)[1]111>diag(1,3,3)
[,1][,2][,3][1,]100[2,]010[3,]001矩陣操作(續(xù))>(m2<-matrix(1:20,4,5))[,1][,2][,3][,4][,5][1,]1591317[2,]26101418[3,]37111519[4,]48121620>lower.tri(m2)[,1][,2][,3][,4][,5][1,]FALSEFALSEFALSEFALSEFALSE[2,]TRUEFALSEFALSEFALSEFALSE[3,]TRUETRUEFALSEFALSEFALSE[4,]TRUETRUETRUEFALSEFALSE>m2[lower.tri(m2)]<-NA>m2
[,1][,2][,3][,4][,5][1,]1591317[2,]NA6101418[3,]NANA111519[4,]NANA
NA1620矩陣加減>(A<-matrix(1:12,3,4))[,1][,2][,3][,4][1,]14710[2,]25811[3,]36912>(B<-matrix(-1:-12,3,4))[,1][,2][,3][,4][1,]-1-4-7-10[2,]-2-5-8-11[3,]-3-6-9-12>(C1<-A+B)[,1][,2][,3][,4][1,]0000[2,]0000[3,]0000>(D1<-A-B)[,1][,2][,3][,4][1,]281420[2,]4101622[3,]6121824矩陣數(shù)乘>(A<-matrix((1:9)^2,3,3))[,1][,2][,3][1,]11649[2,]42564[3,]93681>(AmultipliedByTwo<-2*A)[,1][,2][,3][1,]23298[2,]850128[3,]1872162>(AmultipliedByTwo/2==A)[,1][,2][,3][1,]TRUETRUE
TRUE[2,]TRUETRUETRUE[3,]TRUETRUETRUE矩陣相乘>(E<-matrix(1:9,3,3));(E2<-matrix(1:9,3,3,byrow=T))[,1][,2][,3][1,]147[2,]258[3,]369[,1][,2][,3][1,]123[2,]456[3,]789>E2%*%E[,1][,2][,3][1,]143250[2,]3277122[3,]50122194矩陣逆>(B<-matrix((1:9)^2,3,3))[,1][,2][,3][1,]11649[2,]42564[3,]93681>(Binverse<-solve(B))
[,1][,2][,3][1,]1.291667-2.1666670.9305556[2,]-1.1666671.666667-0.6111111[3,]0.375000-0.5000000.1805556>B%*%Binverse[,1][,2][,3][1,]1.000000e+00-8.881784e-169.714451e-16[2,]-3.552714e-151.000000e+000.000000e+00[3,]6.661338e-16-8.881784e-161.000000e+00>round(B%*%Binverse)[,1][,2][,3][1,]100[2,]010[3,]001矩陣分解因子分解(譜(spectral)分解)>J<-cbind(c(20,3),c(3,18))>j<-eigen(J,symmetric=T)>j$values[1]22.1622815.83772$vectors[,1][,2][1,]-0.81124220.5847103[2,]-0.5847103-0.8112422>C<-j$vec>D<-diag(j$values)>C%*%D%*%t(C)[,1][,2][1,]203[2,]318
C是特征向量構(gòu)成的列矩陣,是特征值構(gòu)成的對角陣。矩陣分解(續(xù))奇異值(singularvalue)分解
>H<-matrix(c(3,2,1,1,1,2),3,2)>(HHt.eigen=eigen(H%*%t(H)))$values[1]1.806226e+011.937742e+002.153937e-15$vectors[,1][,2][,3][1,]0.7286301-0.460386080.5070926[2,]0.5251205-0.09981375-0.8451543[3,]0.43971210.882089490.1690309>U=HHt.eigen$vec[,1:2]>D=diag(sqrt(HHt.eigen$val)[1:2])>V=t(H)%*%U%*%solve(D)>U%*%D%*%t(V)或>svd(H)的特征向量矩陣為,的對角元素是的非零特征值的平方根
矩陣分解(續(xù))Cholesky分解,實(shí)對稱正定矩陣的分解>A<-matrix(c(1,2,2,5),ncol=2)>A>H<-chol(A)>H>t(H)%*%H為非奇異陣,上三角陣列表list可以是任何對象的集合(包括lists)z<-list(1:3,Tom=c(1:2,a=list("R",letters[1:5]),w="hi!"))z[[1]];z[[2]];z$T;z$T$a2;z$T[[3]];z$T$wattributes(z)#屬性!$names[1]"""Tom"attributes(matrix(1:6,2,3))$dim[1]23矩陣,數(shù)組及其維名字
x<-matrix(1:12,nrow=3,dimnames=list(c("I","II","III"),paste("X",1:4,sep="")))
X1X2X3X4I14710II25811III36912y<-array(1:12,c(3,2,2),dimnames=list(c("I","II","III"),paste("X",1:2,sep=""),paste("Y",1:2,sep=""))),,Y1X1X2I14II25III36,,Y2X1X2I710II811III912數(shù)據(jù)框
x<-matrix(1:6,2,3)x<-
as.data.frame(x);xV1V2V311352246x$V2[1]34x$V2[1]34attributes(x)$names[1]"V1""V2""V3"$s[1]"1""2"$class[1]"data.frame"數(shù)據(jù)框(續(xù)
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 單位管理制度呈現(xiàn)合集人事管理篇十篇
- 《證券交易流程》課件
- 《企業(yè)戰(zhàn)略管理》課件
- 新生引航共筑未來
- 學(xué)校三年級班主任工作總結(jié)5篇
- 2023年-2024年新員工入職安全教育培訓(xùn)試題附答案(突破訓(xùn)練)
- 大學(xué)畢業(yè)晚會策劃書合集15篇
- 2023年-2024年新入職員工安全教育培訓(xùn)試題附下載答案可打印
- 2024員工三級安全培訓(xùn)考試題(原創(chuàng)題)
- 保護(hù)環(huán)境的建議書(合集15篇)
- DS12C887電子時(shí)鐘(附程序)
- 新產(chǎn)品試制流程管理辦法
- 王牌電話交換機(jī)說明書
- 列管式換熱器-換熱面積計(jì)算
- 10個(gè)地基基礎(chǔ)工程質(zhì)量通病及防治措施
- 25m預(yù)應(yīng)力混凝土簡支T梁橋設(shè)計(jì)(共30頁)
- 籃球校本課程教案
- 高一學(xué)生文理分班意向表
- 高等傳熱學(xué)部分答案
- 地球物理學(xué)進(jìn)展投稿須知
- 機(jī)床精度檢驗(yàn)標(biāo)準(zhǔn) VDI3441 a ISO230-2
評論
0/150
提交評論