Comments (5)
Hi Priya,
Thanks for your question. In STAAR()
, the eig_sym()
is a c++ function embedded in the Armadillo library that performs eigen-decomposition of the underlying dense symmetric matrix (in this case, the weighted covariance matrix). May I ask what would be the dimension for your eigen-decomposition, i.e. what is the dimension of the G
matrix (you may follow lines 76-102 in the STAAR()
function to obtain this G
matrix from the input genotype
matrix)?
Best,
Xihao
from staar.
Hi, The dimension of the G matrix I was running this function on was 4620 x 285. Thanks, Priya
from staar.
Hi Priya,
Thanks for your reply. The dimension of your G
matrix is manageable for STAAR()
function. In this case, could you please check whether your G
matrix or annotation weights, i.e. w_B
, w_S
, w_A
, contain any missing values? If not, could you please try updating your Rcpp
and RcppArmadillo
packages to the latest version and try it again?
Best,
Xihao
from staar.
Hello, Removing variants with missing values in the annotations fixed the problem for me, thank you so much!! I had one last question: is there a method you would recommend to combine STAAR p-values across multiple studies? Thank you again!
from staar.
Hi Priya,
Thanks so much for letting me know that the problem has been fixed. For your question on combining STAAR p-values, could you please send us an email to have a further discussion, as it may depend on the specific study design of your project?
Thanks again.
Best,
Xihao
from staar.
Related Issues (15)
- error: caught segfault
- Colesky factorization failed HOT 9
- Handling missing annotations HOT 2
- Ranks Used For Generating Phred Scores For Functional Annotations HOT 2
- Using Categorical Annotations with Multiple Categories in STAAR. HOT 2
- Calculation of Phred Scores for Functional/Deleterious Annotations HOT 4
- Question About Assessing Relative Importance of Tests and Annotations HOT 1
- How to generate the genotype file more quickly? HOT 3
- ERROR: compilation failed for package βSTAARβ HOT 3
- Script for calculation of annotation PCs (aPCs) and source file for annotations HOT 1
- docker image of STAAR HOT 1
- How to simply use vcf format file as input HOT 3
- Installation failed HOT 4
- Compilation error HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
π Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. πππ
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google β€οΈ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from staar.