Hola,

Mirando la expresión analítica de la fórmula de Thornthwaite (
http://en.wikipedia.org/wiki/Potential_evaporation):

Thornthwaite <http://en.wikipedia.org/wiki/C._W._Thornthwaite> equation
(1948)

[image: PET = 16
\left(\frac{L}{12}\right)\left(\frac{N}{30}\right)\left(\frac{10\,
T_a}{I}\right)^\alpha]

Where

[image: PET] is the estimated potential evapotranspiration (cm/month)

[image: T_a] is the average daily temperature (degrees Celsius; if this is
negative, use [image: 0]) of the month being calculated

[image: N] is the number of days in the month being calculated

[image: L] is the average day length (hours) of the month being calculated

[image: \alpha = (6.75 \times 10^{-7}) I^3 - (7.71 \times 10^{-5}) I^2 +
(1.792 \times 10^{-2}) I + 0.49239]

[image: I = \sum_{i=1}^{12} \left(\frac{T_{ai}}{5}\right)^{1.514}] is a
heat index which depends on the 12 monthly mean temperatures [image: T_{ai}]
.[1]<http://en.wikipedia.org/wiki/Potential_evaporation#cite_note-thornthwaite48-0>

El valor de PET depende de "Ta" y de "I" que son valores medios, de las
temperaturas diarias y de las mensuales respectivamente.
Estos valores medios variarán si coges un conjunto de datos (Wichita) de
382 entradas (382 meses) frente a tu subconjunto "meu" de los 12 primeros
valores.

Los resultados no difieren mucho porque los meses, salvo alguna anomalía,
suelen tener comportamientos térmicos (de temperatura) equivalentes.
Equivalentes no quiere decir iguales...


Saludos,
Carlos Ortega
www.qualityexcellence.es

El 4 de junio de 2012 18:08, Jaume Tormo <jautorbla@gmail.com> escribió:

> Hola eRReros.
>
> Os lo explico con un ejemplo:
>
> # Cargamos los paquetes y el ejemplo
>    install.packages("SPEI")
>    library("SPEI")
>    data(wichita)
> # los primeros 12 datos
>    head(wichita,12)
>
> # mi subset de los primeros 12 datos
>    meu<-wichita[c(1:12),]
>    meu
>
> # como veis los valores de TMED son iguales en ambos dataframes.
> # ahora viene el problema
>
> # calculamos thornthwaite para toda la matriz
>    tho <- thornthwaite(wichita$TMED ,lat=37.6475)
> # mostramos los 12 primeros
>    head(tho,12)
>
> # calculamos thornthwaite para mi subset (Los 12 primeros)
>    tho.meu <- thornthwaite(meu$TMED ,lat=37.6475)
>    tho.meu
>
> Como veis los valores de tho y tho.meu son diferentes. Lo malo es que no
> son muy diferentes, parece como si hubiera un problema de decimales o algo
> así.
> ¿Se os ocurre que puede estar pasando?
>
> Por si ayuda:
>  SPEI versión 1.2.
>  RStudio 0.96.122
>  R  2.15
>
> Muchas gracias.
>
> jaume.
> --
> Jaume Tormo Blanes
> IMEM Ramon Margalef.
> Carretera San Vicente del Raspeig s/n
> 03690 San Vicente del Raspeig - Alicante
> www.ua.es
> acercad.wordpress.com
>
>        [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

