The data are different sizes. (As I suggested in my first post). The turnpoints function removes "ex-aequos".
Replace the following: ------ minima<-which(tp$pit & data$residual<= 20) ----- With: ------ pits = rep (FALSE, nrow (data) ) pits [tp$pos] = tp$pits minima<-which(pits & data$residual<= -100) ------ Should remove the warning. But there's a second problem: range (data$residual) output: [1] -14.97602 11.53771 There are no "residuals less than -100. So you'll need to fix that too. B. On Sun, Feb 16, 2020 at 10:32 PM Ogbos Okike <giftedlife2...@gmail.com> wrote: > > Dear Abby, > Thank you. I will look at your stimulated data and then run the code with it. > > But since I am dealing with real data and also have volumes of it, I would > like to send my real data to you. > > The OULU05 is attached with dput function. It is labeled Ogbos_dput. I > would be surprised if it was stripped off. Then I will resend it. If, on > the other hand, you don't want it as dput data, then I will email the large > data through your private box. > Warmest regards > Ogbos > > On Sun, Feb 16, 2020, 08:14 Abby Spurdle <spurdl...@gmail.com> wrote: >> >> Note that your post does not contain a minimal reproducible example. >> I, and presumably most other readers, do not have the file "OULU05". >> >> Also, your first post referred to "% CR variation", however, your >> second post referred to "counts". >> >> I created a simple simulated data set: >> >> -------- >> sim.data = function () >> { year = 1:100 >> month = (0:99 %% 12) + 1 >> day = (0:99 %% 28) + 1 >> counts = sample (1:2000, 100) >> data.frame (year, month, day, counts) >> } >> data = sim.data () >> -------- >> >> After replacing the "data" object (as above), everything worked fine. >> (Except that the inequality needed modification based on the value of >> counts). >> >> Maybe the problem is with your dataset...? >> Or maybe there's some step in your code that results in missing >> values, given your input? >> >> Also note that the head and tail functions, are useful for both >> inspecting data, and describing your data to others. >> >> -------- >> head (data) >> tail (data) >> -------- ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.