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

[ASANTOS]

Obrigado Prof. Paulo,

        Vou converter para objetos Spatial do sp e tentar a função 
overlay, estava lendo a função é bem isso que preciso. Já tinha tentando 
quadrat.count() do Rcitrus e quadratcount() do spatstat, mais o problema 
é que estou começando o estudo com simulações e as janelas são 
quadradas, porém quando for para campo e como as áreas possuem contornos 
irregulares, já estou prevendo problemas com as inúmeras parcelas que se 
localizam na borda e não vão formar quadrats de mesmo tamanho, por isso 
preciso estudar abordagem prevendo isso.

Redobrados agradecimentos,

Alexandre

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


Em 16-11-2011 16:34, Leonard de Assis escreveu:
> Prof Paulo
>
> Aproveitando o email, já que tenho uma dúvida no mesmo assunto:
>
> 1. Qual seria uma bibliografia sugerida para se entender essa 
> abordagem "geoestatística" da análise de dados
> 2. Algum curso (presencial ou EAD) previsto para ano que vem?
>
> Pergunto isso pois tenho muitas idéias de aplicação destes conceitos, 
> mas me falta referência (e até incentivo) para pesquisar mais sobre o 
> assunto e formar alguma opinião.
> []s
> Leonard de Assis
> assis<dot>  leonard<at>  gmail<dot>  com
>
> Em 16/11/2011 15:25, Paulo Justiniano escreveu:
>> Alexandre
>>
>> sem fornecer codigo especifico aqui vao algums ideias de por ojnde eu 
>> comecaria
>>
>> algumas alternativas caso nao queira programar isto:
>>
>> 1. represetnar os dados no formato sp. Um como SpatialPoints
>> e outro como SpatialGrid
>> Depois a operaçõa de overlay (over() ) com tabloe vai te fornecer as 
>> contagens
>>
>> 2.
>> este problema é o mesmo do chamado "quadrat counts"
>> que é tratado em pacotes como spatstat (oficial)
>> ou Rcitrus (nao oficial)
>>
>> aqui vai um exemplinho usando o 1o
>>
>> require(spatstat)
>> ## Carregando um conjunto de dados
>> data(swedishpines)
>> X <- swedishpines
>> plot(X)
>> summary(X)
>> ## Contagem de Quadrats
>> Q <- quadratcount(X, nx = 4, ny = 3)
>> Q
>> plot(X)
>> plot(Q, add = TRUE, cex = 2)
>>
>>
>>
>> On Tue, 15 Nov 2011, ASANTOS wrote:
>>
>>> 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
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça 
>>> código mínimo reproduzível.
>>>
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.


-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111116/d412aac0/attachment.html>