Title: DATA 205 Capstone - The Local Food Movement and the Impact on Maryland Farmers Markets and Certified Agricultural Producers
Jannety Mosley Spring 2022 —
Data Source:
Data Montgomery Datasets:
-Farmers Market 2022
-Certified Agricultural Producers
Research Questions:
With the increase of the regional food movement have Certified Agricultural Producer Certification increased in Maryland? Is certification renewal being regulated properly?
What is the location distribution of farmers markets in Maryland in comparison to producer with Agricultural Producer Certification?
Load in the Datasets
setwd ("/cloud/project/")
farmmkts <- read.csv("Farmers_Markets_2022.csv")
agprod <- read.csv("Agricultural_Producer__Certificates.csv")Install libraries
library(stringr)
library(ggplot2)
library(explore)
library(scales)
library(tidyverse)## ── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ──## ✔ tibble 3.1.6 ✔ purrr 0.3.4
## ✔ tidyr 1.2.0 ✔ dplyr 1.0.8
## ✔ readr 2.1.2 ✔ forcats 0.5.1## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ readr::col_factor() masks scales::col_factor()
## ✖ purrr::discard() masks scales::discard()
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()Clean and explore the data variables
Data Summary
summary(agprod)## Application.Type License.Number Processed.Date License.Date
## Length:82 Min. :200026 Length:82 Length:82
## Class :character 1st Qu.:201319 Class :character Class :character
## Mode :character Median :215456 Mode :character Mode :character
## Mean :218330
## 3rd Qu.:224988
## Max. :268783
##
## Issued.Date Last.Renewal.Date Next.Renewal.Date
## Length:82 Length:82 Length:82
## Class :character Class :character Class :character
## Mode :character Mode :character Mode :character
##
##
##
##
## License.Expiration.Date IN_COUNTY OUT_COUNTY
## Length:82 Length:82 Length:82
## Class :character Class :character Class :character
## Mode :character Mode :character Mode :character
##
##
##
##
## PRODUCT_SOLD Farm.name Street.Number Pre.direction
## Length:82 Length:82 Min. : 304 Length:82
## Class :character Class :character 1st Qu.: 4719 Class :character
## Mode :character Mode :character Median :10401 Mode :character
## Mean :12144
## 3rd Qu.:19124
## Max. :28912
## NA's :29
## Street.name.of.work.location Street.Suffix Post.direction
## Length:82 Length:82 Mode:logical
## Class :character Class :character NA's:82
## Mode :character Mode :character
##
##
##
##
## City State ZIP.code Work.Location
## Length:82 Length:82 Min. :19947 Length:82
## Class :character Class :character 1st Qu.:20835 Class :character
## Mode :character Mode :character Median :20866 Mode :character
## Mean :20860
## 3rd Qu.:20876
## Max. :21771
## NA's :27summary(farmmkts)## Organization.s. First.Name Last.Name Alt.Address
## Length:21 Length:21 Length:21 Length:21
## Class :character Class :character Class :character Class :character
## Mode :character Mode :character Mode :character Mode :character
##
##
##
## City State Zip.Code Email
## Length:21 Length:21 Min. :20814 Length:21
## Class :character Class :character 1st Qu.:20850 Class :character
## Mode :character Mode :character Median :20855 Mode :character
## Mean :20866
## 3rd Qu.:20886
## Max. :20912
## Phone website Days.open year.round
## Length:21 Length:21 Length:21 Length:21
## Class :character Class :character Class :character Class :character
## Mode :character Mode :character Mode :character Mode :character
##
##
##
## Season.1 Time.1 Season.2 Time.2
## Length:21 Length:21 Length:21 Length:21
## Class :character Class :character Class :character Class :character
## Mode :character Mode :character Mode :character Mode :character
##
##
##
## Programs.Accepted Georeference
## Length:21 Length:21
## Class :character Class :character
## Mode :character Mode :character
##
##
## Data Description
agprod %>% describe()## # A tibble: 21 × 8
## variable type na na_pct unique min mean max
## <chr> <chr> <int> <dbl> <int> <dbl> <dbl> <dbl>
## 1 Application.Type chr 0 0 1 NA NA NA
## 2 License.Number int 0 0 82 200026 218330. 268783
## 3 Processed.Date chr 0 0 75 NA NA NA
## 4 License.Date chr 0 0 75 NA NA NA
## 5 Issued.Date chr 0 0 64 NA NA NA
## 6 Last.Renewal.Date chr 0 0 27 NA NA NA
## 7 Next.Renewal.Date chr 0 0 69 NA NA NA
## 8 License.Expiration.Date chr 0 0 69 NA NA NA
## 9 IN_COUNTY chr 0 0 2 NA NA NA
## 10 OUT_COUNTY chr 0 0 2 NA NA NA
## # … with 11 more rowsfarmmkts %>% describe()## # A tibble: 18 × 8
## variable type na na_pct unique min mean max
## <chr> <chr> <int> <dbl> <int> <dbl> <dbl> <dbl>
## 1 Organization.s. chr 0 0 21 NA NA NA
## 2 First.Name chr 0 0 17 NA NA NA
## 3 Last.Name chr 0 0 17 NA NA NA
## 4 Alt.Address chr 0 0 21 NA NA NA
## 5 City chr 0 0 15 NA NA NA
## 6 State chr 0 0 1 NA NA NA
## 7 Zip.Code int 0 0 17 20814 20866. 20912
## 8 Email chr 0 0 17 NA NA NA
## 9 Phone chr 0 0 17 NA NA NA
## 10 website chr 0 0 18 NA NA NA
## 11 Days.open chr 0 0 4 NA NA NA
## 12 year.round chr 0 0 3 NA NA NA
## 13 Season.1 chr 0 0 16 NA NA NA
## 14 Time.1 chr 0 0 15 NA NA NA
## 15 Season.2 chr 0 0 5 NA NA NA
## 16 Time.2 chr 0 0 5 NA NA NA
## 17 Programs.Accepted chr 0 0 11 NA NA NA
## 18 Georeference chr 0 0 21 NA NA NABlank Cells to NA
agprod[agprod == ""] <- NA
farmmkts[farmmkts == ""] <- NA Create a new variable indicating whether an Ag Producer was in Montgomery County or not
agprod$IN_COUNTY <- c(Y = "IN_COUNTY", N = "OUT_COUNTY")
agprod$OUT_COUNTY <- c(Y = "OUT_COUNTY", N = "IN_COUNTY")
summary(agprod$IN_COUNTY)## Length Class Mode
## 82 character charactersummary(agprod$OUT_COUNTY)## Length Class Mode
## 82 character characterRecoding Produce Sold into categories
agprod$PRODUCT_SOLD2[agprod$PRODUCT_SOLD == "PRODUCE"] = "1"
agprod$PRODUCT_SOLD2[agprod$PRODUCT_SOLD == "SEAFOOD"] = "2"
agprod$PRODUCT_SOLD2[agprod$PRODUCT_SOLD == "TREES"] = "3"
agprod$PRODUCT_SOLD2[agprod$PRODUCT_SOLD == "FLOWERS AND PRODUCE"] = "4"
agprod$PRODUCT_SOLD2[agprod$PRODUCT_SOLD == "PRODUCE AND TREES"] = "4"
agprod$PRODUCT_SOLD2[agprod$PRODUCT_SOLD == "FLOWERS, PRODUCE AND TREES"] = "5"Extract Year from Processed Date
library(lubridate)##
## Attaching package: 'lubridate'## The following objects are masked from 'package:base':
##
## date, intersect, setdiff, union# extract year from processed date
agprod$ProcessedYear <- year(mdy(agprod$Processed.Date))
# spot check for quality; yes, showing the year
head(agprod)## Application.Type License.Number Processed.Date License.Date
## 1 AGRICULTURAL PRODUCERS CERTIF 208918 07/23/2007 07/23/2007
## 2 AGRICULTURAL PRODUCERS CERTIF 264971 02/13/2020 02/13/2020
## 3 AGRICULTURAL PRODUCERS CERTIF 265612 05/18/2020 05/18/2020
## 4 AGRICULTURAL PRODUCERS CERTIF 265728 06/24/2020 06/24/2020
## 5 AGRICULTURAL PRODUCERS CERTIF 200080 12/01/1999 12/01/1999
## 6 AGRICULTURAL PRODUCERS CERTIF 202287 03/12/2003 03/12/2003
## Issued.Date Last.Renewal.Date Next.Renewal.Date License.Expiration.Date
## 1 08/01/2016 08/10/2020 08/01/2021 08/01/2021
## 2 02/20/2020 <NA> 02/20/2021 02/20/2021
## 3 05/21/2020 <NA> 05/21/2021 05/21/2021
## 4 06/24/2020 <NA> 06/24/2021 06/24/2021
## 5 12/01/1999 <NA> 12/07/2001 12/07/2001
## 6 03/13/2003 03/11/2005 03/11/2006 03/11/2006
## IN_COUNTY OUT_COUNTY PRODUCT_SOLD Farm.name Street.Number
## 1 IN_COUNTY OUT_COUNTY PRODUCE CENTURION FARM 7300
## 2 OUT_COUNTY IN_COUNTY PRODUCE SANDY SPRING GARDENS 304
## 3 IN_COUNTY OUT_COUNTY PRODUCE AIX LA CHAPELLE FARM 19124
## 4 OUT_COUNTY IN_COUNTY PRODUCE THE GARDEN 27455
## 5 IN_COUNTY OUT_COUNTY TREES MUSSER FOREST 6900
## 6 OUT_COUNTY IN_COUNTY SEAFOOD E. GOODWIN & SONS (SEAFOOD) NA
## Pre.direction Street.name.of.work.location Street.Suffix Post.direction
## 1 <NA> MACARTHUR BLVD NA
## 2 <NA> LETHBRIDGE CT NA
## 3 <NA> JERUSALEM RD NA
## 4 <NA> CLARKSBURG RD NA
## 5 <NA> NEW HAMPSHIRE AVE NA
## 6 <NA> <NA> <NA> NA
## City State ZIP.code Work.Location PRODUCT_SOLD2 ProcessedYear
## 1 GLEN ECHO MD 20812 <NA> 1 2007
## 2 ASHTON MD 20861 <NA> 1 2020
## 3 POOLESVILLE MD 20837 <NA> 1 2020
## 4 DAMASCUS MD 20872 <NA> 1 2020
## 5 TAKOMA PARK MD 20912 <NA> 3 1999
## 6 <NA> <NA> NA 29 AND TECH ROAD 2 2003Distribution of Certified Ag Producers by Certification Year
cleanup = theme(panel.grid.major =element_blank(),
panel.grid.minor =element_blank(),
panel.background =element_blank(),
axis.line.x =element_line(color ="black"),
axis.line.y =element_line(color ="black"),
legend.key =element_rect(fill ="white"),
text =element_text(size =15))
#Distribution of price
par(mfrow=c(2,1))
ggplot(agprod) +
cleanup+
geom_histogram(aes(ProcessedYear),fill = 'orange',alpha = 0.85,binwidth = 15) +
theme_minimal(base_size = 13) + xlab("Processed Year") + ylab("Frequency") +
ggtitle("The Distrubition of Certification Year") 
Transformed distribution of Price
ggplot(agprod, aes(ProcessedYear)) +
cleanup+
geom_histogram(bins = 30, aes(y = ..density..), fill = "orange") +
geom_density(alpha = 0.2, fill = "orange") +ggtitle("Transformed distribution of certification year",
subtitle = expression("With" ~'log'[10] ~ "transformation of x-axis")) + scale_x_log10()
## Distribution of Certified Agricultural Producers In Montgomery County
by Year
# Plot of certified agricultural producers broken down by Year
ggplot(agprod, aes(x = ProcessedYear)) +
geom_bar() +
facet_wrap(~ IN_COUNTY)+
theme_bw()+
theme(strip.background =element_rect(fill="orange"))+
theme(strip.text = element_text(colour = 'white'))
Extract Year from License Expiration Date
# extract year from license expiration date
agprod$ExpirationYear <- year(mdy(agprod$License.Expiration.Date))
# spot check for quality; yes, showing the year
head(agprod)## Application.Type License.Number Processed.Date License.Date
## 1 AGRICULTURAL PRODUCERS CERTIF 208918 07/23/2007 07/23/2007
## 2 AGRICULTURAL PRODUCERS CERTIF 264971 02/13/2020 02/13/2020
## 3 AGRICULTURAL PRODUCERS CERTIF 265612 05/18/2020 05/18/2020
## 4 AGRICULTURAL PRODUCERS CERTIF 265728 06/24/2020 06/24/2020
## 5 AGRICULTURAL PRODUCERS CERTIF 200080 12/01/1999 12/01/1999
## 6 AGRICULTURAL PRODUCERS CERTIF 202287 03/12/2003 03/12/2003
## Issued.Date Last.Renewal.Date Next.Renewal.Date License.Expiration.Date
## 1 08/01/2016 08/10/2020 08/01/2021 08/01/2021
## 2 02/20/2020 <NA> 02/20/2021 02/20/2021
## 3 05/21/2020 <NA> 05/21/2021 05/21/2021
## 4 06/24/2020 <NA> 06/24/2021 06/24/2021
## 5 12/01/1999 <NA> 12/07/2001 12/07/2001
## 6 03/13/2003 03/11/2005 03/11/2006 03/11/2006
## IN_COUNTY OUT_COUNTY PRODUCT_SOLD Farm.name Street.Number
## 1 IN_COUNTY OUT_COUNTY PRODUCE CENTURION FARM 7300
## 2 OUT_COUNTY IN_COUNTY PRODUCE SANDY SPRING GARDENS 304
## 3 IN_COUNTY OUT_COUNTY PRODUCE AIX LA CHAPELLE FARM 19124
## 4 OUT_COUNTY IN_COUNTY PRODUCE THE GARDEN 27455
## 5 IN_COUNTY OUT_COUNTY TREES MUSSER FOREST 6900
## 6 OUT_COUNTY IN_COUNTY SEAFOOD E. GOODWIN & SONS (SEAFOOD) NA
## Pre.direction Street.name.of.work.location Street.Suffix Post.direction
## 1 <NA> MACARTHUR BLVD NA
## 2 <NA> LETHBRIDGE CT NA
## 3 <NA> JERUSALEM RD NA
## 4 <NA> CLARKSBURG RD NA
## 5 <NA> NEW HAMPSHIRE AVE NA
## 6 <NA> <NA> <NA> NA
## City State ZIP.code Work.Location PRODUCT_SOLD2 ProcessedYear
## 1 GLEN ECHO MD 20812 <NA> 1 2007
## 2 ASHTON MD 20861 <NA> 1 2020
## 3 POOLESVILLE MD 20837 <NA> 1 2020
## 4 DAMASCUS MD 20872 <NA> 1 2020
## 5 TAKOMA PARK MD 20912 <NA> 3 1999
## 6 <NA> <NA> NA 29 AND TECH ROAD 2 2003
## ExpirationYear
## 1 2021
## 2 2021
## 3 2021
## 4 2021
## 5 2001
## 6 2006Processed Year Table
table(agprod$ProcessedYear)##
## 1999 2000 2001 2002 2003 2004 2006 2007 2008 2009 2011 2012 2013 2014 2015 2016
## 11 10 2 1 3 1 4 2 2 4 4 1 5 2 7 7
## 2017 2018 2019 2020 2021 2022
## 1 5 2 3 4 1Processed Year Filter
agprod %>%
filter(agprod$ProcessedYear %in% c('2020', '2021', '2022'))## Application.Type License.Number Processed.Date License.Date
## 1 AGRICULTURAL PRODUCERS CERTIF 264971 02/13/2020 02/13/2020
## 2 AGRICULTURAL PRODUCERS CERTIF 265612 05/18/2020 05/18/2020
## 3 AGRICULTURAL PRODUCERS CERTIF 265728 06/24/2020 06/24/2020
## 4 AGRICULTURAL PRODUCERS CERTIF 267510 04/27/2021 04/27/2021
## 5 AGRICULTURAL PRODUCERS CERTIF 267835 07/21/2021 07/21/2021
## 6 AGRICULTURAL PRODUCERS CERTIF 267833 07/21/2021 07/21/2021
## 7 AGRICULTURAL PRODUCERS CERTIF 268783 03/01/2022 03/01/2022
## 8 AGRICULTURAL PRODUCERS CERTIF 267355 03/15/2021 03/15/2021
## Issued.Date Last.Renewal.Date Next.Renewal.Date License.Expiration.Date
## 1 02/20/2020 <NA> 02/20/2021 02/20/2021
## 2 05/21/2020 <NA> 05/21/2021 05/21/2021
## 3 06/24/2020 <NA> 06/24/2021 06/24/2021
## 4 05/11/2021 <NA> 05/11/2022 05/11/2022
## 5 07/27/2021 <NA> 07/27/2022 07/27/2022
## 6 07/30/2021 <NA> 07/30/2022 07/30/2022
## 7 03/01/2022 <NA> 03/01/2023 03/01/2023
## 8 03/22/2021 03/30/2022 03/22/2023 03/22/2023
## IN_COUNTY OUT_COUNTY PRODUCT_SOLD Farm.name
## 1 OUT_COUNTY IN_COUNTY PRODUCE SANDY SPRING GARDENS
## 2 IN_COUNTY OUT_COUNTY PRODUCE AIX LA CHAPELLE FARM
## 3 OUT_COUNTY IN_COUNTY PRODUCE THE GARDEN
## 4 OUT_COUNTY IN_COUNTY PRODUCE JAGWOOD FARMS
## 5 IN_COUNTY OUT_COUNTY <NA> BAUGHERS ORCHARD
## 6 OUT_COUNTY IN_COUNTY <NA> BAUGHERS ORCHARD
## 7 IN_COUNTY OUT_COUNTY PRODUCE WATKINS HOUSE PRODUCE LLC
## 8 OUT_COUNTY IN_COUNTY FLOWERS, PRODUCE AND TREES BAUGHERS ORCHARD
## Street.Number Pre.direction Street.name.of.work.location Street.Suffix
## 1 304 <NA> LETHBRIDGE CT
## 2 19124 <NA> JERUSALEM RD
## 3 27455 <NA> CLARKSBURG RD
## 4 5511 <NA> MASSACHUSETTS AVE
## 5 7300 <NA> MAC ARTHUR BLVD
## 6 12300 <NA> FALLS RD
## 7 27130 <NA> RIDGE RD
## 8 9045 <NA> CONGRESSIONAL PKWY
## Post.direction City State ZIP.code Work.Location PRODUCT_SOLD2
## 1 NA ASHTON MD 20861 <NA> 1
## 2 NA POOLESVILLE MD 20837 <NA> 1
## 3 NA DAMASCUS MD 20872 <NA> 1
## 4 NA BETHESDA MD 20816 <NA> 1
## 5 NA GLEN ECHO MD 20812 <NA> <NA>
## 6 NA POTOMAC MD 20854 <NA> <NA>
## 7 NA DAMASCUS MD 20872 <NA> 1
## 8 NA POTOMAC MD 20854 <NA> 5
## ProcessedYear ExpirationYear
## 1 2020 2021
## 2 2020 2021
## 3 2020 2021
## 4 2021 2022
## 5 2021 2022
## 6 2021 2022
## 7 2022 2023
## 8 2021 2023Expiration Year Table
table(agprod$ExpirationYear)##
## 2000 2001 2002 2003 2006 2007 2008 2009 2011 2012 2013 2014 2015 2016 2017 2018
## 2 8 4 4 5 1 3 2 2 2 1 2 2 5 6 1
## 2019 2020 2021 2022 2023
## 6 4 5 3 2Expiration Year Filter
agprod %>%
filter(agprod$ExpirationYear %in% c('2020', '2021', '2022', '2023'))## Application.Type License.Number Processed.Date License.Date
## 1 AGRICULTURAL PRODUCERS CERTIF 208918 07/23/2007 07/23/2007
## 2 AGRICULTURAL PRODUCERS CERTIF 264971 02/13/2020 02/13/2020
## 3 AGRICULTURAL PRODUCERS CERTIF 265612 05/18/2020 05/18/2020
## 4 AGRICULTURAL PRODUCERS CERTIF 265728 06/24/2020 06/24/2020
## 5 AGRICULTURAL PRODUCERS CERTIF 222846 04/28/2015 04/28/2015
## 6 AGRICULTURAL PRODUCERS CERTIF 255274 05/24/2019 05/24/2019
## 7 AGRICULTURAL PRODUCERS CERTIF 226209 07/09/2018 07/09/2018
## 8 AGRICULTURAL PRODUCERS CERTIF 225042 07/12/2016 07/12/2016
## 9 AGRICULTURAL PRODUCERS CERTIF 221765 09/05/2014 09/05/2014
## 10 AGRICULTURAL PRODUCERS CERTIF 267510 04/27/2021 04/27/2021
## 11 AGRICULTURAL PRODUCERS CERTIF 267835 07/21/2021 07/21/2021
## 12 AGRICULTURAL PRODUCERS CERTIF 267833 07/21/2021 07/21/2021
## 13 AGRICULTURAL PRODUCERS CERTIF 268783 03/01/2022 03/01/2022
## 14 AGRICULTURAL PRODUCERS CERTIF 267355 03/15/2021 03/15/2021
## Issued.Date Last.Renewal.Date Next.Renewal.Date License.Expiration.Date
## 1 08/01/2016 08/10/2020 08/01/2021 08/01/2021
## 2 02/20/2020 <NA> 02/20/2021 02/20/2021
## 3 05/21/2020 <NA> 05/21/2021 05/21/2021
## 4 06/24/2020 <NA> 06/24/2021 06/24/2021
## 5 05/04/2015 05/28/2019 05/04/2020 05/04/2020
## 6 05/28/2019 <NA> 05/28/2020 05/28/2020
## 7 07/12/2018 07/17/2019 07/12/2020 07/12/2020
## 8 07/15/2016 08/07/2019 07/15/2020 07/15/2020
## 9 09/08/2014 03/22/2021 09/08/2021 09/08/2021
## 10 05/11/2021 <NA> 05/11/2022 05/11/2022
## 11 07/27/2021 <NA> 07/27/2022 07/27/2022
## 12 07/30/2021 <NA> 07/30/2022 07/30/2022
## 13 03/01/2022 <NA> 03/01/2023 03/01/2023
## 14 03/22/2021 03/30/2022 03/22/2023 03/22/2023
## IN_COUNTY OUT_COUNTY PRODUCT_SOLD Farm.name
## 1 IN_COUNTY OUT_COUNTY PRODUCE CENTURION FARM
## 2 OUT_COUNTY IN_COUNTY PRODUCE SANDY SPRING GARDENS
## 3 IN_COUNTY OUT_COUNTY PRODUCE AIX LA CHAPELLE FARM
## 4 OUT_COUNTY IN_COUNTY PRODUCE THE GARDEN
## 5 OUT_COUNTY IN_COUNTY PRODUCE NORMAN'S FARM
## 6 IN_COUNTY OUT_COUNTY PRODUCE SOMERSET PRODUCE
## 7 OUT_COUNTY IN_COUNTY PRODUCE PARKVIEW FARMS
## 8 IN_COUNTY OUT_COUNTY PRODUCE MOHOLT PRODUCE
## 9 IN_COUNTY OUT_COUNTY PRODUCE BAUGHERS' ENTERPRISES INC
## 10 OUT_COUNTY IN_COUNTY PRODUCE JAGWOOD FARMS
## 11 IN_COUNTY OUT_COUNTY <NA> BAUGHERS ORCHARD
## 12 OUT_COUNTY IN_COUNTY <NA> BAUGHERS ORCHARD
## 13 IN_COUNTY OUT_COUNTY PRODUCE WATKINS HOUSE PRODUCE LLC
## 14 OUT_COUNTY IN_COUNTY FLOWERS, PRODUCE AND TREES BAUGHERS ORCHARD
## Street.Number Pre.direction Street.name.of.work.location Street.Suffix
## 1 7300 <NA> MACARTHUR BLVD
## 2 304 <NA> LETHBRIDGE CT
## 3 19124 <NA> JERUSALEM RD
## 4 27455 <NA> CLARKSBURG RD
## 5 NA <NA> <NA> <NA>
## 6 2600 <NA> BRIGGS CHANEY RD
## 7 17820 <NA> WHITES FERRY RD
## 8 23230 <NA> RIDGE RD
## 9 10900 <NA> SEVEN LOCKS RD
## 10 5511 <NA> MASSACHUSETTS AVE
## 11 7300 <NA> MAC ARTHUR BLVD
## 12 12300 <NA> FALLS RD
## 13 27130 <NA> RIDGE RD
## 14 9045 <NA> CONGRESSIONAL PKWY
## Post.direction City State ZIP.code
## 1 NA GLEN ECHO MD 20812
## 2 NA ASHTON MD 20861
## 3 NA POOLESVILLE MD 20837
## 4 NA DAMASCUS MD 20872
## 5 NA <NA> <NA> NA
## 6 NA SILVER SPRING MD 20905
## 7 NA POOLESVILLE MD 20837
## 8 NA GERMANTOWN MD 20876
## 9 NA POTOMAC MD 20854
## 10 NA BETHESDA MD 20816
## 11 NA GLEN ECHO MD 20812
## 12 NA POTOMAC MD 20854
## 13 NA DAMASCUS MD 20872
## 14 NA POTOMAC MD 20854
## Work.Location
## 1 <NA>
## 2 <NA>
## 3 <NA>
## 4 <NA>
## 5 CORNER OF MASSACHUSETTS AVENUE - ENTRANCE OF WESTLAND MIDDLE SCHOOL; AGRIC
## 6 <NA>
## 7 <NA>
## 8 <NA>
## 9 <NA>
## 10 <NA>
## 11 <NA>
## 12 <NA>
## 13 <NA>
## 14 <NA>
## PRODUCT_SOLD2 ProcessedYear ExpirationYear
## 1 1 2007 2021
## 2 1 2020 2021
## 3 1 2020 2021
## 4 1 2020 2021
## 5 1 2015 2020
## 6 1 2019 2020
## 7 1 2018 2020
## 8 1 2016 2020
## 9 1 2014 2021
## 10 1 2021 2022
## 11 <NA> 2021 2022
## 12 <NA> 2021 2022
## 13 1 2022 2023
## 14 5 2021 2023Distribution of Certified Producers by Expiration Year
cleanup = theme(panel.grid.major =element_blank(),
panel.grid.minor =element_blank(),
panel.background =element_blank(),
axis.line.x =element_line(color ="black"),
axis.line.y =element_line(color ="black"),
legend.key =element_rect(fill ="white"),
text =element_text(size =15))
#Distribution of price
par(mfrow=c(2,1))
ggplot(agprod) +
cleanup+
geom_histogram(aes(ExpirationYear),fill = 'blue',alpha = 0.85,binwidth = 15) +
theme_minimal(base_size = 13) + xlab("Expiration Year") + ylab("Frequency") +
ggtitle("The Distrubition of Expiration Year") ## Warning: Removed 12 rows containing non-finite values (stat_bin).
Filter for License Expiration Date
agprod %>%
filter(ExpirationYear %in% c('2020', '2021', '2022', '2023'))## Application.Type License.Number Processed.Date License.Date
## 1 AGRICULTURAL PRODUCERS CERTIF 208918 07/23/2007 07/23/2007
## 2 AGRICULTURAL PRODUCERS CERTIF 264971 02/13/2020 02/13/2020
## 3 AGRICULTURAL PRODUCERS CERTIF 265612 05/18/2020 05/18/2020
## 4 AGRICULTURAL PRODUCERS CERTIF 265728 06/24/2020 06/24/2020
## 5 AGRICULTURAL PRODUCERS CERTIF 222846 04/28/2015 04/28/2015
## 6 AGRICULTURAL PRODUCERS CERTIF 255274 05/24/2019 05/24/2019
## 7 AGRICULTURAL PRODUCERS CERTIF 226209 07/09/2018 07/09/2018
## 8 AGRICULTURAL PRODUCERS CERTIF 225042 07/12/2016 07/12/2016
## 9 AGRICULTURAL PRODUCERS CERTIF 221765 09/05/2014 09/05/2014
## 10 AGRICULTURAL PRODUCERS CERTIF 267510 04/27/2021 04/27/2021
## 11 AGRICULTURAL PRODUCERS CERTIF 267835 07/21/2021 07/21/2021
## 12 AGRICULTURAL PRODUCERS CERTIF 267833 07/21/2021 07/21/2021
## 13 AGRICULTURAL PRODUCERS CERTIF 268783 03/01/2022 03/01/2022
## 14 AGRICULTURAL PRODUCERS CERTIF 267355 03/15/2021 03/15/2021
## Issued.Date Last.Renewal.Date Next.Renewal.Date License.Expiration.Date
## 1 08/01/2016 08/10/2020 08/01/2021 08/01/2021
## 2 02/20/2020 <NA> 02/20/2021 02/20/2021
## 3 05/21/2020 <NA> 05/21/2021 05/21/2021
## 4 06/24/2020 <NA> 06/24/2021 06/24/2021
## 5 05/04/2015 05/28/2019 05/04/2020 05/04/2020
## 6 05/28/2019 <NA> 05/28/2020 05/28/2020
## 7 07/12/2018 07/17/2019 07/12/2020 07/12/2020
## 8 07/15/2016 08/07/2019 07/15/2020 07/15/2020
## 9 09/08/2014 03/22/2021 09/08/2021 09/08/2021
## 10 05/11/2021 <NA> 05/11/2022 05/11/2022
## 11 07/27/2021 <NA> 07/27/2022 07/27/2022
## 12 07/30/2021 <NA> 07/30/2022 07/30/2022
## 13 03/01/2022 <NA> 03/01/2023 03/01/2023
## 14 03/22/2021 03/30/2022 03/22/2023 03/22/2023
## IN_COUNTY OUT_COUNTY PRODUCT_SOLD Farm.name
## 1 IN_COUNTY OUT_COUNTY PRODUCE CENTURION FARM
## 2 OUT_COUNTY IN_COUNTY PRODUCE SANDY SPRING GARDENS
## 3 IN_COUNTY OUT_COUNTY PRODUCE AIX LA CHAPELLE FARM
## 4 OUT_COUNTY IN_COUNTY PRODUCE THE GARDEN
## 5 OUT_COUNTY IN_COUNTY PRODUCE NORMAN'S FARM
## 6 IN_COUNTY OUT_COUNTY PRODUCE SOMERSET PRODUCE
## 7 OUT_COUNTY IN_COUNTY PRODUCE PARKVIEW FARMS
## 8 IN_COUNTY OUT_COUNTY PRODUCE MOHOLT PRODUCE
## 9 IN_COUNTY OUT_COUNTY PRODUCE BAUGHERS' ENTERPRISES INC
## 10 OUT_COUNTY IN_COUNTY PRODUCE JAGWOOD FARMS
## 11 IN_COUNTY OUT_COUNTY <NA> BAUGHERS ORCHARD
## 12 OUT_COUNTY IN_COUNTY <NA> BAUGHERS ORCHARD
## 13 IN_COUNTY OUT_COUNTY PRODUCE WATKINS HOUSE PRODUCE LLC
## 14 OUT_COUNTY IN_COUNTY FLOWERS, PRODUCE AND TREES BAUGHERS ORCHARD
## Street.Number Pre.direction Street.name.of.work.location Street.Suffix
## 1 7300 <NA> MACARTHUR BLVD
## 2 304 <NA> LETHBRIDGE CT
## 3 19124 <NA> JERUSALEM RD
## 4 27455 <NA> CLARKSBURG RD
## 5 NA <NA> <NA> <NA>
## 6 2600 <NA> BRIGGS CHANEY RD
## 7 17820 <NA> WHITES FERRY RD
## 8 23230 <NA> RIDGE RD
## 9 10900 <NA> SEVEN LOCKS RD
## 10 5511 <NA> MASSACHUSETTS AVE
## 11 7300 <NA> MAC ARTHUR BLVD
## 12 12300 <NA> FALLS RD
## 13 27130 <NA> RIDGE RD
## 14 9045 <NA> CONGRESSIONAL PKWY
## Post.direction City State ZIP.code
## 1 NA GLEN ECHO MD 20812
## 2 NA ASHTON MD 20861
## 3 NA POOLESVILLE MD 20837
## 4 NA DAMASCUS MD 20872
## 5 NA <NA> <NA> NA
## 6 NA SILVER SPRING MD 20905
## 7 NA POOLESVILLE MD 20837
## 8 NA GERMANTOWN MD 20876
## 9 NA POTOMAC MD 20854
## 10 NA BETHESDA MD 20816
## 11 NA GLEN ECHO MD 20812
## 12 NA POTOMAC MD 20854
## 13 NA DAMASCUS MD 20872
## 14 NA POTOMAC MD 20854
## Work.Location
## 1 <NA>
## 2 <NA>
## 3 <NA>
## 4 <NA>
## 5 CORNER OF MASSACHUSETTS AVENUE - ENTRANCE OF WESTLAND MIDDLE SCHOOL; AGRIC
## 6 <NA>
## 7 <NA>
## 8 <NA>
## 9 <NA>
## 10 <NA>
## 11 <NA>
## 12 <NA>
## 13 <NA>
## 14 <NA>
## PRODUCT_SOLD2 ProcessedYear ExpirationYear
## 1 1 2007 2021
## 2 1 2020 2021
## 3 1 2020 2021
## 4 1 2020 2021
## 5 1 2015 2020
## 6 1 2019 2020
## 7 1 2018 2020
## 8 1 2016 2020
## 9 1 2014 2021
## 10 1 2021 2022
## 11 <NA> 2021 2022
## 12 <NA> 2021 2022
## 13 1 2022 2023
## 14 5 2021 2023In_County Recode
agprod$IN_COUNTY2[agprod$IN_COUNTY == "OUT_COUNTY"] = "0"
agprod$IN_COUNTY2[agprod$IN_COUNTY == "IN_COUNTY"] = "1"One-way ANOVA (Expiration Year)
one.way_EY <- aov(agprod$ExpirationYear ~ agprod$PRODUCT_SOLD2, data = agprod)
summary(one.way_EY)## Df Sum Sq Mean Sq F value Pr(>F)
## agprod$PRODUCT_SOLD2 4 1097 274.31 6.315 0.000245 ***
## Residuals 63 2737 43.44
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 14 observations deleted due to missingnessTwo-way ANOVA R Code (Expiration Year)
two.way_EY <- aov(agprod$ExpirationYear ~ agprod$PRODUCT_SOLD2 + agprod$IN_COUNTY2, data = agprod)
summary(two.way_EY)## Df Sum Sq Mean Sq F value Pr(>F)
## agprod$PRODUCT_SOLD2 4 1097.2 274.31 6.217 0.000285 ***
## agprod$IN_COUNTY2 1 0.8 0.77 0.017 0.895550
## Residuals 62 2735.8 44.13
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 14 observations deleted due to missingnessTwo-way ANOVA Interaction (Expiration Year)
interaction_EY <- aov(agprod$ExpirationYear ~ agprod$PRODUCT_SOLD2*agprod$IN_COUNTY2, data = agprod)
summary(interaction_EY)## Df Sum Sq Mean Sq F value Pr(>F)
## agprod$PRODUCT_SOLD2 4 1097.2 274.31 6.005 0.00041 ***
## agprod$IN_COUNTY2 1 0.8 0.77 0.017 0.89736
## agprod$PRODUCT_SOLD2:agprod$IN_COUNTY2 4 86.4 21.61 0.473 0.75532
## Residuals 58 2649.4 45.68
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 14 observations deleted due to missingnessOne-way ANOVA (Processed Year)
one.way_PY <- aov(agprod$ProcessedYear ~ agprod$PRODUCT_SOLD2, data = agprod)
summary(one.way_PY)## Df Sum Sq Mean Sq F value Pr(>F)
## agprod$PRODUCT_SOLD2 4 970 242.42 5.145 0.00102 **
## Residuals 75 3534 47.12
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 2 observations deleted due to missingnessTwo-way ANOVA R Code (Processed Year)
two.way_PY <- aov(agprod$ProcessedYear ~ agprod$PRODUCT_SOLD2 + agprod$IN_COUNTY2, data = agprod)
summary(two.way_PY)## Df Sum Sq Mean Sq F value Pr(>F)
## agprod$PRODUCT_SOLD2 4 970 242.42 5.081 0.00113 **
## agprod$IN_COUNTY2 1 3 2.92 0.061 0.80534
## Residuals 74 3531 47.71
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 2 observations deleted due to missingnessTwo-way ANOVA Interaction (Processed Year)
interaction_PY <- aov(agprod$ProcessedYear ~ agprod$PRODUCT_SOLD2*agprod$IN_COUNTY2, data = agprod)
summary(interaction_PY)## Df Sum Sq Mean Sq F value Pr(>F)
## agprod$PRODUCT_SOLD2 4 970 242.42 5.100 0.00115 **
## agprod$IN_COUNTY2 1 3 2.92 0.061 0.80502
## agprod$PRODUCT_SOLD2:agprod$IN_COUNTY2 4 204 50.91 1.071 0.37739
## Residuals 70 3327 47.53
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 2 observations deleted due to missingnessLeaflet Certified Ag Producers Work Locaation Map
library(zipcodeR)
library(leaflet)
library(leafsync)
Map_A <- agprod %>%
left_join(.,
geocode_zip(agprod$ZIP.code) %>% mutate(zipcode = as.numeric(zipcode)),
by = c("ZIP.code" = "zipcode")) %>%
leaflet() %>%
addTiles() %>%
addMarkers(
~ lng,
~ lat,
popup = ~ as.character(Street.name.of.work.location),
label = ~ as.character(Street.name.of.work.location)
)## Warning in validateCoords(lng, lat, funcName): Data contains 27 rows with either
## missing or invalid lat/lon values and will be ignoredMap_B <- farmmkts %>%
left_join(.,
geocode_zip(farmmkts$Zip.Code) %>% mutate(zipcode = as.numeric(zipcode)),
by = c("Zip.Code" = "zipcode")) %>%
leaflet() %>%
addTiles() %>%
addMarkers(
~ lng,
~ lat,
popup = ~ as.character(Organization.s.),
label = ~ as.character(Organization.s.)
)
sync(Map_A,Map_B)