You need to use p[IND] <- .C(stats:::C_pkstwo, length(x[IND]), p = x[IND], instead of p[IND] <- .C(C_pkstwo, length(x[IND]), p = x[IND], -- View this message in context: http://r.789695.n4.nabble.com/Kolmogorov-distribution-tp4639751p4644945.html Sent from the R help mailing list archive at Nabble.com.