Error following matlab's working with genomes

1 次查看（过去 30 天）
Guy Nir
Guy Nir 2020-9-12
回答： Guy Nir 2020-9-13
I tried following matlab's exapmle of how to memmorry map a genome file ( https://www.mathworks.com/help/bioinfo/examples/working-with-whole-genome-data.html ), but am getting the foollowing error -
Array indices must be positive integers or logical values.
Error in nt2int (line 92)
seq = map((uint8(nt) + 1) - uint8('a'));
I am using 'hg38.fa' from the UCSC genome browser.
Thank you,
Guy
  6 个评论
显示 4更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Guy Nir
Guy Nir 2020-9-13
I added a command to remove any character which is not a letter:
% Remove numbers
numIdx=find(~isletter(charData));
charData(numIdx)='';
% Convert to integers
intData = nt2int(charData);
I think Matlab should consider adjusting some of their functions to whole-genome files (instead of whole chromosmes). Other algorithms use these numbers within the FASTA file for info about these sequences.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Downloads 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by