HSIplotter: Plots habitat suitability index curves

In ecorest: Conducts Analyses Informing Ecosystem Restoration Decisions

Plots habitat suitability index curves

Description

HSIplotter plots all suitability curves.

Usage

HSIplotter(SI, figure.name)

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".

Value

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

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"))

ecorest documentation built on March 31, 2026, 1:07 a.m.