How to eliminate NaN in pearsrnd function in Matlab

Hi,

I have a data set which has

mu=3.82
sigma=1.023
skew=-1.021 
kurt=1.284. 

When I put them in r=pearsrnd(mu,sigma,skew,kurt,39,1) function, I am getting NaN value for 'r'. Could anyone can please help me to eliminate NaN in here?

[Thamasha Samarasinghe]

Since the above combination gives me the NaN value I tried to solve that as follows.

mu=3.82;

sigma=1.023;

x=normrnd(3.82,1.023,39,1);

skew=skewness(x);

kurt=kurtosis(x);

r=pearsrnd(mu,sigma,skew,kurt,39,1);

Could someone please tell me that method is correct or not?

Answer 1

[Shivam Singh]

Some combinations of moments are not valid; in particular, the kurtosis must be greater than the square of the skewness plus 1. For more information, please refer to the "pearsrnd" function .

[Thamasha Samarasinghe]

Thank you. Then do you know any other function that i can use for this combination instead of pearsrnd function ?