#H0: (p,mu,sigma)=(1/3,55,5) psi0=c(1/3,55,5) psihat=MLE[1:3] diff=psihat-psi0 #Wald statistic W=t(diff)%*%solve(Vhat[1:3,1:3])%*%diff W cat("\n p-value is",1-pchisq(W,3))