Я создал остаточные графики, используя функцию plot в R:
# residual diagnostics
par(mfrow = c(2,2))
fitresi<-residuals(fit_crm,type="martingale")
fitpred<-predict(fit_crm)
plot(fitresi,ylab="residuals")
plot(fitpred,fitresi,xlab="prediction",ylab="residuals")
lines(lowess(fitpred,fitresi))
plot(BurnData$inf.time,fitresi,xlab="expression",ylab="residuals")
lines(lowess(BurnData$inf.time,fitresi))
plot(BurnData$inf.time,xlab="TNM",ylab="residuals")
lines(lowess(BurnData$inf.time,fitresi))
Однако я не могу понять, как добавить текстк отдаленным наблюдениям.Например, на верхнем левом графике, кажется, есть выброс около 60-70.Есть также некоторые потенциальные выбросы на других графиках.
Как я могу идентифицировать эти наблюдения? Мои данные имеют следующую структуру:
> head(BurnData)
Treatment Gender Race Surface head buttock trunk up.leg low.leg resp.tract type ex.time excision antib.time antibiotic inf.time infection
1 0 0 0 15 0 0 1 1 0 0 2 12 0 12 0 12 0
2 0 0 1 20 0 0 1 0 0 0 4 9 0 9 0 9 0
3 0 0 1 15 0 0 0 1 1 0 2 13 0 13 0 7 1
4 0 0 0 20 1 0 1 0 0 0 2 11 1 29 0 29 0
5 0 0 1 70 1 1 1 1 0 0 2 28 1 31 0 4 1
6 0 0 1 20 1 0 1 0 0 0 4 11 0 11 0 8 1
ID
1 1
2 2
3 3
4 4
5 5
6 6