Converting variable length vector into variable length input argument for function

3 vues (au cours des 30 derniers jours)
I'm trying to turn a function output vector of variable length into a variable number of arguments to be fed into another function, so each element of the output vector becomes one input argument for the next function below:
function [pd]=fit(xdata,fitname,pdfunc,start)
[params]=mle(xdata, 'pdf',pdfunc, 'start',start);
pd=makedist(fitname,deal(num2cell(params)));
end
xdata is some data array, fitname is something like 'Exponential' or 'GeneralizedPareto', pdfunc is an anonymous function such as pdfgp=@(x,k,sigma,theta) gppdf(x,k,sigma,theta), and start is the initial parameter values required by mle.
Right now this produces a "Not enough input arguments" error. I've also tried variants on using eval, deal and num2cell with no luck. Any help would be appreciated. Thanks!

Réponse acceptée

Walter Roberson
Walter Roberson le 24 Jan 2018
This will not do any good. The output of mle is a numeric parameter estimate scalar or row vector. The parameters of makedist must be first the fitting name, and then name/value pairs.
You can use something like
temp = num2cell(params);
param_names = {'mu', 'sigma'};
pv_pairs = [param_names; temp];
pd = makedist(fitname, pv_pairs{:});
This counts on cell expansion going down columns first.
  1 commentaire
dandan
dandan le 25 Jan 2018
Makedist actually does accept the parameter values without names (at least in my 2016b version of Matlab). For example,
pd = makedist('GeneralizedPareto',1,1,0)
appears to work just as well as
pd = makedist('GeneralizedPareto','k',1,'sigma',1,'theta',0)
Based on your suggestion though, I realized I was just way overcomplicating it, and it turns out
temp=num2cell(params);
pd = makedist(fitname,temp{:});
works perfectly! Thank you!

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Startup and Shutdown dans Help Center et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by