The function implements the Interval Testing Procedure for testing the center of symmetry of a functional population evaluated on a uniform grid. Data are represented by means of the Fourier expansion and the significance of each basis coefficient is tested with an interval-wise control of the Family Wise Error Rate.

1 | ```
ITP1fourier(data, mu = 0, maxfrequency=floor(dim(data)[2]/2), B = 10000)
``` |

`data` |
Pointwise evaluations of the functional data set on a uniform grid. |

`mu` |
The center of symmetry under the null hypothesis: either a constant (in this case, a constant function is used) or a |

`maxfrequency` |
The maximum frequency to be used in the Fourier basis expansion of data. The default is |

`B` |
The number of iterations of the MC algorithm to evaluate the p-values of the permutation tests. The defualt is |

`ITP1fourier`

returns an object of `class`

"`ITP1`

".

An object of class "`ITP1`

" is a list containing at least the following components:

`basis` |
String vector indicating the basis used for the first phase of the algorithm. In this case equal to |

`test` |
String vector indicating the type of test performed. In this case equal to |

`mu` |
Center of symmetry under the null hypothesis (as entered by the user). |

`coeff` |
Matrix of dimensions |

`pval` |
Uncorrected p-values for each frequency. |

`pval.matrix` |
Matrix of dimensions |

`corrected.pval` |
Corrected p-values for each frequency. |

`labels` |
Labels indicating the population membership of each data (in this case always equal to |

`data.eval` |
Evaluation on a fine uniform grid of the functional data obtained through the basis expansion. |

`heatmap.matrix` |
Heatmap matrix of p-values (used only for plots). |

Alessia Pini, Simone Vantini

A. Pini and S. Vantini (2013). The Interval Testing Procedure: Inference for Functional Data Controlling the Family Wise Error Rate on Intervals. *MOX-report 13/2013*, Politecnico di Milano.

See also `ITP1bspline`

, `ITP2bspline`

, `ITP2fourier`

, `ITP2pafourier`

, and `ITPimage`

.

1 2 3 4 5 6 7 8 9 10 11 12 | ```
# Importing the NASA temperatures data set
data(NASAtemp)
# Performing the ITP
ITP.result <- ITP1fourier(NASAtemp$milan,maxfrequency=20,B=1000)
# Plotting the results of the ITP
plot(ITP.result,main='NASA data',xrange=c(1,365),xlab='Day')
# Plotting the p-value heatmap
ITPimage(ITP.result,abscissa.range=c(1,365))
# Selecting the significant coefficients
which(ITP.result$corrected.pval < 0.05)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.