R作业.doc

第二章： 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) #输入2x2阶矩阵 x x= matrix (c(1 ,2 ,3 ,4) ,2,2) #输入2x2阶矩阵 matrix (c(1 ,2 ,3 ,4) ,2,2, byrow = TRUE) #输入2x2阶矩阵，数据按行排列 sqrt(x) #平方根 x^2 #平方 x= rnorm (50) #生成一个个数为50的随机变量(正态分布) y=x+ rnorm (50, mean =50, sd =.1) #mean均值，sd标准差 cor (x,y) #x，y的相关系数 set .seed (1303) #程序运行1303次,都是从同一个种子产生的随机数 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 ) #保存到Figure.pdf中，pdf格式 plot(x,y,col = green ) dev .off () x=seq (1 ,10) #建序列函数，从1到10 x x =1:10 #建序列函数，从1到10 x=seq (-pi ,pi , length =50) #建序列函数，从-pi到pi，长度为50 y=x f= outer (x,y, function (x,y)cos (y) /(1+x^2) )#在绘三维图形时，f并不是简单地关于x与y的某些运算，而是需要function (x,y)cos (y) /(1+x^2)关系下作外积运算形成网络 contour (x,y,f) #二维等值线 contour (x,y,f, nlevels =45, add =T) fa =(f-t(f))/2 contour (x,y,fa , nlevels =15) image (x,y,fa) persp (x,y,fa) #三维网格图 persp (x,y,fa , theta =30) persp (x,y,fa , theta =30, phi =20) persp (x,y,fa , theta =30, phi =70) persp (x,y,fa , theta =30, phi