R语言第六讲 数据的统计分析

基本命令练习

           下面的代码涵盖了一些分析数据常用的一些R语言的命令：

#基本向量、矩阵的一般操作
x <- c(1,3,2,5)
x
x = c(1,6,2)
x
y = c(1,4,3)
length(x)
length(y)
x+y
ls()
rm(x,y)
ls()
rm(list=ls())
?matrix
x=matrix(data=c(1,2,3,4), nrow=2, ncol=2)
x
x=matrix(c(1,2,3,4),2,2)
matrix(c(1,2,3,4),2,2,byrow=TRUE)
sqrt(x)
x^2
x=rnorm(50)
y=x+rnorm(50,mean=50,sd=.1)
cor(x,y)
set.seed(1303)#代码产生完全相同的一组随机数
# 用于设定随机数种子，一个特定的种子可以产生一个特定的伪随机序列，这个函数的主要目的，
# 是让你的模拟能够可重复出现，因为很多时候我们需要取随机数，但这段代码再跑一次的时候，
# 结果就不一样了，如果需要重复出现同样的模拟结果的话，就可以用set.seed()。在调试程序
# 或者做展示的时候，结果的可重复性是很重要的，所以随机数种子也就很有必要。
rnorm(50)
set.seed(3)#里面的参数可以是任意整数
y=rnorm(100)
mean(y)   #计算均值
var(y)    #计算方差
sqrt(var(y)) #计算开方
sd(y)    #计算标准差

#分享一些作图函数
x=rnorm(100)
y=rnorm(100)
plot(x,y)
plot(x,y,xlab="this is the x-axis",ylab="this is the y-axis",main="Plot of X vs Y")
pdf("Figure.pdf")#用jpeg（）函数绘制jpeg图
plot(x,y,col="green")
dev.off()
x=seq(1,10)
x
x=1:10
x
x=seq(-pi,pi,length=50)
# contour()产生一个个等高线图，表示三维数据
y=x
f=outer(x,y,function(x,y)cos(y)/(1+