Return to rogers.exec CVS log

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

Annotation of researchv10dc/cmd/matlab/exec/rogers.exec, revision 1.1.1.1

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

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