Re: [R] Question about KLdiv and large datasets

2010-07-18 Thread Peter Ehlers
I just answered this but realize that I did so off-list. So, for completeness, here's what I said: I think I see the problem. From ?KLdiv, you're getting the modeltools help page. What you need is the flexmix help page for KLdiv. Just get to the flexmix index page (you can do ?flexmix and then cl

Re: [R] Question about KLdiv and large datasets

2010-07-18 Thread Ralf B
Is the 'eps' argument part of KLdiv (was not able to find that in the help pages) or part of a general environment (such as the graphics parameters 'par' ) ? I am asking so that I can read about it what it actually does to resolve the question you already raised about its reliability... Ralf On F

Re: [R] Question about KLdiv and large datasets

2010-07-16 Thread Peter Ehlers
On 2010-07-16 7:56, Ralf B wrote: Hi all, when running KL on a small data set, everything is fine: require("flexmix") n<- 20 a<- rnorm(n) b<- rnorm(n) mydata<- cbind(a,b) KLdiv(mydata) however, when this dataset increases require("flexmix") n<- 1000 a<- rnorm(n) b<- rnorm(n) mydata<- cbin

[R] Question about KLdiv and large datasets

2010-07-16 Thread Ralf B
Hi all, when running KL on a small data set, everything is fine: require("flexmix") n <- 20 a <- rnorm(n) b <- rnorm(n) mydata <- cbind(a,b) KLdiv(mydata) however, when this dataset increases require("flexmix") n <- 1000 a <- rnorm(n) b <- rnorm(n) mydata <- cbind(a,b) KLdiv(mydata) KL se