# NOT RUN {
library(qcr)
m <- 30
tt<-seq(0,1,len=m)
mu<-30 * tt * (1 - tt)^(3/2)
n0 <- 100
set.seed(12345)
mdata<-matrix(NA,ncol=m,nrow=n0)
sigma <- exp(-3*as.matrix(dist(tt))/0.9)
for (i in 1:n0) mdata[i,]<- mu+0.5*mvrnorm(mu = mu,Sigma = sigma )
fdchart <- fdqcd(mdata)
plot.fdqcd(fdchart,type="l",col="gray")
set.seed(1234)
fddep <- fdqcs.depth(fdchart,plot = T)
plot(fddep,title.fdata = "Fdata",title.depth = "Depth")
summary(fddep)
# }
Run the code above in your browser using DataLab