Re: [R] package names in Latex

When putting R package names in Latex, do we use \it or \em, please? As far as I know, the only difference between \it and \em is that \em words are unitalicized in an italics portion of the text, which is desirable in most cases. However, I am not sure that italics is the most common solu

Re: [R] Is profile.mle flexible enough?

Many thanks to Ben and Peter for their help. Even if I am not really happy with it, I used the "environment" solution and got the results I expected. Issue #2 is stickier. I think I must say that the idelology of mle() is that the user passes a likelihood function. If the likelihood function d

[R] Is profile.mle flexible enough?

Hi the list, I am experiencing several issues with profile.mle (and consequently with confint.mle) (stat4 version 2.9.2), and I have to spend a lot of time to find workarounds to what looks like interface bugs. I would be glad to get feedback from experienced users to know if I am really askin

Re: [R] log likelihood and deviance

Hi, If we know the residual of the model, how could we calculate the Log likelihood? It depends on the model (lm? glm? nls?). Why not using directly the logLik function? x <- rnorm(100, 10) y <- rnorm(100, 10) model1 <- lm(y ~ x) logLik(model1) model2 <- glm(y ~ x, family=gaussian) logLik(mode

[R] [R-pkgs] New package: noia

own. Locating the genes requires a QTL detection procedure, for which packages already exist (see Rqtl for instance). Much more information can be found in scientific publications. Two are particularly relevant: * This package and its functions are described in a recent paper: Arnaud L

[R] vectorized mle / optim

Hi the list, I would need some advice on something that looks like a FAQ: the possibility of providing vectors to optim() function. Here is a stupid and short example summarizing the problem: example 1 8< -- library(stats4) data