R语言怎么样定义一个矩阵

Python09

R语言怎么样定义一个矩阵,第1张

>matrix(1:9,3,3,T)

[,1] [,2] [,3]

[1,]123

[2,]456

[3,]789

matrix是定义矩阵的函数，1:9表示1-9连续9个数，第一个3表示3行，第二个3表示3列，T表示转置，TRUE的缩写

运行第一行，后面那是enter后的，有问题可以互相交流

在用 featureCounts 做完表达矩阵的counts值后进行TPM需要注意这个细节问题，在计算TPM时每个基因需要除以各自的基因长度来校正基因长度，每一个样本又要除以它各自的文库大小校正测序深度。

因此，我们的表达矩阵，其实是需要除以两个长度不一的向量，而且方向不一样，一个是按照行来除以（基因长度），一个是按照列来除以（测序深度），在我之前的文章 https://www.jianshu.com/p/cd1c7b4ec6a2 RNA-seq数据分析一：(HISAT2+featureCounts） 中有提到

count 是表达矩阵，kb是不同基因长度向量，而 colSums(count) 是不同样本的长度向量。

矩阵除以向量，是按行的顺序来的，如果需要按列的顺序来(不同样本，一列一列)，就得先转置，再转回来

我们可以看到，一个矩阵除以向量是将向量按照行（第一行第一个，第二行第一个...）的顺序与矩阵中的数据相除。

提供的内容不足,写起来很难给你说明.

我提供一个思路给你:

先将两个矩阵的数据存入二维数组并定义一个新的二维数组用于存入R-L后的结果. 然后这样写:(假设是7*7的矩阵,定义新二维数组为 W(6,6) )

<按键精灵代码>:

For i = 0 To 6//横向递增

For j = 0 to 6 //纵向递增

W(i,j) = R(i,j) - L(i,j) //将计算结果存入W数组中

Next

Next

虽然是按键代码,但是思路都是互通的,自己理解一下.