Dear all,
I get a strange error when I find polychoric correlations with the ML method,
which I have been able to reproduce using randomly-generated data.
What is wrong?
I realize that the data that I generated randomly is a bit strange, but it is
the only way that I duplicate the error message.
> n<-100
> test.x<-rnorm(n, mean=0, sd=1)
> test.c<-test.x + rnorm(n, mean=0, sd=.5)
> thresh.x<-c(-2.5, -1, -.5, .5, 1000)
> thresh.c<-c(-1, 1, 2, 3, 1000)
>
> discrete.x<-discrete.c<-vector(length=n)
>
> for (i in 1:n) {
+ discrete.x[i]<-which.min(thresh.x < test.x[i] )
+ discrete.c[i]<-which.min(thresh.c < test.c[i] )
+ }
> pc<-polychor(discrete.x, discrete.c, std.err=T, ML=T)
Error in optim(c(optimise(f, interval = c(-1, 1))$minimum, rc, cc), f, :
non-finite finite-difference value [1]
In addition: There were 50 or more warnings (use warnings() to see the first 50)
> print(pc)
Error in print(pc) : object "pc" not found
> warnings()
Warning messages:
1: NaNs produced in: log(x)
2: NA/Inf replaced by maximum positive value
3: NaNs produced in: log(x)
---
Thanks,
Janet
--------------------
This email message is for the sole use of the intended recip...{{dropped}}
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
