similar to: logLik for poisson models

Displaying 20 results from an estimated 1000 matches similar to: "logLik for poisson models"

2007 Dec 20
3
mle
Dear all, I'm trying to estimate the parameters of a special case of a poisson model, where the specified equation has an integral and several fixed parameters. I think that the MLE command in STATS4 package could be a good choice, but it's a little complicated. I've got some problems with the offset and I don't understand some of the functions. Do you know where can I find some
2008 Aug 20
3
bug in lme4?
Dear all, I found a problem with 'lme4'. Basically, once you load the package 'aod' (Analysis of Overdispersed Data), the functions 'lmer' and 'glmer' don't work anymore: library(lme4) (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) (gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), family = binomial, data
2008 May 28
1
postscript
Hello, I'm trying to save a graph as a .ps file to include it in a LaTex presentation, but R replies: Error: Invalid font type In addition: Warning messages: 1: font family not found in PostScript font database 2: font family not found in PostScript font database If I load the package grDevices and type "postscript()" then I can save it, but it creates also a text file saying:
2009 Nov 13
1
multivariate meta-analysis with the metafor package
Dear Wolfgang Viechtbauer and R users, I have few questions regarding the development of the package 'metafor. As you suggested , I post to the R-help mailing list. I read you're planning an extension of this method to the multivariate case. I think it would be a useful tool. I'm currently performing some analyses with R on multiple outcomes, using the Stata command mvmeta to get
2008 Jun 08
1
exponential distribution
Dear all, I've tried to solve the Es. 12, cap 4 of "Introduction to GLM" by Annette Dobson. It's about the relationship between survival time of leukemia patients and blood cell count. I tried to fit a model with exponential distribution, first by glm (family gamma and then dispersion parameter fixed to 1) and then with survreg. They gave me the same point estimates but the
2008 Aug 19
1
R vs Stata on generalized linear mixed models: glmer and xtmelogit
Hello, I have compared the potentials of R and Stata about GLMM, analysing the dataset 'ohio' in the package 'faraway' (the same dataset is analysed with GEE in the book 'Extending the linear model with R' by Julian Faraway). Basically, I've tried the 2 commands 'glmmPQL' and 'glmer' of R and the command 'xtmelogit' of Stata. If I'm not
2009 Nov 05
2
new help pages in R 2.10.0
Dear R-users, I've just installed the last version of R (2.10.0) in Windows and found a completely new version of the help pages. Basically, they now appear as a simple (even if slightly embellished) text file. I tried to install it with a customized startup and realized that the old CHM files (the default for Windows up to 2.9.2) are not included anymore as a possible choice. I don't
2009 Jun 01
1
problem with package creation under Windows and R 2.9.0
Dear R users, I found a small problem with the package creation procedure in R 2.9.0 under Windows. Basically, I input in the command prompt; Rcmd BUILD --binary --use-zip namepackage The procedure works, but the titles in the help pages of the created package lose the usual colours and font. Conversely, the result is perfect using R 2.8.1 (changing the path) and the same version of Rtools
2009 Apr 19
1
creating a package in R
Dear R users, I created some functions I want to convert in a package under Windows (Vista or XP). I carefully read both "Writing R extensions" by the R Development Core Team and "Creating R packages: a tutorial" by Friedrich Leisch. I've also installed Rtools as suggested. These tutorials are very useful, but I'm struggling with the last step of the installation
2010 Jan 04
1
glmer (lme4), glmmPQL (MASS) and xtmepoisson (Stata)
Dear R users, I'm trying to specify a generalized linear mixed model in R, basically a Poisson model to describe monthly series of counts in different regions. My aim is to fit subject-specific curves, modelling a non-linear trend for each region through random effects for linear splines components (see Durban et al, Stat Med 2005, or " Semiparametric regression" by Ruppert et al,
2008 Oct 31
1
AIC for quasipoisson link
Dear fellows, I'm trying to extract the AIC statistic from a GLM model with quasipoisson link. The formula I'm referring to is AIC = -2(maximum loglik) + 2df * phi with phi the overdispersion parameter, as reported in: Peng et al., Model choice in time series studies os air pollution and mortality. J R Stat Soc A, 2006; 162: pag 190. Unfortunately, the function logLik
2009 Feb 09
1
gee with auto-regressive correlation structure (AR-M)
Dear all, I need to fit a gee model with an auto-regressive correlation structure and I faced some problems. I attach a simple example: ####################################################### library(gee) library(geepack) # I SIMULATE DATA FROM POISSON DISTRIBUTION, 10 OBS FOR EACH OF 50 GROUPS set.seed(1) y <- rpois(500,50) x <- rnorm(500) id <- rep(1:50,each=10) # EXAMPLES FOR
2009 Jul 06
2
problem with internal functions in Windows
Dear R users, I included 2 internal functions in the package 'dlnm', called 'mkbasis' and 'mklagbasis'. Despite they are not meant to be called by the users, I included them in the namespace in order to make them available, keeping the process more transparent and giving the opportunity to change or improve them. I included an help page 'dlnm-internal.Rd' to
2009 Jul 06
2
problem with internal functions in Windows
Dear R users, I included 2 internal functions in the package 'dlnm', called 'mkbasis' and 'mklagbasis'. Despite they are not meant to be called by the users, I included them in the namespace in order to make them available, keeping the process more transparent and giving the opportunity to change or improve them. I included an help page 'dlnm-internal.Rd' to
2009 Jul 20
0
new package 'dlnm' to run distributed lag non-linear models
Dear R Community, I am pleased to announce the release of a new package called 'dlnm', now available on CRAN (version 0.2.1). The package dlnm provides some facilities to run distributed lag models (DLM's) and their non-linear extension (DLNM's), a modelling framework to describe simultaneously non-linear and delayed effects between predictors and an outcome in time-series
2009 Jul 20
0
new package 'dlnm' to run distributed lag non-linear models
Dear R Community, I am pleased to announce the release of a new package called 'dlnm', now available on CRAN (version 0.2.1). The package dlnm provides some facilities to run distributed lag models (DLM's) and their non-linear extension (DLNM's), a modelling framework to describe simultaneously non-linear and delayed effects between predictors and an outcome in time-series
2009 Aug 03
0
problem with predict(...,type="terms")
Hello R users, I'm trying to understand the meaning of the results returned by the command 'predict' with type="terms". I thought it was the matrix with the fitted values for each term (value of predictor * related coef), but the calculation is different. I include an example below: y <- rnorm(10) x1 <- rnorm(10) x2 <- rnorm(10) model <- lm(y ~ x1 + x2)
2009 May 20
0
function returning the name of an object
Hello R users, I'm trying to solve an apparently simple problem. I include a matrix in a model formula, and I want then to extract the related coefficients. Basically, I want a function func such that: mymatrix <- cbind(rnorm(10),rnorm(10)) y <- b <- rnorm(10) model <- lm(y ~ b + mymatrix) func(mymatrix,model) which returns only the 3rd and 4th coefficients. I thought to
2010 Feb 15
1
get current window (device) size
Dear R-listers, I am writing a function that will generate a simple plot for the user . The size of the window (that contains the plot) may not be what the user wants. Supposing the user adjusts the window to a size of her/his choice, I would like to be able to save the new window size and specify the plot to be of this new window size when the user calls my function again. My question is
2010 Dec 01
2
default arguments and '...' in a function
Dear R-users, I'm trying to work out a way to set default values for arguments in a function which includes the optional argument '...'. In practice, I have a 'plot' method for a function which specifies different types of plots. Every different plot should have different default arguments (for example for 'col', 'ylim' etc), plus the argument '...' to