Re: [R] Dirichlet surface

2011-03-30 Thread Gavin Simpson
On Wed, 2011-03-30 at 11:12 -0700, Kehl Dániel wrote: > Actually, it works for the a=1 case, not for the others. It still gives > the invalid 'zlim argument' error. > I'll try to work it out maybe instead of NULL giving a c which is > dependent on the max(z). Sorry, I misread the helpfile - the

Re: [R] Dirichlet surface

2011-03-30 Thread Kehl Dániel
Actually, it works for the a=1 case, not for the others. It still gives the invalid 'zlim argument' error. I'll try to work it out maybe instead of NULL giving a c which is dependent on the max(z). Daniel 2011-03-30 10:42 keltezéssel, Kehl Dániel írta: It helped a lot indeed, thank you very m

Re: [R] Dirichlet surface

2011-03-30 Thread Kehl Dániel
It helped a lot indeed, thank you very much! Now I understand why it was a problem for persp! Daniel 2011-03-30 10:31 keltezéssel, Gavin Simpson írta: On Wed, 2011-03-30 at 09:55 -0700, Kehl Dániel wrote: Dear David, I think that is a small bug too, maybe because the function is constant? is

Re: [R] Dirichlet surface

2011-03-30 Thread Gavin Simpson
On Wed, 2011-03-30 at 09:55 -0700, Kehl Dániel wrote: > Dear David, > > I think that is a small bug too, maybe because the function is constant? > is there a nice way to put the c(0,2.1) argument optionally, only if all > the parameters are 1? > Should I post the problem somewhere else (developer

Re: [R] Dirichlet surface

2011-03-30 Thread Kehl Dániel
Dear David, I think that is a small bug too, maybe because the function is constant? is there a nice way to put the c(0,2.1) argument optionally, only if all the parameters are 1? Should I post the problem somewhere else (developers maybe?) thanks: Daniel 2011-03-30 04:42 keltezéssel, David W

Re: [R] Dirichlet surface

2011-03-30 Thread David Winsemius
On Mar 29, 2011, at 4:45 PM, Kehl Dániel wrote: Dear list members, I want to draw surfaces of Dirichlet distributions with different parameter settings. My code is the following: # a1 <- a2 <- a3 <- 2 #a2 <- .5 #a3 <- .5 x1 <- x2 <- seq(0.01, .99, by=.01) f <- function(x1, x2){ term1

[R] Dirichlet surface

2011-03-29 Thread Kehl Dániel
Dear list members, I want to draw surfaces of Dirichlet distributions with different parameter settings. My code is the following: # a1 <- a2 <- a3 <- 2 #a2 <- .5 #a3 <- .5 x1 <- x2 <- seq(0.01, .99, by=.01) f <- function(x1, x2){ term1 <- gamma(a1+a2+a3)/(gamma(a1)*gamma(a2)*gamma(a3))