Plot Functional Linear Model Analysis Results of a Continuous Type

Share:

Description

This function produces two plots: An effect of a continuous covariate on activity values by time and the F-test for the effect of the continuous covariate.

Usage

1
2
cont_flm_plot(smoothdata, matchresults, flmresults, xlim, ylim, 
	ftest, nperm, lb, xat, legendx, legendy, L, xlab="Time", ylab="Activity")

Arguments

smoothdata

List output from the fda.smoothdata function.

matchresults

List output from the matchid function.

flmresults

List output from the flm_cate function.

xlim

X axis limits for activity plot.

ylim

Y axis limits for activity plot.

ftest

A logic value indicating whether to implement F test or not. F test will be implement if ftest is TRUE.

nperm

The number of permutations for the F-test.

lb

X-axis labels.

xat

X axis label positions.

legendx

X axis position of the left edge of the legend box.

legendy

Y axis position of the upper edge of the legend box.

L

The length of the time points.

xlab

The label for the x-axis.

ylab

The label for the y-axis.

Value

One plot of the estimated group means and a possible second plot of the F-test results.

Author(s)

William D. Shannon, Tao Li, Hong Xian, Jia Wang, Elena Deych, Carlos Gonzalez

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
	data(act_29pt)
	data(clinic_29pt_ahi)
	
	colnames(act_29pt) <- sub("X", "", colnames(act_29pt)) 
	data <- as.matrix(act_29pt[,-1])
	
	matchid  <- fda.matchid(data, clinic_29pt_ahi, "contin")
	FDcont <- fda.smoothdata(matchid)
	
	L <- nrow(data)
	lb <- c("Midnight", "6AM", "Noon", "6PM", "Midnight") 
	xat <- c(0, L/4, L/2, 3*L/4, L)
	
	geftFDcont <- flm_cate(FDcont)
	predy <- as.vector(geftFDcont$freg$yhatfdobj$y)
	
	xlim <- c(0, L) 
	ylim <-  c(min(predy), max(predy) + 100)
	
	legendx <- 0
	legendy <- max(predy) - 100
	
	cont.flm.results <- cont_flm_plot(FDcont, matchid, geftFDcont, xlim, 
		ylim, TRUE, 10, lb, xat, legendx, legendy, L)