I fit my data to gamma distribution using the gamfit function of MATLAB.

The gamfit function returns MLEs and 95% percent confidence intervals as follows:

• a = 2; b = 4;data = gamrnd(a,b,100,1);

• [p,ci] = gamfit(data)

• p =

• 2.1990 3.7426

• ci =

• 1.6840 2.8298

• 2.7141 4.6554

I am interested in finding the standard error of the coefficients rather than the confidence interval.

I want to know how to compute the standard errors based on the confidence interval??

Moreover, likelihood ratio hypothesis test were used to test hypotheses such as p(1)=0.5, but I don’t know how to obtaine and input the NullLLF and BaseLLF values

Syntax likelihood ratio hypothesis test

[H,pValue,Ratio,CriticalValue]=lratiotest(BaseLLF,NullLLF,DoF,Alpha)