/)[1].gsub(/(.*?)<\/tr>/) do |p| # 対局日の読み込み if ( p =~ / *(.+?)月 *(.+?)日/) then month, day = $1, $2 end # 余分な行の読み捨て row = p.split(/<\/td>/) next if row.size <= 5 # 一局の結果の読み込み case row[4] # 勝敗 when '○', '□' then winLoss = [row[3], row[1]] when '●', '■' then winLoss = [row[1], row[3]] end # テレビ棋戦・女流棋戦・千日手を除外 next if (row[4] =~ /^千?$/ || row[4] =~ '　' || row[5] =~ /女流|レディース|倉敷藤花|鹿島杯/) # 人名表記の「空白」の揺れ、森けい二九段を正規化 winLoss.each do |name| name.sub!(/　雞/, 'けい') name.gsub! (/眞/, '真') name.gsub! (/[ 　]+/, '　') name.gsub! (/^.*アマ$/, 'アマ') end win, loss = *winLoss # 勝敗表の更新 @WinLossHash[win] ||= {'Wins'=>0, 'Losses'=>0} @WinLossHash[win] ['Wins'] += 1 @WinLossHash[loss] ||= {'Wins'=>0, 'Losses'=>0} @WinLossHash[loss] ['Losses'] += 1 if (!@WinLossHash[win] [loss]) then @WinLossHash[win] [loss] = 1 else @WinLossHash[win] [loss] += 1 end end s.close end # # 勝敗表の作成 # @Target = @WinLossHash.keys.reject do |x| @WinLossHash[x]['Wins'] < WinsLimit || @WinLossHash[x]['Losses'] < LossesLimit end @Size = @Target.size @WinLossMatrix = Matrix[* ((0...@Size).collect do |k| ((0...@Size).collect do |i| @WinLossHash[@Target[k]][@Target[i]] || 0 end) end) ] end ####################################### ## レーティング値を計算する ## def rating # レートテーブルの初期値(オール 0) @Rate = Vector[*Array.new(@Size,0)] begin # レーティング・テーブルが予測する勝率行列 @WinRateMatrix = Matrix[* (@Rate.collect do |x| (@Rate.collect do |y| # exp(x) # --------------- # exp(x) + exp(y) 1/(1+Math.exp(y-x)) end) end) ] # 収束補正量の理論値 errorVector = Vector[* ((0...@Size).collect do |k| numerator = 0 # 分子 #--------------------- denominator = 0 # 分母 (0...@Size).each do |i| numerator += @WinLossMatrix[k,i] * @WinRateMatrix[i,k] - @WinRateMatrix[k,i] * @WinLossMatrix[i,k] #------------------------------------------------------ denominator += @WinRateMatrix[i,k] * @WinRateMatrix[k,i] * (@WinLossMatrix[k,i] + @WinLossMatrix[i,k]) end # 0/0 形の場合は、補正なしとする (numerator == 0) ? 0 : numerator / denominator # <=これが問題 end) ] # 収束補正量の平均値を 0 とする(平行移動しても結果は変わらない) errorMean = 0 Array(errorVector).each do |error| errorMean += error end errorVector -= Vector[*Array.new(@Size,errorMean/@Size)] # 収束補正の実行 @Rate += errorVector $stderr.printf "|error| : %5.2e\n", errorVector.r end while (errorVector.r > ErrorLimit) $stderr.printf "\n" self end ####################################### ## 解析結果 ## def result # # 日本将棋連盟の Web のデータを使用した場合 # if (@WinLossHash) then # レーティングの降順インデックスをつくる index = (0...@Size).sort do |x,y| @Rate[y] <=> @Rate[x] end # 結果 printf "以下の#{@Size}棋士を評価しました:\n\n" (0...@Size).each do |k| name = @Target[index[k]] printf("%3d %-14s (%-2s) [%5d] %3d勝%3d敗\n", k+1, name, @@Class[name] , @Rate[index[k]] / Math.log(3) * RateUnit + RateMean, @WinLossHash[name]['Wins'], @WinLossHash[name]['Losses']) end else # # 週刊将棋の集計データを使用した場合 # # 表題 printf "\n\n 【計算結果】注:「対 n 勝数」は勝数の\"予測数(実績数)\"" printf "\n\n No 名前 [ Rate]:" (0...@Size).each do |i| printf(" 対%2d 勝数", i+1) end print "\n" # 内容 (0...@Size).each do |k| # プレーヤー(レート) printf("%2d %-8s[%5d]:", k+1, @Target[k], @Rate[k] / Math.log(3) * RateUnit + RateMean) (0...@Size).each do |i| # 予測(実績) printf("%5.1f(%3d)", @WinRateMatrix[k,i] * (@WinLossMatrix[k,i]+@WinLossMatrix[i,k]), @WinLossMatrix[k,i]) end print "\n" end print "\n" end end ####################################### ## 初期設定 ## def initialize(argv=nil) if (argv) then personalWinLossMatrix(argv) else year2000WinLossMatrix end end end ########################################## ## ## メイン・ルーチン ## ## # 2003年度成績から得るクラス別レーティング RatingTable.new.rating.result # 個人レーティング(暦年2004年1月までの1年のデータに基づくレーティング) RatingTable.new( ['2004/01', '2003/2', '2003/3', '2003/4', '2003/5', '2003/6', '2003/7', '2003/8', '2003/9', '2003/10','2003/11','2003/12']).rating.result