generate random DNA sequence

1 Ansicht (letzte 30 Tage)
Hsinho
Hsinho am 1 Jul. 2016
Beantwortet: Walter Roberson am 1 Jul. 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?

Antworten (1)

Walter Roberson
Walter Roberson am 1 Jul. 2016
You wrote
for i=20
which should have been
for i = 1 : 20

Kategorien

Mehr zu Genomics and Next Generation Sequencing finden Sie in Help Center und File Exchange

Tags

Produkte

Community Treasure Hunt

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

Start Hunting!

Translated by