similar to: LR-based CIs for GLMs

Overview ======== This module implements the Center for Internet Security (CIS) Security Configuration Benchmark for Red Hat Enterprise Linux 6 v.1.1.0 (avilable at http://benchmarks.cisecurity.org). Each scored control has been implemented as a class or a custom fact. Installation ============ Please either: - Clone git repo from https://github.com/arildjensen/cis-puppet - Run "puppet

Dear list, I would like to fit a negative binomial regression model as described in "Byers AL, Allore H, Gill TM, Peduzzi PN., Application of negative binomial modeling for discrete outcomes: a case study in aging research. J Clin Epidemiol. 2003 Jun;56(6):559-64" to my data in which the response is count data. There are also 10 predictors that are count data, and I have also 3

Hi, we try to do a logistic regression with the function glm. But we notice that this function don't give the same results as the SAS proc catmod (differents estimate given). We try to change the contrast on R system with: > options(contrasts=c(unordered="contr.SAS",ordered="contr.poly")) We also try with brlr and logistf functions. Unfortunately, the estimate

G'day list! I'm thinking about programming likelihood and score intervals for generalized linear models in R based on the paper "On the computation of likelihood ratio and score test based confidence intervals in generalized linear models" by Juha Alho (1992) (Statistics in Medicine, 11, 923-930). Being lazy, I thought that I would ask if anyone else on the list has

Version 2.1.1 Platforms: all What is the trellis parameter (or is there a trellis parameter) to set the leading (the gap between lines) when long axis values labels or panel header labels wrap over more than one line? By default, there is a huge gap between lines, and much looking and experimentation has not revealed to me a suitable parameter to adjust this. Tim C

A recent paper from Google Labs, interesting in many respects, not the least the exclusive use of R for data analysis and graphics (alas not cited in the approved manner): http://labs.google.com/papers/disk_failures.pdf Perhaps some of the eminences grises of the R Foundation could prevail upon Google to make some the data reported in the paper available for inclusion in an R library or two, for

List: I'm helping a colleague with some Poisson regression modeling. He uses SAS proc GENMOD and I'm using glm() in R. Note on the SAS and R output below that our estimates, standard errors, and deviances are identical but what we get for likelihoods differs considerably. I'm assuming that these must differ just by some constant but it would be nice to have some confirmation

Hello, I have different results from these two softwares for a simple binomial GLM problem. >From Genmod in SAS: LogLikelihood=-4.75, coeff(intercept)=-3.59, coeff(x)=0.95 >From glm in R: LogLikelihood=-0.94, coeff(intercept)=-3.99, coeff(x)=1.36 Is there anyone tell me what I did wrong? Here are the code and results, 1) SAS Genmod: % r: # of failure % k: size of a risk set data

Dear All, I am analysing a dataset on levels of herbivory in seedlings in an experimental setup in a rainforest. I have seven classes/categories of seedling damage/herbivory that I want to analyse, modelling each separately. There are twenty maternal trees, with eight groups of seedlings around each. Each tree has a TreeID, which I use as the random effect (blocking factor). There are two

Hello, I have to run a general linear mixed model which looks at 2 dependent variables at the same time (var1 divided by var2). I have tryed to search for such a kind of model structure but since I just started using R my search was not successful. Especielly since I only have an old SAS GENMOD code structure from my project supervisor as an indication. My question is no, does there exist a code

I need capabilities, for my data analysis, like the Pinheiro & Bates S-Plus/R package nlme() but with binomial family and logit link. I need multiple crossed, possibly interacting fixed effects (age cohort of twin when entered study, sex of twin, sampling method used to acquire twin pair, and twin zygosity), a couple of random effects other than the cluster variable, and the ability to

Greetings. I'm trying to recreate in R some regression models I've done in SAS, but I'm not getting the same results. My advisor suspects this may be due to differences in precision between R and SAS. Does anyone know where I can find specifications for R's type double? (It doesn't seem to be in the R Language Definition.) Thanks in advance for any help anyone can

Hello Everybody: I am searching for an R package for fitting Generalized Estimating Equations (GEE) with weights (i.e. Weighted Estimating Equations). From the R documentation I found "geese(geepack)" for fitting Generalized Estimating Equations. In this documentation, under the paragraph “weights” it has been written, “an optional vector of weights to be used in the fitting process.

I'm very new to R and modeling but need some help with visualization of glms. I'd like to make a graph of my glms to visualize the different effects of different parameters. I've got a binary response variable (bird sightings) and use binomial glms. The 'main' response variable is a measure of distance to a track and the parameters I'm testing for are vegetation

Dear friends, As we all know, the usual model selection criteria(e.g.deviance,AIC...) in GLMs isn't very good for selecting the best model when overdispersion exist, so we need to adjust the corresponding statistic,see(Fitzmaurice,G.M. (1997) Model selection with overdispersed

hello R-world, according to the strucchange package .pdf, "all procedures in this package are concerned with testing or assessing deviations from stability in the classical linear regression model." i'd like to test/assess deviations from stability in the Poisson model. is there a way to modify the strucchange package to suit my purposes, or should i use be using another

Dear R friends, When fitting GLMs in R, we may need to specify the variance function to do our analysis. I had thought it's the mean value, but it seems not. Could anybody expain the correct meaning of *mu* in the variance function of GLMs? The following content is from the R-hlep. variance for all families other than quasi, the variance function is determined by the family. The quasi

Dear friends, After fitting a generalized linear models ,i hope to get the variance of alpha,variance of beta and their covariance, that is , the variance-covariance matrix/information of alpha and beta , suppose *B* is the object of GLMs, i use attributes(B) to look for the options ,but can't find it, anybody knows how to get it? > attributes(B) $names [1] "coefficients"

V&R describes binomial GLMs with mortality out of 20 budworms. Is it appropriate to use the same approach with mortality out of numbers as low as 3? I feel reticent to do so with data that is not very continuous. There are one continuous and one categorical independent variables. Would it be more appropriate to treat the response as an ordered factor with four levels? If so, what family

Folks, my question is not R-specific, but I've struck out twice on sci.stat.consult, so I'm turning to the R community. Even if it's a silly question, I expect that someone present will probably tell me so... I have been using multiple Poisson GLMs and similar count-re?gression models to analyse forest songbird abundance data. Many of the spe?cies-level models seem to fit the data