Rで正準相関分析

方法

1．データ読み込み

dat＜－read.csv("data1.csv",header=TRUE)
dat＜－data.frame(dat)
datX＜－dat[,c(1,2)]
datY＜－dat[,c(3,4,5)]
datX＜－scale(datX)
datY＜－scale(datY)

２．正準相関分析

result＜－cancor(datX,datY)

その他いろいろ

・正準相関係数の値
result$cor

・正準変数の係数
n＜－dim(dat)[1]
sqrt(n-1)*sd(datX)*result$xcoef
sqrt(n-1)*sd(datY)*result$ycoef

・第１正準相関得点のプロット、第２正準相関得点のプロット
datU＜－as.matrix(datX)%*%(sqrt(n-1)*sd(datX)*result$xcoef)
datV＜－as.matrix(datY)%*%(sqrt(n-1)*sd(datY)*result$ycoef)
plot(datU[,1],datV[,1])
plot(datU[,2],datV[,2])

・正準負荷量
cor(datX,datU)
cor(datY,datV）

・正準寄与率
第一正準変数のＸの寄与率
mean(cor(datX,datU)[,1]^2)

第一正準変数のＹの寄与率
mean(cor(datY,datV)[,1]^2)

第二正準変数のＸの寄与率
mean(cor(datX,datU)[,2]^2)

第二正準変数のＹの寄与率
mean(cor(datY,datV)[,2]^2)


上記＜－は、本当は半角

日立と三菱重、経営統合するの？しないの？

日経の１面とか、ＮＨＫとか

＜日立・三菱重工＞事業統合へ　１３年春にも新会社
http://headlines.yahoo.co.jp/hl?a=20110804-00000013-mai-bus_all

ってことだけど、

日立と三菱重、それぞれ経営統合報道を否定
http://headlines.yahoo.co.jp/hl?a=20110804-00000277-reu-bus_all

ってニュースもある。

日立と三菱重、経営統合するの？しないの？

---

ま、昨日、東芝、携帯電話事業から撤退…富士通に売却へのニュースを取り上げて、

その日のニュースに

日立製作所、テレビの自社生産から撤退を検討　今年度中に
http://headlines.yahoo.co.jp/hl?a=20110803-00000533-san-bus_all


とあり、原発の影響で、東芝も、日立もですか・・・
と思ってたら、今日、これですからね！
本当に原発の影響は大きい

Ｐ．Ｓ

データセンター売却報の東京電力急落、戻り高値から２５％下げ
http://www.asahi.com/business/stock/kabuto/TKZ201107270001.html

なんて話もあったよね。これも、

東京電力、大型データセンター売却報道は当社が公表したものではない
http://news.nifty.com/cs/economy/stockdetail/tech-b-20110727-201107271707/1.htm

って出たけど・・・
（ちなみに、一番初めにこの話題を提供したのは、日経７月２７日号、１面）