初识R语言——用R处理矩阵的运算【随机数矩阵&所有元素求和】

    R语言在处理矩阵运算的时候较java和python来说更加容易,因为R语言本身就是一个统计学使用的软件。再已有的运用R进行的研究中,矩阵是比较常用的一种计算方法,下面对本次编程中矩阵的使用方法做一个总结。

1、产生一个服从泊松分布的随机数矩阵

    首先,产生一组服从泊松分布的随机数:

rpois(M,lambda = N)

    其中,M为产生的数量个数,N为泊松中的均值和方差
    接下来是生成矩阵的代码:

matrix(c(rpois(12,lambda = 2)),nrow=3,ncol=4,byrow=TRUE)

    其中,nrow为矩阵行数,ncol为矩阵列数,byrow为排列方法,“TRUE”为按行排列。【默认为按列排列即“FALSE”】
    之后计算行之和,同样保存在矩阵之中:

S<-matrix(rowSums(C),nrow = 3,ncol = 1,byrow = FALSE)#行之和S

    随机矩阵C:
C
    行之和S:
S

2、矩阵所有元素求和

    方法:求矩阵行和之后(n*1),左乘一个单位矩阵(1*n)

b=C[1,c(3,4)]#选取矩阵C第一行第三列和第四列的两个数
d=C[3,c(3,4)]#选取矩阵C第三行第三列和第四列的两个数
CA<-matrix(c(b,d),nrow=2,ncol=2,byrow=TRUE)#将选出的数组成新的矩阵
CA1<-matrix(rowSums(CA),nrow = 2,ncol = 1,byrow = FALSE)#CA的行之和
y1=matrix(1:1,nrow = 1,ncol = 2)#新建一个一行两列的单位矩阵
Q1=y1%*%CA1#左乘求和
CAA=Q1[1,1]#扫描矩阵所有元素总和CAA

CA:
CA
CA1:
CA1
y1:
y1
Q1:
Q1
    最后将Q[1,1]提取出来就可以了。

债见!
1

猜你喜欢

转载自blog.csdn.net/sinat_41805381/article/details/79626172
今日推荐