writing a script that reads a reference sequence, etc.

2 次查看（过去 30 天）

KayLynn 2013-11-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/105978-writing-a-script-that-reads-a-reference-sequence-etc

I am trying to write a script resolvessamples.m

I have this pseudocode written out based on what the script is to do.

Let rs be the set of reference sequences (structures returned by myfastaread()). Let qs be the set of query sequences (structures returned by myfastaread()). Initialize mouseusage to be a vector of zeros (64 elements). Initialize flyusage to be a vector of zeros (64 elements).

Foreach reference sequence as r: Get codonusage of r.Sequence from getcodonusage(). you only need the first output variable. If r.Header contains "Mus musculus" (use strfind()), then add codonusage to mouseusage. Otherwise add it to flyusage.

Foreach ith query sequence q: Get codon usage qusage of q.Sequence from getcodonusage(). you only need the first output variable. Identify which species qusage is closer to using pickcodonusage(). Print "Sequence1 is from a mouse" or "Sequence1 is from a fruit fly", accordingly. (replace "1" with the current value of "i")

I am not sure how to translate this code to matlab code. Any help is appreciated.