数量遗传学-矩阵

2024-02-24 17:21:11 浏览数 (2)

代码语言:english复制
#矩阵
A<-matrix(c(7,4,-3,4,6,93,0,93,23,5),byrow=TRUE,nrow=5,ncol=2)
dim(A)

B<-diag(c(1,2,3,4,5)) #对角矩阵
I<-diag(c(1),nrow=5,ncol=5) #单位矩阵

矩阵运算

代码语言:english复制
C<-A B #相同位置元素相加或相减
C<-A%*%B#矩阵相乘
C<-A%x%B #直积
tA<-t(A) #转置
trA<-sum(diag(A)) #迹和
r<-qr(D)$rank #秩(rank)
detA<-det(A) # 行列式
Ainv<-solve(A)

#奇异阵的广义逆矩阵
library(MASS)
ginv(D)

#正定矩阵的cholesky分解
T<-chol(A) 

#特征值与特征向量
eigen(A)

0 人点赞