How to obtain a position of last non-zero element Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 23, 2019 at 23:30UTC (7:30pm US/Eastern) Data science time! April 2019 and salary with experience The Ask Question Wizard is Live!How to trim leading and trailing whitespace?How do I replace NA values with zeros in an R dataframe?data.table vs dplyr: can one do something well the other can't or does poorly?Calculate days since last event in REfficient way of taking the max date within groupsTurn off verbose messages when loading tidyverse using library() functionColumn name of last non-NA row per row; using tidyverse solution?Filtering data relative to first and last occurance of an eventdplyr approach to get the last row number with a positive valueA code to imput missing values with linear dependency

How does light 'choose' between wave and particle behaviour?

What order were files/directories output in dir?

How can I set the aperture on my DSLR when it's attached to a telescope instead of a lens?

Sliceness of knots

Why are my pictures showing a dark band on one edge?

If Windows 7 doesn't support WSL, then what is "Subsystem for UNIX-based Applications"?

Why does it sometimes sound good to play a grace note as a lead in to a note in a melody?

A term for a woman complaining about things/begging in a cute/childish way

Electrolysis of water: Which equations to use? (IB Chem)

Intuitive explanation of the rank-nullity theorem

Did Mueller's report provide an evidentiary basis for the claim of Russian govt election interference via social media?

Customizing QGIS plugins

Lagrange four-squares theorem --- deterministic complexity

How to pronounce 伝統色

Where is the Data Import Wizard Error Log

What is an "asse" in Elizabethan English?

Do I really need to have a message in a novel to appeal to readers?

Misunderstanding of Sylow theory

Did any compiler fully use 80-bit floating point?

Crossing US/Canada Border for less than 24 hours

What does Turing mean by this statement?

Drawing spherical mirrors

How long can equipment go unused before powering up runs the risk of damage?

How much damage would a cupful of neutron star matter do to the Earth?



How to obtain a position of last non-zero element



Announcing the arrival of Valued Associate #679: Cesar Manara
Planned maintenance scheduled April 23, 2019 at 23:30UTC (7:30pm US/Eastern)
Data science time! April 2019 and salary with experience
The Ask Question Wizard is Live!How to trim leading and trailing whitespace?How do I replace NA values with zeros in an R dataframe?data.table vs dplyr: can one do something well the other can't or does poorly?Calculate days since last event in REfficient way of taking the max date within groupsTurn off verbose messages when loading tidyverse using library() functionColumn name of last non-NA row per row; using tidyverse solution?Filtering data relative to first and last occurance of an eventdplyr approach to get the last row number with a positive valueA code to imput missing values with linear dependency



.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty height:90px;width:728px;box-sizing:border-box;








18















I've got a binary variable representing if event happened or not:



event <- c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0)


I need to obtain a variable that would indicate the time when the last event happened. The expected output would be:



last_event <- c(0, 0, 0, 0, 5, 5, 5, 5, 5, 5, 5, 5, 13, 13, 13, 13)


How can I obtain that with base R, tidyverse or any other way?










