Hey David, I set it equal to the mean. duplicate="mean"
but still no joy, Thanks On Tue, Dec 10, 2013 at 4:28 PM, David Carlson <dcarl...@tamu.edu> wrote: > The error says you have duplicate points in dat so you need to > set duplicate= to tell interp() how to handle them. > > ?interp > > ------------------------------------- > David L Carlson > Department of Anthropology > Texas A&M University > College Station, TX 77840-4352 > > > -----Original Message----- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Shane Carey > Sent: Tuesday, December 10, 2013 9:36 AM > To: r-help@r-project.org > Subject: [R] 3-D interpretation > > Hi, > > im trying to create a 3-D interpretation of a geological fault > using > the akima package. But my code fails below highlighted in red: > Does anyone have any ideas why this is. > > Thanks > > dat<-read.delim("D:\\fault.txt",header=T,sep=",") > library(rgl) > # data > rgl.spheres(dat$X,dat$Z , dat$Y,1,color="red") > rgl.bbox() > # bivariate linear interpolation > # interp: > akima.li <- interp(dat$X, dat$Y, dat$Z, > xo=seq(min(dat$X), max(dat$X), length = 100), > yo=seq(min(dat$Y), max(dat$Y), length = 100)) > > This is my error: > Error in interp.old(x, y, z, xo = xo, yo = yo, ncp = 0, extrap = > extrap, : > duplicate data points: need to set 'duplicate = ..' > > # interp surface: > rgl.surface(akima.li$X,akima.li$Y,akima.li$Z,color="green",alpha > =c(0.5)) > # interpp: > akima.p <- interpp(dat$X, dat$Y, dat$Z, > runif(200,min(dat$X),max(dat$X)), > runif(200,min(dat$Y),max(dat$Y))) > # interpp points: > rgl.points(akima.p$X,akima.p$Z , > akima.p$Y,size=4,color="yellow") > # bivariate spline interpolation > # data > rgl.spheres(dat$X,dat$Z ,dat$Y,0.5,color="red") > rgl.bbox() > # bivariate cubic spline interpolation > # interp: > akima.si <- interp(dat$X, dat$Y, dat$Z, > xo=seq(min(dat$X), max(dat$X), length = 100), > yo=seq(min(dat$Y), max(dat$Y), length = 100), > linear = FALSE, extrap = TRUE) > > > > > > > -- > Shane > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible > code. > > -- Shane [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
