Hello.
I've just started using R and am trying to figure out if the two codes
snippets below have the same output
gBest<-floor(runif(popsize,min=1,max=top))
velocity[i,j]<<-.4* velocity[i,j] + 1 * runif(1) * (pbestsVar[i,j] -
popVar[i,j]) + 1 * runif(1) * (archiveVar[gBest,j] - popVar[i,j])
and
for (i in 1:popsize) {
for (j in 1:maxvar) {
gBest<-sample(top,size=1)
velocity[i,j]<<-.4* velocity[i,j] + 1 * runif(1) * (pbestsVar[i,j] -
popVar[i,j]) + 1 * runif(1) * (archiveVar[gBest,j] - popVar[i,j])
#velocity[i,j]<-.4* velocity[i,j] + 1 * .3722 * (pbestsVar[i,j] -
popVar[i,j]) + 1 * .3722 * (archiveVar[3,j] - popVar[i,j])
}
many thanks
--
View this message in context:
http://n4.nabble.com/For-loops-in-R-tp1015933p1015933.html
Sent from the R help mailing list archive at Nabble.com.
Hi Definitely not from first I get> gBest<-floor(runif(popsize,min=1,max=top))Error in runif(popsize, min = 1, max = top) : object 'popsize' not found>and from second> for (i in 1:popsize) {+ for (j in 1:maxvar) { + gBest<-sample(top,size=1) + velocity[i,j]<<-.4* velocity[i,j] + 1 * runif(1) * (pbestsVar[i,j] - + popVar[i,j]) + 1 * runif(1) * (archiveVar[gBest,j] - popVar[i,j]) + #velocity[i,j]<-.4* velocity[i,j] + 1 * .3722 * (pbestsVar[i,j] - + popVar[i,j]) + 1 * .3722 * (archiveVar[3,j] - popVar[i,j]) Error: unexpected ')' in: " #velocity[i,j]<-.4* velocity[i,j] + 1 * .3722 * (pbestsVar[i,j] - popVar[i,j])">If you by chance have all data necessary for those codes not to give those stupid errors you could transfer their output to different objects and check those objects if they are equivalent. ?all.equal Regards Petr r-help-bounces at r-project.org napsal dne 17.01.2010 08:17:16:> > Hello. > > I've just started using R and am trying to figure out if the two codes > snippets below have the same output > > gBest<-floor(runif(popsize,min=1,max=top)) > velocity[i,j]<<-.4* velocity[i,j] + 1 * runif(1) * (pbestsVar[i,j] - > popVar[i,j]) + 1 * runif(1) * (archiveVar[gBest,j] - popVar[i,j]) > > and > > for (i in 1:popsize) { > for (j in 1:maxvar) { > gBest<-sample(top,size=1) > velocity[i,j]<<-.4* velocity[i,j] + 1 * runif(1) *(pbestsVar[i,j] -> popVar[i,j]) + 1 * runif(1) * (archiveVar[gBest,j] - popVar[i,j]) > #velocity[i,j]<-.4* velocity[i,j] + 1 * .3722 * (pbestsVar[i,j]-> popVar[i,j]) + 1 * .3722 * (archiveVar[3,j] - popVar[i,j]) > } > > many thanks > > > -- > View this message in context:http://n4.nabble.com/For-loops-in-R-tp1015933p1015933.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 guidehttp://www.R-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code.
Hi r-help-bounces at r-project.org napsal dne 19.01.2010 02:32:45:> > Hello Petr. > > For the random values, I wanted to generate a different random numberfor> each element of my velocity matrix. > > So will this do it? > > rmat <- matrix(runif(1000), 500,2) > rmat2 <- matrix(runif(1000), 500,2) > rindex <- sample(1:500, replace=TRUE) #with repetition > velocity<-0.4 * velocity + rmat * (pbestsVar - popVar) + rmat2 * > (archiveVar[rindex,] - popVar)AFAICS it seems to do what you want. Basically your rmat and rmat2 will have different numbers from uniform distribution. If you wanted different distribution of random numbers you need to use differend random number generator. See ?rnorm. And repetition means that in your index you can have some rows more times and some rows never. Again only you know if this is desired behaviour. If not use replace=FALSE.> > Also, do the apply methods perform better than for loops given the same > function? > sample: > apply(x, fun) > and > for (i in 1:length(x)) fun(x[i])In some quite recent R-News (I believe 2009 or 2008) there is an article about loops. And also R-Inferno by P.Burns is worth reading. My opinion is that simple cycles and apply functions are more or less similar. When I want to scan a file and make same plots for each column I use for cycle and in most other cases I use *apply family. If you have nested cycles which work correctly in reasonable time why not use them. But usually vectorised approach can be far quicker and, when you get used to it, clearer. Regards Petr> > cheers > cjmr > > > -- > View this message in context:http://n4.nabble.com/For-loops-in-R-tp1015933p1017196.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 guidehttp://www.R-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code.