Return to rogers.exec CVS log

Up to [Research Unix] / researchv10dc / cmd / matlab / exec

researchv10 Dan Cross

exec('d.boug');			// reads data
<g,k> = size(p);		// p is matrix of gene frequencies
wv = ncen/sum(ncen);		// ncen contains population sizes
pbar = wv*p;			// weighted average of p
p = p - ones(g,1)*pbar; 	// deviations from mean
p = sqrt(diag(wv)) * p; 	// weight rows of p by sqrt of pop size
h = diag(pbar); h = h*(eye-h);	// diagonal contains binomial variance: p*(1-p)
r = p*inv(h)*p'/k;		// normalized covariance matrix
eig(r)'

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.