similar to: arrows and points for error bars

Can you provide the summary(m2) results? > -----Original Message----- > From: Simon Pickett [mailto:S.Pickett at exeter.ac.uk] > Sent: Wednesday, August 16, 2006 7:14 AM > To: Doran, Harold > Cc: r-help at stat.math.ethz.ch > Subject: [SPAM] - RE: [R] REML with random slopes and random > intercepts giving strange results - Bayesian Filter detected spam > > Hi again,

Hi, I am an R beginner and having problems structuring my REML models. I have a model with y=weight x1=time x2=timesquared id=individual identity I need to structure the model such that in the random effects there is a constant intercept for all individuals but a separate individual slope for both x1 and x2 (a coefficient score for every individual).

Hi R experts, I have been using ReML as follows... model<-lmer(late.growth~mtf+year+treat+hatch.day+hatch.day:year+hatch.day:treat+ mtf:treat+ treat:year+ year:treat:mtf+(1|fybrood), data = A) then I wanted to plot the results of the three way interaction using lines() as follows... tmp<-as.vector(fixef(model)) graph1<-plot(mtf,fitted(f2), xlab=list("Brood Size"),

Hi, I would like to (if possible) set the default width and height for graphs at the start of each session and have each new graphic device overwrite the previous one. I only know how to do this using windows(width=,height=...) which opens up a new plotting device every time, so I end up with lots of graphs all over the place until I get the one I want! Thanks in advance, Simon Simon Pickett

Hi everyone, I have been using REML to derive intercepts and coeficients for each individual in a growth study. So the code is m2 <- lmer(change.wt ~ newwt+(newwt|id), data = grow) Calling coef(model.lmer) gives a matrix with this information which is what I want. However, as a test I looked at each individual on its own and used a simple linear regression to obtain the same information, then

Hi, I have been dabbling with str() to extract values from outputs such as lmer etc and have found it very helpful sometimes. but only seem to manage to extract the values when the output is one simple table, any more complicated and I'm stumped :-( take this example of the extracted coeficients from a lmer analysis... using str(coef(lmer(resp3~b$age+b$size+b$pcfat+(1|sex), data=b)))

Hi all, I have been using princomp() recently, its very useful indeed, but I have a question about how to specify the rows of data you want it to choose. I have a set of variables relating to bird characteristics and I have been using princomp to produce PC scores from these. However since I have multiple duplicate entries per individual (each bird had a varying number of chicks), I only want

Hi R experts, I am interested on the effects of two dietry compunds on the growth of chicks. Rather than extracting linear growth functions for each chick and using these in an analysis I thought using ReML might provide a neater and better way of doing this. (I have read the pdf vignette("MlmSoftRev") and "Fitting linear mixed models in R" by Douglas Bates but I am not

Hi, Thanks for your response, it nearly worked! But it only wrote one coloumn of data and not the three columns I need. Using fixef(m1) doesnt give the same results as coef(m1) when you are using more than one random effect. I need the coefficients for each individual so I use coef(m1) to get this which results in an object of class lmer.coef, 3 columns by 700 rows. as.data.frame() wont work on

Hi all, I have been bumbling around with r for years now and still havent come up with a solution for plotting reliable graphs of relationships from a linear regression. Here is an example illustrating my problem 1.I do a linear regression as follows summary(lm(n.day13~n.day1+ffemale.yell+fmale.yell+fmale.chroma,data=surv)) which gives some nice sig. results Coefficients:

Hi, I am an R newbie and one thing I am having trouble with binding variables that I have created within one data frame into a new data frame when using na.omit(). To illustrate this problem I will give the example I am working on and the approah I have been using:- data.frame1<-filepath.... attach(data.frame1) #create a new variable using a function new.variable<-rep(1,length(weight3))

I have been using lmer in lme4() to analyse the effect of dropping a term from the model as below anova.name<-anova(m1,m2) Where m1 is an original model and m2 has one term removed. I can then create my own type III tables for each variable in the model. However with many variables this becomes exceedingly time consuming! anova.lme in nlme() and Anova() in car dont seem to work for lme

Dear R users, Recent changes to the MuMIn package now means that the model averaging command (model.avg) no longer returns confidence intervals, but instead returns zvalues and corresponding pvalues for fixed effects included in models. Previously I have used this package for model selection/averaging following Greuber et al (2011) where it suggests that one should use confidence intervals from

Hi all, I'm hoping I might be able to get some help with some issues specifying priors for MCMCglmm. I'm trying to fit a gaussian glmm using MCMCglmm to a data set with two (correlated) response variables. The response variables are both logit-transformed proportions (there are a few reasons why I've chosen these with gaussian error over binomal glmm, which I won't go into).

Hello there This is the first time I have used r-help message board so I hope I have got the right address. I am trying to check the residuals of a GLMM model(run using the package lme4). I have been able to check the residiuals of REMLs in lme4 using the following: m1<-lmer(vTotal~Week+fCollar+ (1|fCat), collars) res<-resid(m1) plot(res) qqnorm(res) library(MASS) par(mfrow=c(2,3))

Leroy Tennison wrote: > And I thought it was a Kubuntu (Ubuntu with KDE desktop for those who > aren't familiar) thing! Apparently it's a KDE thing. I haven't > experienced the scrollbar aspect (or maybe I just haven't done what you > do) but my arrows are missing too. I'm thinking this is a KDE Blasted > Ugly Gotcha (BUG). BTW, if you haven't already

Hi, I have a data set with three variables,X Y and Time. X and Y are the coordinates of points, i want to join these points according to the Time sequence using arrows? Demo Example data: > x<-c(1:6) > y<-c(1:6) > time<-c(6:1) > data<-cbind(x,y,time) > data x y time [1,] 1 1 6 [2,] 2 2 5 [3,] 3 3 4 [4,] 4 4 3 [5,] 5 5 2 [6,] 6 6 1 I hope to join

On 10/12/18 8:40 AM, Leroy Tennison wrote: > And I thought it was a Kubuntu (Ubuntu with KDE desktop for those who aren't familiar) thing! Apparently it's a KDE thing. I haven't experienced the scrollbar aspect (or maybe I just haven't done what you do) but my arrows are missing too. I'm thinking this is a KDE Blasted Ugly Gotcha (BUG). BTW, if you haven't already

Leroy Tennison Network Information/Cyber Security Specialist E: leroy at datavoiceint.com 2220 Bush Dr McKinney, Texas 75070 www.datavoiceint.com TThis message has been sent on behalf of a company that is part of the Harris Operating Group of Constellation Software Inc. These companies are listed here . If you prefer not to be contacted by Harris Operating Group please notify us . This message is

On Friday 12 October 2018 12:19:40 Bill Gee wrote: > I agree that issue with the scroll bar jumping all over is really annoying! > It is actually a feature of Gnome and GTK. It can be changed by editing a > file: > > ~/.config/gtk-3.0/settings.ini > > [Settings] > gtk-primary-button-warps-slider=0 > > You may have to create this file. Firefox is a good GTK