How can i compare these two protein sequences?

3 次查看(过去 30 天)
I am trying to compare two cell arrays that contain a protein sequence fragmented up into chunks of two. I am comparing them to see if there are any matches between the sections of two in one sequence to the other, but i keep bumping into this error...
x=length(Prot_H_cell);
y=length(Prot_M_cell);
for ii=Prot_H_cell(1:x);
for jj=Prot_M_cell(1:y);
W(ii,jj)={strcmp(Prot_H_cell(ii),Prot_M_cell(jj))}
end
end
Function 'subsindex' is not defined for values of class 'cell'.- error message

采纳的回答

dbmn
dbmn 2016-12-6
I made a MWE and got a solution. When you see what was wrong you will say "oh man"....
I used the following as an input.
Prot_H_cell = {'AA', 'BB', 'AB', 'AC'};
Prot_M_cell = {'CC', 'AC', 'DA', 'AA'};
This code will fail
x=length(Prot_H_cell);
y=length(Prot_M_cell);
for ii=Prot_H_cell(1:x);
for jj=Prot_M_cell(1:y);
W(ii,jj)={strcmp(Prot_H_cell(ii),Prot_M_cell(jj))}
end
end
And this wont cause the error
x=length(Prot_H_cell);
y=length(Prot_M_cell);
for ii=1:x
for jj=1:y
W(ii,jj)={strcmp(Prot_H_cell(ii),Prot_M_cell(jj))}
end
end
hope this helps.

更多回答(1 个)

KSSV
KSSV 2016-12-6
Try
W{ii,jj}=strcmp(Prot_H_cell(ii),Prot_M_cell(jj)) ;
  5 个评论

请先登录,再进行评论。

Community Treasure Hunt

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

Start Hunting!

Translated by