Leveraging farm production diversity for dietary diversity: evidence from national level panel data

Dietary diversity is the key to improved health and nutrition. Farm production diversity has the potential of enhancing dietary diversity but this interrelationship varies and is ambiguous in many societies. To examine the effect of farm production diversity on household dietary diversity using nationally representative panel data of Bangladesh we have used Bangladesh Integrated Household Survey (BIHS) data collected by International Food Policy Research Institute (IFPRI) in 2011/12, 2015 and 2018/19. Total sample size is 11,720. For assessing dietary diversity we have used different indicators namely household dietary diversity score (HDDS) and food variety score (FVS). We have also used multiple methods for measuring farm production diversity including production diversity score, crop diversity score and Simpson diversification index. Poisson regression model has been used. Results revealed a strong positive association among farm production diversity, income and dietary diversity though the extent of the association is small. The variables such as market orientation, access to market, age and education are also found to influence on household dietary diversity. Our results propose that for increasing dietary diversity efforts should be taken to increase farm production diversity combined with diverse income and market access.

Page 2 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 2014). South Asian countries mainly depend on agriculture for employment and livelihoods (Deb 2014). Bangladesh, a developing and overpopulated nation in South Asia, has achieved commendable growth in reducing poverty and hunger over the last two decades, but several indices of nutrition remain a matter of concern (Belton et al. 2014;Rahman 2010;Rahman and Salim 2013). In rural Bangladesh, one-third of women are malnourished and 36% of children under five are stunted (Ahmed et al. 2012;Osmani et al. 2016). Moreover, the country continues to struggle with micronutrient deficiencies such as zinc, iron, iodine and vitamin A due to inadequate diets that are dominated by rice and lack diversification (Ahmed et al. 2013). Research indicates that the food habit of Bangladeshi people is inadequately diversified and this emerged as a major barrier in acquiring a standard nutritional status. Dietary diversity is a proxy method of nutrient adequacy of an individual's diet (Kennedy et al. 2011). It can be defined as different types of foods available across and within food groups that can ensure adequate intake of vital nutrients for optimal health (WHO/FAO 1996;Ruel 2002). It is linked to adequate vitamin and micronutrient intake and better child nutrition outcomes (Fongar et al. 2019;Kennedy et al. 2007;Muthini et al. 2020;Steyn et al. 2006). Dietary diversity is determined by different factors including agricultural biodiversity and farm production diversity (Jones et al. 2014;Oyarzun et al. 2013), production technology, historical consumption habits, household income (Doan 2014;Drescher et al. 2009) and cultural traits at the regional level. Although dietary diversity is a crucial indicator of nutrition, but for low-income people consumption of nutritious food is quite difficult (FAO 2013;Chegere and Stage 2020). For households engaged in subsistence agriculture, farm production diversity would be a promising approach for achieving nutritional security and therefore nutritionists have given emphasis to harnessing the linkages between farm diversity and diversity in household diets (Malapit et al. 2015). For small landholders, agricultural diversification plays a vital role in dietary diversity as it improves the consumption of own-produced food, as well as purchased food on account of increased income from the selling of farm produce (Chegere and Stage 2020;Herrero et al. 2010;IFAD and UNEP 2013).
Various trails are indicating how agriculture can affect the nutritional outcomes which include increased farm production and consumption of food, increased income, reduction of food prices and women empowerment (Kadiyala et al. 2014;Kanter et al. 2015). Research related to agriculture-nutrition linkages has recently gained increased attention and several studies have been conducted on-farm diversity and nutrition. Different studies used varied methods of assessing the link between production diversity and dietary diversity with varied outcomes. Some studies showed a positive association between farm diversity and household dietary diversity (Chegere and Stage 2020;Jones 2017;Pellegrini and Tasciotti 2014;Saaka et al. 2017). Chinnadurai et al. (2016) showed a positive relationship between household dietary diversity and crop diversity, whereas the association was found negative in case of diversity of vegetable production and dietary diversity. Likewise, Sibhatu and Qaim (2018b) found a mixed result indicating that farm production diversity was positively related to dietary diversity when measured with simple species count, whereas did not find any significant relation with dietary diversity when production diversity was measured in terms of the number of food crops. Production diversity also interrelates with other factors such as household income, market access and commercialization of households. While rising income has proven to be an important indicator of dietary diversity (Dillon et al. 2014;Singh et al. 2020), small effort has been taken to see how dietary diversity and nutrition are affected by diversifying income of the household (Benfica and Kilic 2016). Also in developing countries, farmers are gradually diversifying their sources of income which makes the relationship is critical to understand (Davis et al. 2014(Davis et al. , 2017Winters et al. 2010).
Our present research adds to the literature in various aspects. The current evidence revealed that the association between them is complex and situation-specific (Muthini et al. 2020;Sibhatu et al. 2015). In this strand of literature, one of the drawbacks is that most of the studies were conducted in the African region, where smallholder farmers were mainly subsistence in nature and the green revolution was not very effective (Jones 2017;FAO 2014). Hence, more research is required in the Asian region to determine the association. The intensity of associations may vary by context, but further research from Asian settings is essential to confirm this. Moreover, most of the studies used crosssectional data and one of the main problems of cross-sectional studies is selection bias which violating the true association (Jones 2017). Panel data is needed not only to reduce omitted variable bias but also to increase the accuracy of the estimation procedure (Tae et al. 2007). Also, many studies used small sample sizes and simple econometric methods for assessing the relationship which makes it difficult to draw firm conclusions especially in the policy context. Our research has not only used the large, nationally representative, unique panel data but also used panel econometric procedures. In addition, we have also used alternative measures for measuring both dietary diversity and agricultural production diversity. Hence, it is necessary to address the above shortcomings and gain a deeper knowledge of the relationship between farm production diversity and dietary diversity of households. The present study seeks to answer the research question of whether farm production diversity has any effect on household dietary diversity. Therefore this study's objective is to examine the effect of farm production diversity and dietary diversity using the complete set of (3 rounds) large, unique and nationally representative panel data in Bangladesh.

Conceptual framework
The question of whether a farm produces its food or buys from the market has a significant contribution to its nutrition. The paradigm shows the interplay among farm production diversity, income and dietary diversity. The relationships among these three factors are endogenous implying that they affect each other and the interconnection is ambiguous (Fig. 1). Farm diversity is a direct strategy to promote food diversity in subsistence farm households. In subsistence farming, production decisions are influenced by consumption decisions and vice versa. This means that households must produce what they want to consume and also they have to diversify their diet based on their products only. In the case of commercial farming, households focus on specialization to increase income which has a significant contribution to the diet. Higher agricultural diversity may give farmers greater access to produce diverse food and livestock animal in their farm and more chances to cultivate other crops (i.e., plantation, cash crops) which may also have a contribution on nutrition through improved revenue. Some exogenous Page 4 of 20 Khandoker et al. Agricultural and Food Economics (2022) Fig. 1 Conceptual framework displaying the associations among farm production diversity, income and dietary diversity (Authors' construction based on Bellon et al. 2016) Page 5 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 of rural Bangladesh) and the remaining households fall under a different stratum referred to as the "Feed the Future Zone (FTF)". FTF households have not been used in our study. Besides, we have limited our research to households which are engaged in any of the agricultural production systems (i.e., growing crops, raising livestock/poultry/fisheries, etc.). Therefore, we have used a smaller sample size than the original BIHS data. Of these data, some of the observations are dropped due to probable errors in the demographic and food consumption data. So the final estimation consists of 11,720 samples using 3 round data (Table 1). A slight change in sample size is observed over the year which is mainly due to split-up households.

Estimation strategy
Our main outcome variable is household dietary diversity score (HDDS) which is a count variable that takes the value between 1 and 12 (or between 1 and 9 when includes 9 food groups). We have used Poisson regression with a maximum-likelihood estimation technique. For measuring the association among farm production diversity, income, and dietary diversity we have used the given model: where DD it is dietary diversity. PD it is farm production diversity and HI it is annual household income at time t. γ it is the random error term, and subscripts i and t indicate observations and time, respectively. α 0 , α 1 and α 2 are coefficients to be estimated. Equation (1) indicates that our model has two explanatory variables, i.e., farm production diversity and income. In addition to farm diversity and income there may be other factors, i.e., market orientation, market distance, farm size, participation in social safety net programs and socioeconomic characteristics that affect dietary diversity. So we have included these factors in our extended model and the model takes the following form: where MA it is a vector of variables covering market access, MO it indicates market orientation, SSN it is participation in social safety net programs and HS it household demographic and socioeconomic characteristics.
All the analyses including data coding, merging and reshaping have been done using STATA version 16.1 software.

Dietary diversity
In this study, our primary outcome variable is dietary diversity. It is a crucial element of balanced diet since higher dietary diversity ensures a higher intake of nutrients (Kennedy et al. 2007;Steyn et al. 2006). Food variety score and dietary diversity score (DDS) (1) Page 6 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 are two commonly used measures for determining dietary diversity (FAO 2011). Food variety score (FVS) simply counts the number of food items consumed by the household during the recall period (Drewnowski et al. 1997). Although the food variety score is a suitable measure for assessing nutritional status it is less preferable for cross-country analysis due to differences in dietary habits across countries (FAO 2011;Pellegrini and Tasciotti 2014). In this setting, dietary diversity score is preferred to food variety score. In our study, we have measured dietary diversity using both FVS and HDDS. Household dietary diversity score which is widely accepted measure records the number of food groups eaten by the households over a recall period (Keding et al. 2012;Sibhatu et al. 2015). Although there is no international standard for the number of food groups required to calculate HDDS, FAO proposed of using 12 food groups (FAO 2011). These are Cereals; White tubers and roots; Legumes, nuts and seeds; Vegetables; Fruits; Fish and other seafood; Meat; Eggs; Milk and milk products; Oils and fats; Sweets; and Spices, condiments and beverages (Table 2). Each food group adds one score point toward the HDDS if a food item from that group is consumed by any member of the household in the given period. Thus, HDDS ranges from 0 to 12. For assessing HDDS studies used both 24-h period (e.g. Koppmair et al. 2016;M'Kaibi et al. 2017) and a 7-day recall period (e.g. Chegere and Stage 2020; Jones et al. 2014;Jones 2017;Sibhatu and Qaim 2017). Special days like festival, holidays and fasting is practiced in every community and for this seven-day recall period which captures the daily variation in diets, could be more appropriate than the 24-h (Sariyev et al. 2021). In this study, we have constructed HDDS using a 7-day recall period. Studies indicate that the last three food group's (viz Oils and fats; Sweets; and Spices, condiments and beverages) contribution to dietary diversity is ambiguous and for this we have calculated HDDS both using 12 food groups and 9 food groups. The food variety score (FVS) is our second DD assessment tool, which we created by calculating the amount of various food items ingested by the household during 7-day recall period. The main difference between HDDS and FVS is that HDDS restricts to different food

Spices, condiments and beverages
Page 7 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 groups (not food items) whereas FVS counts each and every food items consumed by the household.

Farm production diversity
Our key independent variable is farm production diversity which can be defined as the variation of crop and animal species produced in a farm (Sibhatu 2016). It seems like a prospective approach for enhancing food and nutritional benefit of smallholders and enhancing their access to the variety of foodstuffs (Burlingame and Dernini 2012; Fanzo et al. 2013). The empirical literature suggests different indicators for measuring farm production diversity. We have used multiple methods for measuring farm production diversity so that we can compare the outcomes and test the reliability of the relationship. These methods are production diversity score, crop diversity score and Simpson diversification index. Farm production diversity score is a simple count of the number of food crops, vegetables, fruits and animal species of the households (Koppmair et al. 2017;Muthini et al. 2020;Sibhatu and Qaim 2018b). Some of the households also grow non-food cash crops like cotton, jute, tobacco, flowers and other minor crops. Although these crops do not directly contribute to dietary diversity but increase the farm income so we have included them in our study. So our second measure is the crop diversity score where we have included all the crops (viz. cash crops, food crops etc.) produced by the household.
Another widely used indicator of farm diversity is Simpson Diversity Index (SID). SID is superior to other methods in the way that it not only captures crop richness but also the crop evenness (Chegere and Stage 2020; Lovo and Veronesi 2019). SID is calculated as: where n = the total number of crops, S i = area proportion of the ith crop in the total cropped area.
Thus, S i denotes the proportion of the area of the ith crop to the total cropped area. SID ranges from zero to one with zero representing no agricultural diversification and one indicating complete diversification. For evaluating the association between farm production diversity and dietary diversity we have used multiple indicators for both outcome variable and main explanatory variable because these allow for a more thorough and accurate result.

Control variables
The relationship between dietary diversity and production diversification may be influenced by socio-economic and demographic characteristics. For this, based on the thorough literature review we have included these factors as control variables in our regression (Annim and Frempong 2018;Babatunde and Qaim 2010;Carletto et al. 2017;Chegere and Stage 2020;Hirvonen and Hoddinott 2017;Jones et al. 2014;Koppmair et al. 2016;Kumar et al. 2015;Muthini et al. 2020;Pellegrini and Tasciotti 2014;Romeo et al. 2016;Sariyev et al. 2021;Sibhatu et al. 2015;Singh et al. 2020;Thorne-Lyman et al. 2009). For assessing whether income has any effect on dietary diversity we have included Page 8 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 the annual income of the households as the independent variable. For getting annual income we have calculated farm and off-farm income. We have taken other variables like market orientation, access to the market, participation in social safety net programs, etc. Individual socioeconomic characteristics such as age, sex and education of the household head, farm size, household size and status of the household's agricultural asset have also been included. Figure 2 at a glance shows the different explanatory variables and outcome variables used for analyzing the relationship between dietary diversity and farm production diversity:

Robustness check
For conducting the robustness check to corroborate the stability of the major findings we have utilized different estimating methods and measures of important variables. The following techniques have been used: Firstly, several factors may affect household dietary diversity other than farm production diversity. Besides, some omitted variables may be interrelated to farm production diversity, thereby biasing the estimated results. For testing omitted variable bias, we have re-estimated regression models by adding household socioeconomic characteristics.
Secondly, to test whether there are some unobserved variables like a special event which occurred from first to the third round or not we have estimated the two-way FE regression model. In this model, we have included year dummy and used 2011-12 as the reference period, 2016 as period 1 and 2018-19 as period 2.
Lastly, as a robustness check we have used different methods both outcome variables and explanatory variables. For measuring dietary diversity we have used HDDS based on 9 food groups and FVS. On the other hand, we have utilized alternate methods of production diversity such as the Simpson Index of diversity and crop diversity score and estimated the regression model for each of these methods with outcome variable dietary diversity (HDDS, HDDS using 9 food groups and FVS).

Descriptive statistics
Descriptive statistics of the variables are shown in Table 3. Households are consumed more than 9 different food groups out of 12 food groups in 7 days. The average HDDS is 9.50 which indicate that the average household has consumed 9.50 food groups out of 12 food groups during the recall period. The average number of food groups eaten by the households is almost constant, especially in the last 2 rounds. HDDS in case of both 12 food groups and 9 food groups has been improved from the first round to the third round. The average HDDS based on healthy food groups (9 food groups) is 7.47. Households consume on an average 32 items in a 7-day recall period. A significant rise is also observed in the case of FVS from the first round to the third. The values of both dietary diversity and FVS are found to be almost constant in the 2nd and 3rd rounds. Household dietary diversity score (HDDS) and production diversity score (PDS) across 7 administrative divisions in Bangladesh are depicted in Fig. 3. The box plot shows HDDS score is almost the same and not the greater variation is found in the HDDS across divisions. Households consume on an average 9 food groups across divisions. However, variation is found in PDS across divisions. Highest PDS is found in Barisal division (6.2) and the lowest is in Chittagong divisions (4.2).
Kernel density distribution of the production diversity score shows that households are producing around five different crops or animal products in the given 12-month period (Fig. 4) The right-skewed distribution curve shows that the average value of farm production diversity is greater than the mode and most of the of household's production diversity is below the average. Households produce on an average 5 different crop species in a year. The value of Simpson index of crop diversification is 0.31. All the indicators of farm production diversity indicate that farm diversity is low in Bangladesh but has increased over time.
The value of market orientation is only 17% indicating that most of the agricultural households experienced subsistence farming and the annual per capita income is comparatively low which is approximately Tk. 1, 38,000 (approximately USD 1665 as at the time of writing). The average distance from home to the nearest market is 1.17 km. The socio-economic characteristics show that around 80% of households are male-headed. The average years of schooling of the household head are very low (only 3 years). The average household size is consisting of 4 people. The average farm size is 75.84 decimal and around 45% of households participate in the different social safety net programs. Table 4 exhibits the results of Poisson regression using pooled, fixed effects and random effects where we have used household dietary diversity score as dependent variable and production diversity score and annual per capita income as explanatory variables. Farm production diversity is positively significantly related to household dietary diversity but the extent of this association is very small. The result of pooled specification indicates that DDS would be increased by 0.7% at the household level if one food group is produced. There is a significant positive association between annual income and HDDS which means that income could improve households' capability to purchase a wider range of foodstuff and thereby increase dietary diversity.

Association between dietary diversity and farm production diversity
Page 10 of 20 Khandoker et al. Agricultural and Food Economics (2022)   Page 11 of 20 Khandoker et al. Agricultural and Food Economics (2022)   Page 12 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 To examine whether market has any effect on the association between diet diversity of households and farm production diversity we have extended the model by using this factor ( Table 5). Because of the significance of market orientation, we have extended the model by adding an interaction term namely market orientation which is recorded by the average percentage of agricultural output sold. For assessing whether distance from home to nearest market has any effect on dietary diversity we have included the variable as market access. Given that households having within walking distance of markets can diversify their eating habits by buying goods from markets, it is reasonable to expect stronger linkages between farm production diversity and household dietary diversity further away from markets. The variable market access is found negative and significant showing that households in more remote areas have less dietary diversity. Simultaneously, the effects of farm diversity on household dietary diversity are still strong. When the estimates are compared, it appears that reducing the walking time to the district market by 1 h would have a greater positive impact on household dietary diversity than growing one more food crop to the farm. Our results noticeably show that market access is important for increasing the dietary diversity of farm households. The variable namely market orientation is found positive and significant suggests that increasing the commercialization may be a better strategy for improving nutrition than encouraging more diverse subsistence agriculture.
The relationship between dietary variety and production diversity also depends on household and socio-economic factors. The extended model of the association among production diversity, income and dietary diversity showed that farm production diversity and annual income have a significant positive association with dietary diversity irrespective of model specifications (Table 6). We have calculated alternative measures of farm diversity (i.e., SI and CDS) and found a significant positive association with dietary diversity. This is common in subsistence farming where the farm household consumes a large portion of the farm's output. The variable market access suggests that dietary diversity can be enhanced by reducing the distance from the market or improving infrastructure.
For testing omitted variable bias, we re-estimated the model using socio-economic variables and found that some of the socioeconomic variables such as age and education of household head and household size including production diversity are positively Table 5 Relationship among production diversity, income, market orientation and market access Authors' analyses using 3 rounds of BIHS panel data. The dependent variable is HDDS ***, and ** indicate significance at the 1% and 5%, levels, respectively; figures in the parentheses indicate standard errors Page 13 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 significantly associated with HDDS in this newly estimated model. This indicates that our primary findings are free from omitted variable bias (Table 6). A positive significant association with age indicates that the households having older heads have revealed higher dietary diversity. The variable education has a significant and positive effect indicating that educated people are more conscious about nutritionally-balanced diets and diversified baskets. Education not only enhances the knowledge of health and nutrition, but it also helps to manage the diversified food at an affordable cost. Household size is positively significantly associated with HDDS indicating higher dietary diversity with more members in the households. Bangladesh Government is conducting different social safety net programs for upgrading the socio-economic status especially for the rural and disadvantaged group such as Stipend for Primary Students, School Feeding Program, Stipend for Dropout Students, Stipend for Secondary and Higher Secondary or Female Student etc. For assessing whether this SSN program has any effect on dietary diversity we have added this variable. We have found a positive significant relation with dietary diversity in the pooled model. Table 6 Association among farm production diversity, income and other confounding factors Authors' analyses using 3 rounds of BIHS panel data. The dependent variable is HDDS ***, and ** indicate significance at the 1% and 5%, levels, respectively; figures in the parentheses are standard errors  Page 14 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15

Discussion
One of the crucial outcomes of this study is that farm production diversity is positively significantly associated with household dietary diversity irrespective of the model specifications and methods used. The extent of the association is relatively smaller and is relevant to previous studies. Most of the earlier studies found the extent of the relationship between 0.05 and 0.20. (Bellon et al. 2016;Koppmair et al. 2017;Sariyev et al. 2021;Sibhatu and Qaim 2018a, b) although Jones et al. (2014) and Chegere and stage (2020) found higher impact using OLS approach. Also this small association is persistent in developing countries like Bangladesh where the majority of farmers are smallholder and subsistent (World Bank 2007). Studies showed that the relationship between farm diversity and dietary diversity depends on the methods used for measuring both. The relationship is positive when both crop and livestock species are taken into account, however, if only crop count is utilized, the relationship become negative (Jones et al. 2014). Sibhatu and Qaim (2018b) showed that in certain cases, other variables such as market access and off-farm income had greater impacts on dietary diversity than farm diversity. However, in our study we have found other factors are also important but not as farm diversity. The variable namely annual income is positively associated with dietary diversity which indicates that increasing income can increase household's ability to purchase diversified food. This result is also consistent with previous studies (Babatunde and Qaim 2010;Singh et al. 2020). The estimated coefficient of market orientation is positive for dietary diversity but the extent of association is smaller than the farm production diversity (0.010 > 0.001). These findings show that although both production diversity and commercialization are crucial tools for increasing dietary diversity but production diversification may be a better approach than commercialization in the long run. However, Koppmair et al. (2016), found that access to markets for buying food and selling farm produce along with modern technology more important for dietary diversity than farm diversity. The variable market distance is negatively associated with dietary diversity which is persistent with previous studies (Muthini et al. 2020;Sibhatu and Qaim 2018b). To know more about the underlying relationship, i.e., whether the relationship is linear or not we have used the interaction term of production diversity squared and estimated the model. A significant positive interaction term (PD 2 ) implies that increased farm production variety is connected with greater dietary diversity, whereas a negatively significant PD 2 suggests that as farm production diversity grows, the magnitude of the association declines (Sibhatu and Qaim 2018b). This study has not found any significant association with dietary diversity after using the interaction term production diversity squared although the sign is negative. Insignificant coefficients do not imply that there is no association or that the interaction component should be deleted but depend on the values of interrelating factors (Ai and Norton 2003;Brambor et al. 2006;Sariyev et al. 2021). It's also likely that production diversity has a significant marginal impact on a range of other parameters that interact with production diversity measures (Sariyev et al. 2021).
We have estimated two way fixed effect regression model to test and control the unobserved bias and found that year dummy is positively significant which indicates that time plays a crucial role and dietary diversity has increased over the years (Table 7).
Finally, all of the regression results using multiple indicators in various situations and specifications show that farm diversity have positive significant relationship with dietary Page 15 of 20 Khandoker et al. Agricultural and Food Economics (2022) 10:15 diversity at the household level (Tables 8, 9). The extent of the association is small but it is nevertheless greater than other coefficients which are significant, indicating the importance of farm diversity on dietary diversity in a rice dominating country.

Conclusions
Dietary diversity is the most commonly used and essential indicator of food and nutrition. At the macro level, this necessitates agricultural diversification. In recent years, nutrition-sensitive agricultural research has gained considerable attention for improving nutritional status and one example of such program is farm diversification. So far different studies have found diverse outcomes. Policy makers are interested to know the exact relationship between farm diversity and dietary diversity. Our study examines the association between farm production diversity and household dietary diversity using a nationally representative and complete (3 rounds) panel dataset of Bangladesh. For doing this we have used alternate measures for assessing both farm diversity and dietary diversity.
We have initiated our analysis with the pooled Poisson regression which assumed that production diversity is exogenous and finally move to the fixed effect model. Our results indicate that both household dietary diversity and farm production diversity increases over the years although production diversity is comparatively low in Bangladesh. Our results provide robust evidence that production diversity is positively significantly associated with dietary diversity although the magnitude is lower. We have also found a positive significant association between income and dietary diversity. Market orientation, access to the market and participation in social safety net programs has also impact on household dietary diversity but the magnitude is lower than production diversity. In Bangladesh as rice is the staple crop, traditionally different government and non-government large-scale support programs have mainly focused on rice cultivation to maintain food security of its growing population, with limited incentive have given to farmers for diversifying production. But in recent years, Government of Bangladesh (GoB) has underlined the necessity of agricultural diversification in different policy papers like National agricultural policy, Country Investment Plan and Five Year Plan. Agriculture nutrition linkages sometimes appear to be too complicated to pursue. In this case our study will help in shaping the agenda of agricultural development strategy that is essential for improving nutrition among farm households. Therefore, from a policy perspective our results suggest that as farm production diversity is an important indicator of household dietary diversity and nutrition agricultural policies and programs should focus on production diversity rather than solely increasing the production of selected staple crops. Furthermore, an attempt can also be taken to encourage income diversification and market access for addressing the nutrition security of SDGs of Bangladesh. One of the potential limitations of this study is that our data is based on three points of time, but dietary diversity of farm households changes throughout the time especially within the years, such as with agricultural seasons. More research using high frequency data is needed to solve the problem of seasonality component of dietary diversity. In addition, we have analyzed effect of farm production diversity on household dietary diversity but not individual dietary diversity. Hence, more research is needed to better understand the intra household dietary diversity especially for women and children.