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....@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 ______________________________________________ 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.