(Difference between revisions)
| Revision as of 12:42, 15 September 2006 Peter (Talk | contribs) (→Step1. Finding the relevant parameters with a principle component analysis) ← Previous diff |
Revision as of 12:43, 15 September 2006 Peter (Talk | contribs) (→Step1. Finding the relevant parameters with a principle component analysis) Next diff → |
||
| Line 43: | Line 43: | ||
| ||[[Image:graph97.jpg|400px|]] | ||[[Image:graph97.jpg|400px|]] | ||
| |} | |} | ||
| + | |||
| + | |||
| + | ===Step 2. Performing the Linear Regression Analysis=== | ||
| + | <br> | ||
| + | |||
| + | {|"class=prettytable" cellpadding="15" border="1" style="border-collapse:collapse" | ||
| + | |width="300"|Activate the ''Tools - Principle Component Analysis'' function. A settings window appears.||[[Image:graph93.jpg|400px|]] | ||
| + | |---- | ||
| + | |width="300"|All parameters are selected.||[[Image:graph94.jpg|graph94.jpg]] | ||
| + | |---- | ||
| + | |width="300"|The output log file gives the most valuable results. The first, second and third component explains over 50% of the variance. The equation will be based on those 3 components.||[[Image:graph96.jpg|400px|]] | ||
| + | |---- | ||
| + | |width="300"|Taking all parameters with correlation above 0.70 or below -0.70 the following list of regression parameters is obtained: | ||
| + | * WRSIfin (final index) | ||
| + | * DEFtot (total deficit) | ||
| + | * DEFflow (total deficit during flowering) | ||
| + | * DEFrip (total deficit during ripening) | ||
| + | * ETAveg (ETA during vegetative stage) | ||
| + | * ETArip (ETA during ripening) | ||
| + | * WEXTot (Total Water Excess) | ||
| + | * WEXveg (Water excess during vegetative stage) | ||
| + | |||
| + | |||
| + | ||[[Image:graph97.jpg|400px|]] | ||
| + | |} | ||
| + | |||
| + | |||
| + | |||
| + | |||
Revision as of 12:43, 15 September 2006
Calibrating Yield
The yield function is a statistically derived function relating the water balance parameters (which constitute the outputs of the water balance model) and the other factors (farm inputs, trend) or NDVI with station yield. Once this function has been established, it can be used for early crop yield forecasting.
Although many different equations are possible, the most widely used one is the outcome of a multiple linear regression procedure:
Y = a + b1X1 + b2X2 + b3X3where b1 to b3 are the corresponding X coefficients.
Using an example for Malawi, first it will be established which water balance parameters are good predictors for yield. Then the multiple linear regression will be performed.
The input data file for Malawi contains multiple lines for stations and years. In every line, besides yield, all possibly relevant water balance output parameters are stored. The file can be downloaded here 
Step1. Finding the relevant parameters with a principle component analysis
Activate the Tools - Principle Component Analysis function. A settings window appears. 
All parameters are selected. 
The output log file gives the most valuable results. The first, second and third component explains over 50% of the variance. The equation will be based on those 3 components. 
Taking all parameters with correlation above 0.70 or below -0.70 the following list of regression parameters is obtained:
- WRSIfin (final index)
- DEFtot (total deficit)
- DEFflow (total deficit during flowering)
- DEFrip (total deficit during ripening)
- ETAveg (ETA during vegetative stage)
- ETArip (ETA during ripening)
- WEXTot (Total Water Excess)
- WEXveg (Water excess during vegetative stage)
Step 2. Performing the Linear Regression Analysis
Activate the Tools - Principle Component Analysis function. A settings window appears. All parameters are selected. The output log file gives the most valuable results. The first, second and third component explains over 50% of the variance. The equation will be based on those 3 components. Taking all parameters with correlation above 0.70 or below -0.70 the following list of regression parameters is obtained:
- WRSIfin (final index)
- DEFtot (total deficit)
- DEFflow (total deficit during flowering)
- DEFrip (total deficit during ripening)
- ETAveg (ETA during vegetative stage)
- ETArip (ETA during ripening)
- WEXTot (Total Water Excess)
- WEXveg (Water excess during vegetative stage)
