27期

Size: px

Start display at page:

Download "27期"

型幼诸
7 years ago
Views:

1 中国农学通报,(7):- Chinese Agricultural Science Bulletin 田间肥效试验数据的频率分析和施肥决策李洪文, 叶和生, 李保华, 李春莲, 施文发, 郭家军, 苏亚芳, 李姝莹 ( 云南省双柏县农业技术推广服务中心 ( 土壤肥料工作站 ), 云南双柏 7; 鄂嘉镇农业技术推广服务中心, 云南双柏鄂嘉镇 77; 双柏县农业局农业执法大队, 云南双柏 7) 摘要 : 为解决农业试验统计分析软件扬州分析器 (.) Excel 进行回归分析建立的部分非典型性肥料效应方程推荐的施肥量出现异常偏大偏低甚至负值, 脱离当地生产实际的问题本文应用概率理论的基本知识和频数的原理和方法, 对油菜田间肥效试验数据应用频率分析法进行统计分析, 得到的油菜优化施肥量产量为 7. kg/hm, 氮磷钾优化组合范围分别为.~. kg/hm.~. kg/hm.~. kg/hm, 与试验地 ( 田 ) 和当地生产实际相吻合本研究是在频率分析的基础上, 用加权平均方法, 确定各种不同产量目标的生产因素组合, 将其中产量高, 出现频率大 ( 稳产 ), 用肥节省的组合作为优化生产措施, 具有增加施肥决策信息量, 减少或避免小概率事件的风险的优点可解决肥效试验数据分析统计汇总, 应用扬州分析器 (.) Excel DPS 等统计分析软件进行回归分析建立的部分非典型肥料效应函数推荐的施肥量和目标产量出现异常值, 与当地生产实际不相符的问题肥料效应函数方程经方差分析不显著的试验, 能否应用该研究方法进行施肥决策, 推荐氮磷钾肥料优化组合和优化施肥量产量, 还需在以后的工作实践中进一步探讨关键词 : 田间肥效试验 ; 肥料效应函数 ; 频率分析 ; 施肥决策中图分类号 :S7. 文献标志码 :A 论文编号 :- Frequency Analysis and Fertilization Decision Data Field Fertilizer Efficiency Test Li Hongwen, Ye Hesheng, Li Baohua, Li Chunlian, Shi Wenfa, Guo Jiajiang, Su Yafang, Li Shuying ( Agricultural Technology Extension and Service Center of Shuangbai County of Yunnan Province (Soil and Fertilizer Station), Shuangbai Yunnan 7; Agricultural Technology Extension and Service Center of EJia Town, EJia, Shuangbai Yunnan 77; Agricultural law enforcement unit of Shuangbai County, Shuangbai Yunnan 7) Abstract: The part of fertilizer effect of atypical equation recommended by agricultural experiments and statistical analysis software, Yangzhou analyzer (.), regression analysis of Excel fertilization appears abnormal low, even negative, from the local production of practical problems. The study adopted the principle and method of basic knowledge and the frequency of using probability theory and carried out statistical analysis on the rape field fertilizer experiment data by frequency analysis method, the optimal fertilization rape yield was 7. kg/hm, nitrogen, phosphorus, potassium optimization combination ranges were.-. kg/hm,.-. kg/hm,.-. kg/hm, coincided with testing (Tian) and local production reality. This research is based on frequency analysis, using the weighted average method with production factors to determine the different target yield, and took the high yield, stable frequency of high yield and optimizing of production measures as the combination of fertilizer saving to increase the advantage of fertilization decision- making, therefore reduce or avoid the risk of low probability. The results of this study could solve the problem of abnormal values caused by Yangzhou analyzer (.), regression analysis of Excel, 基金项目 : 国家测土配方施肥财政补贴项目资金 ( 云财农 []) 第一作者简介 : 李洪文, 男, 年出生, 云南双柏县人, 高级农艺师, 大专, 从事土壤肥料技术研究示范推广工作通信地址 :7 云南双柏县妥甸镇妥甸街号双柏县农业局,Tel:7-777, ynnklhw@.com 收稿日期 :--, 修回日期 :--

