Borrar filtros
Borrar filtros

Converting variable length vector into variable length input argument for function

3 visualizaciones (últimos 30 días)
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!

Respuesta aceptada

Walter Roberson
Walter Roberson el 24 de En. de 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 comentario
dandan
dandan el 25 de En. de 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!

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Fit Postprocessing en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by