rbenttest: test the existence of change point in the bent line...

In Rbent: Robust Bent Line Regression

Description

This function use Wilcoxon score functions for calculating the test statistics and p-value by wild bootstrap.

Usage

rbenttest(y, z, x, NB = 1000, myseed = 1)

Arguments

y	A vector of response
z	A vector of covariates
x	A numeric variable with change point
NB	resampling times
myseed	set seed

Value

A list with the elements

Tn	The statistic based on original data.
Tn.NB	The statistics by wild bootstrap.
p.value	The p-value by wild bootstrap.

Author(s)

Feipeng Zhang

Examples

# for the example of  MRS data
data(data_mrs)
x <- log(data_mrs$mass)
y <- log(data_mrs$speed)
z <- data_mrs$hopper
p.value <- rbenttest(y, cbind(1, z), x, NB = 50)$p.value

# for the example of bedload transport data
data(data_transport)
x <- data_transport$x
y <- data_transport$y
p.value <- rbenttest(y, 1, x, NB = 50)$p.value

Rbent documentation built on May 2, 2019, 8:25 a.m.