fda.usc (version 2.0.1)

fanova.onefactor: One--way anova model for functional data

Description

One--way anova model for k independent samples of functional data. The function contrasts the null hypothesis of equality of mean functions of functional data based on the an asymptotic version of the anova F--test. $$H_0:\, m_1=\ldots=m_k$$

Usage

fanova.onefactor(
  object,
  group,
  nboot = 100,
  plot = FALSE,
  verbose = FALSE,
  ...
)

Arguments

object

functional response data. fdata class object with n curves.

group

a factor specifying the class for each curve.

nboot

number of bootstrap samples.

plot

if TRUE, plot the mean of each factor level and the results of test.

verbose

if TRUE, print intermediate results.

further arguments passed to or from other methods.

Value

Returns:

  • p-value probability of rejecting the null hypothesis H0 at a significance level

  • stat statistic value of test.

  • wm statistic values of bootstrap resamples.

Details

The function returns the p--value of test using one--way anova model over nboot runs.

References

Cuevas, A., Febrero, M., & Fraiman, R. (2004). An anova test for functional data. Computational statistics & data analysis, 47(1), 111-122.

See Also

See Also as: fanova.RPm

Examples

Run this code
# NOT RUN {
data(MCO)
grupo<-MCO$classintact
datos<-MCO$intact
res=fanova.onefactor(datos,grupo,nboot=50,plot=TRUE)
grupo <- MCO$classpermea
datos <- MCO$permea
res=fanova.onefactor(datos,grupo,nboot=50,plot=TRUE)
# }

Run the code above in your browser using DataLab