Description Usage Arguments Details Value See Also Examples
Plots the results of pda.fd, allows the user to group coefficient functions by variable, equation, derivative or combination of them.
1 2 |
x |
an object of class |
whichdim |
which dimension to use as grouping variables
|
npts |
number of points to use for plotting. |
... |
other arguments for 'plot'. |
Produces one plot for each coefficient function in a principle differential analysis.
invisible(NULL)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 | # A pda analysis of the handwriting data
# reduce the size to reduce the compute time for the example
ni <- 281
indx <- seq(1, 1401, length=ni)
fdaarray <- handwrit[indx,,]
fdatime <- seq(0, 2.3, len=ni)
# basis for coordinates
fdarange <- c(0, 2.3)
breaks <- seq(0,2.3,length.out=116)
norder <- 6
fdabasis <- create.bspline.basis(fdarange,norder=norder,breaks=breaks)
# parameter object for coordinates
fdafd0 <- fd(matrix(0,fdabasis$nbasis,1), fdabasis)
fdaPar <- fdPar(fdafd0,int2Lfd(4),1e-8)
# coordinate functions and a list tontaining them
Xfd <- smooth.basis(fdatime, fdaarray[,,1], fdaPar)$fd
Yfd <- smooth.basis(fdatime, fdaarray[,,2], fdaPar)$fd
xfdlist <- list(Xfd, Yfd)
# basis and parameter object for weight functions
fdabasis2 <- create.bspline.basis(fdarange,norder=norder,nbasis=31)
fdafd0 <- fd(matrix(0,fdabasis2$nbasis,1), fdabasis2)
pdaPar <- fdPar(fdafd0,1,1e-8)
pdaParlist <- list(pdaPar, pdaPar)
bwtlist <- list( list(pdaParlist,pdaParlist), list(pdaParlist,pdaParlist) )
# do the second order pda
pdaList <- pda.fd(xfdlist, bwtlist)
# plot the results
plot(pdaList,whichdim=1)
plot(pdaList,whichdim=2)
plot(pdaList,whichdim=3)
plot(pdaList,whichdim=c(1,2))
plot(pdaList,whichdim=c(1,3))
plot(pdaList,whichdim=c(2,3))
plot(pdaList,whichdim=1:3)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.