0

回帰の結果をプロットしようとしていて、推定された係数と 95% の信頼区間をプロットする必要があります (実際には、ベイジアン モデルを当てはめているため、95% の信頼区間がありますが、考え方は同じです)。x 軸には、各変数の名前を入れる必要があります。

ここで試してみましたが、うまくいきませんでした。

library(lattice)
library(Hmisc)

betas1=c(0.4271611, -0.4730473, 0.8428530, -0.4356814,  0.2699920, -0.2308315)
quantiles.beta = matrix(c(-0.09015,-1.30900, 0.43290,  -0.85280, 0.02675, -0.56850, 0.93940, 0.37390, 1.28900, -0.03920, 0.52410, 0.07017), ncol=2)

xYplot(Cbind(betas1,quantiles.beta) ~ seq(0, 125, 25), varwidth = TRUE, ylab="Betas",
xlab="Ano", ylim=c(-1.5, 1.5),  scales=list(cex=1.2, x = list(at=seq(0,125, by=25), 
labels = c("PIB per cap.", "democracia", "ginete_park", "educ","patentes", "FDI" ) ),
abline=c(list(h=0), lty="dotted",
col= "grey69"), main="Betas estimados no modelo 1", 
sub="Inclinação das Covariáveis com respectivos 95% intervalo de credibilidade estimados no modelo1"))

助けてくれてありがとう。

ps .: タイトル、xlab と ylab のフォントとサイズの設定を手伝ってくれると助かります。

4

1 に答える 1

2

申し訳ありませんが、上記のコードにエラーが見つかりました。この場合、私は自分自身の質問に答えますか?

ここで動作するコード...

xYplot(Cbind(betas1,quantiles.beta) ~ seq(0, 125, 25)
, varwidth = TRUE, ylab="Betas",xlab="Ano", ylim=c(-1.5, 1.5),  
scales=list(cex=1.2, x = list(at=seq(0,125, by=25), 
labels = c("PIB per cap.", "democracia", "ginete_park", "educ","patentes", "FDI" ))) ,
abline=c(list(h=0), lty="dotted", col= "grey69"), main="Betas estimados no modelo 1", 
sub="Inclinação das Covariáveis com respectivos 95% intervalo de credibilidade estimados no modelo1")

更新：それはばかげた間違いだったので、結局のところ、自分の答えに賛成票を投じるに値するかどうかはわかりません. 私たちは、人々がばかげた質問をして、賛成票を得るためにそれらに答えるように促したくありません (もちろん、私は意図的にこれを行いませんでした)。いずれにせよ、要求されたとおり、何が間違っていたのか：

「ラベル」の後に括弧がありませんでした。そして、この不足している括弧をコードの最後に挿入しました。だから、括弧を適切な場所に置いただけです。

于 2011-06-01T21:18:35.640 に答える