thierrygosselin / stackr Goto Github PK
View Code? Open in Web Editor NEWstackr: an R package to run stacks software pipeline
Home Page: http://thierrygosselin.github.io/stackr/
stackr: an R package to run stacks software pipeline
Home Page: http://thierrygosselin.github.io/stackr/
Hi -
I'm having trouble using the vcf2dadi() function. First, the example code gives me the error: "Error in as_data_frame(.) : Not a graph object" when I run like bit of code assigned to id.vcf.
Second, it's not clear how to generate the files needed for assigning an outgroup. Do I need to run stacks populations on both ingroup and the outgroup separately to generate a fasta + sumstats file for each? I've tried various renditions of this but am being told there are 0 common markers between my in- and out- group.
Thanks for your help!
When installing the stackr from Rstudio using the #devtools::install_github("thierrygosselin/stackr")# line, I get the following error
I have the last version of R (x64), Rstudio and Rtools.
Attached you can check the full log.
Hi Thierry,
I am using your stackr to do my pictures, but now I have this error, I run stacks in a cluster, and I am trying to do my summarizes in my local Rstudio, can you help me with this?
Thank you for your time
ustacks.summary<- summary_ustacks(
Summarizing information...
| | 0%, ETA NA
Error: Argument 1 must have names
In addition: Warning message:
In .stackr_parallel(X = sample.name, FUN = summarise_ustacks, mc.cores = parallel.core, :
scheduled cores encountered errors in user code
Hi Thierry,
I have problems running the tsv2bam command.
I got the following error:
For progress, look in the log file:
09_log_files/[email protected]
tsv2bam completed
Moving/Renaming stacks tsv2bam log file:
Merging BAM files with SAMtools to generate a catalog.bam file...
Number of bam files to merge: 0
Error: Tibble columns must have compatible sizes.
SPLIT_VEC
.rlang::last_error()
to see where the error occurred.I tried to figure out the problem, do you have any idea how to solve this?
Best
Maria
Hello,
This package is amazing in its completeness. I'm wondering if there is anywhere that I can find anyone's code, who has used the package for a series of analyses, so that I can get a sense for how to implement the workflow that is provided, and how to integrate the various functions (perhaps from the Benestan paper)? Also, I've read through a whole bunch of the help files for a whole lot of functions, but is there a way that I can get all the functions in a single PDF, sort of as documentation?
Many thanks for this great package.
Ella
Hi Thierry,
I was having a problem with reading in fastq files to generate read depth plots. I kept receiving an error message that vroom could not guess the deliminator.
Setting delim = "/n" in the code for both read_depth_plot and clean_fq seemed to do the trick.
Cheers,
Hi,
I would like to know for what version of R is the library available. I have RStudio 1.4 and I cannot install this package.
Cheers
Maria
Hi Thierry,
I have trouble running the cstacks command line with my samples, but it worked with RADproc. I was already using stackr for fasqcr, but I would like to use stackr for the whole workflow instead of switching from command line to R.
However, I have trouble with run_process_radtags. I am doing RADseq with paired end, so my project.info file is expected to have four columns Barcodes, Individuals, Forward and Reverse. Each of my individual is identified by a specific pair of barcodes (one at each end of the sequence). How am I supposed to give both barcodes with only one column?
Thank you in advance,
Teddy Urvois
Dear Thierry,
I am using your package to input a large amount of data to Colony. Our package is the unique way to pass data from stack to colony.
1- First, i use a awk script to filtred out locus without polymorphisms
awk '$2 > 0 {print}' "$src_root"/14_reassignation/input/batch_2.haplotypes.tsv > "$src_root"/14_reassignation/input/TRIM.haplotypes.tsv
2- I use your R package to feed the haplo2colony fonction
Opening & reading data input file: colony2.dat
Marker 2 has the same ID, 169, as marker 1
Errors in DATA. Insufficient data or incorrect format.
Please check DATA and format and then re-run the program
Program stopped in subroutine StopOnDataError
5- After looking into the colony manual user, in the attached file (i modified the extension)
colony2.txt
line 23, the loci name (header) is duplicated... After deleting all duplicates by hand I got a new (and more severe error). :
Opening & reading data input file: colony2.dat
Reading offspring genotype data...
forrtl: Is a directory
forrtl: severe (30): open failure, unit 10, file /home/jean-baptiste/colony/
Image PC Routine Line Source
colony2s.ifort.ou 0000000000633E04 Unknown Unknown Unknown
colony2s.ifort.ou 00000000006493AB Unknown Unknown Unknown
colony2s.ifort.ou 000000000042AE18 Unknown Unknown Unknown
colony2s.ifort.ou 0000000000423E26 Unknown Unknown Unknown
colony2s.ifort.ou 0000000000401EF6 Unknown Unknown Unknown
colony2s.ifort.ou 0000000000401E7E Unknown Unknown Unknown
colony2s.ifort.ou 00000000006E47A4 Unknown Unknown Unknown
Since colony2 inputs are quite plainfull to build-up, i will be very happy to have any inputs.
JB
Hi,
it would be very helpful if you could provide an example project info file as the instructions are a bit unclear about what it should look like if one has multiple plates/lanes in the paired-end scenario. For example: should the file have another column called LANES?
Similarly, it would be very helpful to see how to add multiple barcodes per sample (ie scenario 5 in section 4.1.1 of Stacks manual) - should they be in the same column separated with some delimiter? Or should there be two columns for barcodes? If the latter, what should their names be?
A simple example file would clarify this :-)
This looks like a very exciting way to "tidy up" a typically messy Stacks workflow. Thanks a lot!
from where to downloads file project.info.turtle.tsv
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.