Gaussian fit function question

Emily Stadnicki

2020 10 13

0 个回答

3 次查看（30 天）

0 个投票

I am trying to fit my data to a gaussian curve after cleaning it up for only values that are significant with a p less than 0.05. I have loaded the data and gotten the ANOVA values for all 40 experiments and now must create a gaussian curve but my coade keeps giving me an error at the line when i define the fit. I have attached what I have so far. Any input is helpful! Thank you.

function [outcomes,a,b,c,d] = fit_arm_nerve_practice('folder') 
subdirs = dir('folder');
outcomes = [];
a = []; b = []; c = []; d = [];
for i=1:length(subdirs),
    if subdirs(i).isdir&~strcmp(subdirs(i).name,'.')&~strcmp(subdirs(i).name,'..'),
        filename = ['folder' filesep subdirs(i).name filesep 'mouse_arm_data.txt'];
        if exist(filename),
            disp(['Analyzing file ' filename '.']);
            data = load(filename,'-ascii');
            locations = data(1,:); 
rawdata = data(2:end,:); 
            outcomes(i) = anova1([rawdata;],locations,'off');
            if outcomes<0.05,
% write your own gaussfit function, where x and y
% are your data points
fo = fitoptions('Method','NonLinearLeastSquares','StartPoint',[1 1]);
gauss = fittype('a+b*exp(-((x-c).^2)/((2*d^2)))','options',fo);
gauss = setoptions(gauss,fo);
[gof] = fit(locations',rawdata',gauss);
[a,b,c,d] = gaussfit(x,y);
a(end+1) = a_;
b(end+1) = b_;
c(end+1) = c_;
d(end+1) = d_;
else,
a= NaN;
b= NaN;
c= NaN;
d= NaN;
            
        end;
    end;
end;
end;

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

Asvin Kumar 2021-3-15

What's the error that you're seeing? Error stack would be helpful for anyone trying to answer this question.

Also, this code is hard to work with because we don't have the folders or files that this file is looking for. If possible, attach a MWE which reproduces your issue so the community can help you out.

请先登录，再进行评论。

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

Follow Question