Displaying 3 results from an estimated 3 matches for "hours24".
2010 May 28
1
Comparing and Interpreting GAMMs
Dear R users
I have a question related to the interpretation of results based on GAMMs using Simon Woods package gamm4.
I have repeated measurements (hours24) of subjects (vpnr) and one factor with three levels (pred). The outcome (dv) is binary.
In the first model I'd like to test for differences among factor levels (main effects only):
gamm.11<-gamm4(dv ~ pred +s(hours24), random = ~ (1|vpnr),data=sdata,family=binomial)
In the second model I...
2010 Sep 10
1
modifying axis labels in lattice panels
...idvar = "time", direction="long",
timevar = "PC");
#get descriptive shingle labels
Xlong[,"PC"] <- factor(Xlong[,"PC"], labels = paste("PC",1:2));
xyplot(Comp ~ time | PC ,data = Xlong, pane1l = WeekhourPanel, scales =
list(x=list(at = Hours24-4*3600,
labels=as.character(format(Hours24-4*3600,"%H")))));
WeekhourPanel <- function(x,y,...){
r <- range(x);
#print(r)
Hours8 <- seq(r[1], r[2], by=8*3600);
Hours24 <- seq(r[1]+12*3600, r[2], by=24*3600)
#axis.POSIXct(1, at= Hours8, format="%H");
pane...
2010 Jan 26
1
AIC for comparing GLM(M) with (GAM(M)
Hello
I'm analyzing a dichotomous dependent variable (dv) with more than 100
measurements (within-subjects variable: hours24) per subject and more
than 100 subjects. The high number of measurements allows me to model
more complex temporal trends.
I would like to compare different models using GLM, GLMM, GAM and
GAMM, basically do demonstrate the added value of GAMs/GAMMs relative
to GLMs/GLMMs, by fitting splines...