|
From: | Feldman, Tracy |
Subject: | RE: request for help with fmins or fminbnd |
Date: | Fri, 31 Aug 2007 08:48:57 -0500 |
Dear Michael, Thank you for your response. In this
program (the one that worked in Matlab), I am actually using maximum likelihood
methods to calculate the Poisson confidence interval around the mean of a response
variable in each of several different treatments. I get the impression
that the functions you suggest only give the mle (and AIC, etc.) with standard
deviations on these mle’s. I believe that the fmins or some related
function could work for me if I only knew the syntax... Again, I appreciate
your help. Take care, Please reply to address@hidden or address@hidden From:
address@hidden [mailto:address@hidden On Behalf Of Michael Creel On 8/31/07, To whom it may concern: I am trying to convert a program that had worked well in Matlab, so
that I can use an equivalent program in Octave (I no longer have access to
Matlab). The program calculates the maximum likelihood estimate and the
best-fit parameters for a function (finding poisson confidence limits, in this
case, or in other cases, fitting a logistic function with binomial
error). I have run into difficulties with the fmins function (which
I have used in Octave to replace "fminsearch" in Matlab).
I am trying to call my likelihood function as follows: [maxparam, maxlike] = fmins(@Likelihood_pois1, initparammax, [], [],
pir, boundslowmax, boundsupmax); maxparam is a vector of two values of the maximized function It is received by the Likelihood function as follows: function neglogL = Likelihood_pois1(pp, initparammax, pir,
boundslowmax, boundsupmax); I am not sure what would be the correct syntax, but whatever order I
choose for the variables in parentheses (or if I get rid of one or both of the
empty brackets) I get an error message saying invalid vector index = 4
(or some other #), and that fmins cannot perform any of several binary operations.
I have also tried converting to syntax for fminbnd, hoping it would be
better, but this will not work for me either. Is this more than a simple
syntax issue? I would appreciate any suggestions you might have (or if you would
like more code from me and/or the program files, I can send more--or original
files that worked in Matlab). I really appreciate your time and help. I am not currently
subscribed to the list, so if you can, please respond to me directly at address@hidden. Sincerely, Tracy S. Feldman Postdoctoral Associate The Samuel Roberts Noble Foundation Got a little couch potato?
|
[Prev in Thread] | Current Thread | [Next in Thread] |