share|improve this question




























    18















    I've got a binary variable representing if event happened or not:



    event <- c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0)


    I need to obtain a variable that would indicate the time when the last event happened. The expected output would be:



    last_event <- c(0, 0, 0, 0, 5, 5, 5, 5, 5, 5, 5, 5, 13, 13, 13, 13)


    How can I obtain that with base R, tidyverse or any other way?










    share|improve this question
























      18












      18








      18


      1






      I've got a binary variable representing if event happened or not:



      event <- c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0)


      I need to obtain a variable that would indicate the time when the last event happened. The expected output would be:



      last_event <- c(0, 0, 0, 0, 5, 5, 5, 5, 5, 5, 5, 5, 13, 13, 13, 13)


      How can I obtain that with base R, tidyverse or any other way?










      share|improve this question














      I've got a binary variable representing if event happened or not:



      event <- c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0)


      I need to obtain a variable that would indicate the time when the last event happened. The expected output would be:



      last_event <- c(0, 0, 0, 0, 5, 5, 5, 5, 5, 5, 5, 5, 13, 13, 13, 13)


      How can I obtain that with base R, tidyverse or any other way?







      r tidyverse base






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked Apr 11 at 14:08









      jakesjakes

      481315




      481315






















          4 Answers
          4






          active

          oldest

          votes


















          18














          Taking advantage of the fact that you have a binary vector, the following gives your desired output:



          cummax(seq_along(event) * event)





          share|improve this answer


















          • 6





            Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

            – Konrad Rudolph
            Apr 11 at 14:27






          • 3





            or without multiplication cummax(ifelse(event, seq_along(event), 0))

            – jogo
            Apr 11 at 14:27











          • @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

            – Konrad Rudolph
            Apr 11 at 14:28


















          8














          Whenever you need to fill repetitions with a value, think run-length encoding.



          In this case, you can determine the run lengths and then repeat the indices of count == 0 an according number of times:



          lengths = rle(event == 0)$lengths
          nonzeros = which(event != 0)
          runs = c(0, rep(nonzeros, each = 2))
          result = rep(runs, lengths)


          Alternative, substitute the runs in the RLE and then inverse it:



          rle = rle(event == 0)
          nonzeros = which(event != 0)
          rle$values = c(0, rep(nonzeros, each = 2))
          result = inverse.rle(rle)





          share|improve this answer






























            1














            You can also do somthing like this-



            > zero.locf <- function(x) 
            v <- x!=0
            c(0, x[v])[cumsum(v)+1]


            > zero.locf(1:length(event)*event)

            [1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





            share|improve this answer






























              1














              Another option is to find the index where event == 1 and repeat it based on length.



              rep(c(0, which(event == 1)), tapply(event, cumsum(event == 1), length))
              #[1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





              share|improve this answer























                Your Answer






                StackExchange.ifUsing("editor", function ()
                StackExchange.using("externalEditor", function ()
                StackExchange.using("snippets", function ()
                StackExchange.snippets.init();
                );
                );
                , "code-snippets");

                StackExchange.ready(function()
                var channelOptions =
                tags: "".split(" "),
                id: "1"
                ;
                initTagRenderer("".split(" "), "".split(" "), channelOptions);

                StackExchange.using("externalEditor", function()
                // Have to fire editor after snippets, if snippets enabled
                if (StackExchange.settings.snippets.snippetsEnabled)
                StackExchange.using("snippets", function()
                createEditor();
                );

                else
                createEditor();

                );

                function createEditor()
                StackExchange.prepareEditor(
                heartbeatType: 'answer',
                autoActivateHeartbeat: false,
                convertImagesToLinks: true,
                noModals: true,
                showLowRepImageUploadWarning: true,
                reputationToPostImages: 10,
                bindNavPrevention: true,
                postfix: "",
                imageUploader:
                brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
                contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
                allowUrls: true
                ,
                onDemand: true,
                discardSelector: ".discard-answer"
                ,immediatelyShowMarkdownHelp:true
                );



                );













                draft saved

                draft discarded


















                StackExchange.ready(
                function ()
                StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f55634527%2fhow-to-obtain-a-position-of-last-non-zero-element%23new-answer', 'question_page');

                );

                Post as a guest















                Required, but never shown

























                4 Answers
                4






                active

                oldest

                votes








                4 Answers
                4






                active

                oldest

                votes









                active

                oldest

                votes






                active

                oldest

                votes









                18














                Taking advantage of the fact that you have a binary vector, the following gives your desired output:



                cummax(seq_along(event) * event)





                share|improve this answer


















                • 6





                  Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

                  – Konrad Rudolph
                  Apr 11 at 14:27






                • 3





                  or without multiplication cummax(ifelse(event, seq_along(event), 0))

                  – jogo
                  Apr 11 at 14:27











                • @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

                  – Konrad Rudolph
                  Apr 11 at 14:28















                18














                Taking advantage of the fact that you have a binary vector, the following gives your desired output:



                cummax(seq_along(event) * event)





                share|improve this answer


















                • 6





                  Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

                  – Konrad Rudolph
                  Apr 11 at 14:27






                • 3





                  or without multiplication cummax(ifelse(event, seq_along(event), 0))

                  – jogo
                  Apr 11 at 14:27











                • @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

                  – Konrad Rudolph
                  Apr 11 at 14:28













                18












                18








                18







                Taking advantage of the fact that you have a binary vector, the following gives your desired output:



                cummax(seq_along(event) * event)





                share|improve this answer













                Taking advantage of the fact that you have a binary vector, the following gives your desired output:



                cummax(seq_along(event) * event)






                share|improve this answer












                share|improve this answer



                share|improve this answer










                answered Apr 11 at 14:24









                mgiormentimgiormenti

                444211




                444211







                • 6





                  Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

                  – Konrad Rudolph
                  Apr 11 at 14:27






                • 3





                  or without multiplication cummax(ifelse(event, seq_along(event), 0))

                  – jogo
                  Apr 11 at 14:27











                • @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

                  – Konrad Rudolph
                  Apr 11 at 14:28












                • 6





                  Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

                  – Konrad Rudolph
                  Apr 11 at 14:27






                • 3





                  or without multiplication cummax(ifelse(event, seq_along(event), 0))

                  – jogo
                  Apr 11 at 14:27











                • @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

                  – Konrad Rudolph
                  Apr 11 at 14:28







                6




                6





                Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

                – Konrad Rudolph
                Apr 11 at 14:27





                Yes! So much more elegant than my solution. I was thinking about cumulative sums but I didn’t think of multiplying the indices by the binary vector.

                – Konrad Rudolph
                Apr 11 at 14:27




                3




                3





                or without multiplication cummax(ifelse(event, seq_along(event), 0))

                – jogo
                Apr 11 at 14:27





                or without multiplication cummax(ifelse(event, seq_along(event), 0))

                – jogo
                Apr 11 at 14:27













                @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

                – Konrad Rudolph
                Apr 11 at 14:28





                @jogo That solution makes sense if the type of event is logical. It does work even for a numeric vector due to R’s implicit conversions but … eh.

                – Konrad Rudolph
                Apr 11 at 14:28













                8














                Whenever you need to fill repetitions with a value, think run-length encoding.



                In this case, you can determine the run lengths and then repeat the indices of count == 0 an according number of times:



                lengths = rle(event == 0)$lengths
                nonzeros = which(event != 0)
                runs = c(0, rep(nonzeros, each = 2))
                result = rep(runs, lengths)


                Alternative, substitute the runs in the RLE and then inverse it:



                rle = rle(event == 0)
                nonzeros = which(event != 0)
                rle$values = c(0, rep(nonzeros, each = 2))
                result = inverse.rle(rle)





                share|improve this answer



























                  8














                  Whenever you need to fill repetitions with a value, think run-length encoding.



                  In this case, you can determine the run lengths and then repeat the indices of count == 0 an according number of times:



                  lengths = rle(event == 0)$lengths
                  nonzeros = which(event != 0)
                  runs = c(0, rep(nonzeros, each = 2))
                  result = rep(runs, lengths)


                  Alternative, substitute the runs in the RLE and then inverse it:



                  rle = rle(event == 0)
                  nonzeros = which(event != 0)
                  rle$values = c(0, rep(nonzeros, each = 2))
                  result = inverse.rle(rle)





                  share|improve this answer

























                    8












                    8








                    8







                    Whenever you need to fill repetitions with a value, think run-length encoding.



                    In this case, you can determine the run lengths and then repeat the indices of count == 0 an according number of times:



                    lengths = rle(event == 0)$lengths
                    nonzeros = which(event != 0)
                    runs = c(0, rep(nonzeros, each = 2))
                    result = rep(runs, lengths)


                    Alternative, substitute the runs in the RLE and then inverse it:



                    rle = rle(event == 0)
                    nonzeros = which(event != 0)
                    rle$values = c(0, rep(nonzeros, each = 2))
                    result = inverse.rle(rle)





                    share|improve this answer













                    Whenever you need to fill repetitions with a value, think run-length encoding.



                    In this case, you can determine the run lengths and then repeat the indices of count == 0 an according number of times:



                    lengths = rle(event == 0)$lengths
                    nonzeros = which(event != 0)
                    runs = c(0, rep(nonzeros, each = 2))
                    result = rep(runs, lengths)


                    Alternative, substitute the runs in the RLE and then inverse it:



                    rle = rle(event == 0)
                    nonzeros = which(event != 0)
                    rle$values = c(0, rep(nonzeros, each = 2))
                    result = inverse.rle(rle)






                    share|improve this answer












                    share|improve this answer



                    share|improve this answer










                    answered Apr 11 at 14:23









                    Konrad RudolphKonrad Rudolph

                    404k1017951041




                    404k1017951041





















                        1














                        You can also do somthing like this-



                        > zero.locf <- function(x) 
                        v <- x!=0
                        c(0, x[v])[cumsum(v)+1]


                        > zero.locf(1:length(event)*event)

                        [1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





                        share|improve this answer



























                          1














                          You can also do somthing like this-



                          > zero.locf <- function(x) 
                          v <- x!=0
                          c(0, x[v])[cumsum(v)+1]


                          > zero.locf(1:length(event)*event)

                          [1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





                          share|improve this answer

























                            1












                            1








                            1







                            You can also do somthing like this-



                            > zero.locf <- function(x) 
                            v <- x!=0
                            c(0, x[v])[cumsum(v)+1]


                            > zero.locf(1:length(event)*event)

                            [1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





                            share|improve this answer













                            You can also do somthing like this-



                            > zero.locf <- function(x) 
                            v <- x!=0
                            c(0, x[v])[cumsum(v)+1]


                            > zero.locf(1:length(event)*event)

                            [1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13






                            share|improve this answer












                            share|improve this answer



                            share|improve this answer










                            answered Apr 11 at 14:30









                            RushabhRushabh

                            1,332222




                            1,332222





















                                1














                                Another option is to find the index where event == 1 and repeat it based on length.



                                rep(c(0, which(event == 1)), tapply(event, cumsum(event == 1), length))
                                #[1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





                                share|improve this answer



























                                  1














                                  Another option is to find the index where event == 1 and repeat it based on length.



                                  rep(c(0, which(event == 1)), tapply(event, cumsum(event == 1), length))
                                  #[1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





                                  share|improve this answer

























                                    1












                                    1








                                    1







                                    Another option is to find the index where event == 1 and repeat it based on length.



                                    rep(c(0, which(event == 1)), tapply(event, cumsum(event == 1), length))
                                    #[1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13





                                    share|improve this answer













                                    Another option is to find the index where event == 1 and repeat it based on length.



                                    rep(c(0, which(event == 1)), tapply(event, cumsum(event == 1), length))
                                    #[1] 0 0 0 0 5 5 5 5 5 5 5 5 13 13 13 13






                                    share|improve this answer












                                    share|improve this answer



                                    share|improve this answer










                                    answered Apr 11 at 14:31









                                    Ronak ShahRonak Shah

                                    48.5k104369




                                    48.5k104369



























                                        draft saved

                                        draft discarded
















































                                        Thanks for contributing an answer to Stack Overflow!


                                        • Please be sure to answer the question. Provide details and share your research!

                                        But avoid


                                        • Asking for help, clarification, or responding to other answers.

                                        • Making statements based on opinion; back them up with references or personal experience.

                                        To learn more, see our tips on writing great answers.




                                        draft saved


                                        draft discarded














                                        StackExchange.ready(
                                        function ()
                                        StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f55634527%2fhow-to-obtain-a-position-of-last-non-zero-element%23new-answer', 'question_page');

                                        );

                                        Post as a guest















                                        Required, but never shown





















































                                        Required, but never shown














                                        Required, but never shown












                                        Required, but never shown







                                        Required, but never shown

































                                        Required, but never shown














                                        Required, but never shown












                                        Required, but never shown







                                        Required, but never shown







                                        Popular posts from this blog

                                        រឿង រ៉ូមេអូ និង ហ្ស៊ុយលីយេ សង្ខេបរឿង តួអង្គ បញ្ជីណែនាំ

                                        Crop image to path created in TikZ? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)Crop an inserted image?TikZ pictures does not appear in posterImage behind and beyond crop marks?Tikz picture as large as possible on A4 PageTransparency vs image compression dilemmaHow to crop background from image automatically?Image does not cropTikzexternal capturing crop marks when externalizing pgfplots?How to include image path that contains a dollar signCrop image with left size given

                                        Romeo and Juliet ContentsCharactersSynopsisSourcesDate and textThemes and motifsCriticism and interpretationLegacyScene by sceneSee alsoNotes and referencesSourcesExternal linksNavigation menu"Consumer Price Index (estimate) 1800–"10.2307/28710160037-3222287101610.1093/res/II.5.31910.2307/45967845967810.2307/2869925286992510.1525/jams.1982.35.3.03a00050"Dada Masilo: South African dancer who breaks the rules"10.1093/res/os-XV.57.1610.2307/28680942868094"Sweet Sorrow: Mann-Korman's Romeo and Juliet Closes Sept. 5 at MN's Ordway"the original10.2307/45957745957710.1017/CCOL0521570476.009"Ram Leela box office collections hit massive Rs 100 crore, pulverises prediction"Archived"Broadway Revival of Romeo and Juliet, Starring Orlando Bloom and Condola Rashad, Will Close Dec. 8"Archived10.1075/jhp.7.1.04hon"Wherefore art thou, Romeo? To make us laugh at Navy Pier"the original10.1093/gmo/9781561592630.article.O006772"Ram-leela Review Roundup: Critics Hail Film as Best Adaptation of Romeo and Juliet"Archived10.2307/31946310047-77293194631"Romeo and Juliet get Twitter treatment""Juliet's Nurse by Lois Leveen""Romeo and Juliet: Orlando Bloom's Broadway Debut Released in Theaters for Valentine's Day"Archived"Romeo and Juliet Has No Balcony"10.1093/gmo/9781561592630.article.O00778110.2307/2867423286742310.1076/enst.82.2.115.959510.1080/00138380601042675"A plague o' both your houses: error in GCSE exam paper forces apology""Juliet of the Five O'Clock Shadow, and Other Wonders"10.2307/33912430027-4321339124310.2307/28487440038-7134284874410.2307/29123140149-661129123144728341M"Weekender Guide: Shakespeare on The Drive""balcony"UK public library membership"romeo"UK public library membership10.1017/CCOL9780521844291"Post-Zionist Critique on Israel and the Palestinians Part III: Popular Culture"10.2307/25379071533-86140377-919X2537907"Capulets and Montagues: UK exam board admit mixing names up in Romeo and Juliet paper"Istoria Novellamente Ritrovata di Due Nobili Amanti2027/mdp.390150822329610820-750X"GCSE exam error: Board accidentally rewrites Shakespeare"10.2307/29176390149-66112917639"Exam board apologises after error in English GCSE paper which confused characters in Shakespeare's Romeo and Juliet""From Mariotto and Ganozza to Romeo and Guilietta: Metamorphoses of a Renaissance Tale"10.2307/37323537323510.2307/2867455286745510.2307/28678912867891"10 Questions for Taylor Swift"10.2307/28680922868092"Haymarket Theatre""The Zeffirelli Way: Revealing Talk by Florentine Director""Michael Smuin: 1938-2007 / Prolific dance director had showy career"The Life and Art of Edwin BoothRomeo and JulietRomeo and JulietRomeo and JulietRomeo and JulietEasy Read Romeo and JulietRomeo and Julieteeecb12003684p(data)4099369-3n8211610759dbe00d-a9e2-41a3-b2c1-977dd692899302814385X313670221313670221