data(bambooleaves)
uni.C <- sort( unique(bambooleaves$Code) )
ind <- 1
Data <- bambooleaves[bambooleaves$Code==uni.C[ind], ]
x0 <- Data$x
y0 <- Data$y
dev.new()
plot( x0, y0, asp=1, type="l", cex.lab=1.5, cex.axis=1.5,
xlab=expression(italic("x")), ylab=expression(italic("y")) )
length(x0)
Res1 <- adjdata(x0, y0, ub.np=600, len.pro=1/20)
dev.new()
plot( Res1$x, Res1$y, asp=1, cex.lab=1.5, cex.axis=1.5, type="l",
xlab=expression(italic("x")), ylab=expression(italic("y")) )
graphics.off()
Run the code above in your browser using DataLab