r语言与多元数据.ppt.pdf

数据矩阵及R语言表示 创建一个向量c() x1=c(171,175,159,155,152,158,154,164,168,1 66,159,164)#创建一个向量 x2=c(57,64,41,38,35,44,41,51,57,49,47,46) 向量的长度length() 向量的数据类型mode() length(x1)# 向量的长度 mode(x1)#数据的类型 创建一个矩阵 （1）合并命令 rbind()#按行合并 cbind()#按列合并 rbind(x1,x2)#按行合并 cbind(x1,x2)#按列合并 （2 ）生成矩阵 Matrix(data=NA,rnow=1,ncol=1,byrow=FALSE,dimname=NU LL) Data为数据，rnow为行数，ncol为列数（行与列的乘积等于数 据的个数），byrow控制是否按照行排列，dimname给定行 和列的名称。 matrix(x1,nrow=3,ncol=4)#利用x1数据创建矩阵 matrix(x1,nrow=4,ncol=3)#创建行数列数发生变化的矩阵 matrix(x1,nrow=4,ncol=3,byrow=T)#创建按照行排列的矩阵 A=matrix(1:12,nrow=3,ncol=4)#创建矩阵 矩阵转置 t()或者transpose() t(A)#求矩阵转置 矩阵相加减 用符号+，-表示 A=B=matrix(1:12,nrow=3,ncol=4)#创建两个 相同的矩阵 A+B#矩阵加法 A-B#矩阵减法 矩阵相乘 矩阵的相乘用符号%*%来表示 A=matrix(1:12,nrow=3,ncol=4)#创建矩阵 B=matrix(1:12,nrow=4,ncol=3)#创建矩阵 A%*%B#求矩阵的乘积 矩阵对角元素相关运算 对一个向量应用diag()函数将产生以这个向 量为对角元素的对角矩阵，对一个正整数k 应用diag()将产生k维单位矩阵。 A=matrix(1:16,nrow=4,ncol=4)#创建行列数 相等矩阵 diag(A)#获得矩阵对角线元素 diag(diag(A))#利用对角线元素创建对角矩阵 diag(3)#创建3阶单位矩阵 矩阵求逆 矩阵求逆函数为solve() A=matrix(rnorm(16),4,4)#创建矩阵(用标准正 太分布产生16个随机数) solve(A)#求矩阵的逆 矩阵的特征值与特征向量 A UU  矩阵A 的谱分解为 ,其中 是A 的特征值组成的对角 U 矩阵， 的列为A 的特征值对应的特征向量，得到 和 U  的函数为eign() eigen(x, symmetric, only.values = FALSE, EISPACK = FALSE) x为矩阵，symmetric指定矩阵x是否为对称矩阵，若不指 定，系统将自动检测x是否为对角矩阵。 (A=diag(4)+1)#创建矩阵 (A.e=eigen(A,symmetric=T))#求矩阵的特征值与特征向量 A.e$vectors%*%diag(A.e$values)%*%t(A.e$vectors)#特 征向量矩阵U和特征值矩阵D与原矩阵A 的关系A=UDU 矩阵的Choleskey分解 对于正定矩阵A ，可以进行Choleskey分 解，即A=P’P ，其中P为上三角矩阵， Choleskey分解的矩阵为chol() 。 (A.c=chol(A))#