Actigraphy (version 1.4.0)

cont_flm_plot: Plot Functional Linear Model Analysis Results of a Continuous Type

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

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.

Examples

# NOT RUN {
	data(act_29pt)
	data(clinic_29pt_ahi)
	
	data <- act_29pt
	ahi <- clinic_29pt_ahi
	
	matchid  <- fda.matchid(data, 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)
# }