使用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)