У меня есть два набора переменных; ЛСТ и НДВИ. Значения LST зависят от значений NDVI. Диаграмма рассеяния этого показывает типичную форму. Мне нужно получить уравнения для значений верхнего края и значений нижнего края на графике крестцов.
Я сделал наборы данных верхнего квартиля 90% и данных нижнего квартиля 10%. Но результаты не являются удовлетворительными. Я работаю на платформе R
dat<- read.table("clipboard", header=TRUE)
names(dat)
plot(dat$LST, dat$NDVI, xlab=expression(bold(LST)), ylab=expression(bold(NDVI)))
summary(dat)
quantile(dat$NDVI, 0.9)
quantile(dat$NDVI, 0.1)
a <- dat[dat$NDVI > 0.65,]$NDVI
x <- dat[dat$NDVI > 0.65,]$LST
MO <- lm(a ~ x)
newdata <- seq(10,20, length=100)
lines(newdata, predict(MO, data.frame(x = newdata)), lwd=3,
col="blue")
a1<-dat[dat$NDVI < 0.31,]$NDVI
x1<-dat[dat$NDVI < 0.31,]$LST
MO1 <- lm(a1 ~ x1)
newdata1 <- seq(10,20, length=100)
lines(newdata1, predict(MO1, data.frame(x1 = newdata1)), lwd=3,
col="blue")
Я ожидаю линейного уравнения для этих двух переменных на разных ребрах, то есть ниже и выше. Но не только линейное уравнение верхних значений данных.