ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = NA,
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
setwd("C:/Users/USUARIO/Desktop/OBservData/Datasets_storage")
write_csv(field_level_data, "field_level_data_hevi01.csv")
setwd(dir_ini)
View(field_level_data)
library(tidyverse)
library(openxlsx)
library(sp) #Transforming latitude and longitude
dir_ini <- getwd()
options(digits=14)
datafield <- read.xlsx("Crop_pollination_database_Hevia_et_al_UPDATE.xlsx",
sheet = "field_level_data", startRow = 1)
datafield <- as_tibble(datafield)
datafield <- datafield %>%
rename(field_size = field.size,
fruits_per_plant = mean_fruits_per_plant,
richness_estimator_method = `richness_estimator_&#10;Method`)
datafield$management <- "conventional"
datafield$crop <- "Helianthus annuus"
datafield$study_id <- "hevi01"
chd = substr(datafield$latitude, 3, 3)[1]
chm = substr(datafield$latitude, 6, 6)[1]
chs = substr(datafield$latitude, 11, 11)[1]
cd = char2dms(datafield$latitude,chd=chd,chm=chm,chs=chs)
datafield$latitude <- as.numeric(cd)
chd = substr(datafield$longitude, 2, 2)[1]
chm = substr(datafield$longitude, 5, 5)[1]
chs = substr(datafield$longitude, 10, 10)[1]
cd = char2dms(datafield$longitude,chd = chd,chm = chm,chs = chs)
datafield$longitude <- as.numeric(cd)
field_level_data <- tibble(
study_id = datafield$study_id,
site_id = datafield$site_id,
crop = datafield$crop,
variety = datafield$variety,
management = datafield$management,
country = datafield$country,
latitude = datafield$latitude,
longitude = datafield$longitude,
X_UTM=NA,
Y_UTM=NA,
zone_UTM=NA,
sampling_start_month = datafield$sampling_start_month,
sampling_end_month = datafield$sampling_end_month,
sampling_year = datafield$sampling_year,
field_size = datafield$field_size,
yield=datafield$yield,
yield_units=datafield$yield_units,
yield2=datafield$yield2,
yield2_units=datafield$yield2_units,
yield_treatments_no_pollinators=datafield$yield_treatments_no_pollinators,
yield_treatments_pollen_supplement=NA,
yield_treatments_no_pollinators2=datafield$yield_treatments_no_pollinators2,
yield_treatments_pollen_supplement2=NA,
fruits_per_plant=datafield$fruits_per_plant,
fruit_weight= datafield$fruit_weight,
plant_density=datafield$plant_density,
seeds_per_fruit=datafield$seeds_per_fruit,
seeds_per_plant=datafield$seeds_per_plant,
seed_weight=datafield$seed_weight,
pollinator_richness = datafield$pollinator_richness,
richness_estimator_method = datafield$richness_estimator_method,
abundance = datafield$abundance,
ab_honeybee = datafield$ab_honeybee,
ab_bombus = datafield$ab_bombus,
ab_wildbees = datafield$ab_wildbees,
ab_syrphids = datafield$ab_syrphids,
ab_humbleflies= NA,
ab_other_flies= NA,
ab_beetles=NA,
ab_lepidoptera=NA,
ab_nonbee_hymenoptera=NA,
ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = NA,
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
setwd("C:/Users/USUARIO/Desktop/OBservData/Datasets_storage")
write_csv(field_level_data, "field_level_data_hevi01.csv")
setwd(dir_ini)
View(field_level_data)
field_level_data <- tibble(
study_id = datafield$study_id,
site_id = datafield$site_id,
crop = datafield$crop,
variety = datafield$variety,
management = datafield$management,
country = datafield$country,
latitude = datafield$latitude,
longitude = datafield$longitude,
X_UTM=NA,
Y_UTM=NA,
zone_UTM=NA,
sampling_start_month = datafield$sampling_start_month,
sampling_end_month = datafield$sampling_end_month,
sampling_year = datafield$sampling_year,
field_size = datafield$field_size,
yield=datafield$yield,
yield_units=datafield$yield_units,
yield2=datafield$yield2,
yield2_units=datafield$yield2_units,
yield_treatments_no_pollinators=datafield$yield_treatments_no_pollinators,
yield_treatments_pollen_supplement=NA,
yield_treatments_no_pollinators2=datafield$yield_treatments_no_pollinators2,
yield_treatments_pollen_supplement2=NA,
fruits_per_plant=datafield$fruits_per_plant,
fruit_weight= datafield$fruit_weight,
plant_density=datafield$plant_density,
seeds_per_fruit=datafield$seeds_per_fruit,
seeds_per_plant=datafield$seeds_per_plant,
seed_weight=datafield$seed_weight,
observed_pollinator_richness=NA,
other_pollinator_richness=NA,
other_richness_estimator_method=NA,
abundance = datafield$abundance,
ab_honeybee = datafield$ab_honeybee,
ab_bombus = datafield$ab_bombus,
ab_wildbees = datafield$ab_wildbees,
ab_syrphids = datafield$ab_syrphids,
ab_humbleflies= NA,
ab_other_flies= NA,
ab_beetles=NA,
ab_lepidoptera=NA,
ab_nonbee_hymenoptera=NA,
ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = NA,
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
View(field_level_data)
setwd("C:/Users/USUARIO/Desktop/OBservData/Datasets_storage")
write_csv(field_level_data, "field_level_data_hevi01.csv")
setwd(dir_ini)
library(tidyverse)
library(openxlsx)
library(sp) #Transforming latitude and longitude
dir_ini <- getwd()
options(digits=14)
datafield <- read.xlsx("Crop_pollination_database_Hevia_et_al_UPDATE.xlsx",
sheet = "field_level_data", startRow = 1)
datafield <- as_tibble(datafield)
datafield <- datafield %>%
rename(field_size = field.size,
fruits_per_plant = mean_fruits_per_plant,
richness_estimator_method = `richness_estimator_&#10;Method`)
datafield$management <- "conventional"
datafield$crop <- "Helianthus annuus"
datafield$study_id <- "hevi01"
chd = substr(datafield$latitude, 3, 3)[1]
chm = substr(datafield$latitude, 6, 6)[1]
chs = substr(datafield$latitude, 11, 11)[1]
cd = char2dms(datafield$latitude,chd=chd,chm=chm,chs=chs)
datafield$latitude <- as.numeric(cd)
chd = substr(datafield$longitude, 2, 2)[1]
chm = substr(datafield$longitude, 5, 5)[1]
chs = substr(datafield$longitude, 10, 10)[1]
cd = char2dms(datafield$longitude,chd = chd,chm = chm,chs = chs)
datafield$longitude <- as.numeric(cd)
field_level_data <- tibble(
study_id = datafield$study_id,
site_id = datafield$site_id,
crop = datafield$crop,
variety = datafield$variety,
management = datafield$management,
country = datafield$country,
latitude = datafield$latitude,
longitude = datafield$longitude,
X_UTM=NA,
Y_UTM=NA,
zone_UTM=NA,
sampling_start_month = datafield$sampling_start_month,
sampling_end_month = datafield$sampling_end_month,
sampling_year = datafield$sampling_year,
field_size = datafield$field_size,
yield=datafield$yield,
yield_units=datafield$yield_units,
yield2=datafield$yield2,
yield2_units=datafield$yield2_units,
yield_treatments_no_pollinators=datafield$yield_treatments_no_pollinators,
yield_treatments_pollen_supplement=NA,
yield_treatments_no_pollinators2=datafield$yield_treatments_no_pollinators2,
yield_treatments_pollen_supplement2=NA,
fruits_per_plant=datafield$fruits_per_plant,
fruit_weight= datafield$fruit_weight,
plant_density=datafield$plant_density,
seeds_per_fruit=datafield$seeds_per_fruit,
seeds_per_plant=datafield$seeds_per_plant,
seed_weight=datafield$seed_weight,
observed_pollinator_richness=NA,
other_pollinator_richness=NA,
other_richness_estimator_method=NA,
abundance = datafield$abundance,
ab_honeybee = datafield$ab_honeybee,
ab_bombus = datafield$ab_bombus,
ab_wildbees = datafield$ab_wildbees,
ab_syrphids = datafield$ab_syrphids,
ab_humbleflies= NA,
ab_other_flies= NA,
ab_beetles=NA,
ab_lepidoptera=NA,
ab_nonbee_hymenoptera=NA,
ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = "visit per 100 flowers and hour",
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
setwd("C:/Users/USUARIO/Desktop/OBservData/Datasets_storage")
write_csv(field_level_data, "field_level_data_hevi01.csv")
setwd(dir_ini)
library(tidyverse)
library(openxlsx)
library(sp) #Transforming latitude and longitude
dir_ini <- getwd()
options(digits=14)
datafield <- read.xlsx("Crop_pollination_database_Hevia_et_al_UPDATE.xlsx",
sheet = "field_level_data", startRow = 1)
datafield <- as_tibble(datafield)
datafield <- datafield %>%
rename(field_size = field.size,
fruits_per_plant = mean_fruits_per_plant,
richness_estimator_method = `richness_estimator_&#10;Method`)
datafield$management <- "conventional"
datafield$crop <- "Helianthus annuus"
datafield$study_id <- "hevi01"
chd = substr(datafield$latitude, 3, 3)[1]
chm = substr(datafield$latitude, 6, 6)[1]
chs = substr(datafield$latitude, 11, 11)[1]
cd = char2dms(datafield$latitude,chd=chd,chm=chm,chs=chs)
datafield$latitude <- as.numeric(cd)
chd = substr(datafield$longitude, 2, 2)[1]
chm = substr(datafield$longitude, 5, 5)[1]
chs = substr(datafield$longitude, 10, 10)[1]
cd = char2dms(datafield$longitude,chd = chd,chm = chm,chs = chs)
datafield$longitude <- as.numeric(cd)
field_level_data <- tibble(
study_id = datafield$study_id,
site_id = datafield$site_id,
crop = datafield$crop,
variety = datafield$variety,
management = datafield$management,
country = datafield$country,
latitude = datafield$latitude,
longitude = datafield$longitude,
X_UTM=NA,
Y_UTM=NA,
zone_UTM=NA,
sampling_start_month = datafield$sampling_start_month,
sampling_end_month = datafield$sampling_end_month,
sampling_year = datafield$sampling_year,
field_size = datafield$field_size,
yield=datafield$yield,
yield_units=datafield$yield_units,
yield2=datafield$yield2,
yield2_units=datafield$yield2_units,
yield_treatments_no_pollinators=datafield$yield_treatments_no_pollinators,
yield_treatments_pollen_supplement=NA,
yield_treatments_no_pollinators2=datafield$yield_treatments_no_pollinators2,
yield_treatments_pollen_supplement2=NA,
fruits_per_plant=datafield$fruits_per_plant,
fruit_weight= datafield$fruit_weight,
plant_density=datafield$plant_density,
seeds_per_fruit=datafield$seeds_per_fruit,
seeds_per_plant=datafield$seeds_per_plant,
seed_weight=datafield$seed_weight,
observed_pollinator_richness=NA,
other_pollinator_richness=NA,
other_richness_estimator_method=NA,
abundance = datafield$abundance,
ab_honeybee = datafield$ab_honeybee,
ab_bombus = datafield$ab_bombus,
ab_wildbees = datafield$ab_wildbees,
ab_syrphids = datafield$ab_syrphids,
ab_humbleflies= NA,
ab_other_flies= NA,
ab_beetles=NA,
ab_lepidoptera=NA,
ab_nonbee_hymenoptera=NA,
ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = "visit per 100 flowers and hour",
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
field_level_data <- tibble(
study_id = datafield$study_id,
site_id = datafield$site_id,
crop = datafield$crop,
variety = datafield$variety,
management = datafield$management,
country = datafield$country,
latitude = datafield$latitude,
longitude = datafield$longitude,
X_UTM=NA,
Y_UTM=NA,
zone_UTM=NA,
sampling_start_month = datafield$sampling_start_month,
sampling_end_month = datafield$sampling_end_month,
sampling_year = datafield$sampling_year,
field_size = datafield$field_size,
yield=datafield$yield,
yield_units=datafield$yield_units,
yield2=datafield$yield2,
yield2_units=datafield$yield2_units,
yield_treatments_no_pollinators=datafield$yield_treatments_no_pollinators,
yield_treatments_pollen_supplement=NA,
yield_treatments_no_pollinators2=datafield$yield_treatments_no_pollinators2,
yield_treatments_pollen_supplement2=NA,
fruits_per_plant=datafield$fruits_per_plant,
fruit_weight= datafield$fruit_weight,
plant_density=datafield$plant_density,
seeds_per_fruit=datafield$seeds_per_fruit,
seeds_per_plant=datafield$seeds_per_plant,
seed_weight=datafield$seed_weight,
observed_pollinator_richness=NA,
other_pollinator_richness=NA,
other_richness_estimator_method=NA,
abundance = datafield$abundance,
ab_honeybee = datafield$ab_honeybee,
ab_bombus = datafield$ab_bombus,
ab_wildbees = datafield$ab_wildbees,
ab_syrphids = datafield$ab_syrphids,
ab_humbleflies= NA,
ab_other_flies= NA,
ab_beetles=NA,
ab_lepidoptera=NA,
ab_nonbee_hymenoptera=NA,
ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = "visits per 100 flowers and hour",
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
setwd("C:/Users/USUARIO/Desktop/OBservData/Datasets_storage")
write_csv(field_level_data, "field_level_data_hevi01.csv")
setwd(dir_ini)
library(tidyverse)
library(openxlsx)
library(sp) #Transforming latitude and longitude
dir_ini <- getwd()
options(digits=14)
datafield <- read.xlsx("Crop_pollination_database_Hevia_et_al_UPDATE.xlsx",
sheet = "field_level_data", startRow = 1)
datafield <- as_tibble(datafield)
datafield <- datafield %>%
rename(field_size = field.size,
fruits_per_plant = mean_fruits_per_plant,
richness_estimator_method = `richness_estimator_&#10;Method`)
datafield$management <- "conventional"
datafield$crop <- "Helianthus annuus"
datafield$study_id <- "Violeta_Hevia_Helianthus_annuus_Spain_2017"
chd = substr(datafield$latitude, 3, 3)[1]
chm = substr(datafield$latitude, 6, 6)[1]
chs = substr(datafield$latitude, 11, 11)[1]
cd = char2dms(datafield$latitude,chd=chd,chm=chm,chs=chs)
datafield$latitude <- as.numeric(cd)
chd = substr(datafield$longitude, 2, 2)[1]
chm = substr(datafield$longitude, 5, 5)[1]
chs = substr(datafield$longitude, 10, 10)[1]
cd = char2dms(datafield$longitude,chd = chd,chm = chm,chs = chs)
datafield$longitude <- as.numeric(cd)
field_level_data <- tibble(
study_id = datafield$study_id,
site_id = datafield$site_id,
crop = datafield$crop,
variety = datafield$variety,
management = datafield$management,
country = datafield$country,
latitude = datafield$latitude,
longitude = datafield$longitude,
X_UTM=NA,
Y_UTM=NA,
zone_UTM=NA,
sampling_start_month = datafield$sampling_start_month,
sampling_end_month = datafield$sampling_end_month,
sampling_year = datafield$sampling_year,
field_size = datafield$field_size,
yield=datafield$yield,
yield_units=datafield$yield_units,
yield2=datafield$yield2,
yield2_units=datafield$yield2_units,
yield_treatments_no_pollinators=datafield$yield_treatments_no_pollinators,
yield_treatments_pollen_supplement=NA,
yield_treatments_no_pollinators2=datafield$yield_treatments_no_pollinators2,
yield_treatments_pollen_supplement2=NA,
fruits_per_plant=datafield$fruits_per_plant,
fruit_weight= datafield$fruit_weight,
plant_density=datafield$plant_density,
seeds_per_fruit=datafield$seeds_per_fruit,
seeds_per_plant=datafield$seeds_per_plant,
seed_weight=datafield$seed_weight,
observed_pollinator_richness=NA,
other_pollinator_richness=NA,
other_richness_estimator_method=NA,
abundance = datafield$abundance,
ab_honeybee = datafield$ab_honeybee,
ab_bombus = datafield$ab_bombus,
ab_wildbees = datafield$ab_wildbees,
ab_syrphids = datafield$ab_syrphids,
ab_humbleflies= NA,
ab_other_flies= NA,
ab_beetles=NA,
ab_lepidoptera=NA,
ab_nonbee_hymenoptera=NA,
ab_others = NA,
total_sampled_area = datafield$total_sampled_area,
total_sampled_time = datafield$total_sampled_time,
visitation_rate_units = "visits per 100 flowers and hour",
visitation_rate = datafield$visitation_rate,
visit_honeybee = datafield$visit_honeybee,
visit_bombus = datafield$visit_bombus,
visit_wildbees = datafield$visit_wildbees,
visit_syrphids = datafield$visit_syrphids,
visit_humbleflies = NA,
visit_other_flies = NA,
visit_beetles = NA,
visit_lepidoptera = NA,
visit_nonbee_hymenoptera = NA,
visit_others = NA,
Publication = datafield$Publication,
Credit = datafield$Credit,
Email_contact = datafield$email
)
setwd("C:/Users/USUARIO/Desktop/OBservData/Datasets_storage")
write_csv(field_level_data, "field_level_data_Violeta_Hevia_Helianthus_annuus_Spain_2017.csv")
setwd(dir_ini)
View(datafield)
