如何在matlab中生成随机DNA序列

Question

我是 Matlab 的新手，我正在尝试创建随机 DNA 序列生成器，到目前为止，我只找到了一种以相等概率生成 A、C、G、T 的方法，如何为每个字母分配不同的概率？所有概率加起来必须为 1。

total_bp=10;
%open file
SeqLength=100;
bases = repmat('ACGT', 1, SeqLength/4);
for i=1:total_bp
%random DNA sequence
SeqLength=100;
Seq = bases(randperm(SeqLength));
display(num2str(Seq))
end
fclose('all');

Answer 1

定义一个碱基元胞数组

bases = {'A', 'C', 'G', 'T'};

还有一个概率数组；

probs = [0.1,0.5,0.3,0.1];

然后就可以画一个统一的随机数，装箱了：

bin = discretize(rand,[0,cumsum(probs)]);

然后翻译成一个字母：

seq = bases(bin);

总共：

n = 10;
bases = {'A', 'C', 'G', 'T'};
probs = [0.1,0.5,0.3,0.1];
seq = bases(discretize(rand(1,n),[0,cumsum(probs)]));
disp(strjoin(seq)) %Join to single string.