使用p1 <- prcomp_irlba(beta, n=20)进行PCA计算，一直报错：
Error in (function (A, nv = 5, nu = nv, maxit = 1000, work = nv + 7, reorth = TRUE, : BLAS/LAPACK 例行程序'DLASCL'产生了错误-4信息”
后来发现是beta矩阵中有NA值，去除NA值就可以解决这个问题了。
beta=na.omit(beta)
p1 <- prcomp_irlba(beta, n=20)