2 李洪文等 : 田间肥效试验数据的频率分析和施肥决策 and DPS statistical analysis software. The function of fertilizer effect equation by the analysis of variance test was not significant, and the application of the research methods of fertilization and the recommended optimization of fertilization of nitrogen, phosphorus, potassium still need to practice in further study. Key words: field fertilizer experiment; fertilizer effect equation; frequency analysis; fertilization decision 引言 [] [] 孙克刚等陈建祥等采用三因素二次通用旋转组合回归设计, 对田间试验产量结果经微机分析处理, 建立氮磷钾三元二次肥料效应数学模型, 用微机对回归模型进行模拟寻优, 获得氮磷钾肥优化施肥 [] 量肖志江等对油菜肥效试验产量结果进行回归分析, 建立油菜氮磷钾三元二次肥料效应方程, 依据方程, 按最大边际效应求取偏导数预测最高产量 [] 施肥量和最佳施肥量诸化斌等李大群等 [] 对油菜完全实施的田间肥效试验结果进行回归分析, 分别拟合氮磷钾三元二次二元二次一元二次肥料效应函数模型, 通过求解该函数模型推荐油菜氮磷 [] [7] [] 钾施用量吴志勇等徐卫华等张建平等应用 Excel 的相关功能对田间肥效试验结果进行回归分析和方差分析, 分别建立氮磷钾一元二次二元二次三元二次肥料效应方程, 通过求解方程, 获得氮 [] 磷钾施用量和目标产量张祯勇等应用 R 语言对田间肥效试验结果进行统计分析, 分别建立氮磷钾三元二次二元二次肥料效应回归方程, 求解该 [] 方程, 得到氮磷钾施用量和目标产量戴树荣应用 DPS 数据处理系统对完全实施的田间肥效试验结果进行分析, 分别拟合马铃薯产量与氮磷钾三元二次二元二次一元二次肥料效应函数, 应用计算机对肥料效应函数进行求解, 其中三元二次回归方程进行典型性判别分析, 回归模型为典型式, 采用边际分析求解 ( 回归模型为非典型式的剔除 ), 获得各试验点最高产量施肥量和经济最佳施肥量胡胜勇等 [] [] 肖茂盛等对完全实施的田间肥效试验结果经回归分析和方差, 分别建立氮磷钾一元二次二元二次三元二次肥料效应方程, 由方程求解极值 ( 边际效 [] 应 ) 得出氮磷钾最佳施用量宋晓梅等应用 SPSS 统计分析软件农业试验统计分析软件进行回归分析, 建立三元二次肥料效应方程, 用该方程的交互项的系数进行交互分析获得最佳产量及其施肥量岳素清 [] [] 李洪文使用扬州分析器 (.) 对完全实施的田间肥效试验结果进行回归分析和方差分析, 分别建立氮磷钾一元二次二元二次三元二次肥料效应方程, 得到最高产量及其施肥量和最佳产量施肥量近年来, 全国各省 ( 区, 市 ) 县 ( 旗, 区, 市 ) 在开展的测土配方施肥项目过程中, 实施了大量的田间肥效试验, 取得了大量的试验数据但由于受各种因素的影响, 系统误差较大, 使用 DPS 数据处理系统 SPSS 统计分析软件农业试验统计分析软件等和扬州分析器 (.) 或应用 Excel 的相关功能对田间试验结果数据分析得到的结果达显著极显著的只有试验总量的 % 左右, 其中, 部分试验的统计分析推荐的目标产量肥料施肥量与当地实际相差较大, 出现异常偏大或异常偏低, 甚至负值脱离当地生产实际的现象, 而舍弃这些数据又有些欠妥 [,-7] 到目前为止, 多数研究者刊文报道的是通过回归分析建立氮磷钾三元二次肥料效应方程, 经 F 检验达 F. 显著以上, 方程系数 b b 和 b 均为正值,b b 和 b 均为负值时, 符合肥料报酬递减律, 认为该三元二次方程为典型肥料效应函数, 可进行施肥量推荐, 供生产者用于农业生产 ;F 检验达显著以上, 方程系数未能与 b b b 均为正值和 b b b 均为负值规律相吻合的也可认为拟合成功 [-], 这类函数为非典型肥料效应函数, 可分别建立个两因素的二元二次肥料效应函数和三个单因素的一元二次肥料效应函数, 同样可进行施肥量推荐, 但部分田间试验推荐的施肥量依然存在与当地生产实际不相符的现象因 [] 此, 于年, 何晓滨提出了频率分析法进行施肥决策技术, 认为农业科研工作中, 田间试验影响因素较多, 模拟的肥料效应方程会出现多种情况可能不符合报酬递减律, 肥料效应方程的生产函数不存在极大值, 不能用边际分析方法求解最佳施肥量如果用边际分析方法求解最佳施肥量和最高产量施肥量可能会出现脱离生产实际的推荐用量, 如, 负值, 极大偏高偏低等如, 非典型性水稻肥料效应方程采用边际分析 [] 方法进行施肥决策, 求得的最高产量及其施肥量分别为. kg/hm,=. kg/hm,p O =. kg/hm, K O=-. kg/hm, 显然与生产实际不符 ; 方程 F 检验不显著, 遇到类似这例脱离当地生产实际的统计分析结果时, 不要轻易放弃对田间肥效试验结果数据的进一步分析, 可以用频率分析方法寻优, 挖掘有价值的信息, 提出试验范围内生产因素的优化组合采用频率分析法进行施肥决策, 得到的水稻优化施肥量产量和氮磷钾优化施用量与当地生产实际接近, 解决了测土

