[R] randomForest set.seed()

[Jeff Newmiller]

Please post in plain text, and supply a reproducible example (that 
includes sample data, preferably using the dput function).

The below code obtains repeatable random number sequences in each pass 
through the loop, as confirmed by the first value in each sequence.

##
a <- 0
##
## For loop to test the various combinations of ntree and mtry
for ( j in 1:8 ) {
     for ( i in 1:3 ) {
         b <- 1968 + a
         set.seed( b )
         cat( i, " ", j, " ", runif( 1 ), "\n" )
     }
}


On Thu, 9 Jul 2015, Joao Carreiras wrote:

> Dear forum members,
>
> ?I wrote a piece of code to test various combinations of mtry and ntree, so
> that the best combination (in terms of mse?) could be used. Before each
> randomForest command I included a set.seed() command, so that I can keep
> track of the seed number and replicate the results. However, when I run the
> randomForest command again with the seed number corresponding to the best
> combination I get a different mse. Any suggestions?
> Thanks and best wishes
> Joao
>
> ------------
> ## Define the number of trees to test
> nt <- c(100, 200, 300, 400, 500, 700, 900, 1000)
> ## Define the number of variables to be randomly selected at each node
> mt <- c(1, 2, 3)
> ##
> a <- 0
> ##
> ## For loop to test the various combinations of ntree and mtry
> for (j in 1:8)
>   for (i in 1:3)
> ? ?
> {
>         b <- 1968 + a
>         set.seed(b)
>         HH <- randomForest(HH_MEAN ~ ASF + PALU + FC
> ??
> ,
> ??
>                            data = model.data,
>                            ntree = nt[j],
>                            mtry = mt[i],
>                            replace = T,
> ?                     ?
> ? }?
> ?----------------?
>
> ??
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at 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.

---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnewmil at dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                       Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k