-
Notifications
You must be signed in to change notification settings - Fork 14
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
functions for meta analysis #76
base: devel
Are you sure you want to change the base?
Conversation
R/meta.R
Outdated
BPPARAM=bpparam(), verbose=TRUE) { | ||
|
||
# better way to handle these checks? | ||
if(class(gdsobj) == 'SeqVarBlockIterator' && score.cov){ |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
should be if(is(gdsobj, "SeqVarBlockIterator") && score.cov)
Instead of having a separate score.cov logical argument, select whether to return covariance matrices by input object type. A block iterator is used for single variant tests, so only return Score and Score.SE (no covariance) for SeqVarBlockIterator objects.
…; will need separate method
chr=getChromosome(gdsobj, char=TRUE), | ||
pos=getPosition(gdsobj), | ||
stringsAsFactors=FALSE) | ||
if (alleles) { |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I just made a change in GWASTools so this should work now (allele columns were missing previously)
Names for group IDs in SeqVarListIterator objects may be taken from the names of the GRangesList object defining iterations. If that is NULL, it is not straightforward to construct names since variants in each group are not sequential nor even required to be on the same chromosome. Instead, number groups with sequential integers.
No description provided.