понедельник, 4 января 2010 г.

Графики с планками погрешностей в R


Хороший способ предлагается пакетом gplots. В стандартной поставке его нет, нужно устанавливать дополнительно. После команды
>install.packages("gplots")

R попросит выбрать зеркало и затем сама все скачает и установит.
Загружаем новый пакет:
>library(gplots)

Теперь в R доступна команда plotCI, которая как раз и строит графики "с усами".
>plotCI(x,mean(y), uiw=sd(y), ylim=yy, xlab= (expression("C, "*mu*"M")), ylab="Cell growth, % of control", log="x", bty="l")

По оси X - концентрация в мкМ. Для каждой точки Y есть несколько измерений, поэтому используется среднее - mean(y). uiw это размер усов, т.е. стандартное отклонение sd(y). ylim=yy - жестко задает диапазон значений для шкалы Y (здесь yy вектор 0;120). log(x) указывает, что шкала X - логарифмическая. xlab и ylab подписи осей. bty=l задает тип рамки. L (прямо по форме буквы) обозначает, что рамки будут слева и внизу графика.
Команда построила только точки с усами, если нужно соединить их линией, то достраиваем еще один график, но уже с типом "line":
>plotCI(x,mean(y), ylim=yy, type="line", xlab="", ylab="", log="x", axes=F)

На одном графике можно изобразить множество рядов данных. На рисунке вверху их два.

Комментариев нет:

Отправить комментарий