Расчет стандартного отклонения диаграммы рассеяния в R

Я создал точечную диаграмму двух векторов, используя R, в сочетании с линией (используя abline), которая представляет диагональ x=y. Я хочу рассчитать стандартное отклонение точек от диагонали и раскрасить область, которая находится между первым и третьим квантилями.
Я понятия не имею, как это сделать, и буду признательна за любую помощь!!! Заранее спасибо. Хадж.


r standard-deviation
person HEnav    schedule 23.12.2010    source источник
comment
Попробуйте код на странице 84. ablejec.nib.si/R/I2R/DOC /I2R.pdf   -  person Roman Luštrik    schedule 24.12.2010
comment
@Roman Уверен, что это страница 84? Хотя читать чертовски приятно.   -  person Joris Meys    schedule 24.12.2010
comment
Ты прав. Это 86. Это просто дает вам подсказку о том, как нанести на карту некоторые точки и связать с ними аблайн. Возможно, я недостаточно прочитал вопрос. :)   -  person Roman Luštrik    schedule 24.12.2010

Ответы (1)


arrow_upward
3
arrow_downward

Ну, что вы хотите сделать, так это:

# sample data
x <- rnorm(50,0,2)
y <- x+rnorm(50,0,2)

# construct polygons
div <- quantile(y-x,c(0.25,0.75))
x1 <- min(c(x,y))
x2 <- max(c(x,y))


plot(x,y,type="n")
polygon(x=c(x1,x1,x2,x2),y=c(x1+div,(x2+div)[c(2,1)]),col="grey")
abline(0,1)
points(x,y)

Я бы сделал следующее:

qplot(x,y,geom="point") + stat_smooth(method="lm")

Стандартное отклонение, которое вы хотите рассчитать, равно

sd(y-x)

Правильная мера, которую вы, вероятно, ищете, это:

sd(residuals(lm(y~x)))

Вы должны думать в терминах линейной модели y на x, чтобы получить какой-либо значимый результат, если у вас нет очень веских причин не делать этого. Если отношение между x и y не равно 1 на 1, то предположение, что правильная модель такова, не имеет смысла. И если отношение x к y не равно 1 к 1, y-x не будет нормально распределено, и, следовательно, sd будет трудно осмысленно интерпретировать.

person Joris Meys    schedule 24.12.2010
comment
Благодарю вас! это было очень полезно! - person HEnav; 24.12.2010