writing a script that reads a reference sequence, etc.
4 次查看(过去 30 天)
显示 更早的评论
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.
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Characters and Strings 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!