similar to: Percentage effects in logistic regression

Dear All, I'm trying to use waldtest to test poolability (parameter stability) between two logistic regressions. Because I need to use robust standard errors (using sandwich), I cannot use anova. anova has no problems running the test, but waldtest does, indipendently of specifying vcov or not. waldtest does not appear to see that my models are nested. H0 in my case is the the vector of

Dear All, I would like to analyse the residuals from a generalized linear mixed model (GLMM) that I estimated, with random effects, by means of the command glmmPQL, from the MASS package. It is not very clear to me what the actual residuals to analyse are (Y - Yhat): I obtain two columns of residuals, of which the first are population residuals, and the second refer to the grouping used in the

Dear All, I carry out negative binomial estimations using the glm.nb command from the MASS package. Is there a command or a simple procedure for computing marginal effects from a glm.nb fitted object? If these are the same as for a Poisson fitted object (glm), my question remains how to compute them. Thanks in advance for your help. Roberto Patuelli ******************** Roberto Patuelli, Ph.D.

Dear All, I'm looking to perform an ANOVA between two nested panel fixed effects models. I tried with anova, as well as with waldtest. anova tells me there is no method available for plm objects, while waldtest tells me my models are not nested. I think they are instead. The difference between the two models is that in the second I let the regression coefficients of a given variable variate

Dear All, I've been using already for a year glm.nb() from the MASS package. But today, R gave me an error message when estimating one of my usual models: > depEsf.nb <- glm.nb(depE ~ manuf00E + corps00E + lngdp00E + lngdp00sqE + > lnpop00E + indshE + scishE + mechshE + elecshE + chemshE + drugshE + > urban_dummyE + aggl_dummyE + + eE1 + eE2 + eE3 + eE4 + eE5 + eE6 + eE7 +

Dear Duncan, Thank you a lot for your feedback. I have a problem though. The branch weights counters overflow in some files and thus I get incorrect numbers. Is there any way to find a workaround for that? Is is supposed to be a known bug or is it something that needs configuration on my part? Again, thank you a lot for your reply. Best Regards, Georgios Zacharopoulos 2016-02-03 18:23

Dear Duncan, I am generating branch-weights annotated IR files as described in the documentation of LLVM, using profiling with instrumentation. http://clang.llvm.org/docs/UsersManual.html#profiling-with-instrumentation e.g. llvm-profdata merge -output=$(BENCH).profdata default.profraw > clang -S -emit-llvm -O3 -fprofile-instr-use=$(BENCH).profdata -o > bench.prof.ll bench.c The issue is

Dear All, I was wondering if it is possible to generate a regression summary (it does not matter at this stage if from an lm or for example a glm estimate) in which to obtain the joint significance of a set of regressors? Examples could be looking at the joint significance level of a polynomial, or of a set of exogenous variables of which is of interest the linear combination suggested by the

Estimados usarios de R: Tengo una base de datos madre en formato .sav de SPSS y la quiero modificar usando datos de otras base de datos .sav y otra en .csv a las que llamare hijos. No tengo problema en convertirlas en data.frame. Todos los archivos tienen en común una variable que es única. En aquellos casos que les falta un valor a una variable en el archivo madre lo relleno con el valor del

Hello everybody My collegue and I noticed a strange behaviour of R on different platforms. It's a simple computation, but results are rather different. On Windows XP: > floor(log(8,2)) [1] 3 which is what one should expect. Here's instead the result with Mac OS X (same version, 2.5.0 (2007-04-23)) > floor(log(8,2)) [1] 2 Is it a "bug" in R or in the operating

------------------------------------------------------------------------------------------------ CALL FOR PAPERS 15th International Conference on Managed Languages & Runtimes (ManLang'18) September 10-14, 2018, Linz, Austria http://ssw.jku.at/manlang18/ ------------------------------------------------------------------------------------------------ ManLang (formerly PPPJ) is a

Hi all. I just upgraded my server from 4.7-Stable to 5.0-Release by a full install. Now I hit a problem with the date.. the time goes much faster than normal.. maybe twice as normal... this morning I found the following kernel log in /var/log/messages May 26 17:23:58 jupiter kernel: Timecounter "i8254" frequency 1193182 Hz May 26 17:23:58 jupiter kernel: Timecounter "TSC"

Dear Person That Posted This Email, Please be aware that calls for papers are only on-topic for this email list if they are specifically about LLVM. For example, a call for a workshop on LLVM in Scientific Computing would be on topic whereas a general compiler call for papers (like the one below) is not. That said, I hope you receive many strong submissions. Your conference looks

We apologize for multiple receipts. -------------------------------------------------------------------------------- Eighth Real-Time Linux Workshop October 12-15, 2006 Lanzhou University - SISE Tianshui South Road 222 Lanzhou, Gansu 730000

Hi there. OS version: FreeBSD 4.7-Stable FreeBSD #3 I just added a disk to my box and after playing with fdisk and disklabel with this configuration: fdisk: **************************************************************** g c16383 h16 s63 p 1 165 0 39102336 **************************************************************** disklabel

I''m stumped... I need to return a number of rows with the top & bottom x percentage of rows removed in mysql 4.x e.g. SELECT COUNT(*) value_count, SUM(deals.value) / COUNT(*) average_value FROM deals ORDER BY deals.value Now say that returns the sum of 100 deals, what I actually need to do is return the count & sum of rows 10-90 (i.e. the top & bottom 10% rows removed)

I have been reading a paper whereby the authors took values from Sorensons dissimilarity index and values from a percentage similarity index and applies G-Testing to the table of values. This is carried out to assess the differences in spider faunas (Strattton and Uetz, 1979). I like the method but have been trying to work out what function in R to use to get the percentage similarity. I have

One set of data has censored (less-than detection limits) water chemistry concentrations for 80-100% of all observations. My initial trial-and-error attempts to apply the cenboxplot() method suggests that it has an upper limit to the percentage of censored observations. I do not see this limit in Dennis Helsel's second edition. Has anyone experience plotting censored data and can provide

Hi there, I have some bird flight height data that follows a gamma distribution. The data (x) goes from 0 to 700 meters (n=1055). The calculated parameters calculated from the fitdistr(x) are (shape = 5.1379, rate = 0.017541), and therefore the scale (1/rate) = 57.00929. I would like to calculate the percentage of the function area that occurs below 50 meters (0-50m). Is that

Hi all, Somebody knows how can I assign screen geometry as a percentage of the whole screen like rdesktop does?? For example: "spicec -g 80%". Maybe using xorg.conf specific configuration?? Thanks.