Question: An R error by using xhmmScripts in XHMM
1
gravatar for Shingo
9 months ago by
Shingo10
Shingo10 wrote:

Recently, I used XHMM to call CNV based on WES data. All the parameters were set to as the same as the XHMM workflow showed in the documentation xhmm_tutorial, and also the same as the filename. But I got an R error as follow when I drew CNV plots. Similarity, I even used the demo data set RUN from the documentation.

library(xhmmScripts)

PLOT_PATH="/Users/CJ/Desktop/RUN/"
OB_PREFICES="/Users/CJ/Desktop/RUN/DATA"
OB_TARGETS_TO_GENES="/Users/CJ/Desktop/RUN/annotated_targets.refseq.loci"
SAMPLE_FEATURES=NULL

XHMM_plots(PLOT_PATH, JOB_PREFICES, JOB_TARGETS_TO_GENES, SAMPLE_FEATURES)

Loading data from '/Users/CJ/Desktop/RUN/DATA' XHMM run...
Read 300 records
Read 300 records
Reading 30 x 300 named matrix
awk: calling undefined function gensub
 input record number 1, file 
 source line number 1
awk: calling undefined function gensub
 input record number 1, file 
 source line number 1
Error in strsplit(targets, ":", fixed = TRUE) : non-string argument

How do I draw the plots with the XHMM results? By the way, I used the hg19 for the reference genome. The chromosome ID begins with "chr". Thanks!

software error R genome • 342 views
ADD COMMENTlink modified 9 months ago by zx87548.8k • written 9 months ago by Shingo10
1

Looks like the package is using certain version of awk with a gensub function, inside one their R functions, see: readNamedMatrix.R#L35

And your version of awk doesn't have gensub. Better to contact the developers.

ADD REPLYlink modified 9 months ago • written 9 months ago by zx87548.8k

I never used this tool but are you on Unix, Windows or MacOS ? Because no pathway start with /Users/ in windows nor linux

ADD REPLYlink written 9 months ago by Bastien Hervé4.5k

Thanks for your reply. I ran the XHMM on linux and draw the plot with the results on local Rstudio.

ADD REPLYlink written 9 months ago by Shingo10

Try to install gawk

ADD REPLYlink written 9 months ago by Bastien Hervé4.5k
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 2.3.0
Traffic: 769 users visited in the last hour