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 <ores...@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 <d...@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-boun...@r-project.org] on behalf of Steven LeBlanc
> [ores...@gmail.com]
> > Sent: Wednesday, March 25, 2015 3:43 PM
> > To: r-help@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@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@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@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.
