Hello,
Thank you for paper and your package. I think there is a little mistake in the Example.
The current code for the example is :
n <- 500 X <- cbind(rnorm(n),rnorm(n)) X.NA <- X X.NA[,1] <- ifelse(stats::runif(n)<=0.2, NA, X[,1]) pval <- PKLMtest(X)
I suggest replacing with :
n <- 500 X <- cbind(rnorm(n),rnorm(n)) X.NA <- X X.NA[,1] <- ifelse(stats::runif(n)<=0.2, NA, X[,1]) pval <- PKLMtest(X.NA)
Regards,
Adrien
Hello,
Thank you for paper and your package. I think there is a little mistake in the Example.
The current code for the example is :
n <- 500 X <- cbind(rnorm(n),rnorm(n)) X.NA <- X X.NA[,1] <- ifelse(stats::runif(n)<=0.2, NA, X[,1]) pval <- PKLMtest(X)I suggest replacing with :
n <- 500 X <- cbind(rnorm(n),rnorm(n)) X.NA <- X X.NA[,1] <- ifelse(stats::runif(n)<=0.2, NA, X[,1]) pval <- PKLMtest(X.NA)Regards,
Adrien