Displaying 7 results from an estimated 7 matches for "glm4".
Did you mean:
glm
2004 May 01
0
glm.nb and anova
Hi,
I am trying to fit a negative binomial model with a number of parasite tapeworms as response variable to geographical coordinates
(actually preparing a trend surface before kriging). When I try an anova, I get warnings:
> glm4.nb<-glm.nb(wb~X4+Y4+I(X4^2)+I(Y4^2))
> anova(glm4.nb)
Analysis of Deviance Table
Model: Negative Binomial(0.0463), link: log
Response: wb
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev P(>|Chi|)
NULL 344 225.7
X4 1...
2002 Jan 02
0
comparative rendering of modeling outputs
...presentation in applications contexts.
An extract from the output of 'example(cremo)':
cremo> data(stackloss)
cremo> attach(stackloss)
cremo> glm2 <- glm(stack.loss ~ Air.Flow + Water.Temp + Acid.Conc.)
cremo> glm3 <- glm(stack.loss ~ Air.Flow + Water.Temp)
cremo> glm4 <- glm(stack.loss ~ Air.Flow * Water.Temp)
cremo> try(print(cremo(list(glm2, glm3, glm4), adeqel = c("deviance",
"aic")), quote = F))
var Model 1 Model 2 Model 3
[1,] (Intercept) -39.92 -50.359 22.29
[2,] (11...
2011 Mar 31
1
rank of Matrix
...using the functions/classes from the relatively new
MatrixModels package:
library(MatrixModels)
str(trial <- data.frame(counts=c(18,17,15,20,10,20,25,13,12),
outcome=gl(3,1,9,labels=LETTERS[1:3]),
treatment=gl(3,3,labels=letters[1:3])))
glmS <- glm4(counts ~ 0+outcome + treatment, poisson, trial,
verbose = TRUE, sparse = TRUE)
str(glmS)
str(X <- glmS at pred@X)
str(QR <- qr(X))
I'm not very familiar with matrix decomposition, but is the 'p' slot
from a sparseQR object the same as $pivot from base qr? Ad...
2006 Mar 16
2
DIfference between weights options in lm GLm and gls.
...00
Residual Deviance: 44800 AIC: 11670
> glm3
Call: glm(formula = y ~ x, weights = 1/W)
Coefficients:
(Intercept) x
0.03104 7.31033
Degrees of Freedom: 1242 Total (i.e. Null); 1241 Residual
Null Deviance: 798900
Residual Deviance: 19900 AIC: 5285
> glm4
Call: glm(formula = ys ~ Xs - 1)
Coefficients:
Xs Xsx
2.687 6.528
Degrees of Freedom: 1243 Total (i.e. Null); 1241 Residual
Null Deviance: 4490000
Residual Deviance: 506700 AIC: 11000
With weights, the glm did not give the same results as lm why?
Also for gls, I use varFixe...
2002 Apr 30
1
MemoryProblem in R-1.4.1
...glm3, myfun) #does not works
Error: cannot allocate vector of size 31 Kb
In addition: Warning message:
Reached total allocation of 300Mb: see help(memory.size)
Again I increase the memory.size
> memory.limit(size=320)
NULL
> out3<-lapply(list.glm3, myfun) #works!
> out4<-lapply(list.glm4, myfun) #does not work!
.....
So it seems I have to increase the memory.size each time before applying my
function. This is suprising because I know that returning to the prompt the
memory is fully available again. So being the lists similar, why the same
memory size is not sufficient for every lis...
2011 Mar 27
0
model diagnostics for MatrixModels
...ere facilities to obtain, for example, the rank of
the (sparse) model matrix efficiently? Quite a few tests require rank,
but computing this separately using for example rankMatrix from
package Matrix uses up tonnes of memory, which is exactly what I was
trying to avoid by using sparse matrices and glm4 :-p
Thanks for any pointers,
Carson
> sessionInfo()
R version 2.12.2 (2011-02-25)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_IE.utf8 LC_NUMERIC=C
[3] LC_TIME=en_IE.utf8 LC_COLLATE=en_IE.utf8
[5] LC_MONETARY=C LC_MESSAGES=en_IE.utf8
[7] LC_PAPE...
2003 Apr 24
2
R-1.7.0 build feedback: NetBSD 1.6 (PR#2837)
R-1.7.0 built on NetBSD 1.6, but the validation test suite failed:
Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC)
Remote gcc version: gcc (GCC) 3.2.2
Remote g++ version: g++ (GCC) 3.2.2
Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib
make[5]: Entering directory `/local/build/R-1.7.0/src/library'
>>> Building/Updating