fda.usc (version 0.9.4)

predict.classif.fd: Predicts from a fitted classif.fd object.

Description

Classifier of Functional Data by kernel method using Functional Data Object of class classif.fd.

Usage

predict.classif.fd(object, new.fdataobj,probs.group=FALSE,...)

Arguments

object

Object object estimated by: k nearest neighbors method classif.knn.fd, kernel method classif.kernel.fd or kernel method with fucntional basis representation classif.kernel.fb.

new.fdataobj

New functional explanatory data of fdata class.

probs.group

=TRUE, for each functional data shows the probability of each group membership.

...

Further arguments passed to or from other methods.

Value

• group.pred
  Vector of groups or classes predicted
• prob.group
  For each functional data shows the probability of each group membership.

Details

Returns the groups or classes predicted using a previously trained model.

References

Ferraty, F. and Vieu, P. (2006). Nonparametricc functional data analysis. Springer Series in Statistics, New York. Ferraty, F. and Vieu, P. (2006). NPFDA in practice. Free access on line at http://www.lsp.ups-tlse.fr/staph/npfda/ Ramsay, James O., and Silverman, Bernard W. (2006), Functional Data Analysis, 2nd ed., Springer, New York. Ramsay, James O., and Silverman, Bernard W. (2002), Applied Functional Data Analysis, Springer, New York.

See Also

See also classif.knn.fd, classif.kernel.fd, classif.kernel.fb and summary.classif.fd

Examples

data(phoneme)
mlearn<-phoneme[["learn"]]
glearn<-phoneme[["classlearn"]]
mtest<-phoneme[["test"]]
gtest<-phoneme[["classtest"]]

#	ESTIMATION knn
out1=classif.knn.fd(mlearn,glearn,knn=c(3,5,7))
summary.classif.fd(out1)
#	PREDICTION knn
pred1=predict.classif.fd(out1,mtest,TRUE)
table(pred1$group.pred,gtest)
MisclasPredict <- sum(pred1$group.pred != gtest)/length(gtest)
MisclasPredict

#	ESTIMATION kernel
h=2^(0:5)
out2=classif.kernel.fd(mlearn,glearn,h=h)
summary.classif.fd(out2)
#	PREDICTION kernel
pred2=predict.classif.fd(out2,mtest,TRUE)
table(pred2$group.pred,gtest)
MisclasPredict <- sum(pred2$group.pred != gtest)/length(gtest)
MisclasPredict

#	ESTIMATION kernel functional basis
vbasis=seq(7,37,by=10)
out3=classif.kernel.fb(mlearn,glearn,h=h,type.basis="fourier",
par.basis=list(nbasis=vbasis))
summary.classif.fd(out3)
#	PREDICTION kernel functional basis
pred3=predict.classif.fd(out3,mtest,TRUE)
table(pred3$group.pred,gtest)
MisclasPredict <- sum(pred3$group.pred != gtest)/length(gtest)
MisclasPredict