myscratchpad

データフレームを行列にして、それからtableにする

dat <- data.frame(c(1,2,3), c(4,5,6))
datm <- as.matrix(dat)
datt <- as.table(datm)
dat
class(dat)
datm
class(datm)
datt
class(datt)

とりあえずリンクだけ貼っておく
そのうち調べよう

Imai, Kosuke, Luke Keele, Dustin Tingley, and Teppei Yamamoto. (2010). ``Causal Mediation Analysis Using R,'' in Advances in Social Science Research Using R, ed. H. D. Vinod, New York: Springer, Forthcoming.
Imai, Kosuke, Luke Keele and Dustin Tingley. ``A General Approach to Causal Mediation Analysis."

mediation パッケージ
http://cran.r-project.org/web/packages/mediation/index.html
http://imai.princeton.edu/software/mediation.html

Sobel test はmultilevelパッケージ？
大御所、Kenny先生の解説: http://davidakenny.net/cm/mediate.htm
サンプルデータと結果: http://davidakenny.net/dtt/mediate.htm

RjpwikiのTipsより。感謝。

文字+番号
paste("Sbj_", 1:10, sep = "")
[1] "Sbj_1"  "Sbj_2"  "Sbj_3"  "Sbj_4"  "Sbj_5"  "Sbj_6"  "Sbj_7"  "Sbj_8"  "Sbj_9"  "Sbj_10"

文字列の分解
x <- paste("Sbj_", 1:10, sep = "")
strsplit(x, "_")
# ベクトルにする
unlist(strsplit(x, "_"))
# データフレームにする
data.frame(strsplit(x, "_"))
## 一部を取り出す
data.frame(strsplit(x, "_"))[2,]

全部数値ならdim(...) でもできる

ついでにメモ
rep(1:3, times = 3) # 1,2,3を3回繰り返す。rep(1:3, 3) も同じ
[1] 1 2 3 1 2 3 1 2 3
rep(1:3, each = 3) # それぞれを3回くりかえす
[1] 1 1 1 2 2 2 3 3 3
## 繰り返しは要素ごとに指定もできる
rep(1:3, times = c(3, 2, 5))
 [1] 1 1 1 2 2 3 3 3 3 3

pwrパッケージで検定力分析をしよう、 と思ったらなぜかCRANから削除されていた。
http://cran.r-project.org/web/packages/pwr/index.html
アーカイブのgzというのはLinux用のものだし、困ったと思ってたら
昔のパッケージやらR本体やらは全部補完されていた。便利だのう
http://cran.us.r-project.org/bin/windows/contrib/2.8/pwr_1.1.zip
http://cran.us.r-project.org/bin/windows/contrib/

参考は
http://www.statmethods.net/stats/power.html
川端先生のサイト

構造方程式モデリングについて

ともあれ、cranのsemパッケージ
http://cran.r-project.org/web/packages/sem/index.html

豊田秀樹 (2009). 共分散構造分析 実践編――構造方程式モデリング――　朝倉書店
semパッケージの使い方、コード、自作の関数が記載されている。

しかし、これを購入する前に無料のところで (たとえ英語でも) がんばろう

Rjipwikiの共分散構造分析のまとめ

南山大学・神谷先生のページ

開発者、John Foxのレクチャー: http://socserv.socsci.mcmaster.ca/jfox/Courses/Brazil-2008/index.html

Fox, J. (2006). Structural-Equation Modeling with the sem Package in R. Structural Equation Modeling, 465-486.
pdfがある。http://socserv.mcmaster.ca/jfox/Misc/sem/SEM-paper.pdf
著者のホームページ: http://socserv.mcmaster.ca/jfox/

Revelle, W. (2009). Using the psych package to generate and test structural models. The Comprehensive R Archive Network. <http://cran.r-project.org/web/packages/psych/vignettes/psych_for_sem.pdf>  (Octpber 15, 2009).
著者のWilliam Revelle (psychパッケージの開発者。偉い心理学者) のページ: http://revelle.net/revelle.html

Personality Projectのsemのセクション
http://www.personality-project.org/r/#sem

Bollen, K. A. Latent variables in psychology and the social sciences. Annual Review of Psychology, 2002, 605-634.
これもpdfがある。http://www.unc.edu/%7Ebollen/bollen_2002.pdf
Rのコードはない