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]]
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]]
El 5 de junio de 2012 13:18, Jaume Tormo <jautorbla@gmail.com> escribió:> Hola, > > No acabo de entender la función, pero debe ser eso. Utiliza todos los > meses de la serie (en lugar de solo los meses del año de interés) para > calcular el valor de heat index. Yo estaba convencido de que lo hacia año a > año. Pero claro, se supone que trabajas con series de datos de muchos años. > Muchas gracias Carlos. > > jaume. > > El 5 de junio de 2012 00:45, Carlos Ortega <cof@qualityexcellence.es>escribió: > > 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 >> > > > > -- > 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 >-- 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]]