similar to: lme: extract result-function

Dear R members, I tried to fit an lme model and to use the glht function of multcomp. However, the glht function gives me some errors when using weights=varPower(). The glht error makes sense as glht needs factor levels and the model works fine without weights=. Does anyone know a solution so I do not have to change the lme model? Thanks Sibylle --> works fine

I have some questions about the use of lme(). Below, I constructed a minimal dataset to explain what difficulties I experience: # two participants subj <- factor(c(1, 1, 1, 1, 2, 2, 2, 2)) # within-subjects factor Word Type wtype <- factor(c("nw", "w", "nw", "w", "nw", "w", "nw", "w")) # within-subjects factor

Dear all, I have an experimental setup in which a random variable is nested within a fixed variable; however I have troubles specifying the correct LMM with lme. I have searched the lists but haven't been able to find an example like my setup, which I unfortunately need to get this stuff right. Pinheiro & Bates is great but I still can't figure out how to do it. My

Hello all, I am about to send off a manuscript and, although I am fairly confident I have used the lme function correctly, I want to be 100% sure. Could some kind soul out there put my mind at ease? I am simply interested in whether a predictor (SPI) is related to height. However, there are five different populations, and each may differ in mean level of height as well as the relationship

Dear friends of lme, After so many year with lme, I feel ashamed that I cannot get this to work. Maybe it's a syntax problem, but possibly a lack of understanding. We have growth curves of new dental bone that can well be modeled by a linear growth curve, for two different treatments and several subjects as random parameter. By definition, newbone is zero at t=0, so I tried to force the

Dear all, Suppose I have a nested, repeated measure lme model. Which of the following formulae is correct? (assuming data are sampled from several plots in an agricultural experiment) (1) y~explanatory.variables,random=~time|block/plot/subplot/individual (2) y~explanatory.variables,random=~time|unique.ID.of.every.individual I have read that (2) is the only approach that works. But how could I

Hi, I seem to be facing a strange problem when I use Sweave for creating a LaTeX document of the R lme() output --- The EPS and PDF figure files get created, but are empty. I have attached a reproducible example below (taken from the R lme() help example). ------------------------------------------------------------------------------------ \documentclass[a4paper,10pt]{article}

Dear Madam or Sir I am writing you hoping, that you can help me with a problem concerning the output of regressions done with the function lme in R. I would need the standard deviations for intercepts and predictors, but in the output I can only find those for the intercepts. Could it be, that this is my fault? (I am just a beginner with R and multilevel modeling). I am sorry to annoy you with

Dear R-help I'm working on a large dataset which I have divided into 20 subsets based on similar features. Each subset consists of observations from different locations and I wish to use the location as a random effect. For each group I want to select regressors by a stepwise procedure and include a random effect on the intercept. I use stepAIC() and lme(). (The lmer()-function doesn't

Hi xyplot(y ~ x | subject) plots a separate graph of y against x for each level of subject. But I would like to have an own function for each level. Something like xyplot(y ~ x | subject, panel = function(x,y) { panel.xyplot(x,y) panel.curve(x,y) { # something that dependents on the current subject ... } }) How I get the current

Hi I want to do something like a <- c(10,20,15,43,76,41,25,46) tab <- data.frame(a) name <- "a" for (v in unique(tab[[name]])) { r <- subset(tab, name==v) # this does not work ... } i.e. a "string" on the left side of the select expression (subset). How could I solve this? thx Christof

Hi there are some predefined self-start functions, like SSmicmen, SSbiexp, SSasymp, SSasympOff, SSasympOrig, SSgompertz, SSflp, SSlogis, SSweibull, Quadratic, Qubic, SSexp (nlrwr) Btw, do you know graphic examples for this functions? The SSexpDecay (exponential decay) for y = (y0 - plateau)*exp(-k*x) + plateau from

Hi I like to use "axis.POSIXct" to plot days from 2006 till 2008. But I only have datas for the summer months. Is it possible to get two axis breaks, to have not so long distances without points? thx Christof

I am trying to merge two data frames, but one of the column headings are different in the two frames. How can I rjoin or rbind the tho frames? Johnny # Generate 2 blocks by confounding on abc d1 <- conf.design(c(1,1,1), p=2, block.name="blk", treatment.names = c("A","B","C")) d2 <- conf.design(c(1,1,1), p=2, block.name="blk",

Hi R Users, When I use package nlme for linear model with random effects, there exists errors and I don't know the data structure of lme. Here is my data: Rice<-data.frame(Yield=c(8,7,4,9,7,6,9,8,8,8,7,5,9,9,5,7,7,8,8,8,4,8,6,4,8,8,9), Variety=rep(rep(c("A1","A2","A3"),each=3),3),

Hi guys, Got a few days left and I need to model a random effect of species on the body mass (logM) and temperature (K) slopes. This is what i've done so far that works: model1<-lme(logSSP~logM + K,random=~1|species,data=data1) model2<-lme(logSSP~logM + K,random=~K|species,data=data1) model3<-lme(logSSP~logM + K,random=~logM|species,data=data1) The one I now want is:

Hello, I want to perform an lme on a database with this structure: ID Sequence Temperature Tumour Error 1 5 0 1 8.721872e-08 1 5 0 2 8.695348e-08 1 5 0 3 2.019604e-13 1 5 37 1

Hi all, I have a problem with package installing in Windows, on my PC machine. The end goal is to be able to use the lme() function. Here's what I did so far: > install.packages("lme4") Warning in install.packages("lme4") : argument 'lib' is missing: using 'C:\Users\Angela\Documents/R/win-library/2.9' --- Please select a CRAN mirror for use in this

Dear R-community, I'm analysing some noise using the nlme-package. I'm writing in order to get my usage of lme verified. In practise, a number of samples have been processed by a machine measuring the same signal at four different channels. I want to model the noise. I have taken the noise (the signal is from position 1 to 3500, and after that there is only noise). My data looks like

Hi I would like to shorten mod1 <- nls(ColName2 ~ ColName1, data = table, ...) mod2 <- nls(ColName3 ~ ColName1, data = table, ...) mod3 <- nls(ColName4 ~ ColName1, data = table, ...) ... is there something like cols = c(ColName2,ColName3,ColName4,...) for i in ... mod[i-1] <- nls(ColName[i] ~ ColName1, data = table, ...) I am looking forward to help Christof