R help - Mar 2015 - print and rbind within a for loop

Hi Dave,

Thank you, but that didn't seem to do it. This code produces the same
result:

sim.app.wald<-function(B=0.1,min=10,max=50,alpha=0.05){
	result<-c(fails=0,n=0)
	for(n in min:max){
		x<-seq(1,n-1,1)
		fhat<-x/n
		se<-sqrt((fhat*(1-fhat))/n)
		ci.wald<-fhat+qnorm(1-(alpha/2))*se
		if(sum(ci.wald[ci.wald<=B])){
			result<-rbind(cbind(fails=x[ci.wald<=B],n=n),result)
			print(cbind(X=x[ci.wald<=B],n))
		}
	return(result)
	}
}

Once again, commenting out the return() statement shows the correct output via
the print() statement.

Best Regards,
Steven

On Mar 25, 2015, at 12:51 PM, David R Forrest <drf at vims.edu> wrote:
> From looking at it it looks like nothing is done with the results of
rbind(), and the return() then exits the loop and routine with the unmodified
result.
> 
> Dave
> ________________________________________
> From: R-help [r-help-bounces at r-project.org] on behalf of Steven LeBlanc
[oreslag at gmail.com]
> Sent: Wednesday, March 25, 2015 3:43 PM
> To: r-help at R-project.org
> Subject: [R] print and rbind within a for loop
> 
> Greets,
> 
> I'm trying to iteratively find solutions to a problem given a range of
options. Code is as follows:
> 
> sim.app.wald<-function(B=0.1,min=10,max=50,alpha=0.05){
>        result<-c(fails=0,n=0)
>        for(n in min:max){
>                x<-seq(1,n-1,1)
>                fhat<-x/n
>                se<-sqrt((fhat*(1-fhat))/n)
>                ci.wald<-fhat+qnorm(1-(alpha/2))*se
>                if(sum(ci.wald[ci.wald<=B])){
> #                       rbind(cbind(fails=x[ci.wald<=B],n=n),result)
>                        print(cbind(X=x[ci.wald<=B],n))
>                }
> #       return(result)
>        }
> }
> 
> If you run it like this, you can readily see what I intend to get out;
which is the correct result. If you uncomment the two commented lines, nothing
useful obtains.
> 
> What am I missing here? This seeming codependence of print() and return()
and the failure of rbind() to do anything is quite puzzling.
> 
> Best Regards,
> Steven
>        [[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.

	[[alternative HTML version deleted]]

sim.app.wald<-function(B=0.1,min=10,max=50,alpha=0.05){
        result<-c(fails=0,n=0)
        for(n in min:max){
                x<-seq(1,n-1,1)
                ...
        return(result) # indent!  see what's wrong?
        }
}

Make your indentation (of the return line) correspond to the braces in the
code
and you may see the problem sooner.

Bill Dunlap
TIBCO Software
wdunlap tibco.com

On Wed, Mar 25, 2015 at 1:04 PM, Steven LeBlanc <oreslag at gmail.com>
wrote:
> Hi Dave,
>
> Thank you, but that didn't seem to do it. This code produces the same
> result:
>
> sim.app.wald<-function(B=0.1,min=10,max=50,alpha=0.05){
>         result<-c(fails=0,n=0)
>         for(n in min:max){
>                 x<-seq(1,n-1,1)
>                 fhat<-x/n
>                 se<-sqrt((fhat*(1-fhat))/n)
>                 ci.wald<-fhat+qnorm(1-(alpha/2))*se
>                 if(sum(ci.wald[ci.wald<=B])){
>
> result<-rbind(cbind(fails=x[ci.wald<=B],n=n),result)
>                         print(cbind(X=x[ci.wald<=B],n))
>                 }
>         return(result)
>         }
> }
>
> Once again, commenting out the return() statement shows the correct output
> via the print() statement.
>
> Best Regards,
> Steven
>
> On Mar 25, 2015, at 12:51 PM, David R Forrest <drf at vims.edu>
wrote:
>
> > From looking at it it looks like nothing is done with the results of
> rbind(), and the return() then exits the loop and routine with the
> unmodified result.
> >
> > Dave
> > ________________________________________
> > From: R-help [r-help-bounces at r-project.org] on behalf of Steven
LeBlanc
> [oreslag at gmail.com]
> > Sent: Wednesday, March 25, 2015 3:43 PM
> > To: r-help at R-project.org
> > Subject: [R] print and rbind within a for loop
> >
> > Greets,
> >
> > I'm trying to iteratively find solutions to a problem given a
range of
> options. Code is as follows:
> >
> > sim.app.wald<-function(B=0.1,min=10,max=50,alpha=0.05){
> >        result<-c(fails=0,n=0)
> >        for(n in min:max){
> >                x<-seq(1,n-1,1)
> >                fhat<-x/n
> >                se<-sqrt((fhat*(1-fhat))/n)
> >                ci.wald<-fhat+qnorm(1-(alpha/2))*se
> >                if(sum(ci.wald[ci.wald<=B])){
> > #                      
rbind(cbind(fails=x[ci.wald<=B],n=n),result)
> >                        print(cbind(X=x[ci.wald<=B],n))
> >                }
> > #       return(result)
> >        }
> > }
> >
> > If you run it like this, you can readily see what I intend to get out;
> which is the correct result. If you uncomment the two commented lines,
> nothing useful obtains.
> >
> > What am I missing here? This seeming codependence of print() and
> return() and the failure of rbind() to do anything is quite puzzling.
> >
> > Best Regards,
> > Steven
> >        [[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.
>
>
>         [[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.
>
	[[alternative HTML version deleted]]