Как вручную создать прямоугольные диаграммы в R с двумя категориями по оси X

Question

Я видел очень похожие вопросы, но ни один из них не совсем соответствовал тому, что я пытаюсь сделать. У меня есть следующее RWE:

n1 = 177
avg1 = 7.508192
sd1 = 5.439677
error1 <- qnorm(.975) * sd1/sqrt(n1)
avg1 - error1
avg1 + error1

n2 = 93
avg2 = 6.713011
sd2 = 3.22479
error2 <- qnorm(.975) * sd2/sqrt(n2)
avg2 - error2
avg2 + error2

Я могу видеть из вычисления avg +/- error степень, в которой доверительные интервалы перекрываются или нет;Тем не менее, я хочу построить эти два набора данных рядом с их средними значениями и доверительными интервалами, чтобы показать их в хорошем графике. Я хочу иметь возможность пометить ось X как «data1» и «data2». Я просмотрел функциональность boxplot и не могу понять, как это сделать, когда я не использую данные как таковые, а скорее вручную генерирую доверительные интервалы. Я не уверен, является ли boxplot подходящей функцией для использования, но она находится на уровне того, что я ищу. Любой совет / места, чтобы посмотреть / простые упущения с моей стороны здесь?

Answer 1

Привет, если вы хотите построить 95% -й доверительный интервал, блокпост в базе R может быть не лучшим. Потому что вам придется использовать усы в качестве доверительного интервала? Вы можете использовать geom_point () в сочетании с geom_errorbar (), см. Пример набора данных, который я создал с вашими значениями:

x = data.frame(
x=rep(c("a","b"),each=2),
data=rep(c("A","B"),2),
avg=rep(c(avg1,avg2),2),
lower=rep(c(avg1 - error1,avg2 - error2),2),
upper=rep(c(avg1 + error1,avg2 + error2),2)
)
ggplot(x,aes(x=x,y=avg,col=data,ymin=lower,ymax=upper)) +
geom_point(position=position_dodge(width=0.1)) +
geom_errorbar(width=0.1,position=position_dodge(width=0.1))

Answer 2

boxplot и bwplot требуют фактических данных для построения графика. Но есть функция с именем bxp (в graphics), которая принимает в качестве входных данных сводные статистические данные и создает из них коробочные диаграммы. Он ожидает список параметров, которые будут выводом вызова boxplot. Итак, вы хотите прочитать как ?bxp, так и value раздел ?boxplot. Как минимум, ваш список должен содержать элементы stats и names. Тем не менее, визуализация сводных данных, которые вы показываете выше, в блокпосте, скорее всего, будет сбивать с толку других, потому что вы используете числа, которые отличаются от того, что обычно показывает блокпост - см. Раздел Details в boxplot.stats для часто используемых определений. .

Ваши данные свидетельствуют о том, что вы хотите построить среднее значение и указать 95% доверительный интервал. Это может быть достигнуто с помощью barplot плюс полосы ошибок. Есть много способов сделать это в R - см. Этот пост: Добавление баров ошибок на барплот

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