Please fix this code
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
protein1=input('Please input GenPept code of the first protein: ', 's');
protein2=input('Please input GenPept code of the second protein: ', 's');
disp('Please wait, this can take a while')
for step = 1:50
fprintf('.')
end
seq1 = getgenpept('protein1','SequenceOnly',true);
seq2 = getgenpept('protein2','SequenceOnly',true);
When I input AAA69808 for protein 1, it gives error
3 commentaires
Geoff Hayes
le 5 Jan 2015
Ahmad - if the code is returning an error of some kind, then please copy and paste the full error message to your question. Of particular interest is the line number that generated the error.
Réponses (2)
Image Analyst
le 5 Jan 2015
Modifié(e) : Image Analyst
le 5 Jan 2015
Here's some code you can run to clean things up at the beginning of a script.
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
Try putting that into a quick launch toolbar button called "Clean up" for easy single-click running of it.
2 commentaires
Image Analyst
le 5 Jan 2015
Sorry I can't answer because I don't have the Bioinformatics Toolbox and I don't know what those functions do.
Luuk van Oosten
le 7 Jan 2015
I believe the warning is of no importance to you. Note the difference between different protein annotations, I think that is what is confusing you. It seems that accession numbers are converted to GI nrs.
Have a look here
And here
0 commentaires
Voir également
Catégories
En savoir plus sur Genomics and Next Generation Sequencing dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!