ResSimulLvrRatio(Data,ResH0,Constraint,nsimul,level,eps=1e-04,eps1=1e-04)DataPairCompBradleyEstim corresponding to the result of BradleyEstim() function for T classes (H0)
Constraint=0, the sum of Bradley's scores should be equal to 1. For other values for Constraint, the product of Bradley's scores should be equal to 1 (default Constraint=0).
level=0.05).eps=1e-04).
eps1=1e-04).
LvrRatio with the following components:
nsimul data set is done. We perform estimation of the different parameters for the number of classes defined in the object ResH0 of class BradleyEstim (corresponding to
the null hymothesis) and for one more class corresponding to the alternative hypothesis.We obtain a set of Log Likelihoods under the null and alternative hypothesis on the basis of simulated data and so of the Log Likelihood Ratio Statistic.
We replace the observed value of this statistic for the true data set. And we conclude on the acceptation or not of the null hypothesis (no differences between T and T+1 classes).
data(Cocktail)
ResCock1<-EstimBradley(Cocktail,Constraint=0,Tcla=1,eps=1e-04,eps1=1e-04,TestPi=TRUE)
Res_LvrRatio1<-ResSimulLvrRatio(Cocktail,ResCock1,0,3,level=0.05,eps=0.001,eps1=0.001)
getSimu(Res_LvrRatio1)
getTest(Res_LvrRatio1)
Run the code above in your browser using DataLab