[R-br] juntar base de dados

joão Rodrigo Castro joaorodrigo2005 em gmail.com
Quarta Setembro 2 10:05:03 BRT 2015


Bom dia listeiros,

Estou com um problema ao juntar duas bases de dados, sendo que já tentei o
merge e não ficou bom (encontrei uma solução na lista, em uma postagem de
fevereiro deste ano, e também não deu certo).
Tenho dados de uma estação meteorológica e dados estimados. Preciso juntar
os dados estimados com os dados observados, ou seja, para o mesmo dia e ano
eu preciso dos fatores "Tmax" e "Tmin" (arquivo "dados") e de "Rad". A
dimensão dos arquivos é diferentes pois há dias em que os dados não foram
registrados na estação. Segue parte dos arquivos.

dados = read.table("CANG_NASA.txt", dec=".", h=T ,col.names = c( "ano",
"jday","Rad","Tmax", "Tmin","prec","vento","td","t2m","rh2m"))
> dados1 = read.table("diario_original_CANG.dat", dec=".", h=T ,col.names =
c( "ano", "mes","jday","Rad","Tmax", "Tmin","stopo","trans","dt"))
>
> head(dados)
   ano jday  Rad Tmax Tmin prec vento   td  t2m rh2m
1 2008    1 24.9 20.7 19.2  0.0   8.5 17.8 20.0 87.3
2 2008    2 22.2 21.6 20.8  0.0   8.3 19.1 21.2 87.5
3 2008    3 31.2 22.0 21.2  0.0   7.4 20.1 21.6 91.7
4 2008    4 29.6 22.3 21.7  0.2   7.0 20.3 22.1 89.6
5 2008    5 18.4 22.7 22.1  7.9   8.0 21.0 22.3 92.7
6 2008    6  9.1 22.3 21.9 32.9  11.5 21.3 22.1 95.4
> head(dados1)
   ano mes jday    Rad Tmax Tmin  stopo trans   dt
1 2008   1    1 32.731 32.2 19.2 44.054 0.743 13.0
2 2008   1    2 30.269 30.7 19.6 44.024 0.688 11.1
3 2008   1    3 23.485 32.0 20.1 43.991 0.534 11.9
4 2008   1    4 21.110 26.5 19.6 43.954 0.480  6.9
5 2008   1    5 26.837 27.3 18.1 43.915 0.611  9.2
6 2008   1    6 26.766 26.5 17.4 43.873 0.610  9.1
> dim(dados)
[1] 2557   10
> dim(dados1)
[1] 1911    9


Obrigado.
-- 
João Rodrigo de Castro
Programa de Pós-Graduação em Meteorologia
Bolsista Laboratório de Agrometeorologia - Embrapa Clima Temperado
Universidade Federal de Pelotas
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150902/0ab39477/attachment.html>


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