How to delete repeated values in MCMC sampling and get index of unique values?

7 messages
Open this post in threaded view
|

How to delete repeated values in MCMC sampling and get index of unique values?

 Dear list, How do you delete repeated samples?  In MCMC, when your candidate value has been reject, so you remain on the same point, so you keep that value. Say I have this toy example, > c(1,6,6,6,3,5,4,4,2,3,5) The 6 and 4 are repeated, I only want the index of the non-repeated values. I thought of using which() and unique(), but that does not give you the index of the unique values. Thanks in advance. Mike         [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: How to delete repeated values in MCMC sampling and get index of unique values?

 On 13-01-17 4:50 PM, C W wrote: > Dear list, > How do you delete repeated samples?  In MCMC, when your candidate value has > been reject, so you remain on the same point, so you keep that value. > > Say I have this toy example, > >> c(1,6,6,6,3,5,4,4,2,3,5) > > The 6 and 4 are repeated, I only want the index of the non-repeated values. > > I thought of using which() and unique(), but that does not give you the > index of the unique values. You could use x[!duplicated(x)] or rle(x)\$values, depending on your definition of "repeated".  I hope you're aware that you can't use either for things like quantiles and moments of the limiting distribution.  > x <- c(1,6,6,6,3,5,4,4,2,3,5)  > x[!duplicated(x)] [1] 1 6 3 5 4 2  > rle(x)\$values [1] 1 6 3 5 4 2 3 5 Duncan Murdoch ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: How to delete repeated values in MCMC sampling and get index of unique values?

 Exactly what I am looking for. Thanks a lot! Mike On Thu, Jan 17, 2013 at 4:59 PM, Duncan Murdoch <[hidden email]>wrote: > On 13-01-17 4:50 PM, C W wrote: > >> Dear list, >> How do you delete repeated samples?  In MCMC, when your candidate value >> has >> been reject, so you remain on the same point, so you keep that value. >> >> Say I have this toy example, >> >>  c(1,6,6,6,3,5,4,4,2,3,5) >>> >> >> The 6 and 4 are repeated, I only want the index of the non-repeated >> values. >> >> I thought of using which() and unique(), but that does not give you the >> index of the unique values. >> > > You could use x[!duplicated(x)] or rle(x)\$values, depending on your > definition of "repeated".  I hope you're aware that you can't use either > for things like quantiles and moments of the limiting distribution. > > > x <- c(1,6,6,6,3,5,4,4,2,3,5) > > x[!duplicated(x)] > [1] 1 6 3 5 4 2 > > rle(x)\$values > [1] 1 6 3 5 4 2 3 5 > > Duncan Murdoch > >         [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: How to delete repeated values in MCMC sampling and get index of unique values?

 In reply to this post by tmrsg11 On Jan 17, 2013, at 1:50 PM, C W wrote: > Dear list, > How do you delete repeated samples?  In MCMC, when your candidate value has > been reject, so you remain on the same point, so you keep that value. > > Say I have this toy example, > >> c(1,6,6,6,3,5,4,4,2,3,5) > c(1,6,6,6,3,5,4,4,2,3,5)[!duplicated(c(1,6,6,6,3,5,4,4,2,3,5))] [1] 1 6 3 5 4 2 > x <- c(1,6,6,6,3,5,4,4,2,3,5) > x[!duplicated(x)] [1] 1 6 3 5 4 2 > !duplicated(x)  [1]  TRUE  TRUE FALSE FALSE  TRUE  TRUE  TRUE FALSE  TRUE FALSE FALSE > which(!duplicated(x)) [1] 1 2 5 6 7 9 > > The 6 and 4 are repeated, I only want the index of the non-repeated values. > > I thought of using which() and unique(), but that does not give you the > index of the unique values. > > Thanks in advance. > > Mike > > [[alternative HTML version deleted]] > > ______________________________________________ > [hidden email] 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. David Winsemius Alameda, CA, USA ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: How to delete repeated values in MCMC sampling and get index of unique values?

 In reply to this post by Duncan Murdoch-2 What answer is wanted for c(1,1,1,2,3,1)   ? Note that Duncan's  two suggestions below give different answers for this. -- Bert On Thu, Jan 17, 2013 at 1:59 PM, Duncan Murdoch <[hidden email]> wrote: > On 13-01-17 4:50 PM, C W wrote: >> >> Dear list, >> How do you delete repeated samples?  In MCMC, when your candidate value >> has >> been reject, so you remain on the same point, so you keep that value. >> >> Say I have this toy example, >> >>> c(1,6,6,6,3,5,4,4,2,3,5) >> >> >> The 6 and 4 are repeated, I only want the index of the non-repeated >> values. >> >> I thought of using which() and unique(), but that does not give you the >> index of the unique values. > > > You could use x[!duplicated(x)] or rle(x)\$values, depending on your > definition of "repeated".  I hope you're aware that you can't use either for > things like quantiles and moments of the limiting distribution. > >> x <- c(1,6,6,6,3,5,4,4,2,3,5) >> x[!duplicated(x)] > [1] 1 6 3 5 4 2 >> rle(x)\$values > [1] 1 6 3 5 4 2 3 5 > > Duncan Murdoch > > ______________________________________________ > [hidden email] 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. -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.