How to makedisk for a three parameter weibull distribution

Question

Mustafa Vural 2020-9-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/592708-how-to-makedisk-for-a-three-parameter-weibull-distribution

编辑： Mustafa Vural 2020-9-17

Everytime I run my code with below lines, it displays "Error using makedist (line 101). Invalid parameter name: t0"

t0  = 0.5;
b   = 1:5;
T   = 1:5;
dist = makedist('weibull','a',T,'b',b,'c',t0)
[h,p,adstat,cv] = adtest(data,'Distribution',dist)

How can I integrate a three parameter weibull distribution into "makedist"?

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

请先登录，再进行评论。

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

Answer 1

Jeff Miller 2020-9-12

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/592708-how-to-makedisk-for-a-three-parameter-weibull-distribution#answer_493681

MATLAB's Weibull distribution only allows 2 parameters, so you must handle the offset separately. Try this:

t0  = 0.5; % 
b   = 1.5; % you are using vector 
T   = 1.5; % same here too 
dist = makedist('weibull','a',T,'b',b)
dataAdjusted = data - t0;
[h,p,adstat,cv] = adtest(dataAdjusted,'Distribution',dist)

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

VBBV 2020-9-13

编辑：VBBV 2020-9-13

%if true
  % code
%end
clear all;
n   = 100;
t0  = 0.5;
b = 1:5;
T = 1:5;
for v_T= 1:length(T)
  for v_b= 1:length(b)

Use the for loop as shown above or as below

%if true
  % code
%end
for v_T = [T]
 for v_b = [b]

Mustafa Vural 2020-9-13

编辑：Mustafa Vural 2020-9-17

So like this? Sry I am very noob in matlab.

clear all;
n   = 100;
t0  = 0.5;
b = 1:5;
T = 1:5;
for v_T= 1:length(T)
  for v_b= 1:length(b)
        data(:,v_b,v_T) = wblrnd(v_b,v_T, [n,1]) + t0;
  end
end
 for v_T= 1:length(T)
  for v_b= 1:length(b)
dist = makedist('weibull','a',v_T,'b',v_b)
dataAdjusted = data - t0;
[h,p,adstat,cv] = adtest(dataAdjusted,'Distribution',dist)
  end 
 end

请先登录，再进行评论。

Answer 2

VBBV 2020-9-12

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/592708-how-to-makedisk-for-a-three-parameter-weibull-distribution#answer_493660

编辑：VBBV 2020-9-12

you are using vector of values for 2 parameters and scalar value for other parameter ... try this

%if true
  % code
%end
t0  = 0.5; % 
b   = 1.5; % you are using vector 
T   = 1.5; % same here too 
dist = makedist('weibull','a',T,'b',b,'c',t0)
[h,p,adstat,cv] = adtest(data,'Distribution',dist)

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

VBBV 2020-9-12

编辑：VBBV 2020-9-12

Standard matlab weibull only allows 2 parameter dist as Jeff said But you can use mle to fit a custom distribution Here is the resource below how to do it

http://www.mathworks.com/help/stats/mle.html#btt7ilj

请先登录，再进行评论。