2

我使用plotR 中的函数创建了残差图：

# 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))

在此处输入图像描述

但是，我无法弄清楚如何将文本添加到外围观察中。例如，在左上图中，似乎有一个异常值在 inder 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
4

0 回答 0