similar to: MASS glm.nb: Offset fails

Hello, Has anyone written R functions for applying self-controlled case series methods (http://statistics.open.ac.uk/sccs/). In fact only thing needed is to modify glm function to allow absorption of effect. Eg. in Poisson model individual effect is used as factor, but it is considered as nuisance term where parameter estimates are not needed. Could anyone point how absorbing individual

Dear R people, I am not sure if this is the correct place to tell about problems in evolving programmes, but it seems that the `summary.negbin' function of the excellent `MASS' library is now broken, and gives the following error message: > summary(hm) Error in summary.negbin(hm) : subscript out of bounds `summary.negbin' calls `summary.glm' which seems to work and give the

I upgraded R-base 1.2.0 today, and then I intended to give a finishing touch to some R functions I should send out today. To my surprise, they failed in R 1.2 although they still worked fine in R 1.0.1 (I checked in another machine). The error message I got was like this: > optci(trar.g, water) Error in unique(c("AsIs", class(x))) : Object "x" not found The function

Dear Mr Moderator, please let me through. I want to reply to my own thread. I once subscribed to this list, but then my subscription was not accepted. Now to the business: The final failure came from missing .install.macbinary() function. I grepped R-patched sourcesand the only instance of .install.macbinary() was the call to the function in R- patched/src/library/utils/R/packages2.R (stupid

The R-exts manual says about 'Suggests' field in package DESCRIPTION: "The optional `Suggests' field uses the same syntax as `Depends' and lists packages that are not necessarily needed." However, this seems to be a suggestion you cannot refuse. If you suggest packages: (a line from DESCRIPTION): Suggests: MASS, ellipse, rgl, mgcv, akima, lattice This is what happens:

Message 2 of today: it works now. After re-installing R.app from the same R-2.1.0.dmg file, and fixInNamespace'ing(*) install.packages(), I managed to update 23 outdated packages from CRAN binaries for MacOS X. What I did was to: 1. change call to .install.macbinaries() into call to install.binaries() in install.packages(). install.binaries() is a function defined utils/R/aqua/GUI.R. 2. I

Canberra distance is defined in function `dist' (standard library `mva') as sum(|x_i - y_i| / |x_i + y_i|) Obviously this is undefined for cases where both x_i and y_i are zeros. Since double zeros are common in many data sets, this is a nuisance. In our field (from which the distance is coming), it is customary to remove double zeros: contribution to distance is zero when both x_i

Canberra distance is defined in function `dist' (standard library `mva') as sum(|x_i - y_i| / |x_i + y_i|) Obviously this is undefined for cases where both x_i and y_i are zeros. Since double zeros are common in many data sets, this is a nuisance. In our field (from which the distance is coming), it is customary to remove double zeros: contribution to distance is zero when both x_i

I find it a bit peculiar that a package skeleton created with a utils function package.skeleton() fails subsequent R CMD check. I do understand that the function is intended to produce only a skeleton that should be edited by the package author. I think that it would be justified to say that the skeleton *should* fail the test. However, I have two arguments against intentional failure: * When you

Thanks for the attention paid to my rpoblem. Please find enclosed the matrix with my dissimilarities. This is the only case in which sammon(), from the MASS package, gives me this kind of problems. Domenico > > > -----Messaggio originale----- > > Da: Jari Oksanen [mailto:jarioksa at sun3.oulu.fi] > > Inviato: mercoled?? 20 aprile 2005 11.53 > > A: Domenico Cozzetto >

rolf at math.unb.ca said: > If I say something like > if(x .gt. 42.d0) stop > then indeed everything stops, i.e. R falls over. I'd ***like*** to be > able to print out an informative error message (which I guess could be > done - In Fortran: subroutine foo(..., ier) integer ier ier=0 ... if (x .gt. 42.d0) then ier=1 return endif

After some futile searches, I decided to ask the list to see if any of the sages out there would have an answer: I have a function I wrote a few years ago in S, which calls glim numerous times. I'd like to port it to R, but glm works differently from glim, which takes as part of its input an X design matrix. I probably could write a function to convert glim to glm, but hope this

jari.oksanen at oulu.fi wrote: > Full_Name: Jari Oksanen > Version: 2.6.2 RC (2008-02-07 r44369) > OS: Linux > Submission from: (NULL) (130.231.102.145) > > > Even if function xspline() is called with argument draw=3DFALSE, it req= uires a > graphics device (that it won't use since it was draw=3DFALSE). I run in= to this > because I intended to use xspline within a

Hi all I try to make a weibull survival analysis on R. I know make this on GLIM, and now I try to make the GLIM exercice GLEX8 on R to learning and compare the test. The variables are: time censor group bodymass In GLIM I make: $calc %s=1 $ to fit weibull rather than exponential $input %pcl weibull $ $macro model group*bodymass $endmac$ $use weibull t w %s $ Then, GLIM estimate an alpha for the

Full_Name: Jari Oksanen Version: R 2.1.0 OS: MacOS 10.3.9 Submission from: (NULL) (130.231.102.145) For various reasons (which need not be expanded here) I have tried to update my long neglected R in MacOS X using handy command line tool update.packages() using readily available binaries of contributed packages at CRAN. However, this fails with message saying that packages xxxx_*_tar.gz is not

Dear developers, I just noticed that step() function currently prints the current model using message(), but the resulting model using print(). The relevant commands within the step() body are: if (trace) message("Start: AIC=", format(round(bAIC, 2)), "\n", cut.string(deparse(as.vector(formula(fit)))), "\n") (with example() output:) Start: AIC=190.69

Dear R folks, One of the fantastic new tools in R is `Sweave'. I have tested it so much that I know it works and produces fine documentation, and with (GNU) Emacs/ESS it is nice to work with, too. I started to have a look at including some Swoven (is that a strong verb?) documentation with my R package, but it seems that there is no model to copy among those packages that I have installed in

> On 25 Mar 2016, at 10:08 , Jari Oksanen <jari.oksanen at oulu.fi> wrote: > >> >> On 25 Mar 2016, at 10:41 am, peter dalgaard <pdalgd at gmail.com> wrote: >> >> As I see it, the display showing the first p << n PCs adding up to 100% of the variance is plainly wrong. >> >> I suspect it comes about via a mental short-circuit: If we

plummer at iarc.fr said: > On 24-May-00 Prof Brian D Ripley wrote: > On Wed, 24 May 2000, Christian Posse wrote: > >> >> I just encountered a problem with the ts package: >> >> > > library(ts) >> Error in dyn.load(x, as.logical(local), > as.logical(now)) : >> unable to load shared library >> "/usr/ >

Hi all, In my continuous transition of GLIM to R I try to make a glm with binomial errors. The data file have 3 vectors: h -> the factor that is ajusted (have 3 levels) d -> number of animais alive (the response) n -> total number of animals To test proportion of alive, make d/n. In GLIM: $yvar d$ $error binomial n$ $fit +h$ scale deviance = 25.730 (change = -9.138) at cycle 4