library(ILS)
data(TG)
delta <- seq(from = 40 ,to = 850 ,length.out = 1000 )
curves.fqcd <- lab.fqcd(TG, argvals = delta)
curves.fqcs <- lab.fqcs(curves.fqcd)
summary(curves.fqcs)
names(curves.fqcs)
class(curves.fqcs$mean.i)
xlab <- "Temperature/ C"
ylab <- "Mass/ %"
main <- "Functional Mean Estimation by Laboratory"
p <- dim(curves.fqcd$curves)[3]
col <- 1:p
plot(curves.fqcs$mean.i,main = main, xlab = xlab, ylab = ylab, col = col,legend = FALSE)
legend(45,70,paste("Lab",1:p),
col = col,lty = 1,lwd = 2,cex = 0.7)
Run the code above in your browser using DataLab