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.