3 中国农学通报配方施肥田间肥效试验推荐的施肥量出现极大值和极 [] 小值, 与当地生产实际相差较大的问题但何晓滨的方法是借助 Excel 的相关功能和相关软件来完成其计算, 对于多数县乡镇农业科技人员来说, 由于学历知识办公软件等各方面条件有限, 应用该项技术对田间肥效试验结果数据进行分析计算更有难度为此, [] 笔者参照何晓滨提出的频率分析法施肥决策技术, 将手工分步运算流程进行分述, 供同行探讨参考材料与方法频率分析是考查一组数据中的某个不同的数值出现的频数, 或者是某个数据落入指定区域内的频数, 以了解数据的分布状况的分析方法田间肥效试验数据的频率分析是根据田间肥效试验结果的统计特性, 利用田间肥效试验结果资料, 设定小区产量统计区间, 分析某个小区产量值在统计区间出现的频率 ( 或重现 ) 之间的定量关系其实质是在频率分析的基础上, 用加权平均方法, 确定各种不同产量目标的生产因素组合将其中产量高, 出现频率大 ( 稳产 ), 用肥节省的组合作为优化生产措施显然, 对函数驻点为最大值的肥料效应方程也可以进行频率分析, 只是采用不同的分析方法, 求得的施肥推荐结果就不一致在这种情况下, 频率分析的优点主要是增加施肥决策信息量, 减少或避免小概率事件的风险计算公式为 [-] : p - xi = n j u j n ij (i=,,) () j = n ij: 表示变量值 X ij( 即 u j) 出现的次数,n j: 表示变量 X i 总共出现的次数 S - x i = p n j - (u j -xˉi) n ij () j= - % 置信区间 :( xi - Sxˉi t.(nj - ) ) ~( ) - xi +Sxˉi t.(nj - ), t.(nj - ) 查表可知均数标准差 S x= S () n j. 试验设计采用完全试验设计 ( 表 ), 设 P O K O 肥料个因素, 四个施肥, 其中 : 为不施肥, 为推荐施肥量, 根据历年试验数据和生产实践, 确定 P O K O 肥料推荐施用量分别为 kg/hm kg/hm kg/hm, 为的 /, 为的 / 试验设个处理, 次重复, 个小区随机区组排列小区面积. m. m 中产田和低产田各安排一组试验试验不施任何农家肥小区重复间距. m, 小区间距. m 打塘直播, 塘行距. m. m, 每个小区行, 每行塘, 每塘留苗株, 每个小区株施肥时期及方法, 肥施用总量的 % 和 P O K O 肥料全部作基肥施用, 肥施用总量的 % 在油菜叶期作提苗肥施,% 在抽苔初期施用适时灌水, 加强田间水肥管理试验田油菜成熟收获前和收获后分别选取每个重复的小区中的油菜植株籽粒和土壤样品, 分别制成各处理的植株籽粒和土壤样品检测化验单收单打单晒表油菜测土配方施肥田间肥料效应试验肥料施用量试验编号处理 PO KO 肥料施用纯量 (kg/hm ) PO KO PK PK PK PK PK PK 7 PK PK PK PK PK PK PK PK

