myscratchpad

# オプション状態を表示する
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) # ファイルをダイアログで選ぶ

# クリップボードへ
write.table(dat, "clipboard", sep="\t", row.names=FALSE, quote=FALSE)

# タブ区切りテキスト
write.table(dat, file="dattxt.txt", sep="\t", row.names=FALSE, quote=FALSE)

# タブ区切りcsv
write.csv(dat, file="datcsv.csv")

## fileのところにfile.choose() とすればダイアログで場所を選べる
write.csv(dat, file.choose())

# Excelファイルで保存
library(xlsReadWrite)
write.xls(iris, file.choose())

# クリップボードから
dat <- read.table("clipboard", header=T)

# タブ区切りテキストファイルから
dat <- read.table("dat.txt", header=T)

# タブ区切りcsvファイルから
dat <- read.table("dat.csv", header=T)

# ダイアログで選ぶ。
# file.choose関数はファイルのパスとファイル名を文字列で指定するものだが、日本語は通らない ("デスクトップ") など
dat <- read.table(file.choose(), header=T)
dat <- read.delim(file.choose()) # デフォルトでheader=Tとなってるので便利かもしれない

# Excelのxlsファイルを読む
library("xlsReadWrite")
xls.getshlib() # これを実行しておかないとエラー
dat <- read.xls(file.choose())
## (file.choose, sheet=1) でシートを選べる。デフォルトは1

# Excel2007以降のxlsxファイルを読む
library(xlsx) # javaRuntimeEnvironmentがないとダメみたい。知らない間に入っていたりするけど
file <- system.file("tests", "test_import.xlsx", package = "xlsx") # これを実行しておかないとエラー
dat <- read.xlsx(file.choose(), sheet=1) # シート名を指定しておかないとエラー

## 当然というべきか、Excelファイルを読み込むパッケージは他にもたくさんある。Rコマンダーを使うとシートもダイアログで選択できる。代表的なのはrodbc, gdata, gregmiscなどか。perlのインストールとかめんどそうなのでパスしよう

# webからはクォーテーション内のファイル名をURLにする
dat <- read.table("http://blog.cnobi.jp/v1/blog/user/89d80905c7038b4121822249e9062fba/1255449130",header = T)
library(foreign); dat <- data.frame(read.spss("http://www.indiana.edu/~statmath/stat/all/cfa/values.sav")

# コンマ区切りのときはsep=","とする
dat <- read.table('http://www.ats.ucla.edu/stat/R/faq/hsb2.csv', header=T, sep=",")

# spssのsavファイルから
library(foreign)
dat <- data.frame(read.spss("dat.sav"))
 # 変数ラベルを読み込むときは以下のようにする
x <- read.spss("test2.sav", to.data.frame=T, max.value.labels=99)
head(x)
sapply(x, class)

# webからもいける
dat <- data.frame(read.spss("http://davidakenny.net/dtt/morse_et_al.sav"))


# R以外でspssデータを読んだり出力を開いたりする
## savファイル (データファイル) を開く
## フリーのビューワがある。SPSSのデータ加工のためのもののようで、公式ではない
ViewSav
http://www.asselberghs.dds.nl/stuff.htm
PSPP. SPSSのクローンソフトをつくろうというGNUプロジェクト。まだ開発段階で高度な分析はできないが、とりあえずsavファイルは開ける
http://www.gnu.org/software/pspp/

## spo, spvなどの出力ファイルを開く。
## 公式にフリーのビューワが提供されているのでそれを使う
spoファイル
http://www.spss.co.jp/techsupport/tech-info/spss_statistics_170/000622.html
spvファイル
http://www.spss.co.jp/techsupport/tech-info/spss_statistics_170/000623.html