plot.msa: Plot function for moving subset analysis

In msaFACE: Moving Subset Analysis FACE

Description

Create main or accompanying plot(s) of the moving subset analysis

Usage

## S3 method for class 'MSA_coef'
plot(x, i_var = 1:length(x), main_plot = TRUE, i_acc = 1, 
     axis.param = list(), label.param = list(), ...)

Arguments

x	Object of class MSA_coef. Output of moving_subset_analysis.
i_var	Column number or name of the forcing variable used to define the subset which will be plotted. Default is to plot all available variables.
main_plot	Flag to determine whether to plot the main plot(s) (default) or the accompanying plot(s)
i_acc	If main_plot = FALSE, column number(s) of the experiment-support variable(s) to be plotted as accompanying variable(s)
axis.param	List containing the axis parameters (e.g. lty, col, tck). By default for main plots, first entry represents x-axis (forcing variable), second entry represents first y-axis (regression result), and third entry the second y-axis (dependent variable). See example section.
label.param	List containing the label parameters (e.g. line, cex). For order of the list refer to axis.param. See example section.
...	Parameters passed to generic plot function

Author(s)

Wolfgang A. Obermeier, Lukas W. Lehnert, Joerg Bendix

See Also

moving_subset_analysis, GiFACE, summary.MSA_coef

Examples

data(GiFACE)

### Perform Moving Subset Analysis
msa_data <- moving_subset_analysis(GiFACE$data, "CO2A_Mean", "BYT", 30, group = GiFACE$year)

## Plot all final plots
plot(msa_data)

## Not run: 
## Example to produce similar plots as in article
pdf("Plot_NCC_1.pdf", width = 15, height = 20)
par(mfcol = c(9,4), mar = c(4,4,1,4), las = 0)
for(i in 1:4){
  plot(msa_data, i_var = i, ylim = list(c(0,2.5), c(200,350)), 
       legend = i == 1)
  plot(msa_data, i_var = i, main_plot = F, i_acc = 1:8)
}
dev.off()

## End(Not run)

## Not run: 
## Example to adjust axis and label parameters
par(mar = c(3.2,3.5,1,3), mfrow = c(2,1))
## Plot default for comparison
plot(msa_data, i_var = 2)     

## Change default axis (all parameters are passed to internal function "axis")
## Define x-axis specifications
xaxt_spec  <- list(tck = -0.02, padj = -0.3, cex.axis = 1.5, lwd = 2)
## Define 1st y-axis specifications 
yaxt1_spec <- list(at = c(0,1,2), tck = -0.02, padj = 0.5, cex.axis = 1.5, lwd = 2)
## Define 2nd y-axis specifications 
yaxt2_spec <- list(tck = -0.02, padj = -0.3, col = "gray50", col.axis = "gray50",
                   cex.axis = 1.5, lwd = 2)

## Change default labeling (all parameters are passed to internal function "mtext")
## Define x-axis label specifications 
xlab_spec <- list(side = 1, line = 2, cex = 1.5)
## Define 1st y-axis label specifications
ylab1_spec <- list(side = 2, line = 2, cex = 1.5)
## Define 2nd y-axis label specifications 
ylab2_spec <- list(side = 4, line = 2, cex = 1.5, col = "gray50")
                   
                   
## Define variable to be plotted, y-axis limits
plot(msa_data, i_var = 2, ylim = list(c(0,2.5), c(200,350)),                                     
     axis.param = list(xaxt_spec, yaxt1_spec, yaxt2_spec),
     label.param = list(xlab_spec, ylab1_spec, ylab2_spec),
     ## Define line width for horizontal lines (1) and the dependent variable (2)
     lwd = c(1,2),       
     ## Define point symbol
     pch = 20)                                                                                   
dev.off()

## End(Not run)

Example output

png 
  2 
null device 
          1 

msaFACE documentation built on May 2, 2019, 5:07 a.m.