I am not sure about "overwrite mean() with data".  My purpose was
to generate random numbers that are from a multivariate normal
distribution with the mean vector.

For the var.f function, since I already specify m and J, so the only
variable is really rho, so I wrote it as a function of rho only.

Could you be a little more specific? Thanks a lot again.
2010/5/28 Sarah Goslee <sarah.goslee@gmail.com>

> There are a bunch of problems in your code:
> you overwrite mean() with data, and that could screw things up.
> you have a function var.f that isn't passed all the arguments it needs.
> est.4 is defined several times, each overwriting the previous.
>
> First you need to clean up these sorts of problems since they can
> lead to all kinds of bizarre results.
>
> Then, if you are still getting unexpected results, please send the list
> a minimal example so that we can take a look.
>
> Sarah
>
> On Fri, May 28, 2010 at 10:17 AM, li li <hannah.hlx@gmail.com> wrote:
> > Thanks very much for your reply.  The function is a bit long. I attached
> it.
> >
> > rho.f () is in second text document. It uses rada1.mnorm() function which
> > is contained in the first document.
> >
> > Thank you very very much!!!
> >
>
>  --
> Sarah Goslee
> http://www.functionaldiversity.org
>

	[[alternative HTML version deleted]]

