[next] [tail] [up]

1 Contagens por quadrats - quadrat counts

Um quadrat é uma sub-região definida a partir do número de linhas e número de plantas nas linhas a serem abrangidas. Os dados de Pinta Preta dos Citrus de Itajobi serão usados como exemplo, cujo mapa está na Figura 1.

Carregando os dados:


  > data(gItajobi)


  > plot(gItajobi, pch = 20)

pict

Figura 1: Incidência de pinta preta dos Citrus no talhão de Itajobi.

A função quadrat.count() implementa a contagem por quadrats. Os argumentos dessa função são:


  > args(quadrat.count)

  function (obj, dx, dy = dx, complet = TRUE, random = FALSE, N = NULL,
      p.quadrats = 1, p.quadrats.random = FALSE, criteria = 0.1,
      death = 1, healt = 0)
  NULL

No help pode ser vista uma descrição detalhada de cada argumento. Os principais são obj que indica o objeto que contém os dados, dx que indica o número de linhas abrangido pelo quadrat e dy que indica o número de plantas na linha abrangido pelo quadrat. Por default dy é igual à dx, bastando informar apenas este argumento.

Com esta função podem ser avaliados todos os possíveis quadrats na dimensão especificado pelo número de linhas e plantas na linha. O objeto retornado é uma lista contendo o número de plantas doentes em cada quadrat, o número total de plantas em cada q e o número nominal de plantas em cada quadrat.

Contagem por quadrats 20 × 20:


  > quadrat.count(gItajobi, 20)

  $y
       [,1] [,2] [,3]
  [1,]    1    1    0
  [2,]    8    1    2
  [3,]    4    0    2
  
  $n
       [,1] [,2] [,3]
  [1,]  381  338   76
  [2,]  400  400  196
  [3,]  400  400  305
  
  $dq
  [1] 400

Observa-se no resultado o efeito da irregularidade das bordas.

Contagem por quadrats 20 × 10:


  > quadrat.count(gItajobi, 20, 10)

  $y
       [,1] [,2] [,3] [,4] [,5] [,6]
  [1,]    0    1    0    1    0    0
  [2,]    1    7    0    1    2    0
  [3,]    3    1    0    0    2    0
  
  $n
       [,1] [,2] [,3] [,4] [,5] [,6]
  [1,]  195  186  173  165   76    0
  [2,]  200  200  200  200  188    8
  [3,]  200  200  200  200  200  105
  
  $dq
  [1] 200

Pode-se selecionar os quadrats aleatoriamente:


  > quadrat.count(gItajobi, 20, 10, ran = T)

  $y
   [1] 0 1 1 1 1 0 0 2 0 3 5 0 6 0 2 0 1 5
  
  $n
   [1] 200 200 200 200 200 200 200 200 200 200 200 200 200 200 200 200 200
  [18] 200
  
  $dq
  [1] 200

Pode-se selecionar os quadrats aleatoriamente e sem desconsiderar os quadrats incompletos:


  > quadrat.count(gItajobi, 20, 10, ran = T, comp = F)

  $y
   [1] 1 1 0 2 0 2 0 0 1 1 3 0 3 0 0 0 6 0
  
  $n
   [1] 200 200 195 158 133  97 170 139 200 200 200 200 200   3  56 150 200
  [18] 140
  
  $dq
  [1] 200

[next] [front] [up]