Error in nlmefit function vectorizeSinglePhiGroupX

Question

Erin MacMillan 2013-7-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/81254-error-in-nlmefit-function-vectorizesinglephigroupx

采纳的回答： the cyclist

在 MATLAB Online 中打开

Hi, I am trying to use nlmefit to do a linear random effects model of some data using the following command:

beta = nlmefit(a,b,grp,[],@(PHI,X)(X*PHI),[1 0]);

where

a=[0 1 2 0 1 2 0 1 0 1]';   b=[6.4512 5.9866 4.8052 7.0858 6.2563 5.5702 7.0976 7.7617 5.8512 7.1175]';   grp=['ABC' 'ABC' 'ABC' 'ADW' 'ADW' 'ADW' 'AEC' 'AEC' 'AEL' 'AEL'];

a,b, and grp have been shortened for brevity, they are actually 118x1 vectors.

I get the following error that I do not understand:

In an assignment  A(I) = B, the number of elements in B and I must be the same.
Error in nlmefit>getFcns/vectorizeSinglePhiGroupX (line 1367)
        fval(k) = model(phi(i,:),X(k,:),V(i,:));
Error in nlmefit>getFcns/computeFAllAndCount (line 1450)
    f = modelVectorizedAllGroups(phi, X, VVectorized);
Error in nlmefit>@(Beta,b)computeFAllAndCount(phiAllGroupsFcn(Beta,b)) (line 1459)
computeFAllGroupsFcn = @(Beta,b) computeFAllAndCount(phiAllGroupsFcn(Beta,b));
Error in nlmefit>@(Beta)computeFAllGroupsFcn(Beta,bDummy) (line 478)
    model = @(Beta) computeFAllGroupsFcn(Beta, bDummy);
Error in nlmefit>LMfit (line 1752)
yfit = model(beta);
Error in nlmefit (line 481)
    Beta_hat = LMfit(Y,Beta_hat,LMfitOptions,model,jacobian);

Thanks in advance for any clues!

Erin

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

the cyclist 2013-7-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/81254-error-in-nlmefit-function-vectorizesinglephigroupx#answer_90966

在 MATLAB Online 中打开

The proximate cause of your error is that your function specification

@(PHI,X)(X*PHI)

has only one parameter, PHI, but you have specified a vector with two initial values, [1 0].

The following syntax will execute

beta = nlmefit(a,b,grp,[],@(PHI,X)(X*PHI),[1]);

although I don't know if that is really the function you are trying to solve.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Erin MacMillan 2013-7-8

在 MATLAB Online 中打开

Thank you so much! You are right that the syntax fixed the error but was not actually the function I was trying to solve. If I want my model to be b = Xa + c, where c is some constant, then should the function be

@(PHI,X)(X*PHI(:,1)+PHI(:,2)

The nlmefit function executed without error with this model but I'm not totally sure that this is the correct model. Thanks for your help!

请先登录，再进行评论。