function to find longest string of consecutive ones

Question

SIOUDA abdel hakim 2020-1-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499337-function-to-find-longest-string-of-consecutive-ones

回答： Image Analyst 2020-1-8

hi all, i'm new here. i'v got some difficulties to draft a function to find longest string of consecutive ones (DNA sequence)

for ex: AAAGGTCCATTTTTTTAA

it shows TTTTTTT

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

MaryD 2020-1-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499337-function-to-find-longest-string-of-consecutive-ones#answer_409106

编辑：MaryD 2020-1-8

This is quick made example how it may work. Nbest tells number of letters in the longest string and Lbest tells what letters the string consists of

clear;clc;close all;
str='AAATTTTTCCCCCAATTCCC';
AGTCstr='AGTC';
Strsize=size(str);
AGTCsize=size(AGTCstr);
temp=0;
k=0;
Nbest=0;
for i=1:Strsize(2)
    for j=1:AGTCsize(2)
      h=strcmp(str(i),AGTCstr(j));
      if h==0
          continue
      end
      if h==1
          k=k+1;
           
           if strcmp(temp,str(i))==0
               if k>=Nbest
               Nbest=k;
               Lbest=temp;
              
               end
               k=0;
           end
           temp=str(i);
                      
      end
    end
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Image Analyst 2020-1-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499337-function-to-find-longest-string-of-consecutive-ones#answer_409109

Use findgroups() and regionprops():

s = 'AAAGGTCCATTTTTTTAA' - 'A'
g = findgroups(s)
for k = 1 : max(g)
	props = regionprops(g == k, 'Area');
	largestAreas(k) = max([props.Area]);
	fprintf('The largest stretch of group %d is %d.\n', k, largestAreas(k));
end