Re: Fwd: Junk DNA not so Junky

From: Anders Sandberg (asa@nada.kth.se)
Date: Fri Nov 22 2002 - 16:14:29 MST


On Fri, Nov 22, 2002 at 09:57:23AM -0800, Michael M. Butler wrote:

> I was/am dubious myself. I wonder if the fractality is due to a
> quasifractal distribution of leftovers.

I think the fractal structure can simply be explained by transposomes
rearranging an initially random sequence.

I made the following Matlab program to test it, and it produces
noticeable peaks in the power spectrum. I haven't tested fractality, but
I wouldn't be the least surprised if it showed up. x is the sequence, zz
is the evolution of the power spectrum and xx the evolution of the
sequence.

x=rand(1,1024)<0.5;

zz=[];
xx=[];
for i=1:100
   zz=[zz; abs(fft(x))];
   xx=[xx; x];
   l=100;
   k1=ceil(rand*(1024-l));
   k2=ceil(rand*(1024-l));
   x(k1:(k1+l))=x(k2:(k2+l));
end
subplot(2,1,1)
imagesc(log(zz))
subplot(2,1,2)
imagesc(xx);

-- 
-----------------------------------------------------------------------
Anders Sandberg                                      Towards Ascension!
asa@nada.kth.se                            http://www.nada.kth.se/~asa/
GCS/M/S/O d++ -p+ c++++ !l u+ e++ m++ s+/+ n--- h+/* f+ g+ w++ t+ r+ !y


This archive was generated by hypermail 2.1.5 : Wed Jan 15 2003 - 17:58:19 MST