search for: mnsq

...is code in Splus, the results are the same. Also, if I set the size argument to be the same (or eliminate this line altogether), I get reproducible results.=20=20 =20 I'm using R 2.2.0 =20 Thanks, -Lori =20 test.fun<-function(nsamp, mu,q, seed, nsim=3D2){ set.seed(seed)=20=20=20=20 mnsq<-rep(NA,nsim) for (i in 1:nsim){ y<-rnorm(nsamp,mean=3Dmu) yboot<-sample(y, size=3Dq, replace=3DT) mnsq[i]<-(mean(y)^2) } return(mnsq) } =20 test.fun(50,0,30,12345) test.fun(50,0,10,12345) =20 Lori E. Dodd, Ph.D. National Cancer Institute Biometric Research Branch 6130 Ex...

...,"x5","x6","x7") ltm$loc<-rowSums(ltm)/3 ltm<-round(ltm[order(ltm$loc,decreasing=TRUE),],2) ltm #I get the following order of item parameters:x2,x3,x4,x7,x1,x5,x6 And this is what I got using winsteps.          |           |   INFIT  |  OUTFIT  |  item|   loc   |MNSQ  ZSTD|MNSQ  ZSTD| ------|---------+----------+----------+  x7   |   .49   |     .62  -1.4|1.05    .3|  x2   |   .30       |1.02    .2|1.27    .9|  x3   |   .24   |    1.14    .6|1.05    .3|  x4   |   .00   |    .88   -.4| .80   -.6|  x1   |  -.11       | .91   -.3| .88   -.3|  x5   |  -.28   |    1...