basis: Get the basis functions from a dfrr-object

In dfrr: Dichotomized Functional Response Regression

Get the basis functions from a dfrr-object

Description

Returns the basis functions employed in fitting a dfrr-object.

Usage

basis(object)

Arguments

object

a fitted dfrr-object obtained from invoking the function dfrr.

Value

a basis object used in fitting the functional parameters. The basis object is the one created by the functions create.*.basis of the 'fda' package.

Examples

set.seed(2000)
N<-50;M<-24

X<-rnorm(N,mean=0)
time<-seq(0,1,length.out=M)
Y<-simulate_simple_dfrr(beta0=function(t){cos(pi*t+pi)},
                        beta1=function(t){2*t},
                        X=X,time=time)

#The argument T_E indicates the number of EM algorithm.
#T_E is set to 1 for the demonstration purpose only.
#Remove this argument for the purpose of converging the EM algorithm.
dfrr_fit<-dfrr(Y~X,yind=time,T_E=1)
coefs<-coef(dfrr_fit,return.fourier.coefs=TRUE)

basis<-basis(dfrr_fit)
evaluated_coefs<-coefs%*%t(fda::eval.basis(time,basis))

#Plotting the regression coefficients
oldpar<-par(mfrow=c(1,2))

plot(time,evaluated_coefs[1,],'l',main="Intercept")
plot(time,evaluated_coefs[2,],'l',main="X")

par(oldpar)

dfrr documentation built on May 31, 2023, 5:32 p.m.