Question

我试图绘制回归的结果，我需要绘制估计的系数加上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字体和大小对我来说是一个加分。

Answer 1

对不起伙计们，但我在上面的代码中发现了错误。在这种情况下，我能回答我自己的问题吗？

这里的代码有效......

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

更新：这是一个愚蠢的错误，所以我不知道我是否应该得到我的答案。我们不想激励人们提出愚蠢的问题并回答他们以获得投票（当然，我没有故意这样做）。在任何情况下，根据要求，这里出了什么问题：

“标签”后遗漏了一些药物。我在代码的最后插入了这个缺失的parentesis。所以，我只是将肠胃外放在适当的位置。

使用Hmisc :: xYplot绘制平均值和95％置信区间并调整x轴

1 个答案: