ecorest (version 2.0.2)

HSIplotter: Plots habitat suitability index curves

Description

HSIplotter plots all suitability curves.

Usage

HSIplotter(SI, figure.name)

Value

A multi-panel *.jpeg figure showing all suitability curves.

Arguments

SI

list, matrix, or dataframe of suitability curves ordered as parameter breakpoints and associated suitability indices for each parameter with appropriate column names. Models containing both categorical and continuous parameters must be entered as a dataframe.

figure.name

output figure file name structured as "filename.jpeg".

References

US Fish and Wildlife Service. (1980). Habitat as a basis for environmental assessment. Ecological Services Manual, 101.

US Fish and Wildlife Service. (1980). Habitat Evaluation Procedures (HEP). Ecological Services Manual, 102.

US Fish and Wildlife Service. (1981). Standards for the Development of Habitat Suitability Index Models. Ecological Services Manual, 103.

Examples

#Build and define a matrix of the Barred Owl suitability curves
#Allen A.W. 1982. Habitat Suitability Index Models: Barred owl. FWS/OBS 82/10.143.
#U.S. Fish and Wildlife Service. https://pubs.er.usgs.gov/publication/fwsobs82_10_143.
var1 <- cbind(c(0,2,4,NA), c(0.1,1,1,NA)) #Number of trees > 51cm diameter per 0.4 ha plot
var2 <- cbind(c(0,5,20,NA), c(0,0,1,NA)) #Mean diameter of overstory trees
var3 <- cbind(c(0,20,60,100), c(0,0,1,1)) #Percent canopy cover of overstory trees
barredowl <- cbind(var1, var2, var3)
colnames(barredowl)<- c("tree.num", "tree.num.SIV",
  "avg.dbh.in", "avg.dbh.SIV", "can.cov", "can.cov.SIV")

#Create suitability curve summary plot
HSIplotter(barredowl, tempfile("BarredOwl",fileext=".jpeg"))


#Build and define a matrix of the alewifeSAEL curves
#Pardue, GB. 1983. Habitat suitability index models: alewife and blueback herring. 
#U.S. Dept. Int. Fish Wildl. Serv. FWS/OBS-82/10.58. 22pp.
var1 = data.frame(subs.type.class = c("a", "b", "c", NA, NA, NA), 
                  subs.type.class.SIV = c(1, 0.5, 0.1, NA, NA, NA))
var2 = data.frame("avg.daily.wtr.temp.spwn.C" = c(5, 10, 15, 20, 27, 30), 
                "avg.daily.wtr.temp.spwn.C.SIV" = c(0, 0, 1, 1, 0, 0))
alewifeSAEL = data.frame(var1, var2)

#Create suitability curve summary plot
HSIplotter(alewifeSAEL, tempfile("AlewifeSAEL",fileext=".jpeg"))

#Use HSIplotter with inputs from HSImodels
HSIplotter(HSImodels$barredowl, tempfile("BarredOwl",fileext=".jpeg"))