Thank you very much to both of you, and especially you Phil.

I will tell you if it works.

2008/3/4, Phil Spector <spector@stat.berkeley.edu>:
>
> Pete -
>     As others have told you, outer only works with vectorized
> functions.  An alternative is to use expand.grid to find all
> the combinations of beta and zeta, and then use apply to
> calculate your likelihood for each row.  I believe that this
> will work:
>
> allvals = expand.grid(beta=seq(0.1,0.5,by=10^(-1)),zeta=seq(0.1,0.5
> ,by=10^(-1)))
> answer = cbind(allvals,result =
> apply(allvals,1,function(x)logl(x[1],x[2])))
>
> The columns of answer will be named beta, zeta and result, with
> (hopefully) obvious meanings.
>
>                                         - Phil Spector
>                                          Statistical Computing Facility
>                                          Department of Statistics
>                                          UC Berkeley
>                                          spector@stat.berkeley.edu
>
>
>
>
> On Tue, 4 Mar 2008, Pete Dorothy wrote:
>
> > Hello,
> >
> > I have simulated a set of data which i called "nir" (a vector).
> >
> > I have created a function "logl" which calculates the log-likelihood.
> >
> > logl is a function of 2 real parameters : "beta" and "zeta" (of length
> 1).
> >
> > This function works perfectly well when I try for example "logl(0.1,0.2
> )"
> >
> > Now if I try :
> >
> > "x=seq(0.1,0.5,by=10^(-1))
> > y=seq(0.1,0.5,by=10^(-1))
> > z=outer(x,y,logl)"
> >
> > I get an error.
> >
> > The problem seems to be that inside "logl", the following expression is
> > calculated : "sum( log( beta+(nir-1)*zeta )  )". So it is a vector
> > manipulation. The error tells me that "nir" is not the size of "zeta".
> Yet
> > usually it is no problem since "length(zeta)=1".
> >
> > When I replace "sum( log( beta+(nir-1)*zeta ) )" by a loop, I get no
> > mistake. But I think it slows down the program.
> >
> > Do you have an idea where the problem is ?
> >
> > Thank you very much.
> >
>
> >       [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@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.
> >
>

	[[alternative HTML version deleted]]

