Displaying 20 results from an estimated 100 matches similar to: "influence.measures, cooks.distance, and glm"
2006 Oct 24
1
Cook's Distance in GLM (PR#9316)
Hi Community,
I'm trying to reconcile Cook's Distances computed in glm. The
following snippet of code shows that the Cook's Distances contours on
the plot of Residuals v Leverage do not seem to be the same as the
values produced by cooks.distance() or in the Cook's Distance against
observation number plot.
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
2005 Feb 02
1
anova.glm (PR#7624)
There may be a bug in the anova.glm function.
deathstar[32] R
R : Copyright 2004, The R Foundation for Statistical Computing
Version 2.0.1 (2004-11-15), ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project
2006 May 08
0
Inconsistency in AIC values for glm with family poisson (PR#8841)
This message is in MIME format. The first part should be readable text,
while the remaining parts are likely unreadable without MIME-aware tools.
--27464147-1557463723-1147085467=:8118
Content-Type: TEXT/PLAIN; charset=iso-8859-1; format=flowed
Content-Transfer-Encoding: 8BIT
On Mon, 8 May 2006, x.sole at iconcologia.net wrote:
> Full_Name: Xavier Sol?
> Version: 2.3.0
> OS: Windows
2009 Aug 20
0
Sweave truncation
Peter
Thank you for the information. I accidentally deleted Ken's post without
having read it.
Ken' s thought is great but as you said awful to implement
I thought that capture.output would come in handy some time when I first
saw it on an unrelated reply.
Just thought :- the latex listings package may have alternatives
If I remember correctly it has wrapping and other goodies but I
2013 May 29
1
quick question about glm() example
I don't have a copy of Dobson (1990) from which the glm.D93 example is
taken in example("glm"), but I'm strongly suspecting that these are
made-up data rather than real data; the means of the responses within
each treatment are _identical_ (equal to 16 2/3), so two of the
parameters are estimated as being zero (within machine tolerance). (At
this moment I don't understand
2013 Feb 18
1
nobs() with glm(family="poisson")
Hi!
The nobs() method for glm objects always returns the number of cases
with non-null weights in the data, which does not correspond to the
number of observations for Poisson regression/log-linear models, i.e.
when family="poisson" or family="quasipoisson".
This sounds dangerous since nobs() is, as the documentation states,
primarily aimed at computing the Bayesian
2000 May 09
4
Dispersion in summary.glm() with binomial & poisson link
Following p.206 of "Statistical Models in S", I wish to change
the code for summary.glm() so that it estimates the dispersion
for binomial & poisson models when the parameter dispersion is
set to zero. The following changes [insertion of ||dispersion==0
at one point; and !is.null(dispersion) at another] will do the trick:
"summary.glm" <-
function(object, dispersion =
2005 Oct 15
2
how to import such data to R?
the data file has such structure:
1992 6245 49 . . 20 1
0 0 8.739536 0 . . .
. . . . . "alabama"
. 0 .
1993 7677 58 . . 15 1
0 0
2010 Feb 21
1
tests for measures of influence in regression
influence.measures gives several measures of influence for each
observation (Cook's Distance, etc) and actually flags observations
that it determines are influential by any of the measures. Looks
good! But how does it discriminate between the influential and non-
influential observations by each of the measures? Like does it do a
Bonferroni-corrected t on the residuals identified by
2009 Aug 19
3
Sweave output from print.summary.glm is too wide
Hi all
I am preparing a document using Sweave; a really useful tool. But I am having a problem.
Consider this toy example Sweave file:
\documentclass{article}
\begin{document}
<<echo=TRUE,results=verbatim>>=
options(width=40) # Set width to 40 characters
hide <- capture.output(example(glm)) # Create an example of the problem, but hide the output
summary(glm.D93) #
2005 Apr 06
2
make error in R devel
Dear list,
I just hit an error that stopped my make && make check-devel operation
on my linux box (FC3, i686 P4 2GB RAM). Just to note that I've been
building the development branch(?) for some time and this is the first
hint of a problem.
1) updated the src tree using svn update
2) ran ../configure --with-recommended-package=no from my build directory
3) got:
R is now configured
2012 Jan 12
1
posting for r-help
Hi there
I have a post I would like to put on the "95% confidence intercal with glm"
thread. Thank-you so much!
I am wondering first of all if anyone knows how to calculate confidence
intervals for a GLMM? I use the lme4 library.
Also, I am wondering how to predict a model mean and confidence intervals for a
particular independent variable?
For example in the following example:
1999 Jun 23
1
Influence.measures
I am using rw0641 with Windows 98. To list just the influential
repetitiones that result from "influence.measures", I am using the input
result <- lm(y~x)
and the code from the example in the help for "influence.measures"
INFLM <- function(result){
inflm <- influence.measures(result)
which(apply(inflm$is.inf,1,any))
}
It works fine up to now with the
2008 May 26
1
Sweave does not respect width
Hello,
I'm learning to use Sweave, and I've run into a problem: sometimes, when
entering long lines of input and using long variable names, Sweave will not
insert linebreaks in a way that respects the width setting. This causes
undesirable overflows into the margins in the latex file. For example,
consider the following document (adapted from the GLM example):
\documentclass{article}
2011 May 11
2
New code in R-devel: Rao score test for glm.
I have just committed some code to the r-devel branch to implement the Rao efficient score test. This is asymptotically equivalent to the LRT, but there is some indication that it might have better properties in smaller samples since it is based more directly on the distribution of the sufficient sums under the null hypothesis (e.g., if you have a divergent fit to the model under the alternative,
2009 Feb 23
1
Follow-up to Reply: Overdispersion with binomial distribution
THANKS so very much for your help (previous and future!). I have a two
follow-up questions.
1) You say that dispersion = 1 by definition ....dispersion changes from 1
to 13.5 when I go from binomial to quasibinomial....does this suggest that
I should use the binomial? i.e., is the dispersion factor more important
that the
2) Is there a cutoff for too much overdispersion - mine seems to be
2011 Apr 29
1
logistic regression with glm: cooks distance and dfbetas are different compared to SPSS output
Hi there,
I have the problem, that I'm not able to reproduce the SPSS residual
statistics (dfbeta and cook's distance) with a simple binary logistic
regression model obtained in R via the glm-function.
I tried the following:
fit <- glm(y ~ x1 + x2 + x3, data, family=binomial)
cooks.distance(fit)
dfbetas(fit)
When i compare the returned values with the values that I get in SPSS,
2009 Feb 16
1
Overdispersion with binomial distribution
I am attempting to run a glm with a binomial model to analyze proportion
data.
I have been following Crawley's book closely and am wondering if there is
an accepted standard for how much is too much overdispersion? (e.g. change
in AIC has an accepted standard of 2).
In the example, he fits several models, binomial and quasibinomial and then
accepts the quasibinomial.
The output for residual
2007 Jul 21
1
Gamma MLE
Hello,
I was asked to try the following code on R,
gamma.mles
function (xx,shape0,rate0)
{
n<- length(xx)
xbar<- mean(xx)
logxbar<- mean(log(xx))
theta<-c(shape0,rate0)
repeat {
theta0<- theta
shape<- theta0[1]
rate<- theta0[2]
S<- n*matrix(c(log(rate)-digamma(shape)+logxbar,shape/rate-xbar),ncol=1)
I<- n*matrix(c(trigamma(shape),-1/rate,-1/rate,shape/rate^2),ncol=2)
2012 Apr 09
2
Overall model significance for poisson GLM
Greetings,
I am running glm models for species counts using a poisson link function.
Normal summary functions for this provide summary statistics in the form of
the deviance, AIC, and p-values for individual predictors. I would like to
obtain the p-value for the overall model. So far, I have been using an
analysis of deviance table to check a model against the null model with the
intercept as