4 李洪文等 : 田间肥效试验数据的频率分析和施肥决策. 供试验材料.. 试验地概况该试验位于云南省双柏县妥甸镇中山村委会小龙潭村民小组房子对面普兆禄承包田, 海拔 m, 前作水稻, 紫泥田, 土壤质地中壤, 肥力中等试验田基础土样检测分析结果, 土壤有机质含量. g/kg, 碱解氮 7. mg/kg, 有效磷 7. mg/kg, 速效钾. mg/kg,ph... 供试验肥料 : 云天化金沙江牌尿素 ( 含氮 %),P O : 普通过磷酸钙 ( 勤丰牌含 P O %); K O: 硫酸钾 ( 双狮牌, 含 K O %).. 供试验农作物及其品种供试验农作物为油菜, 品种花油号. 试验经过年月日预整地, 月日施用基肥播种, 月日出苗, 于月日氮素化肥兑水浇施作提苗肥, 月日兑水浇施氮素化肥作蕾苔肥分别于月日月日每公顷用 % 瑞功 7 g 兑水 kg 喷施, 防治蚜虫月日成熟收割结果与分析. 建立三元二次肥料效应数学方程应用扬州分析器 (.) 或者用 Excel 的相关功能对试验数据进行分析, 得到三元二次肥料效应数学表达式 : Y= P-.7K-. -.P -.K +.P+.K+.PK 经方差分析 ( 表 ),F=.>F.=, 方程 F 检验达显著扬州分析器推荐的目标产量肥料施用量 ( 表 ) 与试验田基础土样分析结果和表对照, 综合当时调查了解到的试验田生产力综合分析, 扬州分析器推荐的结果与当地生产实际不符需要使用频率分析方法进行施肥量推荐表方差分析 df SS MS F F. 回归分析. 7..* 残差.7.7 总计 7. 注 :* 表示 F.=>F=., 达显著表扬州分析器推荐的目标产量和肥料施用量方案 /(kg/hm ) PO/(kg/hm ) KO/(kg/hm ) 目标产量 /(kg/hm ) 最大施肥最佳施肥表缺素相对产量相对产量 /% 缺. 缺 P. 缺 K.. 设定统计区间 ( 产量范围 ) 分别统计诸因子 X i 在统计区间内出现的次数, 开展的测土配方施肥田间试验, 设计的是氮磷钾个因素, 四个施肥, 有个处理, 次重复, 共有个小区试验是在水肥栽培管理等条件基本相同的情况下完成的, 个施肥在个小区中出现的频率是多少? 根据田间试验结果 ( 表 ), 设定产量范围 y ~y, 本试验小区实收平均产量最大. kg/ m, 最小. kg/ m, 设定产量范围 ( 统计区间 ) 为. kg<y<. kg/ m, 分别统计诸因子 X i 落入统计区间内出现的次数, 并计算其频数分布变量平均值平均值的标准误 ( S - x i ) % 置信区间和 Y 的个数 K 个组合方案因素及其出现的频率. 统计落入设定产量范围 ( 统计区间 ) 的频率根据表中的小区产量, 分别统计氮磷钾各因素的小区产量落入设定产量范围内出现的频率, 统计结果为表, 表 7, 表. 分别计算氮磷钾平均数标准差均数标准差将表表 7 表统计结果填入表, 并进行相应的分析计算根据式 () 分别得到氮磷钾施用量 (kg/hm ) 的加权平均数 : 氮加权平均数 xˉ =( ) =.;

5 中国农学通报表普兆禄油菜田间试验产量表编号代码肥料施用纯量 /(kg/hm ) PO KO Ⅰ 小区产量 /(kg/ m ) Ⅱ Ⅲ 平均折合 / (kg/hm ) PK..*... PK PK..... PK...7. PK..... PK PK.. *.7. PK..... PK PK. *... PK. *.*.. PK PK....7 PK. *.. 7. 注 : * 表示小区产量值未能落入.<Y<.(kg/ m ) 统计区间表氮各因素的小区产量在设定产量统计区间出现的频率表 7 磷各因素的小区产量在设定产量统计区间出现的频率因素试验设计表中的小区产量落入表中的小区产量 (kg/ m ) 因素试验设计表中的小区产量落入表中的小区产量 (kg/ m ) 频率 / 次统计区间的频率 / 次未能落入统计区间的频率 / 次统计区间的频率 / 次未能落入统计区间的. P. P.... P..... P. 合计合计因素试验设计表中的小区产量落入表中的小区产量 (kg/ m ) K K K K 合计表钾各因素的小区产量在设定产量频率 / 次统计区间出现的频率统计区间的频率 / 次未能落入统计区间的 ; 7. 准差 : 磷加权平均数 xˉ =( ) = 钾加权平均数 xˉ =( ) = 根据式 () 分别得到氮磷钾施用量 (kg/hm ) 的标根据式 () 分别得到氮磷钾施用量 (kg/hm ) 的均数标准差 S x: S = ( - ) +( - ) +( - ) +( - ) - S P = ( - 7) +( - 7) +( - 7) +( - 7) - S K = ( - 7) +( - 7) +( - 7) +( - 7) - 氮的均数标准差 S x =. =.7 =.7 =. =. ;

