This shows you the differences between two versions of the page.
Next revision | Previous revision Next revision Both sides next revision | ||
howtos:calibration:projecting_into_history [2009/12/01 13:51] shona.weldon created |
howtos:calibration:projecting_into_history [2009/12/01 13:57] shona.weldon |
||
---|---|---|---|
Line 8: | Line 8: | ||
<code> | <code> | ||
- | local ratioFamPerPer[cr,th] = familiesTot[cr,th] / CTotPop[cr,th] | + | local ratioFamPerPer[cr,t19782006] = familiesTot[cr,t19782006] / CTotPop[cr,t19782006] |
- | local ratioFamPerPerPH[cr,t19101977] = lintrend (ratioFamPerPer[cr,th]; time=1910) | + | local ratioFamPerPerPH[cr,t19101977] = lintrend (ratioFamPerPer[cr,t19782006]; time=1910) |
- | familiesPH[cr,t19101977] = extract (CEEFpop[cr,t18511990]; time:1910..1977) * ratioFamPerPerPH[cr,t19101977] | + | familiesPH[cr,t19101977] = extract (CTotPopPH[cr,t18511990]; time:1910..1977) * ratioFamPerPerPH[cr,t19101977] |
</code> | </code> | ||
Line 16: | Line 16: | ||
===== log and exponential lintrend===== | ===== log and exponential lintrend===== | ||
- | Another way to use lintrend and avoid the negatives it to use logarithm and exponential | + | Another way is to use lintrend with logarithm and exponential |
<code> | <code> | ||
- | local ratioFamPerPer[cr,th] = familiesTot[cr,th] / CTotPop[cr,th] | + | local ratioFamPerPer[cr,t19782006] = familiesTot[cr,t19782006] / CTotPop[cr,t19782006] |
- | local ratioFamPerPerPH[cr,t19101977] = exp (max (lintrend (loge (ratioFamPerPer[cr,th]); time=1910), 0)) | + | 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] | familiesPH[cr,t19101977] = extract (CEEFpop[cr,t18511990]; time:1910..1977) * ratioFamPerPerPH[cr,t19101977] | ||
</code> | </code> | ||
Line 28: | Line 28: | ||
<code> | <code> | ||
- | local ratioFamPerPer[cr,th] = familiesTot[cr,th] / CTotPop[cr,th] | + | local ratioFamPerPer[cr,t19782006] = familiesTot[cr,t19782006] / CTotPop[cr,t19782006] |
- | local ratioFamPerPerPHx[cr,t19101977] = changeseq (copyshape (ratioFamPerPer[cr,th]); dim=time, start=1910) | + | local ratioFamPerPerPHx[cr,t19101977] = changeseq (copyshape (ratioFamPerPer[cr,t19782006]); dim=time, start=1910) |
- | ratioFamPerPerPHx[cr,t19101977] = insert (ratioFamPerPer[cr,th]; time->time:1978) | + | ratioFamPerPerPHx[cr,t19101977] = insert (ratioFamPerPer[cr,t19782006]; time->time:1978) |
- | ratioFamPerPerPHx[cr,t19101977] = insert (extract (ratioFamPerPer[cr,th]; time:1978, shrink=on) * \ | + | ratioFamPerPerPHx[cr,t19101977] = insert (extract (ratioFamPerPer[cr,t19782006]; time:1978, shrink=on) * \ |
ratioFamPerPer1910to1978[cr]; time:1910) | ratioFamPerPer1910to1978[cr]; time:1910) | ||