[R] why my code is not computing significance level and power correctly

I am trying to compare Levene,bartlett, count5 tests for a random normal sample for the significance level and power, the problem seems bug free but with obvious wrong result, anyone can help, I send the R code as attachment. Thanks http://www.nabble.com/file/p16992744/758%2Bproject 758+project -

[R] help me to debug this part of code?

I am trying to solve the integration equation, for different values of K from 4 to 25, the integration is with respect to u, Here is the equation:gamma(k/2) / ( sqrt(k-1)*gamma((k-1)/2) ) * integrate(f= (1+u^2/k-1)^(-k/2), lower=0, upper= sqrt(a^2*k/(k+1-a^2)) ) = the similar expression as t

[R] use bootstrap to get MLE of lamda from the exp model

the data are 3,5,7,18,43,85,91,98,100,130,230,487, it follows an exp(lamda) model, how to use bootstrap to obtain MLE of lamda, bias, sd. Thanks a lot so far my code is obj<-boot(data,statistic,R=2000) here statistic should call a function to get an estimate of lamda I

[R] how to randomly generate a n by n positive definite matrix in R ?

thanks -- View this message in context: http://www.nabble.com/how-to-randomly-generate-a-n-by-n-positive-definite-matrix-in-R---tp15340155p15340155.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list http

[R] how to calculate chisq value in R

for example, an expression such as chisq(df=1,ncp=0) ? thanks -- View this message in context: http://www.nabble.com/how-to-calculate-chisq-value-in-R-tp15338943p15338943.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-pro