Full_Name: Majid Sarmad Version: 2.8.1 OS: Linux / Windows Submission from: (NULL) (194.225.128.135)
With thanks to Alberto Viglione, in HW.tests function of homtest package, there is the following line V2 <- (sum(ni * ((ti - tauReg)^2 + (t3i - tau3Reg)^2))/sum(ni) )^0.5 which is a mistyping and leads to a miscomputation. It must be V2 <- sum(ni * ((ti - tauReg)^2 + (t3i - tau3Reg)^2) ^0.5) /sum(ni) as it is in help file of the function: V2 = sum[i from 1 to k] ni {(t^(i) - t^R)^2 + (t3^(i) - t3^R)^2}^(1/2) / sum[i from 1 to k] ni Similarly, in V2s[i] <- (sum(ni * ((ti.sim - tauReg.sim)^2 + (t3i.sim - tau3Reg.sim)^2))/sum(ni))^0.5 ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel