FDboost (version 0.2-0)

hmatrix: A S3 class for univariate functional data on a common grid

Description

The hmatrix class represents data for a functional historical effect. The class is basically a matrix containing the time and the id for the observations of the functional response. The functional covariate is contained as attribute.

Usage

hmatrix(time, id, x, argvals = 1:ncol(x), timeLab = "t",
  idLab = "wideIndex", xLab = "x", argvalsLab = "s")

Arguments

time

set of argument values of the response in long format, i.e. at which t the response curve is observed

id

specify to which curve the point belongs to, id from 1, 2, ..., n.

x

matrix of functional covariate, each trajectory is in one row

argvals

set of argument values, i.e., the common gird at which the functional covariate is observed, by default 1:ncol(x)

timeLab

name of the time axis, by default t

idLab

name of the id variable, by default wideIndex

xLab

name of the functional variable, by default NULL

argvalsLab

name of the argument for the covariate by default s

Details

In the hmatrix class the id has to run from i=1, 2, ..., n including all integers from 1 to n. The rows of the functional covariate x correspond to those observations.

See Also

getTime.hmatrix to extract attributes, and ?"[.hmatrix" for the extract method.

Examples

## Example for a hmatrix object
t1 <- rep((1:5)/2, each=3)
id1 <- rep(1:3, 5)
x1 <- matrix(1:15, ncol=5)
s1 <- (1:5)/2
myhmatrix <- hmatrix(time=t1, id=id1, x=x1, argvals=s1, timeLab="t1", argvalsLab="s1", xLab="test")

# extract with [ keeps attributes
# select observations of subjects 2 and 3
myhmatrixSub <- myhmatrix[id1 %in% c(2,3),]
str(myhmatrixSub)
getX(myhmatrixSub)
getX(myhmatrix)

# get time
myhmatrix[,1] # as column matrix as drop=FALSE
getTime(myhmatrix) # as vector

# get id
myhmatrix[,2] # as column matrix as drop=FALSE
getId(myhmatrix) # as vector

# subset hmatrix on the basis of an index, which is defined on the curve level
reweightData(data = list(hmat = myhmatrix), vars = "hmat", index = c(1,1,2))
# this keeps only the unique x values in attr(,'x') but multiplies the corresponding
# ids and times in the time id matrix
# for bhistx baselearner, there may be an additional id variable for the tensor product
newdat <- reweightData(data = list(hmat = myhmatrix,
  repIDx = rep(1:nrow(attr(myhmatrix,'x')), length(attr(myhmatrix,"argvals")))),
  vars = "hmat", index = c(1,1,2), idvars="repIDx")
length(newdat$repIDx)

## use hmatrix within a data.frame
mydat <- data.frame(I(myhmatrix), z=rnorm(3)[id1])
str(mydat)
str(mydat[id1 %in% c(2,3),])
str(myhmatrix[id1 %in% c(2,3),])