This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | |||
howtos:calibration:projecting_into_history [2009/12/01 14:01] shona.weldon |
howtos:calibration:projecting_into_history [2009/12/01 14:02] (current) shona.weldon |
||
---|---|---|---|
Line 14: | Line 14: | ||
The main danger here is that your linear trend variable, ratioFamPerPerPH[cr,t19101977], goes negative | The main danger here is that your linear trend variable, ratioFamPerPerPH[cr,t19101977], goes negative | ||
- | |||
- | ===== log and exponential lintrend===== | ||
- | Another way is to use lintrend with logarithm and exponential | ||
- | |||
- | <code> | ||
- | local ratioFamPerPer[cr,t19782006] = familiesTot[cr,t19782006] / CTotPop[cr,t19782006] | ||
- | local ratioFamPerPerPH[cr,t19101977] = exp (max (lintrend (loge (ratioFamPerPer[cr,t19782006]); time=1910), 0)) | ||
- | familiesPH[cr,t19101977] = extract (CEEFpop[cr,t18511990]; time:1910..1977) * ratioFamPerPerPH[cr,t19101977] | ||
- | </code> | ||
===== Ratio with linint ===== | ===== Ratio with linint ===== | ||
Line 45: | Line 36: | ||
The main problem with this approach is that it is very flat over pre-history time | The main problem with this approach is that it is very flat over pre-history time | ||
+ | |||
+ | ===== log and exponential lintrend===== | ||
+ | Another way is to use lintrend with logarithm and exponential | ||
+ | |||
+ | <code> | ||
+ | local ratioFamPerPer[cr,t19782006] = familiesTot[cr,t19782006] / CTotPop[cr,t19782006] | ||
+ | local ratioFamPerPerPH[cr,t19101977] = exp (max (lintrend (loge (ratioFamPerPer[cr,t19782006]); time=1910), 0)) | ||
+ | familiesPH[cr,t19101977] = extract (CEEFpop[cr,t18511990]; time:1910..1977) * ratioFamPerPerPH[cr,t19101977] | ||
+ | </code> | ||
+ | |||