plot - 在残差图中添加文本以识别R中的潜在异常值

时间:2018-06-07 21:33:09

标签: r plot survival

我使用R:

中的plot函数创建了残差图
# 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))

enter image description here

但是,我无法弄清楚如何在边远观察中添加文字。例如,在左上图中,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

0 个答案:

没有答案