myscratchpad

R プロジェクト検索
http://www.rseek.org/

鹿児島大学の書籍公開
これはいい試みだ
http://www.sci.kagoshima-u.ac.jp/~ebsa/

Fox, J. (2008). Applied regression analysis and generalized linear models (2nd ed.). Los Angeles: Sage.

Aiken, L. S., & West, S. G. (1991).    Multiple Regression: Testing and interpreting interactions . Newbury Park, CA: Sage.

Cohen, J., Cohen, P., West, S. G., & Aiken, L. S. (2003). Applied multiple regression/correlation analysis for the behavioral sciences (3rd ed.). Mahwah, NJ: Erlbaum.

Kline, R. B. (2004). Beyond significance testing: Reforming data analysis methods in behavioral research. Washington, DC: American Psychological Association.
http://www.apa.org/books/resources/kline/homepage.html

足立浩平 (2006). 多変量データ解析法――心理・教育・社会系のための入門――　ナカニシヤ出版

豊田秀樹 (1998). 共分散構造分析入門編――構造方程式モデリング――　朝倉書店

Aguinis, H. (2004). Regression analysis for categorical moderators. New York: Guilford Press.

Aguinis, H., Beaty, J. C., Boik, R. J., & Pierce, C. A. (2005). Effect size and power in assessing moderating effects of categorical variables using multiple regression: A 30-year review. Journal of Applied Psychology, 90, 94-107.

宮本友弘・山際勇一郎・田中　敏 (1991) 要因計画の分散分析 において単純主効果検定に使用する誤差項の選択について 心理学研究，62，207-211．

Howell, D. C. (2002). Statistical methods for psychology (5th ed.). Boston: Duxbury.

Winer, Brown, Michels, 1991. Statistical Principles in Experimental Design (3rd edition). McGraw-Hill, New York.

新村秀一 (1995). パソコンによるデータ解析――統計ソフトを使いこなす――　講談社

太郎丸　博 (2005). 人文・社会科学のためのカテゴリカル・データ解析入門　ナカニシヤ出版

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

ともあれ、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のコードはない

OpenOfficeCalcでRを使う、ROOoというのがあるらしい。
http://rcom.univie.ac.at/
のDownloadから、下にスクロールしていくと見つかる

RExcelをいじってみたときも思ったのだが、
普通にRにコマンドを打つほうがよほど楽な気がするんだけどな…
vbaのバグメッセージとか出されてもわからん

aovとかlmとかの出力は半角スペースがいっぱい入って後から使いづらい。
タブ区切りにしたいと思っていたが、xtableというライブラリでできるようだ。
xtableのヘルプの例を貼り付けておく。

sink関数で作業ディレクトリの"RLOG.html"に出力する。
xtable は出力の中のtable部分だけ適用可能。例えば、summary(glm(...))の最後にあるdevianceやAICは出力されない。

library(xtable)
sink(file = "RLOG.html", split = TRUE, append=TRUE)
data(tli)

## aovで使用
cat("\n=======aov======\n\n")
res.aov <- aov(tlimth ~ sex + ethnicty + grade + disadvg, data=tli)
## 普通の出力
cat("===\n普通の出力\n\n")
print(res.aov)
summary(res.aov)
## htmlで出力
cat("===\nhtmlで出力\n\n")
print(xtable(res.aov), type = "html")
print(xtable(summary(res.aov)), type = "html")

## glmで使用
cat("\n\n=======glm======\n\n")
res.glm <- glm(disadvg ~ ethnicty*grade, data=tli, family=binomial())
## 普通の出力
cat("===\n普通の出力\n\n")
print(res.glm)
summary(res.glm)
## htmlで出力
cat("===\nhtmlで出力\n\n")
print(xtable(res.glm), type = "html")
print(xtable(summary(res.glm)), type = "html")

sink()