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