winners.m

其中提到遺傳學的程式碼與應用提供給次淚相向的研究者參考下載

figure(1);clf
winner=[
     4  9 10 13 32 33 37     % 1
     6 16 28 30 31 35  2     % 2
     7  9 29 34 36 39 16     % 3
    13 25 28 29 30 39 21     % 4
     1  3 11 15 17 39 34     % 5
     2  4 15 16 29 33 39     % 6
    20 25 29 37 38 39 28     % 7
     1  6  7 12 20 42 35     % 8
     7 12 19 20 24 36 26     % 9
     2 10 15 23 30 32 36     %10
     6 20 29 35 38 40 41     %11
    12 20 30 32 35 40 34     %12
     4 10 14 23 29 39 36     %13
     4  5 14 21 27 40 12     %14
     8 12 16 26 34 36  5     %15
     5 10 13 15 29 36  1     %16
     2  5  9 16 25 32  7     %17
     6  9 21 27 31 32  2     %18
     5 18 25 26 35 42 29     %19
    19 21 22 25 31 34 13     %20
    13 20 26 30 32 37  4     %21
     7 10 14 17 27 37 28     %22
     6 20 21 26 36 37 27     %23
    12 14 15 24 36 39  1     %24
    10 12 17 21 22 29 25     %25
     5  6  9 24 27 39 15     %26
     1  2  4  7  8 18  9     %27
    15 19 22 23 33 34 20     %28
    11 23 31 33 36 37 34     %29
     7 15 17 18 28 36  8     %30
    21 25 26 33 37 41 20     %31
     7 15 17 27 30 42  5     %32
     2  9 18 20 27 36  4     %33
     2 14 21 23 27 32 35     %34
     3  4 12 19 32 35 29     %35
     6 10 19 22 27 34  7     %36
    18 19 27 34 39 41 40     %37
     8 11 14 24 31 36 23     %38
    10 15 16 26 27 42  3     %39
  ];
 
 % Calculate mean and std
 [M,N]=size(winner);
 seq=1:M;seq=seq';
 wmean=mean(winner')';
 wstd =std( winner')';
 winner=[seq winner];
 % Find average of mean
 haaa=1:length(wmean);haaa=haaa';
 ahhh=cumsum(wmean)./haaa;
 aha =cumsum(wstd)./haaa;
 % TTmean and TTstd will be used to cancel out those numbers outside.
 TTmean=mean(ahhh); TTstd=mean(aha);
 
 figure(1);subplot(211);plot(seq,wmean,seq,wmean,'x',seq,ahhh,seq,TTmean*ones(size(seq)));grid
 figure(1);subplot(212);plot(seq,wstd ,seq,wstd ,'x',seq,aha, seq,TTstd*ones(size(seq)));grid
 
 % Plot winner curves
 Y=winner(:,2:7);Y=Y';X=1:6;X=X';
 figure(2);plot(X,Y)
 figure(1)