generate random DNA sequence

3 vues (au cours des 30 derniers jours)
Hsinho
Hsinho le 1 Juil 2016
I would like to randomly generate a set of 7-bases DNA sequences with a weighted probability. I want to have 20 sequences.
I code it as
for i=20
R_1{i}= randsample('ACGT',1,true,[1 0 0 0]);
R_2{i}= randsample('ACGT',1,true,[0.5 0 0.5 0]);
R_3{i}= randsample('ACGT',1,true,[0.25 0.25 0.25 0.25]);
R_4{i}= randsample('ACGT',1,true,[0.5 0 0 0.5]);
R_5{i}= randsample('ACGT',1,true,[0 0 0 1]);
R_6{i}= randsample('ACGT',1,true,[0.5 0 0 0.5]);
R_7{i}= randsample('ACGT',1,true,[0.5 0 0.5 0]);
end
But, it turns out that each cell has no string.
How can I code it correctly? In addition, how can I concatenate them into a single string as a DNA sequence?

Réponses (1)

Walter Roberson
Walter Roberson le 1 Juil 2016
You wrote
for i=20
which should have been
for i = 1 : 20

Catégories

En savoir plus sur Genomics and Next Generation Sequencing dans Help Center et File Exchange

Tags

Produits

Community Treasure Hunt

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

Start Hunting!

Translated by