ch1-R軟件基本操作_第1頁
ch1-R軟件基本操作_第2頁
ch1-R軟件基本操作_第3頁
ch1-R軟件基本操作_第4頁
ch1-R軟件基本操作_第5頁
已閱讀5頁,還剩32頁未讀 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(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)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論