data(starfish)
uni.C <- sort( unique(starfish$Code) )
ind <- 2
Data <- starfish[starfish$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=400, times=1.2, len.pro=1/20)
x1 <- Res1$x
y1 <- Res1$y
dev.new()
plot( x1, y1, asp=1, type="l", cex.lab=1.5, cex.axis=1.5,
xlab=expression(italic("x")), ylab=expression(italic("y")) )
length(x1)
graphics.off()
Run the code above in your browser using DataLab