[Rd] Lattice: Drawing a single contour line with a negative value fails

2009-08-19 Thread Thorn Thaler
Hi everybody, I want to add a single contourline to a levelplot. While everything works fine if the value at which the line should be drawn is positive, there is an error if the value is negative: library(lattice) my.panel <- function(..., at, contour=FALSE, labels=NULL) { panel.levelplot(

[Rd] Extended Quasi Likelihood Library

2009-07-13 Thread Thorn Thaler
Hi everybody, for my master thesis I'm currently working on a package providing a possibility to estimate the parameter(s) of a family of variance functions in a GLM. My approach is based on the EQL function (cf Nelder & Pregibon: An extended quasi-likelihood function. Biometrika 74 (1987), 2

Re: [Rd] Speed up code, profiling, optimization, lapply vs. loops

2009-07-07 Thread Thorn Thaler
Thanks for ur good advices. I've run Rprof and figured out that the bottleneck lies within the deviance function, which uses "integrate" to determine the deviance numerically. Since I'm not aware of a closed form solution for the integral \int 1/(t^k*(1-t)^l) dt, \forall k,l \in R , I've to

[Rd] Speed up code, profiling, optimization, lapply vs. loops

2009-07-06 Thread Thorn Thaler
High everybody, currently I'm writinig a package that, for a given family of variance functions depending on a parameter theta, say, computes the extended quasi likelihood (eql) function for different values of theta. The computation involves a couple of calls of the 'glm' routine. What I'm

Re: [Rd] Rd: pdf manual: package information not on first page

2009-06-26 Thread Thorn Thaler
Duncan Murdoch schrieb: How did you produce the manual? When I do R CMD Rd2dvi on a package, I get the package man page first. There are other ways to produce the manual, and in some cases we may not control the sort order. Well, I see. There is a misunderstanding regarding the term "package

Re: [Rd] Rd: pdf manual: package information not on first page

2009-06-26 Thread Thorn Thaler
Duncan Murdoch schrieb: So my question is, if there is any possibility to get the package documentation on the first page? That's where it is supposed to be. Could you look for errors in the one where it comes last? If you can't spot anything, you could send me the tar.gz version of it and

[Rd] Rd: pdf manual: package information not on first page

2009-06-25 Thread Thorn Thaler
unction in both packages starting with a letter prior to the first letter of the package name. So my question is, if there is any possibility to get the package documentation on the first page? Best regards, Thorn Thaler __ R-devel@r-project.org m