How can I define a foor loop in which there is a bootci function?

1 回表示 (過去 30 日間)
Ashkan Rigi
Ashkan Rigi 2021 年 10 月 30 日
コメント済み: Ashkan Rigi 2021 年 10 月 30 日
Hello. Firstly I estimate cubic splines coefficients and then I would like to implement bootstrap function on each interval of cubic splines. This is my code and I got the error " Not enough input arguments ".
t=0:.1:pi;
y=sin(t);
N=length(t); %number of points
n=N-1; %number of subintervals
h=(t(N)-t(1))/n; %step size
% h=2;
Trid=diag(4*ones(1,n-1))+diag(ones(1,n-2),1)+diag(ones(1,n-2),-1);
for i=1:n-1
z(i)=6/h^2*(y(i+2)-2*y(i+1)+y(i));
end
z=z';
w=inv(Trid)*z;
sigma=[0;w;0];
for i=1:n
d(i)=y(i);
b(i)=sigma(i)/2;
a(i)=(sigma(i+1)-sigma(i))/(6*h);
c(i)=(y(i+1)-y(i))/h-h/6*(2*sigma(i)+sigma(i+1));
end
r=1; %number of subintervals
hh=h/r; %step size of subintervals
x=t(1):hh:t(N);
for i=1:n
for j=r*(i-1)+1:r*i
s(j)=a(i)*(x(j)-t(i))^3+b(i)*(x(j)-t(i))^2+c(i)*(x(j)-t(i))+d(i);
end
end
s(r*n+1)=y(N);
for i=1:n
s_i{i}=@(x) a(i)*(x-t(i)).^3+b(i)*(x-t(i)).^2+c(i)*(x-t(i))+d(i);
end
k=1;
for i=1:n
B = [-.161 0 1 0];
beta0(i,:)=[a(i) b(i) c(i) d(i)];
x=t';
s_i{i}=@(x) a(i)*(x-t(i)).^3+b(i)*(x-t(i)).^2+c(i)*(x-t(i))+d(i);
modelfunc(i)=s_i{i};
beta(i) = @(predictor,response)nlinfit(predictor,response,modelfunc(i),beta0(i,:))
c(:,k) = bootci(1000,beta(i),t',y,'Type','norm');
k=k+2;
end

採用された回答

Alberto Cuadra Lara
Alberto Cuadra Lara 2021 年 10 月 30 日
I have not checked all the code, but that error is because you are not declaring the inputs well. Use { } inside bottci function.
bootci(1000,{beta(i),t',y},'Type','norm');
  1 件のコメント
Ashkan Rigi
Ashkan Rigi 2021 年 10 月 30 日
I did the justification that you mentioned but I took again the error text.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeResampling Techniques についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by