Summary of the invention
The purpose of the present invention is to provide it is a kind of reproduction initial stage prediction tomato whether the method for the underproduction, this method can be
Tomato reproduction initial stage predict tomato whether the underproduction.
To achieve the above object, the technical solution of the present invention is as follows:
A method of reproduction initial stage prediction tomato whether the underproduction, comprising the following steps:
1) as independent variable and variable combination is arranged in the physical signs of the tomato leaf of flower bud differentiation period, according to different
The combination of variable and dependent variable establish Logic Regression Models respectively;The physical signs includes that object element content, organic matter contain
Amount;The object element includes metallic element and nonmetalloid, the metallic element be calcium, iron, potassium, magnesium, sodium, in zinc extremely
Few one kind;The nonmetalloid is nitrogen and/or phosphorus;The dependent variable be whether the underproduction;
2) using the physical signs of the tomato leaf in training set sample as independent variable, using whether the underproduction is used as because becoming
Amount, is trained the Logic Regression Models, determines the logistic regression coefficient of model;The training set sample includes normal group
Sample and nutritional deficiency group sample;
3) physical signs of the tomato leaf in verifying collection sample is brought into as independent variable and logic has been determined in step 2)
The model of regression coefficient, obtains prediction result;The actual result of the prediction result of different models and verifying collection sample is compared, screening
The highest model of prediction accuracy out;The verifying collection sample includes normal group sample and nutritional deficiency group sample;
4) it is predicted using the model that step 3) filters out whether tomato to be predicted the underproduction occurs.
The present invention screens the physical signs of tomato leaf in tomato reproductive development " flower bud differentiation period " stage early stage,
Selected section physical signs, then establish index and whether the model of the underproduction, predicted using model tomato whether the underproduction.This hair
It is bright that for tomato, whether the underproduction using the physical signs of tomato itself as marker establishes the prediction model of superelevation accuracy rate.
The present invention just establishes Production Forecast Models at tomato reproduction initial stage, facilitate peasant household or enterprise can also intervene yield when
Phase carries out effective operation, avoids the generation of the underproduction, be more of practical significance.Whether of the invention predicts tomato at reproduction initial stage
The method of the underproduction is that the following tomato underproduction early warning and high yield optimal conditions have all established theoretical basis.
The variable combination includes that entire variable combines, and the entire variable combination is by all indexs in the physical signs
It is obtained as independent variable.Entire variable combination is all calculated all factors relevant to result, because the factor of consideration is more,
The result theoretically obtained can be more accurate.
The variable combination includes AIC optimum combination, and the AIC optimum combination is that the AIC value of corresponding model is the smallest
Variable combination.The smallest set of variables of AIC value build jointly vertical model for prediction tomato whether the underproduction accuracy rate it is higher.
The variable combination includes that significant difference combines, and the significant difference combination includes normal group sample and nutritional deficiency
The variable of group sample physical signs significant difference.Influence of the physical signs variable for yield with significant difference may be compared with
Greatly, therefore by significant difference combination it takes into account.
The tomato leaf is plant top the 4th leaf down.Since first and second leaves at the top of plant are new
Leave piece can only guarantee that the locating growth phase of plant is the same and cannot be guaranteed the development of all plant in practical planting process
Time phase is just the same, therefore cannot select newborn blade.In addition, also there are no complete for the third piece leaf at the top of plant down
Sizing, length and width are also inconsistent.And from the 4th leaf of plant down for the tomato of contemporaneity, form and state are
Compare unification, the physiological phenomenon of aging will not occur.
Embodiment 1
The present embodiment reproduction initial stage prediction tomato whether the underproduction method the following steps are included:
(1) tomato is cultivated
To tomato wild type " beautiful spring " seed bought in scientia Agricultura Sinica research institute (Beijing Haidian) seed retail department
(Lycopersicon esculentum Mill ev.Lichun) is cultivated, after seedling length to one heart stage of two leaves, field planting
It is 20cm to diameter, in the big flowerpot of a height of 20cm, plants 152 basins (strain) altogether, be then placed in temperature-controlling green house.When experiment
Between from July to October, greenhouse daytime temperature be 25~30 DEG C, nocturnal temperature be 18~20 DEG C, relative humidity be 60~
80%, daily illumination 16h.
Matrix in flowerpot is the turfy soil that mass ratio is 1:1 and vermiculite head.In 152 plants of plantation, wherein 76 plants long-term
Watering, 76 plants of pouring 1/4MS nutrient solutions are (in 1/4MS nutrient solution: calcium nitrate 945mg/L, potassium nitrate 607mg/L, ammonium phosphate
115mg/L, magnesium sulfate 493mg/L, iron salt solutions 2.5mg/L, microelement 5mg/L, pH=6.0)
(2) measurement of the content of physical signs
Plant grows to flower bud differentiation period, chooses at the top of plant the 4th leaf down, in an oven by leaf at a temperature of 180 DEG C
Sheet material is dried completely, and material is crushed to obtain to dry powder completely with grinding, then carries out element determination.
1) metallic element is tested
It weighs 0.1g dry powder to be added in 100mL kelvin bottle, addend, which drips, keeps sample wet, and it is (close that the 3mL concentrated sulfuric acid is then added
Degree is 1.84g/mL, is analyzed pure) and 10 drop perchloric acid (mass fraction is 60~70%, is analyzed pure), it shakes up.It is small that one is put in bottleneck
Funnel, heating, which disappears, on adjustable electric cooking stove boils, until solution colour turns white and transparent in kelvin bottle, is further continued for boiling 20min, cooling
The boil liquid that disappears afterwards is carefully washed in 100mL volumetric flask from kelvin bottle with water, and repeated flushing kelvin bottle keeps mixture in bottle complete
Portion washes in volumetric flask, is then diluted with water to scale, shakes up, and stands clarification, and carefully Aspirate supernatant carries out subsequent phosphorus member
The measurement of cellulose content.Make reagent blank control group with deionized water simultaneously.Calcium is used respectively, and iron, potassium, magnesium, sodium, Zn-ef ficiency standard specimen matches
Standard curve processed.Then successively liquid containing is entered in ICP, carries out constituent content mensure.
2) measurement of organic matter
It takes 0.1g dry powder to be put into Boiling tube (outer diameter 25mm, length 100mm), is added 5mL 0.8000mol/L's with suction pipe
K2Cr2O7Then the 5mL concentrated sulfuric acid (density 1.84g/mL, chemistry are pure) is added with syringe in standard solution, and carefully rotation is shaken
It is even.Oil bath pan (built-in solid paraffin or blade oil) is heated to 180~190 DEG C in advance, Boiling tube is inserted into wire mesh cage, then
It puts it into oil bath pan and heats, kettle temperature should be controlled at this time at 170~180 DEG C, and solution is made to keep boiling 5min, then
Wire mesh cage is taken out, after test tube is slightly cold, cleaning the oil liquid outside test tube with clean paper, (solution after such as boiling is indicated in green
Potassium bichromate dosage is insufficient, should weigh blade dry powder again and reform;Solution after boiling such as solution is in yellow or yellow green, then continues
It is tested).Mixture in test tube is washed in 250mL conical flask after cooling, makes bottle inner volume in 60~80mL or so, adds 3
~4 drop N- phenylanthranilic acid indicator, are titrated, colourshifting process is by brownish red through purple with 0.2mol/L ferrous sulfate solution
It is terminal to blue-green.Record ferrous sulfate dosage.
When measuring the content of organic matter, it is necessary to make 2~3 blank calibration.Blade dry powder is not added in blank calibration, but is added
0.1~0.5g quartz sand, it is identical when other steps are with measurement blade dry powder, ferrous sulfate dosage is recorded, then asks it flat
Mean value.
0.8000mol/L K used2Cr2O7Standard solution is prepared by following methods: 39.2245g potassium bichromate adds 400mL
Water, heating make to dissolve, and are settled to 1L with water after cooling.0.2mol/L ferrous sulfate solution used is prepared by following methods: 56.0g
Ferrous sulfate or 80.0g iron ammonium sulfate, are dissolved in water, add the 15mL concentrated sulfuric acid, are settled to 1L with water.N- phenylanthranilic acid
Indicator is prepared by following methods: molten 0.2g indicator in 100mL mass fraction be 0.2% sodium carbonate liquor in, slightly it is hot simultaneously
It is stirred continuously, promotes the indicator dissolution for floating on surface.
The method for calculating the content of organic matter in blade according to the dosage of ferrous sulfate are as follows:
Organic matter (%)=organic carbon (%) × 1.724 (2)
In formula:The concentration of standard solution, mol/L;
The volume of standard solution, mL;
V0--- the average value of the volume for the ferrous sulfate solution that blank is spent when demarcating, mL;
V --- the volume of the ferrous sulfate solution spent when titration plant sample, mL;
0.003 --- the molal weight of 1/4 carbon atom, g/mol;
1.1 --- oxidation adjusting coefficient;
1.724 --- organic carbon is converted into the coefficient of organic matter;
m1--- the quality of blade dry powder, g;
k2--- the moisture conversion coefficient for being converted into drying will be air-dried, is 1 in the present embodiment.
3) measurement of nitrogen content
It weighs 0.1g dry powder and is sent into dry kelvin bottom of bottle portion, after adding a small amount of deionized water to soak, 2g accelerator is added
(100g potassium sulfate, 10g Salzburg vitriol and 1g selenium powder are finely ground in mortar, are sufficiently mixed and are uniformly made.) and the 5mL concentrated sulfuric acid
(density 1.84g/mL is analyzed pure), shakes up.Kelvin bottle is placed in disappear to boil to disappear on furnace and is boiled, boil liquid to be disappeared all become it is greyish white slightly
After band green, it is further continued for disappearing and boils 1h.Disappear to boil and finish, it is cooling, it is to be distilled.Disappear boil while, two parts of blank determinations are done, except being not added
Outside blade dry powder, other operations are identical.
It is measured on the directly upper azotometer of the well-done sample that will disappear (if being full nitrogen with full-automatic azotometer final result
Percentage composition).The calculation of measurement result are as follows:
In formula: V --- the volume of sour standard solution used, mL when titration test solution;
V0--- the average external volume of sour standard solution used, mL when titration blank;
CH--- the concentration of sour standard solution, mol/L;
0.014 --- mM quality of nitrogen-atoms, g/mmol;
M --- dry powder quality, g.
When measuring nitrogen content using azotometer, when alkalization alkali used be 10mol/L sodium hydroxide solution (sodium hydroxide is
It is industrial or chemical pure).Boric acid used in absorbing ammonia is boric acid-indicator mixed solution.Wherein boric acid is to analyze pure, mass fraction
For 2% boric acid solution, indicator is that a small amount of 95% ethyl alcohol is added in 0.5g bromocresol green and 0.1g methyl red in agate mortar,
Grinding is until reagent all after dissolution plus 95% ethyl alcohol is formulated to 100mL.Before use, 20mL is added in every liter of boric acid solution
Mixed indicator, and reddish violet is adjusted to diluted alkaline.This liquid standing time is unsuitable too long, and such as pH value has change in use
Change, should be adjusted at any time with diluted acid or diluted alkaline.Using sour standard solution titrate when it is used acid be 0.005mol/L sulfuric acid or
0.01mol/L hydrochloric acid standard solution.
4) measurement of phosphorus content
It takes the 5~10mL of supernatant drawn in metal element content continuous mode to be placed in 50mL volumetric flask, adds water to 15
~20mL, adds 1 drop 2,4- dinitrophenol dinitrophenolate indicator, with 4mol/L sodium hydroxide solution adjust pH to solution be just in it is faint yellow, add
The anti-indicator of 5mL molybdenum antimony, is settled to scale with water, shakes up.After 30min on spectrophotometer, with 2cm cuvette, 700nm is selected
Wavelength colorimetric, blank test solution are reference liquid, adjust absorption value to zero, then survey the absorption value of each prepare liquid developing solution.In work
Make the phosphorus ppm number that developing solution is found on curve.
The calculating process of phosphorus content are as follows:
Full phosphorus pentoxide (P2O5, %) and=full phosphorus (P, %) × 2.29 (2)
In formula: the phosphorus ppm number (test value) for the developing solution that C-is looked into from working curve;
V-developing solution volume, 50mL;
Ts-points takes multiple,
(100 milliliters of prepare liquid volume, draw 5 milliliters of testing liquid)
M-blade dry powder quality, g;
106- the divisor for being converted into microgram gram;
2.29-are converted into phosphorus the coefficient of phosphorus pentoxide.
Molybdenum antimony it is anti-than method require developing solution in sulfuric acid concentration be 0.23~0.33mol/L.If acidity is less than 0.23mol/
It is shorter to stablize the time although colour developing is accelerated by L;If acidity is greater than 0.33mol/L, develop the color relatively slow.Molybdenum antimony is anti-than method requirement
Colour temp is 15 DEG C or more, if room temperature is lower than 15 DEG C, can be placed in 30~40 DEG C of insulating box and keeps 30min, is taken out
Colorimetric after cooling.
(3) Production rate
After tomato grows to " picking time ", all fruits are won, carry out every plant of weighing, and record.
(4) statistical analysis and model construction
1) it statisticallys analyze
Using the plant of pouring 1/4MS nutrient solution as normal group of nutrition (YZ), watering group is as nutritional deficiency group (YQ).It is right
The content of the physical signs of measurement take statistics credit analysis, as shown in table 1.
1 results of statistical analysis of table
As shown in Table 1, nutrition, which is normally organized, generates apparent yield difference with nutritional deficiency group tomato.Meanwhile in bud point
Change interim, N, P, K, there is also significant differences for the content of tetra- kinds of elements of Mg.Its medium multiple is approximately equal to YQ average value and YZ is average
The ratio of value, it is therefore intended that judge normal group of the nutrition relationship between each physical signs of nutritional deficiency group whether with reduction
Yield has correlation.
2) model construction
Using whether the underproduction as final result variable, is labeled as TC;Exposure variable includes N content % (labeled as R1), P content
(being labeled as R2), content of organic matter % (are labeled as R3), Ca content (ppm is labeled as R4), Fe content (ppm is labeled as R5), K
Content (ppm is labeled as R6), Mg content (ppm is labeled as R7), Na content (ppm is labeled as R8), Zn content (ppm, label
For R9).
It is combined by variable, establish Logic Regression Models using R language glm function: group is unified: R1+R2+R3+R4+R5+
R6+R7+R8+R9, combination two: R1+R2+R6+R7, combination three: R1+R3+R7+R9.
It randomly selects 50% inside the 50% and YZ data inside YQ data and is used as training set construction model;It is remaining
Data verify model as verifying collection.
(I) according to the data of variable each in training set, the unification of calculating group, combination two and combine three logistic regression coefficient and
AIC value.Calculated result is as shown in 2~table of table 4.
The logistic regression coefficient and AIC value of 2 groups of unifications of table
The logistic regression coefficient and AIC value of the combination of table 3 two
The logistic regression coefficient and AIC value of the combination of table 4 three
According to the AIC value for the model that group unification~combination three is established, top-priority model should combine three moulds established
Type.
(II) model evaluation is carried out with verifying collection
What model (respectively Model I, modelⅱ and model III) the evaluation verifying established with group unification~combination three was concentrated
Plant yield, the actual production then concentrated with verifying compare, by establishing the accuracy of prediction orthogonal matrix judgment models, in advance
It is as shown in table 5 to survey orthogonal matrix.
Table 5 predicts orthogonal matrix
In prediction orthogonal matrix, 0 is indicated normal (i.e. the not underproduction), and 1 indicates to lack (i.e. the underproduction), accuracy rate=(tp+
tn)/(tp+tn+fp+fn).The underproduction is thought when the yield of single plant fruit is subnormal 10% in the present embodiment.
The Model I established using group unification~combination three~the results are shown in Table 6 for model III judgement.
6 Model Is of table~model III judging result
According to table 6 it is found that the accuracy rate of Model I is 91%, the accuracy rate of modelⅱ is 92%, and the accuracy rate of model III is
95%.Top-priority according to accuracy rate should be that model III combines three.
(III) ROC curve evaluation model is utilized
ROC curve is made to Model I~model III of foundation, as shown in Figure 1.Model I~model III is determined according to ROC curve
AUC value be respectively 0.949,0.984,0.959.According to AUC value it is found that model III is the model preferably considered.
(5) forecast production
Characteristic value in conjunction with three models is as shown in table 7, confirms model III, that is, N content % (R1), content of organic matter %
(R3), the parameter combination of Mg content (ppm, R7), Zn content (ppm, R9) can achieve under the prediction of Logic Regression Models
The yield of the prediction tomato of accuracy rate 95% is very suitable to future implementation and application in the production of tomato.
7 three model feature values of table are summarized
In other embodiments of the invention, Logic Regression Models also are established to the conjunction of other set of variables, by verifying,
The model that the predictablity rate that his set of variables builds vertical Logic Regression Models jointly is established without the combination three in embodiment 1
Accuracy rate is high.