Comments (15)
Hi Madison. Sorry for my delayed reply. It looks like something is wrong with the loci.file. Can you show me what the contents of the spoofed loci.file look like? Either by pasting it into a code block here, or sending me the file?
Also, are there supposed to be only 2 loci?
from clinehelpr.
I should clarify that if you set loci.file=NULL
then it should make a spoofed one. It should be called "loci_map.txt"
from clinehelpr.
No worries, I appreciate your help. In my files, there are many more than 2 loci. But I'm only really wanting to output the Phi plots, so I thought the loci.file was optional. The spoofed loci.file the program creates the output below.
CHROM POS 0 0 -1 -1
(It doesn't really look like it in the code block here, but there are two columns in the file with three rows: the headers of CHROM and POS, the 0s, and the -1s.)
I have it running in a job script so that it continues on into the creation of the Phi plots, but it always fails at this step, and the loci.file error is the only one I get at that time. Thanks!
from clinehelpr.
I think something might be wrong in a previous step. The spoofed loci.file gets made by looking at the df.list
object, so something might be wrong with df.list. Can you show me what the contents of df.list
looks like? Or if you want, you can send me your BGC output and I can try to see what is going wrong.
from clinehelpr.
I output the contents of df.list (which I input as the object 'bgc.genes') into the text file below. I would send the bgc output, but those files are fairly large for the most part and stored on my HPCC. Thank you again for all your help!
from clinehelpr.
Hi there, I hope you are doing well. I was just checking in to see if you've had a chance to look at the output I sent above as of yet. Thank you!
from clinehelpr.
Hi Madison, Sorry for the slow responses -- Bradley and I both just started new positions, so have not had much time lately :)
I think I have identified the issue and have pushed a potential fix, just testing it now and will get back to you as soon as I have confirmed that it is working.
from clinehelpr.
Ok, I was able to complete an analysis using the spoof loci file with ccd77e4
df.list<-combine_bgc_output("~/bgc_test/moch_bgc/",
prefix="moch",
thin=50)
bgc.outliers <- get_bgc_outliers(df.list,
admix.pop="Sympatric",
popmap="~/bgc_test/moch.popmap",
loci.file=NULL)
alphaBetaPlot(bgc.outliers, showPLOTS = TRUE, saveToFile = FALSE)
from clinehelpr.
Congratulations on your new positions! No worries, thank you for responding. I really appreciate it. I'm going to try my run again, I'll let you know if it works. Thanks again!
from clinehelpr.
Hi again! So I attempted to rerun the program, and am now getting a new error.
Error in 0:((nrow(df.list[[2]][, 1]) - 1)) : argument of length 0
Calls: get_bgc_outliers -> data.frame
Execution halted
Any ideas? Thanks!
from clinehelpr.
Can you save your bgc.genes object with saveRDS() and send it to me? If you don't want to share the data publicly you can email it to [email protected] and I will look into it.
from clinehelpr.
Hi, just saw your email. You can either reduce the file size by re-running combine_bgc_outputs with higher thinning/ burn-in (which might be good anyways so test plotting will run faster), or upload the big file to google drive and send me a link via email.
from clinehelpr.
Hi, I don't mean to open an old thread, but have you figured out the reason behind this error?
Error in 0:((nrow(df.list[[2]][, 1]) - 1)) : argument of length 0
Calls: get_bgc_outliers -> data.frame
Execution halted
Thanks!
from clinehelpr.
Hi Allyson,
My collaborator @tkchafin was mostly working on this bug last time it seems. I don't know if he had fixed the bug or figured out what the cause was. Tyler, did you push an update or create a pull request for this?
Sorry for the issue. We will try to fix it ASAP. Per the comments above, we might need you to send us some files that we can debug it with.
-Bradley
from clinehelpr.
Thank you for the quick response! I will happily send over files for debugging - just let me know what files you need from me and where to send them.
from clinehelpr.
Related Issues (17)
- Hybrid Index HOT 1
- Error with parsing alleles depth HOT 7
- Error HOT 9
- Rename snps$crazy.a and snps$crazy.b to something more intuitive
- Rebuild docker image with updates to vcf2bgc.py HOT 1
- Fix run_bgc.sh estpost LnL file naming HOT 1
- run bgc error in docker HOT 3
- ClineHelpR R package conda install? HOT 1
- jinja2 HOT 2
- R function: combine_bgc_output HOT 2
- GSL error. HOT 2
- phiPlot visualization issue HOT 1
- Assessing the convergence of runs HOT 3
- Add vcf2bgc script
- vcf2bgc.py file HOT 6
- pafscaff HOT 6
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 clinehelpr.