5

r で lme フィット オブジェクトをプロットし、外れ値 ID (studyID) をグラフに表示していますが、プロット オブジェクトでそれらを検索して、これらの ID に自動的にアクセスしたいと考えています。これを行う方法がわかりません。私は多くの分析を行っているため、外れ値の ID 番号について各グラフを実際に調べるのではなく、これを自動的に実行できると役立ちます。

これが私がやっていることの簡単な例です:

fit <- lme(dv ~ studyID + Gender + Group * DOP, random=~1|studyID, cor=corSymm(), na.action="na.omit", method="ML", data=x$data)

require (car)

plotObject <- plot(fit, resid(., type = "p") ~ fitted(.) | Group*DOP, abline = 0, id=.05)

私ができるようにしたいのは、plot ステートメントから生じるグラフ内の外れ値を識別するために使用される ID 番号を格納する plotObject の属性にアクセスすることです。

ありがとうございました。

4

1 に答える 1

7

この情報が実際にプロット オブジェクトに格納されているかどうかはわかりませんが、自分で計算するのは簡単です。から?plot.lme:

 id: an optional numeric value, or one-sided formula. If given as
          a value, it is used as a significance level for a two-sided
          outlier test for the standardized, or normalized residuals.
          Observations with absolute standardized (normalized)
          residuals greater than the 1 - value/2 quantile of the
          standard normal distribution are identified in the plot using
          ‘idLabels’.

だから私は次のようなことを言うでしょう

library(nlme)
fm1 <- lme(distance ~ age, data = Orthodont) # random is ~ age
which(abs(residuals(fm1,type="normalized"))>qnorm(0.975))
## M09 M09 M13 
##  34  35  49 
plot(fm1,id=.05)  ## for comparison

トリックを行うようです。

于 2012-12-11T21:05:56.130 に答える