logo Use CA10RAM to get 10%* Discount.
Order Nowlogo
(5/5)

Choose two species that you think interact in some way, either as competitors, predator-prey

INSTRUCTIONS TO CANDIDATES
ANSWER ALL QUESTIONS

TASK:

Choose two species that you think interact in some way, either as competitors, predator-prey, 

mutualists or host-parasite. You should choose species that have locality data available on GBIF. If 

your species don’t have data, pick another pair that do. Then complete the following tasks. For each 

task, provide a brief legend (no more than 150 words) summarising the analysis used to generate each the plot(s) and describing what the results show.

1. Run linear models to predict the present-day distribution of species 1 using climate variables and 

use them to present a map of its current distribution. Which set of climatic variables best explain the 

current distribution of the species?

2. Run linear models to predict the present-day distribution of species 2 using climate variables and 

use them to present maps of its current distribution. Which set of climatic variables best explain the 

current distribution of the species?

3. Plot the overlap in distribution of the two species. Devise and calculate a metric for the degree of 

overlap between their ranges, explaining how you calculated it.

4. Use a linear model to test whether the distribution of species 1 at the present time depends on the distribution of species 2, while also taking account of the effects of climatic variables. 

5. Predict the future distribution of each species using CMIP6 data for future climate and predict how  the degree of overlap in ranges change will change in that time. Do you expect the two species to co-occur more often or less often than at the present?

 

MY CODE SO FAR (for task 1, 2 and beginning of 3)

 

#choose two interacting species

 

#Pyrrhura molinae - parakeet

#fig tree - Ficus calyptroceras

 

#run the packages needed

library('dismo')

library('rworldmap')

library('sf')

library('geodata')

library('ggplot2')

library('rmarkdown')

 

wrld_simpl<-getMap(resolution = "coarse")

plot(wrld_simpl)

 

##TASK 1

#linear model to predict present day distribution of species 1 

#using climate variables

#use them to present a map of its current distribution

 

bird <- gbif("pyrrhura", "molinae", geo=TRUE)

bird.coords<-cbind(bird$lon,bird$lat)

bird.coords<-na.omit(bird.coords)

bird.coords<-data.frame(bird.coords)

colnames(bird.coords)<-c("lon","lat")

 

plot(wrld_simpl, xlim=range(bird.coords$lon), ylim=range(bird.coords$lat), axes=TRUE, col="light blue")

points(bird.coords, col='orange', cex=0.75)

 

#can see the distribution for the species is mainly in South America

#we can zoom into South America to look at its distribution more closely 

 

trim.coords<-function (x,latmin,latmax,lonmin,lonmax) {

  if (sum(x$lon < lonmin)>0) {

    tmp<-which(x$lon < lonmin)

    x<-x[-tmp,]}

  if (sum(x$lon > lonmax)>0) {

    tmp<-which(x$lon > lonmax)

    x<-x[-tmp,]}

  if (sum(x$lat < latmin)>0) {

    tmp<-which(x$lat < latmin)

    x<-x[-tmp,]}

  if (sum(x$lat > latmax)>0) {

    tmp<-which(x$lat > latmax)

    x<-x[-tmp,]}

  return(x) }

 

bird.coords.trim<-trim.coords(bird.coords,latmin=-50,latmax=25,lonmin=-75,lonmax=-20)

plot(wrld_simpl, xlim=range(bird.coords.trim$lon), ylim=range(bird.coords.trim$lat), axes=TRUE, col="light blue")

points(bird.coords.trim, col='orange', cex=0.75)

 

#extracting climatic variables 

bio.values <- extract(bio.data, bird.coords)[,-1]

 

bio.data<-worldclim_global(var="bio",res=10,path=getwd())

names(bio.data)<-paste0("bio",1:19)

species.data<-cbind(bird.coords,bio.values)

species.data<-na.omit(species.data)

write.csv(species.data,file="species.data.csv")

#extract mean values, min and max

rbind(mean=colMeans(species.data),

      min=apply(species.data, 2, min),

      max=apply(species.data, 2, max))

 

#use these to model current distribution

#make pseudo-absence points

 

#create raster

ext <- extent(wrld_simpl)

xy <- abs(apply(as.matrix(bbox(ext)), 1, diff))

n <- 5

r <- raster(ext, ncol=xy[1]*n, nrow=xy[2]*n)

mask <-rasterize(wrld_simpl, r)

 

#setting a plotting window with least grey space 

plot(mask, axes = TRUE, col = "grey",

     xlim = c(-80, -30),

     ylim = c(-40,-5)); box()

points(bird.coords,col = "orange", pch = 20, cex = 1)

 

 

#set an area of extent and create background points 

e <- extent(-80, -30, -40,-5)

plot(e, add=TRUE, col='black')

bg <- randomPoints(mask, 500, ext=e,extf=1)

points(bg,col="black", pch=20)

colnames(bg)<-c("lon","lat")

 

(5/5)
Attachments:

Related Questions

. Introgramming & Unix Fall 2018, CRN 44882, Oakland University Homework Assignment 6 - Using Arrays and Functions in C

DescriptionIn this final assignment, the students will demonstrate their ability to apply two ma

. The standard path finding involves finding the (shortest) path from an origin to a destination, typically on a map. This is an

Path finding involves finding a path from A to B. Typically we want the path to have certain properties,such as being the shortest or to avoid going t

. Develop a program to emulate a purchase transaction at a retail store. This program will have two classes, a LineItem class and a Transaction class. The LineItem class will represent an individual

Develop a program to emulate a purchase transaction at a retail store. Thisprogram will have two classes, a LineItem class and a Transaction class. Th

. SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of Sea Ports. Here are the classes and their instance variables we wish to define:

1 Project 1 Introduction - the SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of

. Project 2 Introduction - the SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of Sea Ports. Here are the classes and their instance variables we wish to define:

1 Project 2 Introduction - the SeaPort Project series For this set of projects for the course, we wish to simulate some of the aspects of a number of

Ask This Question To Be Solved By Our ExpertsGet A+ Grade Solution Guaranteed

expert
Um e HaniScience

612 Answers

Hire Me
expert
Muhammad Ali HaiderFinance

992 Answers

Hire Me
expert
Husnain SaeedComputer science

528 Answers

Hire Me
expert
Atharva PatilComputer science

633 Answers

Hire Me
April
January
February
March
April
May
June
July
August
September
October
November
December
2025
1950
1951
1952
1953
1954
1955
1956
1957
1958
1959
1960
1961
1962
1963
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017
2018
2019
2020
2021
2022
2023
2024
2025
2026
2027
2028
2029
2030
2031
2032
2033
2034
2035
2036
2037
2038
2039
2040
2041
2042
2043
2044
2045
2046
2047
2048
2049
2050
SunMonTueWedThuFriSat
30
31
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
1
2
3
00:00
00:30
01:00
01:30
02:00
02:30
03:00
03:30
04:00
04:30
05:00
05:30
06:00
06:30
07:00
07:30
08:00
08:30
09:00
09:30
10:00
10:30
11:00
11:30
12:00
12:30
13:00
13:30
14:00
14:30
15:00
15:30
16:00
16:30
17:00
17:30
18:00
18:30
19:00
19:30
20:00
20:30
21:00
21:30
22:00
22:30
23:00
23:30