《统计建模与R软件》第二章课后作业解题答案

x<-c(1,2,3);y<-c(4,5,6);e<-c(1,1,1)
z=2*x+y+e;z1=crossprod(x,y);z2=tcrossprod(x,y)
z;z1;z2

A<-matrix(1:20,c(4,5));A
B<-matrix(1:20,nrow=4,byrow=T);B
C=A+B;C
#无法实现D的AB相乘的写法
E=A*B;E
F<-A[1:3,1:3];F
#G有两种解题答案
H<-matrix(c(1,2,4,5),nrow=1);H
G<-B[,H];G
G<-B[,-3];G

#两种解题答案
c(rep(1,5),rep(2,3),rep(3,4),rep(4,2));c(rep(1:4,c(5,3,4,2)))

#第一道题第一种解题答案
i<-matrix(c(rep(1:5,5)),nrow=5,byrow=T);i
j<-matrix(c(rep(1:5,5)),nrow=5);j
H<-1/(i+j-1);H
#第一道题第二种解题答案
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
#第一道题第三种解题答案
library(Matrix)
H<-Hilbert(5);H
#第二题和第三题解题答案
N<-solve(H);N
E<-eigen(H);E

#输入数据
s<-data.frame("姓名"=c("张三","李四","王五","赵六","丁一"),
"性别"=c("女","男","女","男","女"),
"年龄"=c(14,15,16,14,15),
"身高(cm)"=c(156,165,157,162,159),
"体重(kg)"=c(42.0,49.0,41.5,52.0,45.5)
)
#把数据框s在工作目录里输出,输出文件名s.txt
write.table(s,file="s.txt")
#读取文档s.txt
S<-read.table("s.txt");S
#把数据框s在工作目录中输出，输出文件名为s.csv，用excel工具打开
write.table(s,file="s.csv")

Fun<-function(n){
if(n<=0)
list("请输入一个正整数")
else{
repeat{
if(n==1)  break
else if(n%%2==0){n<-n/2}
else n<-3*n+1
}
list("运算成功")
}
}
#输入数字进行验证
Fun(2);Fun(-2)