myscratchpad

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()