1日坊主回避

今日は統計赤本と、統計のための行列代数上が届くので後で受け取りに行く。なお財布の中身が足りない模様。

【Rによる統計的検定と推定】
・パッケージのインスコ：utils:::menuInstallPkgs()

・箱ひげ図：boxplot(y ~ x, data)

・2つの平均の検定(t検定)：t.test(data, paired, alternative)
・3つ以上の平均の検定(一元配置分散分析)：oneway.test(y ~ x, data, var.equal)
・2因子以上の平均の検定(二元配置分散分析)：
・2つの分散の検定：var.test(data, alternative)
・3つ以上の分散の検定：bartlett.test(y ~ x, data)

【ネットワーク】
以前ダイクストラ法復習として記事書いたんですけど、ネットワークが分断されてるものに対応してなかったので、対応させた。
とりあえず隣接行列の方のみ。

import numpy as np
inf=np.float('inf')

def dykstra_linkmatrix(link, start):
 end=np.zeros(link.shape[0])+inf
 pre_end=end.copy()
 end[start]=0
 while np.sum(end==pre_end)!=end.shape[0]:
  pre_end=end.copy()
  length=link+end.reshape(link.shape[0], 1)
  length[:, end!=inf]=inf
  if np.min(length)!=inf:
   end[np.argmin(length)%link.shape[0]]=np.min(length)
 return end