原创力文档- 1、本文档共62页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
2. An introduction to R (II);2.8 Matrices and arrays;矩阵的维数 dim( );拼接矩阵cbind() and rbind();取矩阵的上下三角部分lower.tri() upper.tri();upper.tri()的运算结果;提取上下三角矩阵;生成对角阵和矩阵提取对角元素 diag();矩阵的拉直as.vector( ); a - array(1:(2*4*3),dim=c(2,4,3));a
, , 1
[,1] [,2] [,3] [,4]
[1,] 1 3 5 7
[2,] 2 4 6 8
, , 2
[,1] [,2] [,3] [,4]
[1,] 9 11 13 15
[2,] 10 12 14 16
, , 3
[,1] [,2] [,3] [,4]
[1,] 17 19 21 23
[2,] 18 20 22 24;广义转置(Array Permutation);The recycling and arithmetic rules;矩阵运算;矩阵运算;矩阵运算;矩阵运算;矩阵行列运算 ;Array的维运算 ; 与tapply对比;Array和矩阵/向量/array之间的运算;Array和矩阵/向量/array之间的运算;Array和矩阵/向量/array之间的运算;矩阵运算;矩阵运算;验证;警告:计算机中的0是什么?;矩阵的Choleskey分解 chol()
对于正定矩阵A，可对其进行Choleskey分解，即：A=PP，其中P为上三角矩阵
A
[,1] [,2] [,3] [,4]
[1,] 2 1 1 1
[2,] 1 2 1 1
[3,] 1 1 2 1
[4,] 1 1 1 2
P=chol(A); P
[,1] [,2] [,3] [,4]
[1,] 1.414214 0.7071068 0.7071068 0.7071068
[2,] 0.000000 1.2247449 0.4082483 0.4082483
[3,] 0.000000 0.0000000 1.1547005 0.2886751
[4,] 0.000000 0.0000000 0.0000000 1.1180340;应用;矩阵运算 Singular Value Decomposition;矩阵的QR分解qr()
Am×n=QR,其中,Qm×n是正交阵,Rn×n是上三角阵
A=matrix(1:12,3,4)
B=qr(A);B
$qr
[,1] [,2] [,3] [,4]
[1,] -3.7416574 -8.552360 -1.336306e+01 -1.817376e+01
[2,] 0.5345225 1.963961 3.927922e+00 5.891883e+00
[3,] 0.8017837 0.988693 3.443426e-16 -3.439089e-16
注:当 mn时, R有很多元素是零,故QR分解可以大大减少存储空间和计算时间.;访问R阵和Q阵;还原;最小二乘拟合 Least Squares Fit;使用QR分解;两种方法结果比较;2.9 Lists; li - list(pi=pi, e = exp(1))
li
$pi
[1] 3.141593
$e
[1] 2.718282
li[1]
$pi
[1] 3.141593
li[[1]]
[1] 3.141593;递回索引（recursive indexing into lists ）;属性——名字;给阵列各维命名 ;列表的修改;列表的修改;列表的串联Concatenation;展为向量unflatten all data;2.10 Data frames;创建数据框;as.data.frame( );命名（类似矩阵情形）;修改变量名;查询数据框;attach and detach;增加变量;合并两个数据框merge();;;拉直/拆分数据框或列表;fix(my.dataset);2.11 Time series;创建时间序列 ts( );;2.12 Expressions;
文档评论（0）