data(kp)
uni.C <- sort( unique(kp$Code) )
ind <- 1
Data <- kp[kp$Code==uni.C[ind], ]
x0 <- Data$x
y0 <- Data$y
Res1 <- adjdata(x0, y0, ub.np=200, 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")) )
# \donttest{
x1 <- Res1$x
y1 <- Res1$y
x0.ini <- mean( x1 )
y0.ini <- mean( y1 )
theta.ini <- pi
a.ini <- 0.9
n1.ini <- c(1, 4)
n2.ini <- 5
n3.ini <- c(5, 10, 15)
ini.val <- list(x0.ini, y0.ini, theta.ini,
a.ini, n1.ini, n2.ini, n3.ini)
Res2 <- fitGE( GE, x=x1, y=y1, ini.val=ini.val,
m=3, simpver=5, nval=1, unit="cm",
par.list=FALSE, fig.opt=TRUE, angle=NULL,
control=list(reltol=1e-20, maxit=20000),
np=2000 )
# }
graphics.off()
Run the code above in your browser using DataLab