R软件及统计分析因子列表数据框及输入输出PPT.ppt

实验目的;因子（factor）和有序因子 ;因子是一种特殊的字符型向量，其中每一个元素取一组离散值中的一个，而因子对象有一个特殊属性levels表示这组离散值（用字符串表示）。例如： x - c(男, 女, 男, 男， 女) y - factor(x) y [1] 男 女 男 男 女 Levels: 男 女 ;可以自行指定各离散取值水平（levels），不指定时由x的不同值来求得。 labels可以用来指定各水平的标签，不指定时用各离散取值的对应字符串。 exclude参数用来指定要转换为缺失值（NA）的元素值集合。如果指定了levels，则因子的第i个元素当它等于水平中第j个时元素值取“j”，如果它的值没有出现在levels中则对应因子元素值取NA。 ordered取真值时表示因子水平是有次序的（按编码次序）。 可以用is.factor()检验对象是否因子，用as.factor()把一个向量转换成一个因子。; x=c(1,0,1,1,0) y=factor(x,levels=sort(unique(x),decreasing=T),labels=c(男, 女),exclude=NA,order=F) y [1] 男 女 男 男 女 Levels: 男 女 ;列表（list）定义; rec = list(name=李明, age=30, scores=c(85, 76, 90)) rec $name [1] 李明 $age [1] 30 $scores [1] 85 76 90 ;在定义列表时如果指定了元素的名字（如rec中的name，age，scores），则引用列表元素还可以用它的名字作为下标，格式为“列表名[[元素名]]”，如： rec[[age]] [1] 30 另一种格式是“列表名$元素名”，如： rec$age [1] 30 其中“元素名”可以简写到与其它元素名能够区分的最短程度，比如“rec$s”可以代表 “rec$score”。 这种写法方便了交互运行，编写程序时一般不用简写，以免降低程序的可读性。;使用元素名的引用方法可以让我们不必记住某一个下标代表那一个元素，而直接用易记的元素名来引用元素。事实上，已知向量和矩阵都可以指定元素名、行名、列名。 定义列表使用list()函数，每一个自变量变成列表的一个元素，自变量可以用“名字＝值”的方式给出，即给出列表元素名。 自变量的值被复制到列表元素中，自变量如果是变量并不会与该列表元素建立关系（改变该列表元素不会改变自变量的值）。(例见后);修改列表 ; rec$sex - 男 rec[[6]] - 161 rec ;第五号元素因为没有定义所有其值是“NULL”，这是空对象的记号。 如果rec是一个向量，则其空元素为“NA”，这是缺失值的记号。从这里我们也可以体会“NULL”与“NA”的区别。 几个列表可以用连接函数c()连接起来，结果仍为一个列表，其元素为各自变量的列表元素。如： list.ABC - c(list.A, list.B, list.C) （注意在S中句点是名字的合法部分，一般没有特殊意义。） ;几个返回列表的函数 ;一、特征值和特征向量 函数eigen(x)对对称矩阵x计算其特征值和特征向量，返回结果为一个列表，列表的两个成员（元素）为values和vectors。例如： ev=eigen((1:3)%o%(1:3)) ev $values [1] 1.400000e+01 9.176554e-16 -1.023391e-15 $vectors [,1] [,2] [,3] [1,] -0.2672612 0.6067775 0.7485937 [2,] -0.5345225 0.5530299 -0.6390960 [3,] -0.8017837 -0.5709457 0.1765328 特征向量按矩阵存放，每一列为一个特征向量。;二、奇异值分解及行列式 svd()函数进行奇异值分解 自学 ;三、最小二乘拟合与QR分解 关于最小二乘拟合还可参见ls.diag() 函数（查看帮助）。 函数qr(x)返回x的QR分解结果。矩阵X的QR分解为X=Q*R ，Q为对角线元素都等于1的下三角阵，R为上三角阵。 函数结果为一个列表，成员qr 为一个矩阵，其上三角部分（包括对角线）分解的R，其下三角部分（不包括对角线）为分解的Q。其它成员为一些辅助信息。; x=rbind(c(1,2,3),c(4,5,6),c(7,8,9))