Re: [R] Design and analysis of mixture experiments

2008-04-21 Thread k . jewell
work in R (or even contribute a package!) if only I had the necessary skills and time. Never mind, I have other tools which will do the job. It's just a little iconoclastic finding needs that R can't meet :-O Regards to all. Keith Jewell mailto:[EMAIL PROTECTED] telephone (direct) +4

[R] Design and analysis of mixture experiments

2008-04-17 Thread k . jewell
such experiments, but I haven't been able to find them (I have looked quite hard!). Can anyone point me in the right direction? Thanks in advance, Keith Jewell mailto:[EMAIL PROTECTED] telephone (direct) +44 (0)1386 842055

[R] Solved (??) Behaviour of integrate (was 'Poisson-lognormal probab ility calculations')

2008-02-18 Thread k . jewell
integration as shown here: -- function (x, meanlog = 0, sdlog = 1, sdlim=6, rel.tol = .Machine$double.eps^0.5,...) { #+ # K. Jewell Feb 2008 # Based on VGAM::dpolono, modified to cover wider range of x # Some simplification, but major change is to avoid

[R] Behaviour of integrate (was 'Poisson-lognormal probability calcul ations')

2008-02-15 Thread k . jewell
Hi again, Adding further information to my own query, this function gets to the core of the problem, which I think lies in the behaviour of 'integrate'. - function (x, meanlog = 0, sdlog = 1, ...) { require(stats) integrand <- function(t, x, meanlog, sdl

[R] Poisson-lognormal probability calculations

2008-02-15 Thread k . jewell
Hi, just for the record, although I don't think it's relevant (!) - > sessionInfo() R version 2.6.0 (2007-10-03) i386-pc-mingw32 locale: LC_COLLATE=English_United Kingdom.1252;LC_CTYPE=English_United Kingdom.1252;LC_MONETARY=English_United Kingdom.1252;LC_NUM