Я создал диаграмму рассеяния (несколько групп GRP) с IV=time
, DV=concentration
. Я хотел бы добавить кривые регрессии квантилей (0.025,0.05,0.5,0.95,0.975)
к моему сюжету.
И, кстати, это то, что я сделал для создания рассеянного графика:
attach(E) ## E is the name I gave to my data
## Change Group to factor so that may work with levels in the legend
Group<-as.character(Group)
Group<-as.factor(Group)
## Make the colored scatter-plot
mycolors = c('red','orange','green','cornflowerblue')
plot(Time,Concentration,main="Template",xlab="Time",ylab="Concentration",pch=18,col=mycolors[Group])
## This also works identically
## with(E,plot(Time,Concentration,col=mycolors[Group],main="Template",xlab="Time",ylab="Concentration",pch=18))
## Use identify to identify each point by group number (to check)
## identify(Time,Concentration,col=mycolors[Group],labels=Group)
## Press Esc or press Stop to stop identify function
## Create legend
## Use locator(n=1,type="o") to find the point to align top left of legend box
legend('topright',legend=levels(Group),col=mycolors,pch=18,title='Group')
Поскольку данные, которые я создал здесь, являются небольшим подмножеством моих больших данных, может показаться, что он может быть аппроксимирован как прямоугольная гипербола. Но пока я не хочу называть математическую связь между моими независимыми и зависимыми переменными.
Я думаю, что nlrq
из пакета quantreg
может быть ответом, но я не понимаю, как использовать эту функцию, когда я не знаю отношения между моими переменными.
Я нахожу этот график из научной статьи, и я хочу сделать точно такой же график:
Снова, спасибо за вашу помощь!
Обновление
Test.csv Мне было указано, что мои выборочные данные не воспроизводятся. Вот пример моих данных.
library(evd)
qcbvnonpar(p=c(0.025,0.05,0.5,0.95,0.975),cbind(TAD,DV),epmar=T,plot=F,add=T)
Я также пробовал qcbvnonpar:: evd, но кривая не кажется очень гладкой.