Rで判別分析（１）－分散が等しい場合：線形判別分析－

内容

グループ１（Ｇ１）、グループ２（Ｇ２）とあって、それぞれのデータがある。
このとき、あるデータＡがあったとき、Ａは、グループ１に属するか、
グループ２に属するかを見たい。


方法
０．インクルード(MASSを読み込む）

library(MASS)

１．データを作成・読み込む

＜＜グループ１、グループ２のデータが１項目しかないとき＞＞

例：ある試験　受験者１０人（うち、５人合格）
　　グループ１（合格者）：100,75,80,90,95
　　グループ２（不合格）：68,65,50,40,55

Rでの入力：
label＜-c(rep(1,5),rep(2,5))
g1＜-c(100,75,80,90,95)
g2＜-c(68,65,50,40,55)
dat＜-c(g1,g2)
dat＜-data.frame(dat)


＜＜グループ１、グループ２のデータが２項目のとき＞＞

例：国語と英語の試験　受験者１０人（うち、５人合格）
　　グループ１（合格者）の国語：100,75,80,90,95
　　グループ１（合格者）の英語：100,75,80,90,95

　　グループ２（不合格）の国語：68,65,50,40,55
　　グループ２（不合格）の英語：68,65,50,40,55

Rでの入力：
label＜-c(rep(1,5),rep(2,5))
kokugo＜-c(100,75,80,90,95,68,65,50,40,55)
eigo＜-c(90,80,76,74,78,68,55,30,20,45)
dat＜-cbind(kokugo,eigo)
dat＜-data.frame(dat)


２．判別分析の実行（以降、１項目でも２項目でも同じ）
result＜-lda(label~.,data=dat)

３．判別分析の結果表示
係数を出す
result$scaling

定数項を出す
-mean(result$means %*% result$scaling)

判別得点の表
cbind(label,predict(result)$class,predict(result)$x)

判別結果
table(label,predict(result)$class)

４．新しいデータの判別

＜＜１項目の場合＞＞
x＜-c(70)
x %*% result$scaling-mean(result$means %*% result$scaling)

＜＜２項目の場合＞＞

x＜-c(87,70)
x %*% result$scaling-mean(result$means %*% result$scaling)

ひろゆき「フジテレビと日本テレビの放送免許が危ない」…外国人の株保有比率が20%以上に

ここの痛いニュース

ひろゆき「フジテレビと日本テレビの放送免許が危ない」…外国人の株保有比率が20%以上に
http://blog.livedoor.jp/dqnplus/archives/1647182.html

(以下斜体は上記サイトより引用）

日本には、電波法という法律があって、携帯電話の電波やテレビの電波などの管理をしています。んで、限られた電波資源を外国が乗っ取ったりしたら大変なので、テレビ局などは、外国人株主の比率をは20%以下にしなければいけないと決められています。（電波法 第五条）

んで、外国人直接保有比率は、証券保管振替機構で見られるんですが、
2011/07/26の数字を見ると、、、

日本テレビ 22.66%
TBS　7.19%
フジテレビ 28.59%
テレビ朝日 14.57%
テレビ東京 1.26%

フジテレビは、外国人直接保有比率が28.59%と、20%を大きく超えちゃってます。


フジ、これはひどすぎでは（＾＾；）
だから、日本じゃない人向けと思われる番組内容になったりする？
（韓日戦とか）

どうにかしないと・・・
やっぱ、第三者割り当て増資なんでしょうか？