私のデータのロジスティック回帰曲線をプロットしたいが、プロットしようとすると複数の曲線が得られる。ここに私の最後の試みの絵です:プロットロジスティック回帰曲線in R
ここで私が使用しています関連するコードは次のとおりです。
fit = glm(output ~ maxhr, data=heart, family=binomial)
predicted = predict(fit, newdata=heart, type="response")
plot(output~maxhr, data=heart, col="red4")
lines(heart$maxhr, predicted, col="green4", lwd=2)
私の教授は、次のコードを使用していますが、私はそれを実行しようとすると、私はエラーを取得します最後の行で、xとyの長さが一致しないと言っています:
# fit logistic regression model
fit = glm(output ~ maxhr, data=heart, family=binomial)
# plot the result
hr = data.frame(maxhr=seq(80,200,10))
probs = predict(fit, newdata=dat, type="response")
plot(output ~ maxhr, data=heart, col="red4", xlab ="max HR", ylab="P(heart disease)")
lines(hr$maxhr, probs, col="green4", lwd=2)
助けてください。
編集: mtcarsデータセットを使用して、再現可能なコードを要求したよう
:
fit = glm(vs ~ hp, data=mtcars, family=binomial)
predicted= predict(fit, newdata=mtcars, type="response")
plot(vs~hp, data=mtcars, col="red4")
lines(mtcars$hp, predicted, col="green4", lwd=2)
SO [ロジスティック回帰曲線をプロット上の質問の多くを(持っているのhttp ://stackoverflow.com/search?q = [r] +プロット+ロジスティック)。それらのいずれかが助けてくれますか? – eipi10
リンクeipiを参照してください、またはあなたの例を再現可能にしてください。すでに提供しているコードに適合するデータをシミュレートします。 –
私は最初にそれを検索しようとしましたが、質問のほとんどは私の頭の上にあったものか、私が持っている問題に対処しなかったものでした。問題を再現できるように、組み込みのmtcarsデータセットを使用するコードをいくつか掲載しました。 – cafemolecular