User Tools

Site Tools


howtos:calibration:projecting_into_history

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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>​
 +
  
  
  
howtos/calibration/projecting_into_history.1259676116.txt.gz ยท Last modified: 2009/12/01 14:01 by shona.weldon