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

[Leonard de Assis]

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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111116/0d2c4fd3/attachment.html>