seq_microRNA = fasta read (' mir156.fa');
n = length (seq_mircroRAN);
k = 1;
for i = 1:n;
k=k + 1;
[A,B,D] = rnafold (seq_mirna (i), sequence);
L = length (seq_mirna(i), sequence);
AMFE = B/L;
C = dimercount (seq_mirrna(i), sequence);
GC = C.GC/(L-1);
MFEI(i) = AMFE/GC;
end
阅读(2201) | 评论(0) | 转发(0) |
