统计建模与R软件第二章习题…
发布日期:2021-10-16 07:12:20 浏览次数:32 分类:技术文章

本文共 1996 字,大约阅读时间需要 6 分钟。

原文地址:
作者:
   R,从零水平开始。国内真的没有一本像样的R教科书啊!勉强用用薛毅编的《统计建模与R软件》吧,找不出更好的了……

  工作环境仍是linux。
  第二章答案:
Ex2.1

x<-c(1,2,3)

y<-c(4,5,6)

e<-c(1,1,1)

z=2*x+y+e

z1=crossprod(x,y)#z1为x1与x2的内积 或者 x%*%y

z2=tcrossprod(x,y)#z1为x1与x2的外积 或者 x%o%y

z;z1;z2

要点:基本的列表赋值方法,内积和外积概念。内积为标量,外积为矩阵。


Ex2.2

A<-matrix(1:20,c(4,5));A

B<-matrix(1:20,nrow=4,byrow=TRUE);B

C=A+B;C

#不存在AB这种写法

E=A*B;E

F<-A[1:3,1:3];F

H<-matrix(c(1,2,4,5),nrow=1);H

#H起过渡作用,不规则的数组下标

G<-B[,H];G

要点:矩阵赋值方法。默认是byrow=FALSE,数据按列放置。

      取出部分数据的方法。可以用数组作为数组的下标取出数组元素。

 

Ex2.3

x<-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.4

n <- 5; H<-array(0,dim=c(n,n))

for (i in 1:n){for (j in 1:n){H[i,j]<-1/(i+j-1)}};H

G <- solve(H);G #求H的逆矩阵

ev <- eigen(H);ev #求H的特征值和特征向量

要点:数组初始化;for循环的使用

待解决:如何将很长的命令(如for循环)用几行打出来再执行?每次想换行的时候一按回车就执行了还没打完的命令...


Ex2.5

StudentData<-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(StudentData,file="studentdata.txt")

#把数据框StudentData在工作目录里输出,输出的文件名为studentdata.txt.


StudentData_a<-read.table("studentdata.txt");StudentData_a

#以数据框的形式读取文档studentdata.txt,存入数据框StudentData_a中。


write.csv(StudentData_a,"studentdata.csv")

#把数据框StudentData_a在工作目录里输出,输出的文件名为studentdata.csv,可用Excel打开.

要点:读写文件。read.table("file")

                write.table(Rdata,"file")

                read.csv("file")

                write.csv(Rdata,"file")

      外部文件,不论是待读入或是要写出的,命令中都得加双引号。


Ex2.7

Fun<-function(n){

if(n <= 0)

 list(fail="please input a integer above 0!")

else{

repeat{

  if(n==1) break

  else if(n%%2==0){n<-n/2}

  else n<- 3*n+1

}

list("sucess!")

}

}

在linux下新建一个R文件,输入上述代码,保存为"2.7.R"

然后在当前目录下进入R环境,输入source("2.7.R"),即打开了这个程序脚本。

然后就可以执行函数了。

输入Fun(67),显示

"sucess!"

输入Fun(-1),显示

$fail

[1] "please input a integer above 0!"
待解决:source("*.R")是可以理解为载入这个R文件吧?如何在R环境下关闭R文件呢?


OK,自己写的第一个R程序~~

转载地址:https://blog.csdn.net/sjpljr/article/details/70168871 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!

上一篇:统计建模与R软件第四章习题…
下一篇:统计建模与R软件第三章习题…

发表评论

最新留言

初次前来,多多关照!
[***.217.46.12]2024年04月13日 05时34分09秒

关于作者

    喝酒易醉,品茶养心,人生如梦,品茶悟道,何以解忧?唯有杜康!
-- 愿君每日到此一游!

推荐文章