similar to: pdIdnot / logLik in glmmPQL

Hi all, I've been attempting to fit a logistic glmm using glmmPQL in order to estimate variance components for a score test, where the model is of the form logit(mu) = X*a+ Z1*b1 + Z2*b2. Z1 and Z2 are actually reduced rank square root matrices of the assumed covariance structure (up to a constant) of random effects c1 and c2, respectively, such that b1 ~ N(0,sig.1^2*I) and c1 ~

Hi, all I can not figure this out, please have a look and help me out. thank you! Note: this is in SPLUS, not R. I have following code *********************************** modfit<-function(yir,yew, ft) { n<-length(yew) yew<-yew[1:(n-1)] yy<-yir-ft xx<-yew-ft n<-length(xx) xx0<-xx[2:n] yy0 <-yy [2:n] xx1<-xx[1:(n-1)] fit <- garch(yy0~xx0 + xx1+var.in.mean,

Hello, I am a bit confused about gamm in mgcv. Consulting Wood (2006) or Ruppert et al. (2003) hasn't taken away my confusion. In this code from the gamm help file: b2<-gamm(y~s(x0)+s(x1)+s(x2)+s(x3),family=poisson,random=list(fac=~1)) Am I correct in assuming that we have a random intercept here....but that the amount of smoothing is also changing per level of the

Hi all - R2.0.1, OS X Perhaps while there is some discussion of lme going on..... I am trying to execute a glmm using glmmPQL from the MASS libray, using the example data set from McCullagh and Nelder's (1989, p442) table 14.4 (it happens to be the glmm example for GENSTAT as well). The data are binary, representing mating success (1,0) for crosses between males and females from two

I'm a beginer programming C and I have the following problem: I have the following C-code file #include <stdlib.h> void gen(int *n, int *a, int *c, int *m, int *x0, int *x); main(){ int nn = 31; int aa = 3; int cc = 0; int mm = 31; int xx0 = 9; int xx[nn]; int i; gen(&nn,&aa,&cc,&mm,&xx0,xx); for (i = 0; i <= nn-1; i++) printf("%d

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,

Hello, I have a question about the add1() function and quasilikelihoods for GLMs. I am fitting quasi-Poisson models using glm(, family = quasipoisson). Technically, with the quasilikelihood approach the deviance does not have the interpretation as a likelihood-based measure of sample information. Functions such as stepAIC() cannot be used. The function add1() returns the change in the scaled

Dear R users, I'm attempting to fit a GLM with random effects using the tweedie family for the error structure. I'm getting the error: iteration 1 Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3) I'm running V2.1.0 I notice from searching the lists that the same error was reported in May 2004 by Spencer Graves, but no-one was able to

Hi to all, i have dovecot 2.2.5 and i have implemented lda rejection through quota full this is my dovecot conf protocol lda { mail_fsync = optimized auth_socket_path = /usr/local/var/run/dovecot/auth-userdb deliver_log_format = msgid=<%f>-<%s>-%m: %$ hostname = mail.cgilfe.it info_log_path = /var/log/dovecot/dovecot-deliver.log log_path =

Hello, I have the following r-codes for solving a quasilikelihood estimating equation: >library(geepack) >fit<-geese(y~x1+x2+x3,jack=TRUE,id=id,scale.fix=TRUE,data=dat,mean.link = "logit", corstr="independence") Now my question is how can I calculate the confidence interval of the parameters of the above model "fit"? [[alternative HTML version deleted]]

Is it possible to define an external special function (SPECIALSXP)? I'm trying to do some language-level work, and don't want my arguments evaluated before they hit C. It looks like the only way to define a SPECIALSXP is by using XX0 in the `eval' field of R_FunTab; is there any way to make this applicable to externally defined functions?

Hello, Why would I get an error message with the following code for gamm? I want to fit the a gam with different variances per stratum. library(mgcv) library(nlme) Y<-rnorm(100) X<-rnorm(100,sd=2) Z<-rep(c(T,F),each=50) test<-gamm(Y~s(X),weights=varIdent(form=~1|Z)) summary(test$lme) #ok summary(test$gam) Gives an error message: Error in inherits(x, "data.frame")

Dear R-helpers, I'd like to understand how to test the statistical significance of a random effect in gamm. I am using gamm because I want to test a model with an AR(1) error structure, and it is my understanding neither gam nor gamm4 will do the latter. The data set includes nine short interrupted time series (single case designs in education, sometimes called N-of-1 trials in medicine)

Hello, I have the following r-codes for solving a quasilikelihood estimating equation: >library(geepack) >fit<-geese(y~x1+x2+x3,jack=TRUE,scale.fix=TRUE,data=dat,mean.link = "logit", corstr="independence") Now my question is how can I calculate the confidence interval of the parameters of the above model "fit"? [[alternative HTML version deleted]]

Hi list, I try to connect to a GW which have one domain eg sip.mydomain.com and have few IPs related to this domain. I register * to this domain with host=sip.mydomain.com and type=user. So DNS will decide on which IP of my domain I will register (or redirection on the GW side). If an incoming call arrive, I would guess that, as type=user, it will not try to match the IP from INVITE as I

Hi R helpers, I'm trying to fit a rather complex model to some simulated data using lme and am not getting the correct results. It seems there might be some identifiability issues that could possibly be dealt with by specifying starting parameters - but I can't see how to do this. I'm comparing results from R to those got when using GenStat... The raw data are available on the

I have been using lme from nlme to do a 3-way anova with all the effects treated as random. I was wondering if someone could direct me to an example of how to do this using lmer from lme4. I have 3 main effects, tim, trt, ctr, and all the interaction effects tim*trt*ctr. The response variable is ge. Here is my lme code: dat <-

Douglas Bates wrote: > > The GLMM function in the lme4 package allows you to specify crossed > random effects within the random argument without the need for the > pdBlocked and pdIdent constructions. Simply ensure that your grouping > factors are defined in such a way that each distinct group has a > different level in the grouping factor (this is usually not a problem

Hi all, I have an experiment where 5 raters assessed the quality of 24 web sites. (each rater rated each site once). I want to come up with a measure of reliability of the ratings for the web sites ie to what extent does each rater give the same (or similar) rating to each web site. My idea was to fit a random effects model using lme and from that, calculate the intraclass correlation as a

Consider the Assay data where block, sample within block and dilut within block is random. This model can be fitted with (where I define Assay2 to get an ordinary data frame rather than a grouped data object): Assay2 <- as.data.frame(Assay) fm2<-lme(logDens~sample*dilut, data=Assay2, random=list(Block = pdBlocked(list(pdIdent(~1), pdIdent(~sample-1),pdIdent(~dilut-1))) )) Now, block