[R] Random Forest prediction questions

Liaw, Andy andy_liaw at merck.com
Mon Mar 1 16:09:29 CET 2010


From: Dror
> 
> Hi,
> I need help with the randomForest prediction. i run the folowing code:
>  
> > iris.rf <- randomForest(Species ~ ., data=iris,
> > importance=TRUE,keep.forest=TRUE, proximity=TRUE)
> > pr<-predict(iris.rf,iris,predict.all=T)
> > iris.rf$votes[53,]
>     setosa versicolor  virginica 
>  0.0000000  0.8074866  0.1925134 
> > table(pr$individual[53,])/500
> 
> versicolor  virginica 
>      0.928      0.072 
> > 
> 
> why the voting is not the same for the same data? what do i do wrong?

It's because the $votes components reflects the OOB predictions, whereas
predict() gives you predictions based on all of the trees in the forest.
 
> another 2 questions:
> 1. i tries to debug another problem in which the individual vector was
> smaller the tree number in the forest.
>    i noticed that in this row of code:
> 
>    treepred <- matrix(object$classes[t1$treepred], nrow = 
> length(keep), 
>                 dimnames = list(rn[keep], NULL))
> the t1$treepred has values of 0 (i have 2 classes) and they 
> droped from the
> results
> what does this 0 mean?

Not sure why you're debugging that portion of the code.  That is just to
dimension the array passed back from C into a matrix.  What is "t1"?

> 2. how can i drop a tree from the forest?

Look at the $forest component of the randomForest object, and subset the
dimension that correspond to ntree in all of its components.  Change
$ntree accordingly.

Andy
 
> Thanks,
> Dror
> -- 
> View this message in context: 
> http://n4.nabble.com/Random-Forest-prediction-questions-tp1573
> 530p1573530.html
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> R-help at r-project.org 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.
> 
Notice:  This e-mail message, together with any attachme...{{dropped:10}}



More information about the R-help mailing list