386314_en-US

liliconcom

How is Std Error for Beta Binomial Distribution Analysis calculated? - differences between R and JMP

Hi All!

 

We are calculating beta binomial distribution parameters in R and wanted to compare the results and the std error.

Typically, the values for location and dispersion are up to 1% higher in JMP and there is also some discrepancy in the std error.

 

We would like to match here our R-code. 

May anybody provide a hint how the standard error for mu and rho are calculated? Fisher Matrix?

 

Thanks

AnalyticsRe: How is Std Error for Beta Binomial Distribution Analysis calculated? - differences between R and JMP

Thank you very much. This is what I thought.

Re: How is Std Error for Beta Binomial Distribution Analysis calculated? - differences between R and JMP

JMP uses the inverse of the information matrix to estimate the covariance matrix (and standard errors). https://en.wikipedia.org/wiki/Fisher_information

Re: How is Std Error for Beta Binomial Distribution Analysis calculated? - differences between R and JMP

No, that quantity is the variance of the random variable when assuming a beta binomial distribution.

Re: How is Std Error for Beta Binomial Distribution Analysis calculated? - differences between R and JMP

I found this in your documentation:

Var(x) = np(1-p)[1+(n-1)δ]

Is this the Std Error shown in the table?

 

 

Re: How is Std Error for Beta Binomial Distribution Analysis calculated? - differences between R and JMP

Have you seen the JMP Help about this topic? Is it possible that the R function or model uses a different parameterization?