source('../../Data/get_food_import_data.R') DF = get_food_import_data() ## With Mexico: ## m = lm(Number_of_shipments_refused ~ Total_value_food_imports, data=DF) print(summary(m)) ## Without Mexico: ## m2 = lm(Number_of_shipments_refused ~ Total_value_food_imports, data=DF[c(-2), ]) print(summary(m2)) ## Plot both models: ## plot(DF$Total_value_food_imports, DF$Number_of_shipments_refused, 'p', pch=19, xlab='Total Value Food Imports (Billions of Dollars)', ylab='Number of Shipments Refused (Count)') grid() abline(m, col='red') abline(m2, col='green')