2013-05-01から1日間の記事一覧

R 練習用データをSASサンプルデータより取得する

#こちらでSASのサンプルデータを取得可能 http://ftp.sas.com/samples/A58519 #Rで以下の通り変換すれば、SASのデータをRで使用可能 install.packages("sas7bdat") library(sas7bdat) setwd("ダウンロードしたディレクトリ") data <- read.sas7bdat("拡張子…

Linux Tips (ファイルの移動、PATHの変更)

#ファイル検索してコピー for i in $(find ./ -name 'XXXXX') do echo "copy $i" cp $i /directory/ done #directory, XXXXはその場に応じて変える #ファイルを探す locate file名 #パスの変更の仕方の例 cd /directory/R-3.0.0/.bashrc vi .bashrc #最後の…

R Discriminant Analysis (判別分析)

#線形判別関数はlda(),2次判別の関数はqda() ###################1次判別の例################### x1<-data.frame(hikki=c(68,85,50,54,66),mensetsu=c(65,80,95,70,75)) x2<-data.frame(hikki=c(35,56,25,43,70),mensetsu=c(55,65,75,50,40)) xm1 <- c(mean…

R Linear Regression (線形回帰)

#線形回帰 yi = b0 + b1xi +εiにおいて、誤差項εiの仮定は次の4つ。 #1.εiの期待値は0である(i=1,2,…,n) #2.εiの各分散はすべて等しい(i=1,2,…,n) #3.εi,εuは無相関である(i≠u; i,u=1,…,n) #4.εiは正規分布にしたがう(i=1,…,n) #最小自乗法による推定の良さ…