6 李洪文等 : 田间肥效试验数据的频率分析和施肥决策 7 表普兆禄油菜频率分析 PO KO 因素施用量 / (kg/hm ) 频次频率 /% 施用量 /(kg/hm ) 频次频率 /% 施用 /(kg/hm ) 频次频率 /% 小计平均数 xˉ 7 7 标准差 S..7.7 均数标准差 Sxˉ = S n... % 置信区间 ( xˉ ±. Sxˉ ).~. kg/hm.~. kg/hm.~. kg/hm 注 : 频次指因素在田间试验中出现的次数 ; 频率 (%) 等于频次除以频次小计 ;T=,k= 磷的均数标准差 S x =.7 =. ; 钾的均数标准差 S x =.7 =.. 计算氮磷钾优化组合范围应用 % 置信区间 ( xˉ ±. Sxˉ ), 分别推荐氮磷钾施用量 % 置信区间其中, 氮推荐施用量 % 置信区间 ( xˉ ±. Sxˉ )=.±..=.~. kg/hm ; 磷推荐施用量 % 置信区间 ( xˉ ±. Sxˉ )=7.±..=.~. kg/hm ; 钾推荐施用量 % 置信区间 ( xˉ ±. Sxˉ ) =7. ±..=.~. kg/hm. 计算优化施肥量产量将经频率分析得到的 P O K O 肥料施用量平均数分别为 kg/hm 代入油菜肥料效应数学方程 : Y= P-.7K-. -.P -.K +.P+.K+.PK 通过计算得到油菜优化施肥量产量 y=7. kg/hm 氮磷钾优化组合范围 :=.~. kg/hm, P O =.~. kg/hm,k O=.~. kg/hm, 氮磷钾肥料比例 :P O :K O=:.~.7:.~.7 结论与讨论本研究经回归分析建立的油菜肥料效应方程为非典型肥料效应函数, 不符合报酬递减律,F 检验达 F. 显著 ( 表 ), 用边际分析方法求解最佳施肥量和最高产量施肥量出现脱离生产实际的推荐施用量, 油菜最大施肥量 ( 表 ) P O K O 施用量分别为 kg/hm, 目标产量 7. kg/hm ; 最佳施肥 P O K O 施用量分别为 kg/hm, 目标产量 77. kg/hm 经频率分析得到的油菜优化施肥量产量 7. kg/hm, 肥料优化组合范围, 氮.~. kg/hm, 磷.~. kg/hm, 钾.~. kg/hm, 氮磷钾肥料比例 :P O :K O= :.~.7:.~.7 与当地生产实际相吻合多数研究者报道的是使用 DPS 数据处理系统 SPSS 统计分析软件农业试验统计分析软件和扬州分析器 (.) 应用 Excel 的相关功能对田间试验结果数据进行回归分析, 建立的方程为典型肥料效应函数, 求解该方程, 得到氮磷钾施用量和目标产量方差分析达 F. 显著, 非典型肥料效应函数的试验, 用边际分析方法求解最佳施肥量和最高产量施肥量出现脱离生产实际的推荐施用量本研究应用频率分析方法对完全实施的田间肥效试验结果进行分析, 是在频率分析的基础上, 用加权平均方法, 确定各种不同产量目标的生产因素组合, 将其中产量高, 出现频率大 ( 稳产 ), 用肥节省的组合作为优化生产措施, 具有增加施肥决策信息量, 减少或避免小概率事件的风险的优点本项试验研究在田间试验期间, 正逢百年不遇的持续冬春干旱, 受其影响, 试验小区产量偏低, 从而对频率分析推荐的优化施肥量产量和氮磷钾肥料优化

