A very limited set of attachment types are allowed on this list.. yours was not one of them. Reading the Posting Guide will inform you about many useful things.
Self-contained R code is the expected mode of communication here. >From your description, you might consider monotonic spline fitting the diff() >of your data and integrating the fitted result... but without your sample >plots I could be way off. --------------------------------------------------------------------------- Jeff Newmiller The ..... ..... Go Live... DCN:<jdnew...@dcn.davis.ca.us> Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer (Solar/Batteries O.O#. #.O#. with /Software/Embedded Controllers) .OO#. .OO#. rocks...1k --------------------------------------------------------------------------- Sent from my phone. Please excuse my brevity. On October 23, 2014 9:35:28 AM PDT, Ingrid Charvet <ingrid.char...@rms.com> wrote: >In the figure attached "test", I would like to interpolate the empty - >"NA" section to connect the two curves. > >1/ I tried using "na.approx" and "na.spline", without success. > >na.approx(test,na.rm=FALSE) >na.spline(test,na.rm=FALSE) > >In the first case it just draws a straight line between the end and >beginning points of each curve, respectively, which doesn't look >smooth. In the second case the cubic spline interpolation results >overshoot above 1.0 which I would like to avoid (see attached test.na) >. > >2/ I have also used splinefun for a monotonic increasing spline, which >is the best result so far >Yfun <- splinefun(X,test,method="monoH.FC") >test.spline <- Yfun(X) > >however the resulting curve displays a "kink" at the limit of the >interpolation zone (see attached figure test.splinefun). > >is there a way to constrain the spline in na.spline in 1/, for example, >to force it to be monotonic strictly increasing? >Alternatively can we alter the smoothing performed by splinefun somehow >to avoid this "kink"? > >Any other suggestions welcome! >Thank you >IC > > >________________________________ >This message and any attachments contain information that may be RMS >Inc. confidential and/or privileged. If you are not the intended >recipient (or authorized to receive for the intended recipient), and >have received this message in error, any use, disclosure or >distribution is strictly prohibited. If you have received this message >in error, please notify the sender immediately by replying to the >e-mail and permanently deleting the message from your computer and/or >storage system. > > > > >------------------------------------------------------------------------ > >______________________________________________ >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. ______________________________________________ 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.