[R] Finding which seed resulted in a specific sample

Thierry Onkelinx thierry.onkelinx at inbo.be
Mon Apr 13 12:12:38 CEST 2015

you need set.seed(seed) instead of set.seed = seed

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2015-04-13 11:12 GMT+02:00 Isabel Natario <icn op fct.unl.pt>:

> Hello!
> I was wanting to find out which random seed could generate the characters
> in the word "love", for example, when I sample with replacement from the
> vector of the letters. So I've written the code below:
> seed=0
> set.seed=0
> x<-sample(letters,4,replace=T)
> while (sum(x==c("l","o","v","e"))<4){
> seed<-seed+1
> set.seed=seed
> x<-sample(letters,4,replace=T)
> }
> When I run this code I always get that the final x vector is
> c("l","o","v","e"), but afterwards, when I try to set.seed iqual to the
> value stored in object "seed" I can never reproduce this vector
> c("l","o","v","e"), by doing sample(letters,4,replace=T). Why is that?
> Thank you very much,
> Isabel Natario
> --
> Dep. Matemática
> Faculdade de Ciências e Tecnologia
> Universidade Nova de Lisboa
> 2829-516 Caparica, Portugal
> icn op fct.unl.pt
>         [[alternative HTML version deleted]]
> ______________________________________________
> R-help op 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.

	[[alternative HTML version deleted]]

More information about the R-help mailing list