[R-br] Frequência dentro de parcelas selecionadas com table()

ASANTOS alexandresantosbr em yahoo.com.br
Terça Novembro 15 14:04:22 BRST 2011


Boa tarde pessoal,

     Estou tentando calcular a frequência de pontos dentro de parcelas 
selecionadas em uma área através da função table() e estou tendo 
dificuldades,
bom consegui calcular para todas as parcelas da área mais não estou 
conseguindo utilizar o comando para selecionar algumas parcelas de 
interesse fiz:

require(splancs)
#Geração dos pontos aleatórios
x <- runif(n=500,min=0, max=20)
y <- runif(n=500,min=0, max=20)

plot(x=x,y=y,axes=F, xlab="", ylab="",
      xlim=c(-1,20),ylim=c(0,21))##Plota os pontos

##Divide em parcelas
segments(x0=rep(0,20),y0=0:20,
          x1=rep(20,20),y1=0:20, col="gray", lty=2)
segments(x0=0:20,y0=rep(0,20),
          x1=0:20,y1=rep(20,20),col="gray", lty=2)

pt0 <- c(0.5,0.5)##Define o centro da primeira parcela

N_S=seq(pt0[1], by=1, length.out=20)##Define o centro das parcelas
L_O=seq(pt0[2],by=1, length.out=20)
grid0 <- expand.grid(N_S=seq(pt0[1], by=1, length.out=20),
L_O=seq(pt0[2],by=1, length.out=20))

##Cria o limite da área
limx<-c(2.5,2.5,18,18)
limy<-c(2.5,18,18,2.5)
lim=cbind(limx,limy)
lim<-as.matrix(lim)
polygon(lim)##Desenha o limite

grid0<-as.matrix(grid0)
m.int<-pip(grid0,lim,out=FALSE)##Seleciona as parcelas dentro do limite
points(m.int,col="blue", pch=10)

##Frequencia em cada parcela
x.parcela <- cut(x,breaks=0:20)
y.parcela <- cut(y,breaks=0:20)
contag <- data.frame(table(x.parcela,y.parcela))

Aqui esta minha dificuldade, se uso o cut(), só consigo dividir a área 
toda, porém não estou conseguindo  selecionar as parcelas localizadas 
dentro do limite,
alguém poderia me dar uma luz,
Obrigado,

-- 
Alexandre DOS SANTOS
Engenheiro Florestal, Msc.
Laboratório de Pragas Florestais
Departamento de Entomologia
Universidade Federal de Lavras
Caixa Postal 3037
37200-000 - Lavras - Minas Gerais - Brasil
Tel: +55 35 92230304



Mais detalhes sobre a lista de discussão R-br