myscratchpad

R のブログを集めたサイト
http://www.r-bloggers.com/

まじめなのからこんなのやこんなのもあって面白い

options(repos="http://cran.md.tsukuba.ac.jp")
install.packages("aspace")      # 三角関数
install.packages("BiplotGUI")   # バイプロット
install.packages("ca")          # 対応分析
install.packages("car")         # 分散分析と回帰分析
install.packages("ctv")         # CRAN Task View
install.packages("doBy")        # カテゴリ別集計
install.packages("foreign")     # ファイル読み書き
install.packages("FactoMineR")  # 探索的解析いろいろ
install.packages("GPArotation") # 因子分析の回転
install.packages("Hmisc")       # カテゴリカルデータ解析 (?)
install.packages("homals")      # 等質性分析
install.packages("lavaan")      # 構造方程式モデル、潜在変数解析
install.packages("lme4")        # マルチレベルモデル
install.packages("MASS")        # 相関の生成とか
install.packages("mvnormtest")  # 多変量正規分布とか
install.packages("multilevel")  # 媒介分析
install.packages("nlme")        # マルチレベルモデル
install.packages("polychor")    # ポリコリック相関
install.packages("psych")       # 因子分析他心理学の解析
install.packages("psychometric")# 心理計量
install.packages("pwr")         # 検定力
install.packages("Rcmdr")       # Rコマンダー
install.packages("reshape2")    # データ整形
install.packages("rgl")         # 3Dプロット
install.packages("sem")         # 構造方程式モデリング
install.packages("vcd")         # カテゴリカルデータ解析


# 読み込まれているパッケージ
search()

# インストールされているパッケージ
library()

# 有用なパッケージをまとめてインストール (CRAN Task View)
library(ctv)
install.views("Psychometrics")
install.views("SocialSciences")
install.views("Multivariate")

# パッケージの更新
options(repos="http://cran.md.tsukuba.ac.jp")
update.packages(ask=F)

# ヘルプ表示法の変更
options(help_type = "text") # テキスト
options(help_type = "html") # html

# Rソースコードの読み込み
 # C:\Program Files\R\R-2.12.0\etc　の　Rprofile.site　を書き換える
source("C:/Documents/all.R") # など

# オプション状態を表示する
options()
# オプションで設定可能な項目
names(options())
? options

# これまでの作業で使用した命令一覧
history(max.show=Inf)
history() # 直近25行
# 命令の保存と読み込み
savehistory(file="mycode.Rhistory")
savehistory(file="mycode.Rhistory")

# コンソールの内容 (命令、出力） を全て保存する
ファイル -> ファイルの保存
## スクロールして見れるものしか保存できない。コマンドでやるにはどうするのかな…

# 作業スペースの保存
save.image("wspace.Rdata)

# オブジェクトの保存と読み込み
save(dat, file="dataobj.Rdata")
load(dat, file="dataobj.Rdata")

# 作業スペースからオブジェクトを削除
rm(dat)

# コンソールへの出力を保存する
sink("Csl.txt", append=T, split=T) # append=Tは追加モード。Fだと上書き, split=Tはコンソールに表示しつつ保存
# 日付・時間が入ったファイル名で保存していく。最初に実行しておくとログを記録し続ける
sink(file = paste(Sys.Date(), gsub(":", "-", format(Sys.time(), "%X")), ".txt", sep = "_"), split = TRUE, append=TRUE)

# 作業環境のオブジェクトをclassつきでリストにする
objlist <- sapply(ls(), function(x) class(get(x)))## classが各1つずつのときはベクトル、複数もつものがあるときはリストで格納
(objnames <- names(objlist)) # オブジェクト一覧
(clstbl <- table(unlist(objlist))) # class一覧

# classが"data.frame"のオブジェクトだけまとめる
dfobjname <- names(objlist)[c(grep("data.frame", objlist))] # "data.frame"オブジェクトの名前をとりだす。"lm", "aov"などとすれば、それぞれのオブジェクトを列挙できる
dflist <- list()
for(i in 1:length(dfobjname)) {
  dflist[[i]] <- get(dfobjname[i]) # リスト化
}
names(dflist) <- dfobjname # リストの要素名をデータフレーム名にする
dflist

# データフレームオブジェクトを調べ、全てcsvで保存する
for ( i in 1:length(dfobjname)) {
write.csv(get(dfobjname[i]), file=paste(dfobjname[i], "_.csv",
sep=""), quote=FALSE)
}



# 出力をhtmlファイルに保存していく。多分あんま使わない
http://www.statmethods.net/interface/output.html

# サンプルデータ
dat <- iris

# 作業環境のオブジェクト一覧
ls()
# 組み込みデータセット一覧
data()


# データ変数名の確認
names(dat)

# データ変数の型を調べる
sapply(dat, class)

# 行数、列数を確認する
dim(dat)
nrow(dat)
ncol(dat)

# データフレームの変数と型の一覧をつくる
data.frame(no=1:ncol(dat), cls=sapply(dat, class)) # classを複数もつ変数にはエラー
# -> 変数の型ごとにデータフレームを分ける

# 因子変数の水準を調べる
levels(dat[,5])
# 因子変数の水準数
nlevels(dat[,5])

# 順序尺度の順序を調べる
x <- ordered(dat[,5]) # irisのSpeciesを順序づける
x["Levels"]

source("http://blog.cnobi.jp/v1/blog/user/89d80905c7038b4121822249e9062fba/1269058643")

# クリップボードからの入力
dat <- read.c()
# クリップボードへ出力
write.c(dat)
# csvから入力
dat <- rdcsv() # ファイルをダイアログで選ぶ

# ファイルから入力
dat <- read.t("filename.txt")
" ファイルへ出力
write.t(dat, "filename.csv")
# csvで出力
wrcsv(dat) # ファイルをダイアログで選ぶ