График среднего и 95% доверительный интервал с Hmisc :: xYplot и корректировка оси x - PullRequest
Новая
       5

График среднего и 95% доверительный интервал с Hmisc :: xYplot и корректировка оси x

0 голосов
/ 02 июня 2011

Я пытаюсь построить результаты регрессии, и мне нужно построить расчетные коэффициенты плюс 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"))

Заранее спасибо за любую помощь.

пс .: помощь в установке шрифта и размера заголовка, xlab и ylab будет для меня плюсом.

1 Ответ

2 голосов
/ 02 июня 2011

Извините, ребята, но я нашел ошибки в коде выше.В этом случае я отвечу на свой вопрос?

Вот код, который работает ... 

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

обновление: Это была глупая ошибка, поэтому я не знаю, заслуживаю ли я в конечном итоге за мой ответ.Мы не хотим побуждать людей задавать глупые вопросы и отвечать на них, чтобы получить право голоса (я не сделал этого нарочно, конечно).В любом случае, как и просили, вот что было не так:

Произошел пропущенный парентез после «ярлыков».И я вставил этот отсутствующий парентез в конец кода.Итак, я просто поставил парентез в нужном месте.

...