Tente usar esse script de uma "moving window" feito pelo Prof. Dr.
Jon Grahan.
(http://www.math.umt.edu/graham/stat544/movewin.r). Não sei se o
link ainda está ativo.
Nesse script vc. define o tamanho (m) da janela móvel (nas direções
x e y), se haverá sobreposição entre janelas adjacentes (overlap) e
a variável em foco (velocidade?)
Ele retorna o número de observações em cada janela, a média/desvio
padrão das observações ali encontradas.
Obs.
1. as coordenadas de seus dados deverão estar em UTM (coordenadas
métricas).
2. as classes que vc. possui deverão ser "traduzidas" num
espaçamento regular fixo em x e em y (o tamanho da janela)
FCosta
Em 22/4/2016 16:52, Rodrigo Campos
escreveu:
Olá Leonardo
Boa tarde. Obrigado pelo retorno.
Acredito que possa compreender melhor a partir deste email.
Abaixo, meu arquivo data.frame ("full_data.txt") com dim()
= 298995 observações e 6 variáveis:
cod_barco
lat_dec
long_dec
mes
ano
velocidade
1
892
-23.0086
-44.3131
8
2015
2.4
2
445
-23.0086
-44.3131
8
2015
3.7
3
67
-23.0086
-44.3131
8
2015
3.1
4
223
-23.0086
-44.3131
8
2015
8.9
...
...
...
...
...
...
...
...
...
...
...
...
...
...
298995
678
-23.0234
-44.4543
12
2015
4.5
Gostaria de obter a frequência de observações (número de
amostras) por classe de latitude e longitude.
As classes estão pré-estabelecidas em outro arquivo txt
independente ("lat_long_classes.txt");
no qual tenho 1350 classes de pares de coordenadas
(lat e long), conforme demostrado abaixo (dim()=1350
observações e 2 variáveis):
Lat_class
Long_class
1
-44.4251
-23.6433
2
-44.4201
-23.6433
3
-44.4152
-23.6432
4
-44.43
-23.6388
5
-44.4251
-23.6388
6
-44.4202
-23.6387
...
...
...
...
...
...
1350
-44.4349
-23.6343
Alguém saberia uma dica de como posso gerar a frequência
(n) de observações do meu dataset "full_data.txt" conforme as
classes de pares de latitude e longitude determinadas em outro
arquivo ("lat_long_classes.txt")?
Ou seja, preciso agregar 298995 amostras dentro de 1350
classes estabelecidas em outro data.frame!
Sua explicação não é
detalhada o suficiente para a gente dar uma dica que
funcione sem adaptações. O cut() pode ser usado para
categorizar variáveis numéricas, e o merge (), para
mesclar duas data.frames.
Em Qui 14 abr. 2016, às 09:30, Rodrigo Campos
escreveu:
Amigos usuários do
R,
Bom dia.
Tenho um data.frame ("full_data.txt") com
cerca de 300 mil observações de posições
geográficas (em duas colunas distintas com as
variáveis latitude (y) e longitude (x) em
decimal de grau).
Tenho um segundo arquivo
("lat_long_classes.txt") apenas com classes
de pares de coordenadas (também em
decimal de grau). São 1350 classes de lat long
(como centroides de um grid! Posições
equidistantes!).
Alguém poderia me dar uma dica de como
posso obter as frequências de observações
(constantes em "full_data.txt") por classe de
coordenadas (conforme "lat_long_classes.txt")?