7 中国农学通报组合将受一定的影响本项研究方法, 可以应用于当前在全国各地开展的测土配方施肥项目田间肥效试验结果统计分析与汇总, 可解决近年来在全国各地开展的测土配方施肥田间试验, 应用扬州分析器 (.) Excel DPS 数据处理系统 SPSS 统计分析软件农业试验统计分析软件统计分析, 非典型肥料效应函数推荐的施肥量和目标产量出现负值, 或极大值极小值, 与当地生产实际不符合的问题回归分析建立的肥料效应函数方程, 经方差分析未达 F. 显著的田间试验, 能否应用该研究方法进行施肥决策, 推荐氮磷钾肥料优化组合和优化施肥量产量, 还需在以后的工作实践中进一步探讨参考文献 [] 孙克刚, 杨稚娟, 王守刚, 等. 油菜目标产量施肥数学模型及氮磷钾肥推荐措施研究 [J]. 土壤肥料,():-. [] 陈建祥, 张安平, 周文娟, 等. 油菜田间肥效试验研究 [J]. 科园月刊,():-. [] 肖志江, 贺峰, 蒋代铭, 等. 油菜施肥效果试验 [J]. 湖北农业科学,():-7. [] 诸化斌, 魏彬彬. 油菜肥效田间试验报告 [J]. 现代农业科技, ():-,. [] 李大群, 雷廷英, 汪森富. 油菜肥料效应试验研究 [J]. 现代农业科技, ():-. [] 吴志勇, 杜之虎, 徐鸿, 等. 利用 Excel 进行肥效统计分析 [J]. 新疆农垦科技,7():-. [7] 徐卫华, 吴雪芬, 张远照. 基于 Excel 的田间肥效试验模型的应用研究 [J]. 现代农业科技, (7):-. [] 张建平, 李爱群, 叶茂华. 甜玉米田间肥效试验 [J]. 现代农业科技,():-. [] 张祯勇, 高明文, 肖启银, 等. 基于 R 语言的肥效试验的统计分析 [J]. 中国农学通报,7(7):7-. [] 戴树荣. 应用试验设计建立二次肥料效应函数寻求马铃薯氮磷钾适宜施肥量的研究 [J]. 中国农学通报,():-. [] 胡胜勇, 王龙, 吴恩华, 等. 铜梁县玉米测土配方施肥试验及模型拟合研究 [J]. 西南师范大学学报 : 自然科学版,():-. [] 肖茂盛, 潘礼斌, 赵吉胜, 等. 水稻肥效田间试验报告 [J]. 现代农业科技, ():-+. [] 宋晓梅, 蒋太明, 刘洪武, 等. 基于田间试验的水稻施肥模型的研究 - 以重庆市合川区为例 [J]. 西南师范大学学报 : 自然科学版, ():-. [] 岳素清, 金喜强. 玉米肥料效应函数建立与应用的研究 [J]. 内蒙古农业科技,():-. [] 李洪文, 李保华, 李春莲, 等. 紫砂泥田水稻肥料效应田间试验 [J]. 现代农业科技,():-. [] 程秋华, 陆萍, 冯建忠, 等. 水稻氮磷钾肥效试验初报 [J]. 上海农业学报,,():7-. [7] 宋朝玉, 高峻岭, 张清霞, 等. 肥料试验结果统计分析方法的探讨 [J]. 山东农业科学,():-. [] 朱显海. 概率论与数理统计 [M]. 长春 : 东北师范大学出版社,: -. [] 方萍. 实用农业试验设计与统计分析指南 [M]. 北京 : 中国农业出版社,. [] 何晓滨, 熊艳, 孔令郁. 田间试验统计分析方法 [C].. 致谢 : 此课题在试验数据分析整理过程中, 承蒙双柏县职业高级中学数学教师李琼英同志的悉心指导, 谨此致谢

Microsoft Word 李强.doc

Microsoft Word 李强.doc 土壤 (Soils), 2017, 49(3): 565 571 DOI: 10.13758/j.cnki.tr.2017.03.020 1 曲靖烤烟镁含量特征及其影响因素 1,2 2 2 2* 3 4 4 3 4 (1 410128 2 410128 3 655000 4 ( ) 650231) 3 506 4.10 g/kg (0.36 ~ 17.70 g/kg) 46.36%