Great, thanks. Eric's suggestion is the most simple. Here's a resulting problem. In the call to ame.bopa in a loop, I like inputs in the call to ame.bopa to be bop1, bop2, bop3,... Thanks.
for (im in 1:m) { ame<-ame.bopa(bop,y1.level=y1value,y2.level=y2value,jindex=jindex1,vb.method="invH",joint12=TRUE, printing=FALSE,testing=TRUE) } On 6/24/2024 7:57 PM, Eric Berger wrote: > assign(paste0("bop",im),boprobit(eqs,mydata,wt=weight,method="BHHH",tol=0,reltol=0,gradtol=1e-5,Fisher=TRUE)) > > > On Mon, Jun 24, 2024 at 2:56 PM Steven Yen <st...@ntu.edu.tw> wrote: > > Thanks Eric. I am not following your suggested line. Would you > just edit my line 4? Thanks. > > On 6/24/2024 7:51 PM, Eric Berger wrote: >> How about >> >> assign(paste0("bop",im), boprobit( etc )) >> >> >> >> On Mon, Jun 24, 2024 at 2:41 PM Steven Yen <st...@ntu.edu.tw> wrote: >> >> I would like a loop to >> >> (1) read data files 2010midata1,2010midata2,2010midata3; and >> >> (2) name OUTPUT bop1,bop2,bop3. >> >> I succeeded in line 3 of the code below, >> >> BUT not line 4. The error message says: >> >> Error in paste0("bop", im) <- boprobit(eqs, mydata, wt = >> weight, method >> = "NR", : target of assignment expands to non-language object >> Please >> help. Thanks. >> >> m<-3 >> for (im in 1:m) { >> mydata<-read.csv(paste0("2010midata",im,".csv")) >> >> paste0("bop",im)<-boprobit(eqs,mydata,wt=weight,method="BHHH",tol=0,reltol=0,gradtol=1e-5,Fisher=TRUE) >> } >> >> >> >> [[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 >> <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.