npregfast (version 1.5.2)

criticaldiff: Differences between the critical points for two factor's levels

Description

Differences between the estimation of critical for two factor's levels.

Usage

criticaldiff(model, level1 = NULL, level2 = NULL, der = NULL)

Value

An object is returned with the following elements:

critical.diff

a table with a couple of factor's level where it is used to calculate the differences between the critical points, and their 95% confidence interval (for the estimation, first and second derivative).

Arguments

model

Parametric or nonparametric regression model obtained by frfast function.

level1

First factor's level at which to perform the differences between critical points.

level2

Second factor's level at which to perform the differences between critical points.

der

Number which determines any inference process. By default der is NULL. If this term is 0, the calculate of the differences for the critical point is for the estimate. If it is 1 or 2, it is designed for the first or second derivative, respectively.

Author

Marta Sestelo, Nora M. Villanueva and Javier Roca-Pardinas.

Details

Differences are calculated by subtracting a factor relative to another (\(level2 - level1\)). By default level2 and level1 are NULL, so the differences calculated are for all possible combinations between two factors. Additionally, it is obtained the 95% confidence interval for this difference which let us to make inference about them.

References

Sestelo, M. (2013). Development and computational implementation of estimation and inference methods in flexible regression models. Applications in Biology, Engineering and Environment. PhD Thesis, Department of Statistics and O.R. University of Vigo.

Sestelo, M., Villanueva, N.M., Meira-Machado, L., Roca-Pardinas, J. (2017). npregfast: An R Package for Nonparametric Estimation and Inference in Life Sciences. Journal of Statistical Software, 82(12), 1-27.

Examples

library(npregfast)
data(barnacle)
fit2 <- frfast(DW ~ RC : F, data = barnacle, seed = 130853, nboot = 100) # with interactions
criticaldiff(fit2)
criticaldiff(fit2, der = 1)
criticaldiff(fit2, der = 1, level1 = "lens", level2 = "barca")