Comments (12)
Hello Joon,
It's unexpected that MiTCR output do not have the first line with metadata. The MiTCR parser function tries to skip first line and then starts to read column names and rows. Have you removed somehow the very first line from the MiTCR output and then tried to parse it?
Anyway, you could use this function to parse MiTCR files without skipping the first line:
https://gist.github.com/vadimnazarov/1b438863748fdc3bccf1
from tcr.
First of all, thank you so much for your quick reply.
No. I ran MiTCR with one of the output format options '-level=2'.
It seems that the first line with metadata is added to the output when -level=3. Can these parse functions work with level 3 MiTCR outputs?
from tcr.
I just realised that this file structure could be due to the MiTCR run parameters. If you run MiTCR with -level 3
parameter then parse.mitcr
should work. I will try to fix this issue with the first line in the couple of days.
from tcr.
Yes!
from tcr.
I will post in this topic when I fix this issue to let you know about this update. Please let me know about any issues with parse.mitcr
with -level 3
and/or other bugs, I would greatly appreciate this.
from tcr.
It works well with level 3 data set. Thank you so much for your support.
from tcr.
Dear Joon,
I've just added the feature to parse different levels of MiTCR's output to the parse.mitcr
function. If you would like to use it, please let me know if it works for you.
from tcr.
Thanks a lot. I'll try later and let you know. ^^
from tcr.
Dear Vadim,
I have the same problem as Joon, except that using level 3 output does not fix my problem:
Error in [.data.frame
(df, , make.names(c(.barcodes, .umi.prop, .reads, :
undefined columns selected
My first line contains metadata:
MiTCRFullExportV1.1 HomoSapiens TRB /media/sf_VMdata/mitcr/file1.fastq.gz 1.01.0.3TRBHomoSapiens25both-41123both-14122310010both0.153oneMismatch0.33
Read count Percentage CDR3 nucleotide sequence CDR3 nucleotide quality Min quality CDR3 amino acid sequence V alleles V segments J alleles J segments D alleles D segments Last V nucleotide position First D nucleotide position Last D nucleotide position First J nucleotide position VD insertions DJ insertions Total insertions
5888 0.014855894293045634 TGTGCCAGCAGTGATTCGCGCGTCCAAGGGGCCAACGTCCTGACTTTC HHHHHHHHGHHIHIIHIHGGGGGGGGHHGGGGGGHGGHHGGHHHHHHH 38 CASSDSRVQGANVLTF TRBV6-1_01 TRBV6-1 TRBJ2-6_01 TRBJ2-6 TRBD1_01, TRBD2_02, TRBD2_01 TRBD1, TRBD2 13 21 24 27 7 2 9
4048 0.010213427326468873 TGTGCCAGCAACCCTGCGGGACAGTCAAAGGGACAGCAGTTCTTC HHHHHHHHGHHHHGHHHHGGGGHGHHHHHHHGHGGGHHHHHHHHH 38 CASNPAGQSKGQQFF TRBV6-1_01 TRBV6-1 TRBJ2-1_01 TRBJ2-1 TRBD1_01 TRBD1 9 17 23 34 7 10 17
Do you have any idea what might cause this?
Thanks in advance!
from tcr.
Hi!
Could you please update to the latest development version via
install.packages("devtools")
devtools::install_github("imminfo/tcR", build_vignettes = FALSE, ref = "develop")
and try to parse files again?
from tcr.
Hi, thanks for the quick response!
The installation failed with this message:
library(devtools)
devtools::install_github("imminfo/tcR", build_vignettes = FALSE, ref = "develop")
Downloading github repo imminfo/tcR@develop
Installing tcR
"C:/PROGRA1/R/R-321.1/bin/x64/R" --no-site-file --no-environ --no-save --no-restore CMD INSTALL
"C:/Users/WSUser/AppData/Local/Temp/RtmpS4Ig9x/devtools1e445fbb3a1a/imminfo-tcr-3a98516"
--library="C:/Users/WSUser/Documents/R/win-library/3.2" --install-tests
- installing source package 'tcR' ...
** libs
*** arch - i386
Warnung: Ausführung von Kommando 'make -f "Makevars" -f "C:/PROGRA1/R/R-321.1/etc/i386/Makeconf" -f "C:/PROGRA1/R/R-321.1/share/make/winshlib.mk" CXX='$(CXX1X)
ERROR: compilation failed for package 'tcR'
- removing 'C:/Users/WSUser/Documents/R/win-library/3.2/tcR'
- restoring previous 'C:/Users/WSUser/Documents/R/win-library/3.2/tcR'
Error: Command failed (1)
from tcr.
Hi,
could you please install Rtools (http://cran.r-project.org/bin/windows/Rtools/) for Windows, that is needed to build packages?
from tcr.
Related Issues (20)
- Problem parsing immunoseq for tcR package -- Error: Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : scan() expected 'an integer', got 'nodata' HOT 1
- parsing mixcr 3.0.3 file HOT 4
- Error when trying to parse data HOT 1
- Issue with name fixes HOT 1
- tcR conflict with other R packages?
- New immunoSEQ files won't parse. HOT 4
- pca.segments.2D
- plot "vis.gene.usage" by "Read.proportion" HOT 1
- V-segments usage plot only shows fixed 50 V?
- Differential abundance calculation? HOT 1
- Skipping files while parsing HOT 1
- Chao1 tcR
- signif.digits in vis.heatmap does not work as espected HOT 1
- importing mixcr dataset HOT 1
- TRBV6-8 and TRBV6-9 absence form the HUMAN_TRBV set HOT 1
- Names of gene segments in tcR HOT 1
- vis.rarefaction Warning: Ignoring unknown aesthetics: Group HOT 1
- Sample names in mixcr input HOT 1
- Cosine similiarity HOT 1
- Parsing immunoseq data for tcR package
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 tcr.