diff --git a/01_Mesobromion.R b/01_Mesobromion.R index 816b33edbcf1604dfb1732938142b3d27c589d05..98aa7d63326229340d3d0fa544ff650623a47e19 100644 --- a/01_Mesobromion.R +++ b/01_Mesobromion.R @@ -902,14 +902,17 @@ ggsave("_data/Mesobromion/PCoA_fuzzy.png", width=6, height=4, dpi=300, last_plot library(rgdal) library(sp) library(sf) +coords <- read_delim(file = "_data/Mesobromion/Mes2_somecoordinates.csv", delim = ";") + env <- read_delim(file = "_data/Mesobromion/env.10perc.txt", delim = "\t") env0 <- read_delim("_data/Mesobromion/GVRD_MES2_site.csv", delim = ",") -header <- "/data/sPlot/users/Francesco/Project_11/Germany/_data/tvhabita.dbf" +#header <- "/data/sPlot/users/Francesco/Project_11/Germany/_data/tvhabita.dbf" env.all <- env0 %>% - left_join(foreign::read.dbf(header) %>% - as.data.frame() %>% - dplyr::select(RELEVE_NR, LAT, LON), + left_join(coords %>% + dplyr::select(RELEVE_NR=MES2_releve_nr., + LON=DD_X_, + LAT=DD_Y_), by="RELEVE_NR") %>% filter(!is.na(LAT)) %>% filter(!(LAT==0 | LON==0))