2012-12-13 4 views
0

Уважаемые эксперты по стековому потоку,решетки, включают в себя погрешности

Я строю бархаты с lattice. У меня есть две отдельные таблицы, одна со средствами и другими со стандартными ошибками. Как включить стандартные ошибки в барах?

table.df<-means 
table.std<-table_error 
library(lattice) 
library(gridExtra) 
graph1<-barchart(value1~Var1|Var2+Var3,data=table.df, ylab=NULL) 
graph2<-barchart(value2~Var1|Var2+Var3,data=table.df, ylab=NULL) 
grid.arrange(graph1,graph2, nrow=2, ncol=2) 

Любое предложение? Это возможно?

barchart example

+1

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

+1

Это, похоже, также связано с созданием «динамитных сюжетов». Этот дизайн многими считается «плохой формой». Вот это Фрэнк Харрелл взял на себя: http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/DynamitePlots. Поэтому совет будет заключаться в использовании лучшего дизайна сюжета, такого как ящики или скрипки. –

+1

Здесь есть два совершенно разных вопроса, и было бы лучше разделить их на две отдельные должности. –

ответ

0

Хотя я согласен с комментариями, что это, вероятно, плохая идея, вот это решение в любом случае:

x <- rnorm(100, 10, 1) 
ind <- sample(c("one", "two"), 100, TRUE) 

dd <- data.frame(
    m = tapply(x, ind, mean), 
    ind = unique(ind), 
    se = tapply(x, ind, function(x) {sd(x)/sqrt(length(x))}) 
) 

lattice::barchart(ind ~ m, data = dd, xlim = c(8, 12), 
        se = dd$se, 
        panel = function(x, y, ...) { 
        panel.barchart(x, y, ...) 
        panel.segments(y0 = y, x0 = x - dd$se, 
            y1 = y, x1 = x + dd$se, ...) 
        }) 

Imgur

Смежные вопросы