"Lapack routine dgesv: system is exactly singular: U[5,5] = 0"
Calls: commonfactorGWAS -> %dopar% -> <Anonymous>
Here is my code:
AF_factor <- commonfactorGWAS(covstruc = LDSCoutput, SNPs = AF_MGWAS_sumstats , estimation = "DWLS", cores = 15, toler = 1e-60 , SNPSE = FALSE, parallel = TRUE,GC="standard")
here is the head of LDscoutput
$V
[,1] [,2] [,3] [,4] [,5]
[1,] 2.031495e-04 2.199144e-05 -5.543132e-06 1.345035e-05 1.619604e-05
[2,] 2.199144e-05 6.194667e-05 9.888094e-06 2.138668e-06 1.043033e-05
[3,] -5.543132e-06 9.888094e-06 2.537404e-04 -2.656666e-06 -1.785127e-06
[4,] 1.345035e-05 2.138668e-06 -2.656666e-06 1.996487e-05 4.346950e-06
[5,] 1.619604e-05 1.043033e-05 -1.785127e-06 4.346950e-06 1.558995e-05
[6,] 1.222046e-05 2.238214e-05 2.051274e-05 5.457345e-06 5.443163e-06
[7,] -2.144103e-05 3.450537e-06 5.315150e-05 4.470802e-06 4.603180e-06
[8,] 6.777807e-06 3.250754e-06 2.442156e-06 5.488622e-06 3.181726e-06
[9,] 5.057405e-06 8.183235e-06 5.704989e-06 3.786424e-06 7.641225e-06
[10,] 1.398211e-05 2.353221e-05 1.925616e-04 -9.741097e-06 -8.569595e-06
[11,] 7.281276e-06 6.307137e-07 1.089621e-05 9.320047e-06 6.610598e-07
[12,] -7.893921e-06 3.513355e-06 2.615673e-05 6.529312e-07 9.366111e-07
[13,] 3.366773e-06 9.578546e-07 -1.238012e-06 4.608643e-06 2.202693e-06
[14,] 3.125222e-06 1.306617e-06 4.998588e-07 2.830650e-06 2.877961e-06
[15,] 5.490855e-06 3.225999e-06 2.702072e-06 1.414718e-06 6.787058e-06
[,6] [,7] [,8] [,9] [,10]
[1,] 1.222046e-05 -2.144103e-05 6.777807e-06 5.057405e-06 1.398211e-05
[2,] 2.238214e-05 3.450537e-06 3.250754e-06 8.183235e-06 2.353221e-05
[3,] 2.051274e-05 5.315150e-05 2.442156e-06 5.704989e-06 1.925616e-04
[4,] 5.457345e-06 4.470802e-06 5.488622e-06 3.786424e-06 -9.741097e-06
[5,] 5.443163e-06 4.603180e-06 3.181726e-06 7.641225e-06 -8.569595e-06
[6,] 8.376462e-05 9.735481e-06 7.078647e-06 1.870274e-05 1.926272e-05
[7,] 9.735481e-06 1.681985e-04 4.081684e-07 2.183360e-06 1.067483e-04
[8,] 7.078647e-06 4.081684e-07 1.219203e-05 5.847097e-06 6.888067e-06
[9,] 1.870274e-05 2.183360e-06 5.847097e-06 1.799373e-05 7.548939e-06
[10,] 1.926272e-05 1.067483e-04 6.888067e-06 7.548939e-06 1.976499e-03
[11,] 9.084638e-07 1.870238e-05 5.336604e-06 1.878276e-06 1.063224e-05
[12,] 5.075382e-06 3.036960e-05 -4.578612e-07 2.516906e-06 3.263163e-05
[13,] 4.015436e-06 1.462215e-06 3.490226e-06 2.597928e-06 8.458484e-06
[14,] 3.946939e-06 1.497038e-06 3.456198e-06 3.526601e-06 7.533071e-06
[15,] 1.204805e-05 4.903666e-06 3.745328e-06 1.020440e-05 9.716050e-07
[,11] [,12] [,13] [,14] [,15]
[1,] 7.281276e-06 -7.893921e-06 3.366773e-06 3.125222e-06 5.490855e-06
[2,] 6.307137e-07 3.513355e-06 9.578546e-07 1.306617e-06 3.225999e-06
[3,] 1.089621e-05 2.615673e-05 -1.238012e-06 4.998588e-07 2.702072e-06
[4,] 9.320047e-06 6.529312e-07 4.608643e-06 2.830650e-06 1.414718e-06
[5,] 6.610598e-07 9.366111e-07 2.202693e-06 2.877961e-06 6.787058e-06
[6,] 9.084638e-07 5.075382e-06 4.015436e-06 3.946939e-06 1.204805e-05
[7,] 1.870238e-05 3.036960e-05 1.462215e-06 1.497038e-06 4.903666e-06
[8,] 5.336604e-06 -4.578612e-07 3.490226e-06 3.456198e-06 3.745328e-06
[9,] 1.878276e-06 2.516906e-06 2.597928e-06 3.526601e-06 1.020440e-05
[10,] 1.063224e-05 3.263163e-05 8.458484e-06 7.533071e-06 9.716050e-07
[11,] 4.708489e-05 1.654007e-05 4.655501e-06 2.956237e-06 6.488007e-07
[12,] 1.654007e-05 4.461739e-05 2.343577e-08 -3.819908e-07 1.210018e-06
[13,] 4.655501e-06 2.343577e-08 8.741990e-06 3.314809e-06 2.685256e-06
[14,] 2.956237e-06 -3.819908e-07 3.314809e-06 4.864370e-06 4.376309e-06
[15,] 6.488007e-07 1.210018e-06 2.685256e-06 4.376309e-06 2.338487e-05
$S
HF CARDMYO PERICARD PAROXTAC AF
[1,] 0.10573958 0.05862415 0.04589085 0.03256111 0.04449670
[2,] 0.05862415 0.07481068 0.03640735 0.03375220 0.04592359
[3,] 0.04589085 0.03640735 0.01491200 0.01898525 0.02707625
[4,] 0.03256111 0.03375220 0.01898525 0.03585887 0.03168534
[5,] 0.04449670 0.04592359 0.02707625 0.03168534 0.07513127
$I
[,1] [,2] [,3] [,4] [,5]
[1,] 1.04776847 0.156393797 0.042357034 0.016483349 0.12526071
[2,] 0.15639380 1.030435868 0.041963538 -0.002700575 0.20654119
[3,] 0.04235703 0.041963538 1.004571015 -0.003469138 0.03657651
[4,] 0.01648335 -0.002700575 -0.003469138 1.008450766 0.19945965
[5,] 0.12526071 0.206541195 0.036576509 0.199459646 1.10944217
$N
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,] 23209.93 30176.08 10202.76 55203.64 81224.31 35903.94 12689.86 68660.03
[,9] [,10] [,11] [,12] [,13] [,14] [,15]
[1,] 101023.9 4104.781 23215.59 34157.62 124732 187739.4 280205.7
$m
[1] 1173569
And I also looked through previous similarly question, and I check there is no 0 for any b and se column.
I'm wondering is that because we have a sumstats dataset only include 4104 individuals? if it's not, could you help me to fix this problem? Thanks in advance!
[1] "The common factor initially failed to converge. A lower bound of 0 on residual variances has been added to try and troubleshoot this."
[1] "The common factor model failed to converge on a solution. Please try specifying an alternative model using the usermodel function."
Then I tried 2-factor and 3- factor models
Call:
factanal(factors = 2, covmat = Ssmooth, rotation = "promax")
Uniquenesses:
AF_GRCh37 Cardiomyopathy
0.442 0.005
Heart_valve_diseases_(aortic_stenosis) HeartFailure
0.850 0.251
Paroxysmal_tachycardia Pericarditis
0.308 0.005
Loadings:
Factor1 Factor2
AF_GRCh37 0.134 0.664
Cardiomyopathy 1.011
Heart_valve_diseases_(aortic_stenosis) -0.240 0.466
HeartFailure 0.908
Paroxysmal_tachycardia 0.338 0.595
Pericarditis 1.004
Factor1 Factor2
SS loadings 2.227 1.838
Proportion Var 0.371 0.306
Cumulative Var 0.371 0.677
Factor Correlations:
Factor1 Factor2
Factor1 1.000 0.555
Factor2 0.555 1.000
And the corresponding CFA model is like this:
Call:
factanal(factors = 3, covmat = Ssmooth, rotation = "promax")
Uniquenesses:
AF_GRCh37 Cardiomyopathy
0.422 0.005
Heart_valve_diseases_(aortic_stenosis) HeartFailure
0.005 0.308
Paroxysmal_tachycardia Pericarditis
0.264 0.005
Loadings:
Factor1 Factor2 Factor3
AF_GRCh37 0.750
Cardiomyopathy 0.992
Heart_valve_diseases_(aortic_stenosis) 0.994
HeartFailure -0.135 0.877
Paroxysmal_tachycardia 0.198 0.737
Pericarditis 0.981
Factor1 Factor2 Factor3
SS loadings 2.005 1.875 1.004
Proportion Var 0.334 0.313 0.167
Cumulative Var 0.334 0.647 0.814
Factor Correlations:
Factor1 Factor2 Factor3
Factor1 1.0000 -0.0415 0.610
Factor2 -0.0415 1.0000 0.269
Factor3 0.6099 0.2686 1.000
But it shows error:
Error in eigen(V_LD) : 0 x 0 matrix
I didn't find any answer about this error, thanks for your time !
LDSCoutput
$V
[,1] [,2] [,3] [,4] [,5]
[1,] 6.634380e-04 5.922035e-05 6.105079e-05 2.452001e-05 1.219305e-04
[2,] 5.922035e-05 6.192872e-05 -1.713810e-05 4.796429e-06 2.585079e-05
[3,] 6.105079e-05 -1.713810e-05 9.566928e-04 4.561080e-06 6.410901e-05
[4,] 2.452001e-05 4.796429e-06 4.561080e-06 5.219829e-05 4.651135e-06
[5,] 1.219305e-04 2.585079e-05 6.410901e-05 4.651135e-06 8.368324e-05
[6,] 2.026113e-04 2.190326e-04 -6.466438e-05 1.702795e-05 8.990660e-05
[7,] 1.124142e-05 2.148922e-05 -2.966638e-05 4.472015e-06 1.356114e-05
[8,] 2.229469e-05 -8.367692e-06 3.494793e-04 7.467856e-06 1.360570e-05
[9,] -2.804225e-06 5.919124e-06 1.890440e-05 8.293495e-06 1.424664e-06
[10,] 1.293769e-06 1.284126e-05 -1.066283e-05 -3.923771e-06 1.255282e-05
[11,] 2.696624e-05 8.506095e-05 -1.143686e-04 1.901643e-05 3.849050e-05
[12,] 3.118345e-04 1.183130e-04 -3.986870e-04 8.661613e-05 5.143988e-05
[13,] -3.111370e-05 -7.582109e-06 7.189813e-05 7.258930e-06 2.720809e-07
[14,] 3.392233e-05 -3.769431e-06 1.708443e-04 -4.221119e-06 1.021996e-06
[15,] 8.148975e-05 -3.334130e-05 1.285410e-03 2.680681e-05 5.117130e-05
[16,] 1.046518e-05 4.156791e-06 1.562696e-05 4.198132e-06 4.263664e-06
[17,] -4.929012e-07 4.162904e-06 -9.542051e-06 6.022992e-06 7.521804e-06
[18,] -1.157290e-05 2.060207e-05 6.765845e-05 2.932552e-05 4.395275e-06
[19,] 5.671485e-05 6.876690e-06 2.432035e-05 -8.332178e-07 3.432684e-05
[20,] 2.728035e-06 4.574546e-05 -3.519962e-05 -1.380948e-05 4.475861e-05
[21,] 1.262746e-04 2.724515e-04 -3.714117e-04 5.139420e-05 1.685144e-04
[,6] [,7] [,8] [,9] [,10]
[1,] 2.026113e-04 1.124142e-05 2.229469e-05 -2.804225e-06 1.293769e-06
[2,] 2.190326e-04 2.148922e-05 -8.367692e-06 5.919124e-06 1.284126e-05
[3,] -6.466438e-05 -2.966638e-05 3.494793e-04 1.890440e-05 -1.066283e-05
[4,] 1.702795e-05 4.472015e-06 7.467856e-06 8.293495e-06 -3.923771e-06
[5,] 8.990660e-05 1.356114e-05 1.360570e-05 1.424664e-06 1.255282e-05
[6,] 7.768122e-04 7.642568e-05 -3.408236e-05 2.080900e-05 4.637237e-05
[7,] 7.642568e-05 8.776124e-05 -3.976552e-05 4.198629e-06 6.547519e-06
[8,] -3.408236e-05 -3.976552e-05 6.179345e-04 1.094728e-05 1.130748e-06
[9,] 2.080900e-05 4.198629e-06 1.094728e-05 2.120082e-05 -9.727880e-08
[10,] 4.637237e-05 6.547519e-06 1.130748e-06 -9.727880e-08 3.372584e-05
[11,] 3.017917e-04 2.783681e-04 -9.721988e-05 1.588716e-05 3.541168e-05
[12,] 4.377122e-04 9.250781e-05 -1.042353e-03 -1.758833e-05 -8.279628e-05
[13,] -2.461815e-05 4.899213e-06 -9.594649e-06 7.205405e-06 1.545697e-06
[14,] -1.210695e-05 -1.445649e-05 1.035925e-04 -1.950742e-06 -9.501641e-06
[15,] -1.348041e-04 -1.409700e-04 2.251624e-03 3.864469e-05 5.655547e-06
[16,] 1.491790e-05 2.395895e-06 4.877120e-06 1.444903e-06 1.633178e-06
[17,] 1.491519e-05 5.879769e-06 -5.122866e-06 4.138441e-06 2.107078e-06
[18,] 7.265659e-05 1.516162e-05 3.810204e-05 7.454811e-05 -1.339967e-07
[19,] 2.393177e-05 6.496462e-06 2.678122e-05 6.475219e-07 1.655893e-05
[20,] 1.656262e-04 2.262409e-05 6.540637e-06 -2.883009e-07 1.194850e-04
[21,] 9.752449e-04 1.099707e-03 -5.263886e-04 4.887686e-05 8.464518e-05
[,11] [,12] [,13] [,14] [,15]
[1,] 2.696624e-05 3.118345e-04 -3.111370e-05 3.392233e-05 8.148975e-05
[2,] 8.506095e-05 1.183130e-04 -7.582109e-06 -3.769431e-06 -3.334130e-05
[3,] -1.143686e-04 -3.986870e-04 7.189813e-05 1.708443e-04 1.285410e-03
[4,] 1.901643e-05 8.661613e-05 7.258930e-06 -4.221119e-06 2.680681e-05
[5,] 3.849050e-05 5.143988e-05 2.720809e-07 1.021996e-06 5.117130e-05
[6,] 3.017917e-04 4.377122e-04 -2.461815e-05 -1.210695e-05 -1.348041e-04
[7,] 2.783681e-04 9.250781e-05 4.899213e-06 -1.445649e-05 -1.409700e-04
[8,] -9.721988e-05 -1.042353e-03 -9.594649e-06 1.035925e-04 2.251624e-03
[9,] 1.588716e-05 -1.758833e-05 7.205405e-06 -1.950742e-06 3.864469e-05
[10,] 3.541168e-05 -8.279628e-05 1.545697e-06 -9.501641e-06 5.655547e-06
[11,] 1.203079e-03 2.397177e-04 2.250930e-05 -4.862213e-05 -3.444602e-04
[12,] 2.397177e-04 1.573092e-02 3.998182e-04 1.646396e-05 -3.911817e-03
[13,] 2.250930e-05 3.998182e-04 1.959533e-04 3.012789e-05 -3.901054e-05
[14,] -4.862213e-05 1.646396e-05 3.012789e-05 3.349500e-04 3.731003e-04
[15,] -3.444602e-04 -3.911817e-03 -3.901054e-05 3.731003e-04 8.215347e-03
[16,] 1.180938e-05 3.600969e-05 5.159459e-06 3.952247e-06 1.847447e-05
[17,] 1.228824e-05 5.307077e-06 5.834441e-06 -9.994379e-07 -1.879787e-05
[18,] 5.741482e-05 -6.236128e-05 2.519566e-05 -6.241884e-06 1.347104e-04
[19,] 2.437687e-05 3.616198e-05 -3.065632e-07 -1.065065e-05 9.997289e-05
[20,] 1.229302e-04 -2.993341e-04 4.703682e-06 -3.360864e-05 2.923468e-05
[21,] 3.496527e-03 1.331281e-03 7.259531e-05 -1.618936e-04 -1.873581e-03
[,16] [,17] [,18] [,19] [,20]
[1,] 1.046518e-05 -4.929012e-07 -1.157290e-05 5.671485e-05 2.728035e-06
[2,] 4.156791e-06 4.162904e-06 2.060207e-05 6.876690e-06 4.574546e-05
[3,] 1.562696e-05 -9.542051e-06 6.765845e-05 2.432035e-05 -3.519962e-05
[4,] 4.198132e-06 6.022992e-06 2.932552e-05 -8.332178e-07 -1.380948e-05
[5,] 4.263664e-06 7.521804e-06 4.395275e-06 3.432684e-05 4.475861e-05
[6,] 1.491790e-05 1.491519e-05 7.265659e-05 2.393177e-05 1.656262e-04
[7,] 2.395895e-06 5.879769e-06 1.516162e-05 6.496462e-06 2.262409e-05
[8,] 4.877120e-06 -5.122866e-06 3.810204e-05 2.678122e-05 6.540637e-06
[9,] 1.444903e-06 4.138441e-06 7.454811e-05 6.475219e-07 -2.883009e-07
[10,] 1.633178e-06 2.107078e-06 -1.339967e-07 1.655893e-05 1.194850e-04
[11,] 1.180938e-05 1.228824e-05 5.741482e-05 2.437687e-05 1.229302e-04
[12,] 3.600969e-05 5.307077e-06 -6.236128e-05 3.616198e-05 -2.993341e-04
[13,] 5.159459e-06 5.834441e-06 2.519566e-05 -3.065632e-07 4.703682e-06
[14,] 3.952247e-06 -9.994379e-07 -6.241884e-06 -1.065065e-05 -3.360864e-05
[15,] 1.847447e-05 -1.879787e-05 1.347104e-04 9.997289e-05 2.923468e-05
[16,] 8.809876e-06 1.230295e-06 5.411681e-06 4.567398e-06 5.781078e-06
[17,] 1.230295e-06 1.696943e-05 1.451396e-05 -3.088723e-06 7.591733e-06
[18,] 5.411681e-06 1.451396e-05 2.624528e-04 2.492117e-06 -1.968992e-07
[19,] 4.567398e-06 -3.088723e-06 2.492117e-06 6.577462e-05 5.869771e-05
[20,] 5.781078e-06 7.591733e-06 -1.968992e-07 5.869771e-05 4.237385e-04
[21,] 3.188145e-05 7.488160e-05 1.777220e-04 7.602458e-05 2.939419e-04
[,21]
[1,] 1.262746e-04
[2,] 2.724515e-04
[3,] -3.714117e-04
[4,] 5.139420e-05
[5,] 1.685144e-04
[6,] 9.752449e-04
[7,] 1.099707e-03
[8,] -5.263886e-04
[9,] 4.887686e-05
[10,] 8.464518e-05
[11,] 3.496527e-03
[12,] 1.331281e-03
[13,] 7.259531e-05
[14,] -1.618936e-04
[15,] -1.873581e-03
[16,] 3.188145e-05
[17,] 7.488160e-05
[18,] 1.777220e-04
[19,] 7.602458e-05
[20,] 2.939419e-04
[21,] 1.384528e-02
$S
AF_GRCh37 Cardiomyopathy Heart_valve_diseases_(aortic_stenosis)
[1,] 0.5006282 0.091405496 0.154347289
[2,] 0.0914055 0.067310004 0.004267099
[3,] 0.1543473 0.004267099 1.146664084
[4,] 0.0827238 0.020425839 0.070998694
[5,] 0.1170491 0.044840690 0.037167858
[6,] 0.3268329 0.244734845 0.011726896
HeartFailure Paroxysmal_tachycardia Pericarditis
[1,] 0.08272380 0.11704914 0.32683288
[2,] 0.02042584 0.04484069 0.24473485
[3,] 0.07099869 0.03716786 0.01172690
[4,] 0.03613277 0.03212762 0.07319881
[5,] 0.03212762 0.06690655 0.15889559
[6,] 0.07319881 0.15889559 0.85456863
$I
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1.10527062 0.129914322 0.0343249836 0.207421816 0.2055565524 0.129827078
[2,] 0.12991432 1.049058888 0.0021870492 0.032163612 0.1604980339 1.013417009
[3,] 0.03432498 0.002187049 1.0343552317 0.068662865 0.0001504762 0.002425674
[4,] 0.20742182 0.032163612 0.0686628646 1.010830695 -0.0040220051 0.032114943
[5,] 0.20555655 0.160498034 0.0001504762 -0.004022005 1.0281760643 0.160713166
[6,] 0.12982708 1.013417009 0.0024256741 0.032114943 0.1607131660 1.049198112
$N
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,] 44008.5 31679.94 14409.16 73188.62 39477.84 14367.18 22904.1 10397.85
[,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16]
[1,] 52028.28 28252.99 10281.42 4720.47 24131.34 12955.93 4714.237 123831.4
[,17] [,18] [,19] [,20] [,21]
[1,] 64850.56 23606.46 35561.04 12813.83 4709.425
$m
[1] 1173569