统计建模与R软件 习题

2.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的内积 z2=tcrossprod(x,y)#z1为x1与x2的外积 z;z1;z2

2.2

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

B<-matrix(1:20,nrow=4,byrow=TRUE);B C=A+B;C #D不存在 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

2.3

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

2.4

n<-5;

H<-array(0,c(n,n))#array不能用matrix代替，原因？ for(i in 1:n) {

for(j in 1:n)

{H[i,j]<-1/(i+j-1) }

} H

2.5

T2.3<-read.table(\"表2.3 学生数据.txt\

2.6

write.table(T2.3,file=\"123.txt\")

#把数据框T2.3在工作目录里输出，输出的文件名为123.txt. T456<-read.table(\"123.txt\

#以数据框的形式读取文档123.txt，该数据框的名字为T456.

write.csv(T456,\"7.csv\")

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

2.7

ST<-function(n){

if(n<=0) list(\"Please enter a positive integer.\") else{repeat{

if(n==1) break

else if(n/2==as.integer(n/2)) n<-n/2 else n<-3*n+1 }

list(\"Successful operation\") }

}

ST(10);ST(-4);ST(11)