统计建模与R软件课后习题答案2-5章.docVIP

PAGE/NUMPAGES

?R，从零水平开始。国内真的没有一本像样的R教科书啊！勉强用用薛毅编的《统计建模与R软件》吧，找不出更好的了……??工作环境仍是linux。??第二章答案：Ex2.1x-c(1,2,3)y-c(4,5,6)e-c(1,1,1)z=2*x+y+ez1=crossprod(x,y)#z1为x1与x2的内积或者x%*%yz2=tcrossprod(x,y)#z1为x1与x2的外积或者x%o%yz;z1;z2要点：基本的列表赋值方法，内积和外积概念。内积为标量，外积为矩阵。Ex2.2A-matrix(1:20,c(4,5));AB-matrix(1:20,nrow=4,byrow=TRUE);BC=A+B;C#不存在AB这种写法E=A*B;EF-A[1:3,1:3];FH-matrix(c(1,2,4,5),nrow=1);H#H起过渡作用，不规则的数组下标G-B[,H];G要点：矩阵赋值方法。默认是byrow=FALSE,数据按列放置。??????取出部分数据的方法。可以用数组作为数组的下标取出数组元素。?Ex2.3x-c(rep(1,times=5),rep(2,times=3),rep(3,times=4),rep(4,times=2));x#或者省略times=，如下面的形式x-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2));x要点：rep（）的使用方法。rep（a,b）即将a重复b次?Ex2.4n-5;H-array(0,dim=c(n,n))for(iin1:n){for(jin1:n){H[i,j]-1/(i+j-1)}};HG-solve(H);G#求H的逆矩阵ev-eigen(H);ev#求H的特征值和特征向量要点：数组初始化；for循环的使用待解决：如何将很长的命令（如for循环）用几行打出来再执行？每次想换行的时候一按回车就执行了还没打完的命令...Ex2.5StudentData-data.frame(name=c(zhangsan,lisi,wangwu,zhaoliu,dingyi),sex=c(F,M,F,M,F),age=c(14,15,16,14,15),height=c(156,165,157,162,159),weight=c(42,49,41.5,52,45.5));StudentData要点：数据框的使用待解决：SSH登陆linux服务器中文显示乱码。此处用英文代替。Ex2.6?write.table(Rdata,file)????????????????read.csv(file)????????????????write.csv(Rdata,file)??????外部文件，不论是待读入或是要写出的，命令中都得加双引号。Ex2.7Fun-function(n){if(n=0)?list(fail=pleaseinputaintegerabove0!)else{repeat{??if(n==1)break??elseif(n%%2==0){n-n/2}??待解决：source(*.R)是可以理解为载入这个R文件吧？如何在R环境下关闭R文件呢？OK,自己写的第一个R程序~~

第二章答案：Ex2.1x-c(1,2,3)y-c(4,5,6)e-c(1,1,1)z=2*x+y+ez1=crossprod(x,y)#z1为x1与x2的内积或者x%*%yz2=tcrossprod(x,y)#z1为x1与x2的外积或者x%o%yz;z1;z2要点：基本的列表赋值方法，内积和外积概念。内积为标量，外积为矩阵。Ex2.2A-matrix(1:20,c(4,5));AB-matrix(1:20,nrow=4,byrow=TRUE);BC=A+B;C#不存在AB这种写法E=A*B;EF-A[1:3,1:3];FH-matrix(c(1,2,4,5),nrow=1);H#H起过渡作用，不规则的数组下标G-B[,H];G要点：矩阵赋值方法。默认是byrow=FALSE,数据按列放置。取出部分数据的方法。可以用数组作为数组的下标取出数组元素。Ex2.3x-c(rep(1,times=5),rep(2,times=3),rep(3,times=4),rep(4,times=2));x#或者省略times=，如下面的形式x-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2));x要点：rep（）的使用方法。rep（a,b）即将a重复b次Ex2.4n-5;H-array(0,dim=c(n,n))for(iin1:n){for(jin1:n){H[i