oncoSNP stops run giving MATLAB function error
0
0
Entering edit mode
4.0 years ago
SJ Basu ▴ 50

Hello everyone,

Hope this message find everyone in good health.

Our lab has generated SNP array data of tumor samples using Illumina CytoSNP 850K beadchip array. Now I am trying to to set up a pipeline to analyze the data and after reading through a few papers, found PennCNV and QuantiSNP/OncoSNP to be popular tools. So I extracted the LRR values from the the data files using bedtools plugins and illumina's EGT normalization files. now the data looks like this:

SNP Name        Chr     Position        Log R Ratio     B Allele Freq
rs4477212       chr1    82154   0.314346        0.00559425
rs3131972       chr1    752721  -0.270539       0.514414
rs12562034      chr1    768448  -0.305087       0.942383
rs12124819      chr1    776546  0.925362        0.00254118

This was given as input to OncoSNP v1.4 (2014). The tool then starts running then stops mid-way giving the following error:

OncoSNP v1.4 (2014). Internal version v3.1
OncoSNP: Performing unpaired tumour sample analysis.
OncoSNP: Running in single sample analysis mode.
OncoSNP: Reading: /home/sangramjit/Desktop/SNParray_analysis/oncosnp_output/txt_4m_gtc/204240030107_R08C01.txt
OncoSNP: Doing GC content correction: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 
Warning: Window size must be less or equal as the number of elements.
> In nanmoving_average at 83
  In initialiseModelParams at 36
  In oncosnp_main at 42
  In oncosnp at 506 
OncoSNP: Loading system options ...
OncoSNP: EM Iteration: 1
Undefined function or variable "loglik".

Error in calcLikelihoodQuick (line 298)

Error in runEMQuick (line 23)

Error in oncosnp_main (line 59)

Error in oncosnp (line 506)


MATLAB:UndefinedFunction

I have installed OncoSNP from this site and corresponding MCR library(R2013b) installer too.

Please help.

Thanks in advance

SNP snp-array array oncosnp matlab • 697 views
ADD COMMENT

Login before adding your answer.

Traffic: 3155 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6