вращая сетку для построения горизонтальных баров ошибок с помощью Hmisc :: xYplot в R - PullRequest
Новая
       16

вращая сетку для построения горизонтальных баров ошибок с помощью Hmisc :: xYplot в R

3 голосов
/ 18 июня 2011

Я использую xYplot для построения результатов регрессии с ошибками. Тем не менее, xYplot только строит горизонтальные полосы ошибок, и мне нужны вертикальные полосы ошибок. Ища решение, я нашел эту ветку , где кто-то задал примерно такой же вопрос. После некоторых сообщений пользователь, который задал вопрос, говорит, что «я только что обнаружил, что используя xYplot (Hmisc) и вращая сетку окно просмотра (и ярлыки и т. д.) дали мне именно то, что мне нужно ".

Итак, я посмотрел вокруг, как вращать сетку, и обнаружил, что используя библиотеку сетки, pushviewport и т. Д. Вы можете вращать сетку. Тем не менее, мой код не работает. Вот что я попробовал до сих пор: 

estimate=structure(list(coefi = c(-5.08608456607495, -4.17906898422091, 
-2.85696514398422, -3.06968196245069, -2.73660002514793, -1.0017403629931, 
-1.66291850690335, 0.431265159072978, -0.472895611533531, 0.845421348865878, 
-0.437434919008876, 0.269041451577909, -0.233066564595661, 0.0137190330621302, 
-2.94808159763314, 1.9166875739645), lower = c(-8.1895, -6.8485, 
-5.214125, -5.532875, -5.106625, -3.271625, -3.97375, -0.09773, 
-1.340625, 0.415125, -0.86615, 0.02665125, -0.5861, -2.079, -5.626625, 
0.8115125), upper = c(-2.11475, -1.611125, -0.5602375, -0.7309625, 
-0.3721375, 1.259875, 0.7167875, 0.9672875, 0.39035, 1.30025, 
-0.05634125, 0.5115, 0.07237875, 2.14275, -0.3653, 4.202625), 
x = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 
16)), .Names = c("coefi", "lower", "upper", "x"), row.names = c("alpha.1.", 
"alpha.2.", "alpha.3.", "alpha.4.", "alpha.5.", "alpha.6.", "alpha.7.", 
"b.income", "b.democracy", "b.ginete", "b.educ", "b.patent", 
"b.fdi", "b.0", "mu.alpha", "sigma.alpha"), class = "data.frame")

legenda=c(as.character(seq(1970,2000,5)),"PIB_pc", "democ",  "legis", "educ", "patent", "FDI",  "b.0", "mu.ano", "var.ano" )

grid.newpage()
pushViewport(viewport(angle = 90, name = "VP"))
upViewport() 
xYplot(Cbind(coefi,lower, upper) ~x, data=estimate, , varwidth = TRUE, ylab="Betas",
xlab="Inclinação das Covariáveis com respectivos 95% intervalos de credibilidade \n N=409",
ylim=with(estimate,c(min(lower)-.5, max(upper)+.5)),  scales=list(cex=1.2, x = list(at=seq(1,16,     by=1), labels = legenda)) ,abline=c(list(h=0), lty="dotted", col= "grey69"), xlab.top="Adesão ao Tratado de Cooperação de Patentes, 1970-2000", draw.in = "VP")

Буду признателен за любую помощь.

Обновление: были комментарии, указывающие, что код был правильным. Поэтому мне интересно, не могу ли я сообщить, что я хочу, или это ошибка ... Итак, я сейчас выложу изображение вывода моего кода, и вы скажете мне, есть ли код на вашем компьютере дает тот же вывод или другой:

enter image description here

1 Ответ

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

Вы должны напрямую вызвать print для объекта решетки (это упоминается в документации по сетке RShowDoc("grid", package="grid") - "Добавление решетки в сетку"): 

require(grid)
grid.newpage()
pushViewport(viewport(angle = 90, name = "VP"))

print(
    xYplot(Cbind(coefi,lower, upper)~x, data=estimate, , varwidth = TRUE,
        ylab="Betas", xlab="Inclinaçao das Covariáveis com respectivos 95% intervalos de credibilidade \n N=409",
        ylim=with(estimate,c(min(lower)-.5, max(upper)+.5)),
        scales=list(cex=1.2, x = list(at=seq(1,16,by=1), labels = legenda)),
        abline=c(list(h=0), lty="dotted", col= "grey69"),
        xlab.top="Adesao ao Tratado de Cooperaçao de Patentes, 1970-2000",
        draw.in = "VP"
    ),
    newpage=FALSE
)

Rotated xYplot

...