similar to: formula method with "special" characters

Hi, Currently the `gaussian' family's initialization code signals an error if any response data are zero or negative and a log link is used. Given that zero or negative response data are perfectly legitimate under the GLM fitted using `gaussian("log")', this seems a bit unsatisfactory. Might it be worth changing it? The current offending code from `gaussian' is:

G'day all, I had a look at the GLM code of R (1.4.1) and I believe that there are problems with the function "glm.fit" that may bite in rare circumstances. Note, I have no data set with which I ran into trouble. This report is solely based on having a look at the code. Below I append a listing of the glm.fit function as produced by my system. I have added line numbers so that I

Is there an easy way to compile a packages's fortran so that it has the information for ddd or gdb to follow the lines of the code? I seem to remember seeing something about this, but I cannot find it in Writing R Extensions. (I see about Valgrind, debugging a segfault, and some other tricks, but what I really want is to follow the fortran logic.) Thanks for any pointers, Paul

I encountered a segfault running glm() and wonder if it could have something to do with the way memory is handled in a dual core system (which I just set up). I'm running R-base-2.4.0-1, installed from the SuSE 10.1 x86_64 rpm (obtained from CRAN). (My processor is an AMD Athlon 64 x2 4800+). The error and traceback are *** caught segfault *** address 0x8001326f2b, cause 'memory not

With R-rc_2010-10-11_r53297.tar.gz I am getting ... ** building package indices ... ** testing if installed package can be loaded Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared object '/home/mfa/gilp/toolchain/R/src/R-rc/library/cluster/libs/cluster.so': /home/mfa/gilp/toolchain/R/src/R-rc/library/cluster/libs/cluster.so: undefined symbol: cl_daisy_

Is anyone else having trouble connecting to R-forge svn? (Perhaps I missed an outage announcement, or is it bad weather?) Paul ==================================================================================== La version française suit le texte anglais. ------------------------------------------------------------------------------------ This email may contain privileged and/or confidential

In one package (pk1) I define some S3 generics and methods. In pk2 I extend the methods to some S4 objects following the description in ?Methods, with setMethod() and setGeneric() either explicit or implicit (tried both). I export the S3 generic and methods from pk1 NAMESPACE and the S4 from pk2. This works, except when I try to use a function defined in pk1, Within the function the S4 method

I am experimenting with putting vignettes in a vignettes/ directory as suggested in the new R release news. If I understand correctly, it is now suggested that the vignette source (.Stex in my case) should be in the vignettes/ directory and the pdf should be put in the inst/doc/ directory, for the case when it is not being built by R CMD build. I am getting an error (R-2.14.0): (cd inst/doc ;

I have recently upgraded samba to 3.6.23-9.el5_11 and have noticed all new users created can not be seen with using the samba config tool within CentOS or within the samba webmin module. I can manually log in and look the account information within the smbpasswd file and can query it using pdbedit. I need to be able to edit this information from within the webmin module as there are a number

Is CRAN having trouble or is it just me? (The web interface is very slow and install.packages() is timing out and giving me "service not available"). Paul ==================================================================================== La version française suit le texte anglais. ------------------------------------------------------------------------------------ This email

I am trying to extend an S3 method to work with an S4 object as well as the S3 objects it works with, but UseMethod does not seem to recognize the S4 class and dispatches to the default method. Is this to be expected or should I be looking for an error in my code? If it is not an error in my code, is there an easy way to do this, or do I have to convert the generic to S4 and then make those

I am trying to check my packages with R-rc and latex is failing to find utf8.def. I presume my latex installation is defective or too old, or is this file distributed with R and I just don't have a path set correctly? Paul ______________ * checking PDF version of manual ... WARNING LaTeX errors when creating PDF version. This typically indicates Rd problems. LaTeX errors found: !

Bug summary: glm() causes a segfault if the argument 'data' is a data frame with more than 16384 rows. Bug demonstration: -------input --------------- N <- 16400 df <- data.frame(x=runif(N, min=1,max=2),y=rpois(N, 2)) glm(y ~ x, family=poisson, data=df) ------ output --------------- *** caught segfault *** address (nil),

Details in documentation: "All of ?weights?, ?subset?, ?offset?, ?etastart? and ?mustart? are evaluated in the same way as variables in ?formula?, that is first in ?data? and then in the environment of ?formula?." In fact, `data` is usually not an environment, and I have not seen arguments evaluated in `environment(formula)` when `data` is provided. (Information in `environment(formula)`

I am having the following problem building R-beta_2009-10-12_r50045.tar.gz on Ubuntu 9.04 - the Jaunty Jackalope Paul ________ ... make[2]: Entering directory `/home/paul/toolchain/R/R-beta/src/library/Recommended' begin installing recommended package boot gzip: invalid option -- 'x' Try `gzip --help' for more information. tar: This does not look like a tar archive

Is it possible in R to call a fortran routine that sets variables in a common block and expect the values to persist when a call is made from R to a second routine that uses the common block? If not (as I suspect), is it possible to use a common block in a group of routines that are used together, for as long the routines do not return to R? Paul

Hi, We're running Samba 3.0.26a on AIX 5.3 using security=ads. Users get access denied when ADS user token size gets over an unknow limit even though the ACL settings on JFS2 files/directories should give them access. To demonstrate and based on "id userid" output, user1 AD user token size is over 9k in size since it belongs to about 400 AD groups while user2 AD user token size

I am trying to understand the missing vs default value for an argument to an S4 method. (I'm not sure if this is a bug or my confusion.) In S3 I can specify a default value for an argument and then both check if the argument was missing in the call, and use it because it takes on the default value. I do not seem to be able to do that in S4 as shown by this simplified example:

Hi everyone: I am just wondering how can I use weights with MCMClogit function (in MCMCpack package). For example, in case of glm function as given below, there is weights option in the arguments. Aparently there is no option of using weights in MCMClogit. glm(formula, family = gaussian, data, weights, subset, na.action, start = NULL, etastart, mustart, offset, control =

Dear users, I'm working on binary classification problem using Support Vector Machines (SVM). My objective is to train a series of SVM models on a grid of hyperparameters and then select those that maximize the AUC based on an independent validation sample. My attempted code is shown below. It runs well on "small" data sets but when I use it on a slightly larger sample (e.g., my