On Jun 7, 2012, at 12:55 , Gabor Grothendieck wrote: > On Wed, Jun 6, 2012 at 4:40 PM, May Katharina > <may.kathar...@googlemail.com> wrote: >> On Jun 6, 2012, at 10:19 , Gabor Grothendieck wrote: >> >>> On Wed, Jun 6, 2012 at 3:55 PM, May Katharina >>> <may.kathar...@googlemail.com> wrote: >>>> Hello, >>>> >>>> I'm trying to use na.spline (package zoo) to fill some missing data in a >>>> time series. >>>> this works fine, however, if I apply the 'maxgap' argument, I always get >>>> the error: >>>> <------ >>>> Error in na.spline.vec(x., coredata(object.), xout = xout., ...) : >>>> attempt to apply non-function >>>> ------> >>>> >>>> I couldn't find a similar error for this case in the mailing lists and zoo >>>> vignette, I hope I didn't >>>> miss any basic thing, but as I don't know how to proceed, any help on why >>>> this isn't working >>>> is very much appreciated (I guess I don't have to mention I'm new to time >>>> series & zoo and that I'm terrible sorry if this is >>>> an absolute dull question). >>>> >>>> The data I use is rather large, but the error can be reconstructed in a >>>> small example (s. below). >>>> If of any interest: I use R 2.15.0 GUI 1.51 Leopard build 32-bit (6148). >>>> >>>> Many thanks and best regards, Katharina >>>> >>>> > > Yes, it seems that that bug hid a second bug and I have now hopefully > removed both. Try this (or download and build the zoo package from > R-Forge): > >> library(zoo) >> >> # source revised code >> source("http://r-forge.r-project.org/scm/viewvc.php/*checkout*/pkg/zoo/R/na.spline.R?root=zoo") >> source("http://r-forge.r-project.org/scm/viewvc.php/*checkout*/pkg/zoo/R/na.approx.R?root=zoo") >> na.spline(z, maxgap = 2) > 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 > 1915 > 1 2 3 4 NA NA NA 8 9 10 11 12 13 14 NA NA > 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 > NA 18 19 20 21 22 23 24 25 26 27 28 29 30 >>
Many thanks, it works like a dream now! Great package and great support... Best regards, Katharina ______________________________________________ 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.