PlotCov_Fn: Plots impact of covariate on surfaces

In aaronmberger/Geo_dGLMM_habitat: Package for conducting spatial estimation of CPUE index standardization (fisheries, ecology, etc).

Usage

PlotCov_Fn(Report, NN_Extrap, X_xj, FileName, ControlList = list(Width = 5 * 3, Height = 2 * 3, Res = 200, Units = "in"))

Arguments

Report
NN_Extrap
X_xj
FileName
ControlList

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (Report, NN_Extrap, X_xj, FileName, ControlList = list(Width = 5 * 
    3, Height = 2 * 3, Res = 200, Units = "in")) 
{
    D_it = Report$D_xt[NN_Extrap$nn.idx, ]
    R1_it = Report$R1_xt[NN_Extrap$nn.idx, ]
    R2_it = Report$R2_xt[NN_Extrap$nn.idx, ]
    for (CovI in 1:ncol(X_xj)) {
        png(file = paste0(FileName, colnames(X_xj)[CovI], ".png"), 
            width = ControlList$Width, height = ControlList$Height, 
            res = ControlList$Res, units = ControlList$Units)
        par(mfrow = c(2, 6), oma = c(0, 0, 2, 0), mar = c(2, 
            2, 0, 0), mgp = c(1.5, 0.5, 0), tck = -0.02)
        for (t in 1:length(Year_Set)) {
            plot(x = X_xj[NN_Extrap$nn.idx, CovI], y = R2_it[, 
                t], type = "p", col = "red", ylim = c(0, max(R2_it)))
            points(x = X_xj[NN_Extrap$nn.idx, CovI], y = R1_it[, 
                t] * max(R2_it), type = "p", col = "blue")
            points(x = X_xj[NN_Extrap$nn.idx, CovI], y = D_it[, 
                t], type = "p", col = "black")
        }
        mtext(side = 3, outer = TRUE, colnames(X_xj)[CovI])
        dev.off()
    }
  }

aaronmberger/Geo_dGLMM_habitat documentation built on May 10, 2019, 3:20 a.m.