Reanalysis of the probabilities of stop codons for each pair of aminoacids

PROBABILITIES IN THE GENES:

>> genes=genoma.gene.sec’;
>> genes=genes(genes~=0);
>> [probstop,casosfavorables,casostotales]=gen2probstop(genes,codigo);
>> subplot(1,2,1)
>> imagesc(sum(probstop{1}(:,:,:),3))
>> colorbar
>> subplot(1,2,2)
>> imagesc(sum(probstop{2}(:,:,:),3))
>> colorbar

THEORETICAL PROBABILITIES:

>> probcodones=hist(codones,1:65);
>> probcodones=probcodones(1:64);
>> probcodones=probcodones/sum(probcodones);
>> probstop_teor=paresaa2probs(codigo,probcodones);
>> clf
>> subplot(1,2,1)
>> imagesc(sum(probstop_teor{1}(:,:,:),3))
>> codones=gen2codones(genes);
>> colorbar
>> subplot(1,2,2)
>> imagesc(sum(probstop_teor{2}(:,:,:),3))
>> colorbar

RATIO:

>> subplot(1,2,1)
>> imagesc(sum(probstop{1}(:,:,:),3)./sum(probstop_teor{1}(:,:,:),3))
>> colorbar
>> subplot(1,2,2)
>> imagesc(sum(probstop{2}(:,:,:),3)./sum(probstop_teor{2}(:,:,:),3))
>> colorbar

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

%d bloggers like this: