## Not run:
#
# # example for fdata class of 1 dimension (curve)
# a1<-seq(0,1,by=.01)
# a2=rnorm(length(a1),sd=0.2)
# f1<-(sin(2*pi*a1))+rnorm(length(a1),sd=0.2)
# nc<-10
# np<-length(f1)
# tt=seq(0,1,len=101)
# mdata<-matrix(NA,ncol=np,nrow=nc)
# for (i in 1:nc) mdata[i,]<- (sin(2*pi*a1))+rnorm(length(a1),sd=0.2)
# fdataobj<-fdata(mdata,tt)
# res=plot.fdata(fdataobj,type="l",col=gray(1:nrow(mdata)/nrow(mdata)))
# lines(func.mean(fdataobj),col=3,lwd=2) #original curve
#
# # example for fdata2d class of 2 dimension (surface)
# t1 <- seq(0, 1, length= 51)
# t2 <- seq(0, 1, length= 31)
# z<-array(NA,dim=c(4,51,31))
# for (i in 1:4) z[i,,] <- outer(t1, t2, function(a, b) (i*a)*(b)^i)
# z.fdata<-fdata(z,list(t1,t2))
# plot(z.fdata,time=2)
# plot(z.fdata,mfrow=c(2,2),type="persp",theta=30)
# ## End(Not run)
Run the code above in your browser using DataLab