Question: Error in solveUserSEW0 when using get.targets in TEQC
0
gravatar for Johanna Hasmats
8.5 years ago by
Johanna Hasmats20 wrote:
Dear all, I am trying to load the exome bed file "2.1M_Human_Exome.bed" provided by NimbleGen (http://www.nimblegen.com/downloads/annotation/seqcap_exome/2.1M_Human _Exome _Annotation.zip) as a target by using get.targets in the TEQC package. However, I get this error message: Error in solveUserSEW0(start = start, end = end, width = width) : solving row 175279: range cannot be determined from the supplied arguments (too many NAs) and I cannot seem to come around it. Tried if conditions as well as na.omit/exclude but without results. > traceback() 4: .Call("solve_user_SEW0", start, end, width, PACKAGE = "IRanges") 3: solveUserSEW0(start = start, end = end, width = width) 2: IRanges(start = dat[, startcol], end = dat[, endcol]) 1: get.targets("2.1M_Human_Exome.bed", chrcol = 1, startcol = 2, endcol = 3) > sessionInfo() R version 2.12.2 (2011-02-25) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] TEQC_0.99.2 IRanges_1.8.9 loaded via a namespace (and not attached): [1] Biobase_2.10.0 tools_2.12.2 Could somebody help me out with this? Thanks in advance, Johanna ############################################################ Kindly ? Johanna Hasmats, M.Sc, PhD Student Division of Gene Technology, School of Biotechnology, KTH SciLifeLab Phone: +46(0)73 625 14 60 ? Postal address: Box 1031 171 21 Solna ? Delivery address: KISP (Karolinska Institutet Science Park) Tomtebodav?gen 23B 171 65 Solna Sweden
teqc • 709 views
ADD COMMENTlink modified 8.5 years ago by Hervé Pagès ♦♦ 14k • written 8.5 years ago by Johanna Hasmats20
Answer: Error in solveUserSEW0 when using get.targets in TEQC
0
gravatar for Hervé Pagès
8.5 years ago by
Hervé Pagès ♦♦ 14k
United States
Hervé Pagès ♦♦ 14k wrote:
Hi Johanna, As you can see, line 175279 in the 2.1M_Human_Exome.bed file is different from the other lines: > read.delim("2.1M_Human_Exome.bed", as.is=TRUE, skip=175277, nrows=5, header=FALSE) V1 V2 V3 1 chrY 26179378 26179600 2 chrY 26179986 26180066 3 track name=tiled_region description=NimbleGen Tiled Regions NA NA 4 chr1 58924 59145 5 chr1 59087 59167 I think that this how tracks are separated in a BED file. I don't know the details of how the TEQC::get.targets() function should be used but it seems that it is currently unable to load a BED file that contains more than 1 track. You could still use the 'skip' and 'nrows' argument to load only the track you want though. For example, to load the first track: > aa <- get.targets("2.1M_Human_Exome.bed", chrcol = 1, startcol = 2, endcol = 3, nrows=175278) [1] "read 175278 (non-overlapping) target regions" But that requires that you have a way to find out the right values to use for 'skip' and 'nrows'. A more convenient way to achieve the same result is to load the file with the import() function defined in the rtracklayer package: > library(rtracklayer) > bb <- import("2.1M_Human_Exome.bed", format="bed") > bb RangedDataList of length 2 names(2): target_region tiled_region Then if you want the first track: > as(bb[[1L]], "RangedData") This will give you the same thing as 'aa'. Finally note that TEQC has not been released yet. It's in BioC devel (BioC 2.8) and you are using BioC release (BioC 2.7). This is very likely to bring you troubles. Cheers, H. On 03/03/2011 08:08 AM, Johanna Hasmats wrote: > Dear all, > > I am trying to load the exome bed file "2.1M_Human_Exome.bed" provided by > NimbleGen > (http://www.nimblegen.com/downloads/annotation/seqcap_exome/2.1M_Hum an_Exome > _Annotation.zip) as a target by using get.targets in the TEQC package. > > However, I get this error message: > > Error in solveUserSEW0(start = start, end = end, width = width) : > solving row 175279: range cannot be determined from the supplied arguments > (too many NAs) > > and I cannot seem to come around it. Tried if conditions as well as > na.omit/exclude but without results. > >> traceback() > 4: .Call("solve_user_SEW0", start, end, width, PACKAGE = "IRanges") > 3: solveUserSEW0(start = start, end = end, width = width) > 2: IRanges(start = dat[, startcol], end = dat[, endcol]) > 1: get.targets("2.1M_Human_Exome.bed", chrcol = 1, startcol = 2, > endcol = 3) > >> sessionInfo() > R version 2.12.2 (2011-02-25) > Platform: i386-pc-mingw32/i386 (32-bit) > > locale: > [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United > States.1252 > [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C > > [5] LC_TIME=English_United States.1252 > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] TEQC_0.99.2 IRanges_1.8.9 > > loaded via a namespace (and not attached): > [1] Biobase_2.10.0 tools_2.12.2 > > Could somebody help me out with this? > > Thanks in advance, > > Johanna > > > ############################################################ > Kindly > > Johanna Hasmats, M.Sc, PhD Student > Division of Gene Technology, School of Biotechnology, KTH > SciLifeLab > Phone: +46(0)73 625 14 60 > > Postal address: > Box 1031 > 171 21 Solna > > Delivery address: > KISP (Karolinska Institutet Science Park) > Tomtebodav?gen 23B > 171 65 Solna > Sweden > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor -- Hervé Pagès Program in Computational Biology Division of Public Health Sciences Fred Hutchinson Cancer Research Center 1100 Fairview Ave. N, M2-B876 P.O. Box 19024 Seattle, WA 98109-1024 E-mail: hpages at fhcrc.org Phone: (206) 667-5791 Fax: (206) 667-1319
ADD COMMENTlink written 8.5 years ago by Hervé Pagès ♦♦ 14k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 165 users visited in the last hour