如何识别使用splom制作的散点图矩阵中的异常值?

时间:2012-12-29 09:16:03

标签: r lattice

我有关于2001 - 2009年期间印度各地的气温,降雨量等天气变量的数据。

要查看这些变量之间的关系,我使用splom创建一个散点图矩阵,并带有以下代码。

   library(lattice)
   library(RColorBrewer)
   splom(~aod_corr[c(5,8,11,14)], lower.panel = panel.splom,
  upper.panel = function(x, y, ...) {
      panel.fill(col = brewer.pal(9, "RdBu")[ round(cor(x, y) * 4 + 5)])
      panel.text(mean(x), mean(y), round(cor(x, y),2), font=2)
  },

 scales = list(x = list( draw = TRUE, cex=0.1)), type = c("g", "p", "smooth"),layout =     c(1, 1))

我现在希望确定每个面板中可以追溯到数据的异常值。我怎么能这样做?

2 个答案:

答案 0 :(得分:5)

似乎这是函数panel.link.splom的工作。看这里 http://procomun.wordpress.com/2011/03/18/splomr/

答案 1 :(得分:4)

根据@Karsten给出的链接,我给你一个与splom对象互动的场景:

  1. 我策划我的splom
  2. 我确定了小组
  3. 我使用link splom
  4. 选择一些积分
  5. 我逃跑以获得结果
  6. 细节:

    1. 我使用格子splom创建一个splom
    2. 相互作用:

      ## I get the focus on the major panel (red color)
      trellis.focus('panel', 1, 1)
      ## I link my panel to the splom
      idx <- panel.link.splom(pch=13, cex=0.6, col='green')
      
    3. 我用鼠标选择一些点。它们将以绿色着色

    4. 逃避以获得焦点

    5. 检查结果

      str(idx)   num [1:5] 111 95 132 96 15