统计建模与R软件课后参考答案

第二章

2.1

> x<-c(1,2,3);y<-c(4,5,6) > e<-c(1,1,1) > z<-2*x+y+e;z [1] 7 10 13

> z1<-crossprod(x,y);z1 [,1] [1,] 32 > z2<-outer(x,y);z2 [,1] [,2] [,3] [1,] 4 5 6 [2,] 8 10 12 [3,] 12 15 18 2.2

(1) > A<-matrix(1:20,nrow=4);B<-matrix(1:20,nrow=4,byrow=T) > C<-A+B;C (2) > D<-A%*%B;D (3) > E<-A*B;E (4) > F<-A[1:3,1:3] (5) > G<-B[,-3]

1 / 43

2.3

> x<-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2));x 2.4

> H<-matrix(nrow=5,ncol=5) > for (i in 1:5) + for(j in 1:5) + H[i,j]<-1/(i+j-1) (1)> det(H) (2)> solve(H) (3)> eigen(H) 2.5

> studentdata<-data.frame(姓名=c('张三','李四','王五','赵六','丁一') + ,性别=c('女','男','女','男','女'),年龄=c('14','15','16','14','15'),

+ 身高=c('156','165','157','162','159'),体重=c('42','49','41.5','52','45.5')) 2.6

> write.table(studentdata,file='student.txt') > write.csv(studentdata,file='student.csv') 2.7

count<-function(n) { if (n<=0)

print('要求输入一个正整数')

2 / 43

else{ repeat{ if (n%%2==0) n<-n/2 else n<-(3*n+1) if(n==1)break }

print('运算成功')} }

第三章

3.1

首先将数据录入为x。利用data_outline函数。如下 > data_outline(x) 3.2

> hist(x,freq=F)

> lines(density(x),col='red') > y<-min(x):max(x)

> lines(y,dnorm(y,73.668,3.9389),col='blue') > plot(ecdf(x),verticals=T,do.p=F) > lines(y,pnorm(y,73.668,3.9389))

3 / 43

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4