R的操作入门熟悉
产生向量:
a=c(1,2,3) //产生 1,2,3向量
a=1:10-1 //产生 0 - 9数字
a=seq(5,20,by=2) // 以2增长
a=seq(5,120,length=10) //5到20,总数为10,长度系统自己计算得到
矩阵:
生成矩阵:A=matrix(1:16,nrow,ncol,byrow=T)
取矩阵列:temp<-A[,c(1,3)]
PS:为什么c前面有逗号?矩阵二维,逗号分割行、列索引,[, c(....)]按列subgroup,[c(....), ]按行subgroup。
which函数
which.max(a) //向量a中,最大值的下标
which(a>1 & a<6) // 1<a<6之间的a下标
t(a) //a的转置
a%*%b //矩阵相乘
diag(a) //求对角元素 diag(4) :4阶对角矩阵
solve(a,b) //ax=b 方程组解
a.e=eigen(a,symmetric=T) //求出特征向量,特征根。 若需要引用特征根: a.e$values 特征向量 a.e$vectors
数据框
x=data.frame(x1,x2) //可理解为每一行为观察值,x1,x2长度要相等
plot(x) //作出x的图
读文件:
工作目录下放置文本。
x=read.table("abc.txt")
excel 先另存为prn文本,w<-read.table("test.prn",header=T)
for循环
for(i in 1:59) (a[i]=i*2+3)
while(a[i]<121) (i=i+1; a[i]=a[i-1]+2)
R脚本
source("D:\\h.r") //其中print() 函数用来显示
分布函数
rnorm(n,mean=0,sd=1) 高斯(正态)
rexp(n,rate=1) 指数分布
rgamma(n,shape,scale=1) 伽玛分布
rpois(n,lambda) Poisson分布
rweibull(n,shape,scale=1) Weibull(韦伯)分布
rcauchy(n,location=0,scale=1) Cauchy分布
rbeta(n,shape1,shape2) 贝塔分布
rt(n,df) t 分布
rf(n,df1,df2) F分布
rchisq(n,df) 卡方分布
rbinom(n,size,prob) 二项分布
rgeom(n,prob) 几何分布
rhyper(nn,m,n,k) 超几何分布
rlogis(n,location=0,scale=1) logisitic分布
rlnorm(n,meanlog=0,sdlog=1) 对数分布
rnbinom(n,size,prob) 负二项分布
runif(n,min=0,max=1) 均匀分布
rwilcox(nn,m,n) ,rsignrank(nn,n) Wilcoxon分布
写文件
x=data.frame(num,x1,x2,x3)
write.table(x,file="d:\\mark.txt",col.name=F,row.name=F sep=" ")
colMeans(x)[c("x1")] //计算平均值,只显示x1列
apply(x,2,mean) //在列的方向求平均 2表示列
apply(x[("x1","x2")],1,sum] //对每一行求和