search for: longdata

Fellow R-users, I have a longitudinal data set with missing values in it. I would like to produce a residual plot for each time using panel.xyplot function but I get an error message. Here's a simple example, library(nlme) set.seed(1544) longdata <- data.frame(ID=gl(10,1,50), y=rnorm(50), time = as.numeric(gl(5,10,50)), x = rnorm(50)) longdata$y[5] <- NA longdata$x[35] <- NA modlme<- lme(fixed=y ~ time + x,                 random= ~ 1 | ID,                 na.action=na.exclude, data=longdata) plot( modlme, abs( resid(., type...

...I build a package. They was a mistake in it, but R CMD check did not find it. Is that normal ? Here is what Kurt gets (which is right, I did this mistake): --- 8< ---------------- * checking for code/documentation mismatches ... WARNING S4 class codoc mismatches from documentation object 'LongData-class': Slots for class 'LongData' Code: id other time traj varName Docs: id time traj varName --- 8< ---------------- But here is what I get : --- 8< ---------------- R CMD check longitudinalData * checking for working pdflatex ... OK * using log directory 'C:/Documents...

...lways the same, I tell people to remove the package Matrix, and then kml works fine. I wonder what my code is not compatible with the Matrix package. Here is a (simplified) version of the kml code : --- 8< ------------ library(longitudinalData) setClass(Class="CLD", contains=c("LongData","ListPartition")) myCld <- new("CLD", idAll=as.character(c(100,102,103,109,115,123)), idFewNA=as.character(c(100,102,103,109,115,123)), time=c(1,2,4,8,15), varNames="P", traj=matrix(c(1,2,3,1,4, 3,6,1,8,10, 1,2,1,3,2, 4,2,5,6,3...

...takes hours to do the calculation on a fast notebook under Windows. I suspect by(.) to be a poor implementation for doing individual analysis on subgroups. Is there an alternative and more efficient way for doing by-group processing within lme(.). Here some code to give you a glimpse gfit <- by(longdata, gen, function(x) lme(fixed=response ~ dye + C(treat, base = 4 ), data=x,random =~ 1 | slide) ) Thanks in advance & regards Gerhard Krennrich ------------------------------------------------------------------------ BASF Aktiengesellschaft GVC/S - B009 Gerhard Krennrich Tel.: +49...

...I build a package. They was a mistake in it, but R CMD check did not find it. Is that normal ? Here is what Kurt gets (which is right, I did this mistake): --- 8< ---------------- * checking for code/documentation mismatches ... WARNING S4 class codoc mismatches from documentation object 'LongData-class': Slots for class 'LongData' Code: id other time traj varName Docs: id time traj varName --- 8< ---------------- But here is what I get on my computer: --- 8< ---------------- C:\Documents and Settings\Christophe\Mes documents\Recherche\Trajectoires\kmeal> R CMD chec...

...that? Thanks for your help Christophe --- 8< ----------------- Code for reproductible example ------------------- library(kml) dn <- gald(1) dn["traj"] library(lme4) dn["traj"] setMethod( ### Simplified version ### "[", signature=signature(x="ClusterLongData", i="character", j="ANY",drop="ANY"), definition=function (x, i, j="missing", ..., drop = TRUE){ x <- as(x, "LongData") return(x[i, j]) } ) dn["traj"] dn["traj"] --- 8< ----------------- Execution of...

...all the package and I try an example: --- 8< -------------- library(kml) dn <- as.cld(gald()) kml(dn) # XXX ~ Fast KmL ~ # Erreur dans as.vector(x, mode) : argument 'mode' incorrect --- 8< -------------- So I make some verifications: --- 8< ---- class(dn) # [1] "ClusterizLongData" # attr(,"package") # [1] "kml" getMethod("kml","ClusterizLongData") # Method Definition: # # function (Object, nbClusters = 2:6, nbRedrawing = 20, saveFreq = 100, # maxIt = 200, trajMinSize = 2, print.cal = FALSE, print.traj = FALSE, # imputationM...