ch1-R软件基本操作

R软件基本操作R软件和软件包安装 数学运算 向量运算 矩阵运算 列表、数据框 概率函数

软件安装(Windows)

当前版本3.2.2。至http://www.77cn.com.cn, 点 击左边download下链接CRAN ,任意选择一个 镜像站点，比如USA下 http://www.77cn.com.cnr.Berkeley.edu 。在Download and Install R中选择Download R for Windows ,点击 base,选择文件Download R 3.2.2 for Windows , 点击鼠标右键另存为保存到指定目录中。下载好 后安装R。

软件包安装 比如安装“BRugs”,可以使用命令 >install.packages("BRugs",repos="http://www.77cn.com.cnr.B erkeley.edu",dep=TRUE) 调用软件包 >require(BRugs) # >library(BRugs) 升级软件包 >update.packages(repos="http://www.77cn.com.cnr.Berkeley. edu",dep=TRUE) 删除软件包 >remove.packages( "BRugs",lib = file.path("path", "to", "library"))

命令的执行 主窗口 R Editor窗口 (File-New script) require(Rcmdr)

常用的小技巧 查看工作目录： > getwd()[1] "D:/Program Files/R/R-2.11.1” 设定工作目录： > setwd("d:/Program Files/R") > getwd() [1] "d:/Program Files/R” 取得帮助： > help.start() #打开帮助文档 > help.search('panel data') > RSiteSearch('Kalman Filter') > ?lm 查看对象类型： > typeof() > class() 单引号和双引号的区别： > c("single ' qutoe", 'double " quote', "double \" quote", 'single \' quote') [1] "single ' qutoe" "double \" quote" "double \" quote" "single ' quote" 更改选项设定 > options(prompt='R>', continue='R')

运算符号 +-*/^ < <= > >= == != %%余数 %/%商 %in%判断 > 17%in%1:100 [1] TRUE

基本数学运算 加减乘除乘方运算 > 5+6*7/12^(1/3) [1] 23.34519 指数对数运算 > exp(2)#e^2 [1] 7.389056 > log(2,10)#lg(2) [1] 0.30103 > log(exp(1))#ln(e) [1] 15 6 3 7 12

求导 > D(expression(x^2), "x") 2*x > D(expression(sin(x^2)), "x") f<-function(x){1/sqrt(2*pi)*exp(-x^2/2)} cos(x^2) * (2 * x) integrate(f,-1.96,1.96)$value 积分：利用函数integrate() > integrate(dnorm,-1.96,1.96)$value [1] 0.9500042 > integrate(dnorm,-Inf,Inf)$value [1] 1 #对于多重积分计算，利用adapt库中函数adapt,不允许用无 穷界限。

基本数学运算(续)

复数运算 ‘Re’-实部 ‘Im’-虚部 ‘Mod’-模 ‘Arg’-角度 ‘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] x^2+2x-3=0 [1,] 2 3 [2,] 1 1 2x+3y=1 x+y=-1 > B<-c(1,-1) > solve(A,B) [1] -4 3

集合运算 并 > union(x,y) 交 > intersect(x,y) 补 > setdiff(x,y)> x<-1:20 > y<-3:15 > intersect(x,y) [1] 3 4 5 6 7 8 9 10 11 12 13 14 15 > union(x,y) [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 > match(y,x) [1] 3 4 5 6 7 8 9 10 11 12 13 14 15 > z<-c(23,15,12) > match(y,z) [

1] NA NA NA NA NA NA NA NA NA 3 NA NA 2 > setdiff(x,y) #从x中消去y [1] 1 2 16 17 18 19 20

排列组合 choose(n,k)-组合数 combn(n,k)-列出所有组合 factorial(n)-阶乘 require(e1071) > permutations(3) [,1] [,2] [,3] [1,] 1 2 3 [2,] 2 1 3 [3,] 2 3 1 [4,] 1 3 2 [5,] 3 1 2 [6,] 3 2 1

序列和向量 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)

Error in rep(3, 1:3) : invalid 'times' argumentx <- rep(3:5,1:3) >x [1] 3 4 4 5 5 5 x <- rep(c(1,10),c(4,5)) w <- c(1,3,x,z);w[3]

向量运算 x <- rep(0,10);z <- 1:3;x+z [1] 1 2 3 1 2 3 1 2 3 1 Warning longer object length is not a multiple object length in: x + z x*z [1] 0 0 0 0 0 0 0 0 0 0 Warning longer object length is not a multiple object length in: x * z rev(x) z <- c("no cat","has","nine","tails") z[1]=="no cat" [1] TRUE message: of shorter message: of shorter

向量名和元素附加 x <- 1:3;names(x) <- LETTERS[1:3] x ABC 123 set.seed(10);append(x,runif(3),after=2) A B C 1.0000 2.0000 0.5075 0.3068 0.4269 3.0000

向量赋值 z <- 1:5 z[7] <- 8;z [1] 1 2 3 4 5 NA 8 z <- NULL z[c(1,3,5)] <- 1:3; z [1] 1 NA 2 NA 3 rnorm(10)[c(2,5)] z[-c(1,3)] #去掉第1、3元素. z[(length(z)-4):length(z)] #最后五个元素.

向量的大小次序 set.seed(10);z<-sample(1:100,10);z #比较sample(1:100,10,rep=T) [1] 51 31 42 68 9 22 26 94 57 40 order(z)#给出数据从小到大的位置，即最小的是第5个 [1] 5 6 7 2 10 3 1 9 4 8 z[order(z)] [1] 9 22 26 31 40 42 51 57 68 94 sort(z) [1] 9 22 26 31 40 42 51 57 68 94 which(z==max(z))#给出下标 [1] 8

矩阵操作 > (A<-matrix(1:12,3,4)) [,1] [,2] [,3] [,4] [1,] 1 4 7 10 [2,] 2 5 8 11 [3,] 3 6 9 12 > (A<-matrix(1:12,3,4,byrow=T)) [,1] [,2] [,3] [,4] [1,] 1 2 3 4 [2,] 5 6 7 8 [3,] 9 10 11 12 > dim(A) [1] 3 4

矩阵操作(续) > nrow(A) [1] 3 > dim(A)[1] [1] 3 > ncol(A) [1] 4 > dim(A)[2] [1] 4 > A[3,] [1] 9 10 11 12 > A[,3] [1] 3 7 11 > A[3,3] [1] 11

本文来源：https://www.bwwdw.com/article/gvbi.html