Learn R Programming

⚠️There's a newer version (2.5) of this package.Take me there.

itsadug (version 2.2)

Interpreting Time Series and Autocorrelated Data Using GAMMs

Description

GAMM (Generalized Additive Mixed Modeling; Lin & Zhang, 1999) as implemented in the R package 'mgcv' (Wood, S.N., 2006; 2011) is a nonlinear regression analysis which is particularly useful for time course data such as EEG, pupil dilation, gaze data (eye tracking), and articulography recordings, but also for behavioral data such as reaction times and response data. As time course measures are sensitive to autocorrelation problems, GAMMs implements methods to reduce the autocorrelation problems. This package includes functions for the evaluation of GAMM models (e.g., model comparisons, determining regions of significance, inspection of autocorrelational structure in residuals) and interpreting of GAMMs (e.g., visualization of complex interactions, and contrasts).

Copy Link

Version

Install

install.packages('itsadug')

Monthly Downloads

3,643

Version

2.2

License

GPL (>= 2)

Maintainer

Jacolien van Rij

Last Published

June 13th, 2016

Functions in itsadug (2.2)

get_difference

Get model predictions for differences between conditions.
fvisgam

Visualization of nonlinear interactions, summed effects.
get_coefs

Get coefficients for the parametric terms (intercepts and random slopes).
get_predictions

Get model predictions for specific conditions.
get_modelterm

Get estimated for selected model terms.
gamtabs

Convert model summary into Latex/HTML table for knitr/R Markdown reports.
fadeRug

Fade out the areas in a surface without data.
get_fitted

Get model all fitted values.
check_resid

Inspect residuals of regression models.
diff_terms

Compare the formulas of two models and return the difference(s).
acf_plot

Generate an ACF plot of an aggregated time series.
derive_timeseries

Derive the time series used in the AR1 model.
compareML

Function for comparing two GAMM models.
get_pca_predictions

Return PCA predictions.
inspect_random

Inspection and interpretation of random factor smooths.
find_difference

Return the regions in which the smooth is significantly different from zero.
missing_est

Return indices of data that were not fitted by the model.
info

Information on how to cite this package
infoMessages

Turn on or off information messages.
plot_data

Visualization of the model fit for time series data.
plot_diff2

Plot difference surface based on model predictions.
get_random

Get coefficients for the random intercepts and random slopes.
itsadug

Interpreting Time Series, Autocorrelated Data Using GAMMs (itsadug)
plot_diff

Plot difference curve based on model predictions.
plot_modelfit

Visualization of the model fit for time series data.
plot_topo

Visualization of EEG topo maps.
rug_model

Add rug to plot, based on model.
report_stats

Returns a description of the statistics of the smooth terms for reporting.
simdat

Simulated time series data.
plot_pca_surface

Visualization of the effect predictors in nonlinear interactions with principled components.
print_summary

plot_smooth

Visualization of smooths.
pvisgam

Visualization of partial nonlinear interactions.
plot_parametric

Visualization of group estimates.
resid_gam

Extract model residuals and remove the autocorrelation accounted for.
summary_data

Print a descriptive summary of a data frame.
wald_gam

Function for post-hoc comparison of the contrasts in a single GAMM model.
timeBins

Label timestamps as timebins of a given binsize.
start_value_rho

Extract the Lag 1 value from the ACF of the residuals of a gam, bam, lm, lmer model, ...
start_event

Determine the starting point for each time series.
diagnostics

Visualization of the model fit for time series data.
acf_resid

Generate an ACF plot of model residuals. Works for lm, lmer, gam, bam, ....
convertNonAlphanumeric

Prepare string for regular expressions (backslash for all non-letter and non-digit characters)
acf_n_plots

Generate N ACF plots of individual or aggregated time series.
eeg

Raw EEG data, single trial, 50Hz.