similar to: Likelihood ratio test for comparing non-nested cox models

When running a = runif(10) class(a) = "foo" Math.foo = function(x, ...) { NextMethod(.Generic) } signif(a, 3) cumsum(a) I don't understand why cumsum strips the class, but signif does not. Both claim in the documentation that "These are generic functions: methods can be defined for them individually or via the ?Math? group generic." -- Edzer Pebesma Institute for

Dear All, I am using the function vuong from pscl package to compare 2 non nested models NB1 (negative binomial I ) and Zero-inflated model. NB1 <-  glm(, , family = quasipoisson), it is an object of class: "glm" "lm" zinb <- zeroinfl( dist = "negbin") is an object of class: "zeroinfl"   when applying vuong function I get the following: vuong(NB1,

Dear R folks While preparing figures of 'envfit' plots with vegan for publication, I ran into a layout problem, that I found no solution for in the literature or the help archive: The labels of the vectors that indicate correlations of environental variables sometimes overlap with each other or with other content of the figure. Hence, I would like to rearrange them. Is there a

Dear All, I am using the function vuong of the package pscl to compare 2 non nested glm models with a numeric response. I did the following m1<-glm(y ~x ,data=xxx) m2<-glm(y ~z , data=xxx) When calling the vuong function I get the following message: > vuong(m1,m2) Error in predprob.glm(m1) : your object of class glm is unsupported by predprob.glmyour object of class lm is unsupported

I want to compare the poisson and the zero-inflated poisson distribution on describing the data. So, after using the GLM and the ZEROINFL function, I used the Voung function to compare them. Here is my code: library(pscl) glm1 <- glm(nmer9_1[, 1] ~ 1, family = poisson) zip <- zeroinfl(nmer9_1[, 1] ~ 1) vuong(glm1, zip) However, R returns the following error message: "Error: cannot

Dear All, I am using the function vuong of the package pscl to compare 2 non nested glm models with a numeric response. I did the following m1<-glm(y ~x ,data=xxx) m2<-glm(y ~z , data=xxx) When calling the vuong function I get the following message: > vuong(m1,m2) Error in predprob.glm(m1) : your object of class glm is unsupported by predprob.glmyour object of class lm is unsupported

Dear All, I want to plot a heat map with annotated columns. Both functions heatmap (stats) and heatmap.2 (gplots) can plot a horizontal side bar that can be used to visualize a categorical variable. In addition to a categorical variable, I would like to visualize a continuous variable. This could be done by small bars, a curve or simply numbers above the columns. (The Sample names are already

Hi! I noticed that `ssh -D /tmp/socket $myhost` will actually provide a socks server listing on the unix domain socket /tmp/socket (this does not appear to be documented, though it is semi-documented for -L and -R). Do you know of any "socksification" script to go along with it, in the spirit of tsocks of socksify (from dante)? Those do not seem to support socks servers running on

Am Di, 20. Feb 2018, 23:13:16 -0800 schrieb Dan Kaminsky: > Date: Tue, 20 Feb 2018 23:13:16 -0800 > From: Dan Kaminsky <dan at doxpara.com> > To: J? Fahlke <jorrit at jorrit.de> > Cc: openssh-unix-dev at mindrot.org > Subject: Re: Is there socksify script for dynamics forwardings to unix > domain sockets? > > Whoa. That's pretty cool. > >

Dear all, I want to build a model in R based on animal collection data, that look like the following Nr Village District Site Survey Species Count 1 AX A F Dry B 0 2 AY A V Wet A 5 3 BX B F Wet B 1 4 BY B V Dry B 0 Each data point shows one collection unit in a certain Village, District, Site, and Survey for a certain Species. 'Count' is the number of animals collected in that

Dear list, I'm using the windows version of R. When plotting a curve or a line for time series with annual data , e. g. GDP growth 1991-2003, the line seems to exist of a lot of smaller lines. Printing the results the curves and lines seems to be "unclean" (because of using small resolution bitmaps?). Comparing the result of R with the same results of Excel the lines in excel seems

Hi Dominik, from the output, the master process could not "listen" on the port where it expects a connection from the worker. We need to find out why. I'd recommend first to create a minimal reproducible example (and one that does not use future, only parallel, and a minimal number of threads, ideally just 2). Then I'd recommend to check if the problem still exists with

Hi all, The call parallel::makeCluster(1L) hangs infinitely on my MacOS machine which seems to be already reported by some people (e.g., https://stat.ethz.ch/pipermail/r-devel/2018-February/075565.html). However, the solutions posted on SO, GH or R-devel do not work in my case. So far, I unsuccessfully tested ? 1. Couple of reboots 2. Adding 192.0.0.1 to /etc/hosts 3. Using R.app

List, In order to get rid of some old, unreadable S3 code in package sp, I'm trying to rewrite things using S4 methods. Somewhere I fail, and I cannot sort out why. In order to isolate the problem, I created two functions, doNothing<- and dosth, and both should do nothing. The issue is that in most cases they do nothing, but in some cases dosth(obj) changes the class of obj and breaks with

A while ago our institute was approached by a large company, say XXX, who asked whether we would like to become an Xxx development center, as we are active in the same area. After some negotation we did, believing there will be mutual benefits. Some time later, after putting up the Xxx development center logo on the institute home page, my feeling grew that we are also, and probably to a larger

Hi Thomas, thanks for your reply (and thanks for your patience...). I am now using the following minimal reprex: > library(parallel) > cl <- makeCluster(2L) I freshly started the machine and did not open any other app. Just R.app (3.6.1). After executing the second line of code, R seems to hang infinitely and does not respond. The R process itself uses almost no CPU. Unfortunately,

Hello everybody, for serveral hours we had "Had Sum mismatch" errors for the following packages from the mirror of U of M?nster -- r-project works fine as a temp. fallback: r-cran-cluster 2.0.6-2~stretchcran.0 r-cran-foreign 0.8.69-1~stretchcran.0 r-cran-kernsmooth 2.23-15-3~stretchcran.0 r-cran-class 7.3-14-2~stretchcran.0 r-cran-nnet 7.3-12-2~stretchcran.0 r-cran-spatial

Dear list, are there any possibilities to fill a polygon with a point pattern or with a symbol pattern like '+' oder '-' instead of shading lines? Thanks in advance Dr. Michael Wolf Bezirksregierung M??nster Dezernat 61 Domplatz 1-3 48161 M??nster Tel.: ++ 49 (02 51) / 4 11 - 17 95 Fax.: ++ 49 (02 51) / 4 11 - 8 17 95 E-Mail: michael.wolf at bezreg-muenster.nrw.de

I would like to test in R what regression fits my data best. My dependent variable is a count, and has a lot of zeros. And I would need some help to determine what model and family to use (poisson or quasipoisson, or zero-inflated poisson regression), and how to test the assumptions. 1) Poisson Regression: as far as I understand, the strong assumption is that dependent variable mean = variance.

Hi List, I just updated an old sendmail installation to a modern postfix/dovecot system with mysql backend. Everything works fine but SMTP-AUTH fails with some clients. The problem is, that with PLAIN (old) clients send the base64 coded form: USERNAME\0USERNAME\0PASSWORD which, as i understand dovecot interprets as MASTERUSER\0USERNAME\PASSWORD. New clients send \0USERNAME\0PASSWORD which works