ランク

• 連立方程式ができたら、それを、行列を用いてと表せば、線形独立かどうかは、のランクを計算するとわかる

m<-matrix(c(1,1,1,0,0,0,
0,0,0,1,1,1,
1,0,0,1,0,0,
0,1,0,0,1,0,
0,0,1,0,0,1,
1,1,1,1,1,1, # 以上６行が周辺度数の制約条件に対応する
1,0,0,0,0,0, # (1,1)セルに値を与える式に対応する行
0,1,0,0,0,0), # (1,2)セルに値を与える式に対応する行ncol=6,byrow=TRUE)

# qr() はランクを求める関数
> 6-qr(m[1:6,])$rank # 第1-6行のみを用いると自由度は2になる
[1] 2
> 6-qr(m[1:8,])$rank # 第1-8行のすべてを用いると自由度は0になる
[1] 0
>