[R] Linear mixed effect model doubt:

2013-07-02 Thread Paqui GG
We have a doubt on whether we are applying a lme model correctly. Our experimental/sampling design is as follows: We are studying the relationship between two quantitative and continuous variables at different sampling stations along a transect in the ocean. At each station, we take three w

Re: [R] Linear mixed-effect models and model selection

2012-07-24 Thread David_txertudi
Not sure about the first question without knowing more about your model and research aims. As for the second, there are a number of methods that can be used to select models--F tests and other variance comparisons being among the most common. Given your examiner's comment about parsimony, I'm t

[R] Linear mixed-effect models and model selection

2012-07-23 Thread fariba moslih
Hi, I am looking at the effect of allelochemicals produced by two freshwater macrophyte species on two different algal species at different days. I am comparing the effect of each macrophyte on each algae at each day. I received help from someone doing the LMEM (Linear mixed-effect models) and

Re: [R] LINEAR MIXED EFFECT

2009-11-13 Thread tlumley
On Fri, 13 Nov 2009, milton ruser wrote: Hi Ana, you did again :-) require(fortunes) fortune("dog") Firstly, don't call your matrix 'matrix'. Would you call your dog 'dog'? Anyway, it might clash with the function 'matrix'. -- Barry Rowlingson R-help (October 2004) if you call your dat

Re: [R] LINEAR MIXED EFFECT

2009-11-13 Thread milton ruser
Hi Ana, you did again :-) > require(fortunes) > fortune("dog") Firstly, don't call your matrix 'matrix'. Would you call your dog 'dog'? Anyway, it might clash with the function 'matrix'. -- Barry Rowlingson R-help (October 2004) if you call your data "data" you will crach data() function.

Re: [R] LINEAR MIXED EFFECT

2009-11-13 Thread ANARPCG
Douglas Bates-2 wrote: > > On Thu, Nov 12, 2009 at 10:14 AM, milton ruser > wrote: >> Hi Ana, >> >> I am not quite sure if it is the problem, but if you call your data.frame >> as >> exp, >> you will crash exp() function... try use another name for your >> data.frame. >> >> By the way, I sugge

Re: [R] LINEAR MIXED EFFECT

2009-11-12 Thread Douglas Bates
On Thu, Nov 12, 2009 at 10:14 AM, milton ruser wrote: > Hi Ana, > > I am not quite sure if it is the problem, but if you call your data.frame as > exp, > you will crash exp() function... try use another name for your data.frame. > > By the way, I suggest you not use attach(). > > Try something lik

Re: [R] LINEAR MIXED EFFECT

2009-11-12 Thread milton ruser
Hi Ana, I am not quite sure if it is the problem, but if you call your data.frame as exp, you will crash exp() function... try use another name for your data.frame. By the way, I suggest you not use attach(). Try something like model<-lme(weight~date*diet, data=yourdataframe) good luck milton

Re: [R] LINEAR MIXED EFFECT

2009-11-12 Thread ANARPCG
Milton's point is dead-on, and I would highly encourage you to give the posting guide a look. That said... you might try "na.action = na.omit" in your call to... actually, we don't know what function you are using (see first point). Regardless, sounds like you have missing data and na.action

Re: [R] LINEAR MIXED EFFECT

2009-11-11 Thread Dave Atkins
Milton's point is dead-on, and I would highly encourage you to give the posting guide a look. That said... you might try "na.action = na.omit" in your call to... actually, we don't know what function you are using (see first point). Regardless, sounds like you have missing data and na.action is

Re: [R] LINEAR MIXED EFFECT

2009-11-11 Thread milton ruser
Dear Ana Golveia, It is completelly impossible someone realise what kind or help you need or what is happening. I suggest you give a look on the posting guide, mainly that part about a minimum reproducible code with self explaining information, etc. Cheers milton On Wed, Nov 11, 2009 at 7:22 AM

[R] LINEAR MIXED EFFECT

2009-11-11 Thread ANARPCG
CAN ANYONE PLEASE HELP ME WITH THIS i HAVE TO DO A MIXED EFFECT LINEAR MODEL WITH MY DATA DUE TO THE FACT THAT I have pseudoreplication! Although after reading and trying it for several times can get around due to "Error in na.fail.default(list(date = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, : missing v