统计建模与R软件第二章作业

    xiaoxiao2021-03-25  76

    1、

    x<-matrix(c(1,2,3),3,1) y<-matrix(c(4,5,6),3,1) e<-matrix(c(1,1,1),3,1) z<-2*x+y+e crossprod(x,y)###x,y的内积///x%*%y tcrossprod(x,y)###x,y的外积///x%o%y

    2、

    A<-matrix(c(1:20),4,5,byrow=FALSE) B<-matrix(c(1:20),4,5,byrow=TRUE) c<-A+B D<-tcrossprod(A,B) E<-A*B F<-A[c(-4,-5),c(-4,-5)]###保留A的前三行,前三列。 G<-B[,-3]###不包含B的第三列。

    3、

    X<-c(rep(1,5),rep(2,3),rep(3,4),rep(4,2)) 结果:[1] 1 1 1 1 1 2 2 2 3 3 3 3 4 4

    4、生成一个5阶的HilbertJ矩阵n<-5

    H<-matrix(0,5,5) for(i in 1:n){ for(j in 1:n){ H[i,j] = 1/(i+j-1) } } det(H)###求行列式 solve(H)###求矩阵逆 ev<-eigen(H)###求特征值和特征向量 错误形式: H<-matrix(0,5,5) i=1 j=1 if(j<=5){ if(i<=5){ H[i,j]<-1/(i+j-1) i=i+1 } j=j+1 }##############*不是循环语句,循环只能用while,for,loop,repeat。*

    5、建立一个数据框:

    df<-data.frame(Num=c(1,2,3,4,5),Name=c(“张三”,”李四”,”王五”,”赵六”,”丁一”),Gender=c(“f”,”m”,”f”,”m”,”f”),Age=c(14,15,16,14,15),Height=c(156,165,157,162,159),Weight=c(42.0,49.0,41.5,52.0,45.5))

    6、读取文件 rt<-read.table(“foo.txt”) write.csv(df,file=”foo.csv”) 注意:读取文件夹中的数据时,读取格式如下:

    read.table(“E:/研究生文件/R语言与统计分析/R自测/R老虎机.R”) 不是 \ , 而是 /。

    7、编写一个函数:

    myfun<-function(n){ if(n<=0) {print(“要求输入一个正整数”) + } while (n!=1){ + if(n%%2==0) + {n<-(n/2)} +else n<-3*n+1 +} +print(“运算成功”) +}

    用myfun(0),myfun(8)验证!

    转载请注明原文地址: https://ju.6miu.com/read-33062.html

    最新回复(0)