ご存じかもしれませんが、正確なロジスティック回帰と呼ばれる手法があります。UCLA の idre は、優れた分析例を提供します。( http://www.ats.ucla.edu/stat/r/dae/exlogit.htm ) 分析例を実行すると、knitr (+ lyx) のエラーに直面します。
Web ページの構文リストの途中で、次の構文を使用するとエラーが発生します。Knitr はチャンクをうまく解釈しているようですが、lyx は pdf 出力に関連するエラーを生成します。
<<>>=
## model with female predictor only'
m.female <- elrm(formula = admit/ntrials ~ female, interest = ~female, iter = 22000, dataset = cdat, burnIn = 2000)
@
lyx は次のようなエラー メッセージを生成します。
Running: pdflatex "logistics.regression.tex" > /dev/null
11:50:28.071: Error while exporting format: pdf2../../../src/support/Systemcall.cpp(270): Systemcall: 'pdflatex "logistics.regression.tex"' finished with exit code 1
一方、lyx のエラー ボックスには、次のメッセージが表示されます。
\end{verbatim}
You need to provide a definition with \DeclareInputText or \DeclareInputMath before using this key
ただし、上記のチャンクの構文は、R ターミナルで問題なく機能します。チャンクに Knitr オプションを追加する必要がありますか? では、それは何ですか？前もって感謝します。
ジョンファ
require(knitr)
require(elrm)
dat <- read.table(text = "
female apcalc admit num
0 0 0 7
0 0 1 1
0 1 0 3
0 1 1 7
1 0 0 5
1 0 1 1
1 1 0 0
1 1 1 6",
header = TRUE)
dat
summary(dat)
dat <- dat[rep(1:nrow(dat), dat$num), -4]
summary(dat)
xtabs(~ female + apcalc, data=dat)
xtabs(~ female + admit, data=dat)
xtabs(~ apcalc + admit, data=dat)
xtabs(~ female + apcalc + admit, data=dat)
x <- xtabs(~admit + interaction(female, apcalc), data = dat)
x # view cross tabs
cdat <- data.frame(female = rep(1:0, 2), apcalc = rep(1:0, each = 2), admit = x[1, ], ntrials = colSums(x))
cdat # view collapsed data set
## model with female predictor only
m.female <- elrm(formula = admit/ntrials ~ female, interest = ~female, iter = 22000, dataset = cdat, burnIn = 2000)
####### the last syntax causes an error in lyx (+ knitr)