Skip to content
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

Error using assocTestAggregate #105

Open
pjuge opened this issue Nov 20, 2023 · 1 comment
Open

Error using assocTestAggregate #105

pjuge opened this issue Nov 20, 2023 · 1 comment

Comments

@pjuge
Copy link

pjuge commented Nov 20, 2023

Hi,

I am using the GENESIS package to perform a SKAT-O analysis on every chromosome, variants being aggregated by genes.
I am using the same code for all of the chromosomes, and it worked perfectly until the chromosome 19.

Running the code

assoc <- assocTestAggregate(iterator,
+                             mod_1,
+                             test="SKATO",
+                             AF.max=0.01,
+                             weight.beta=c(1,25),
+                             verbose = FALSE)

I have the following message:
Error in seqBlockApply(gdsobj, "genotype", function(x) { :
Invalid position in CIndex.

I can provide more details if needed.
Thanks!

@smgogarten
Copy link
Collaborator

It sounds like the code is looking for a position that doesn't exist in the genotype file. I can't tell you more than that without seeing the actual data file, but I suggest creating an iterator object just for chrom 19 with the same genes you're using in your analysis and trying to do a simple task with it (maybe printing out the number of variants?), so you can see on what region it's failing.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants