similar to: different results from lme() and lmer()

Hi, I have this: lme(y~x1+x2,random=~1|x1/x2) How to make this random effect using lmer? I try this: lmer(y~x1+x2+(1|x1/x2) But it dont work. Any idea? Thanks Ronaldo -- System halted! -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia ] | / \ [36570-000 Vi??osa -

Sorry for bothering you - but I have a little problem. I`m completely new in R and trying to use lmer. That works, but I m not sure if my code is right, so I would just like to get an opinion from the experts. I have a Growthrate of plants (RGR Fixed Factors: Fertilizer, Status (alien or native plants) and their interaction Random Factors: Age is a covariable; than I have plantfamily and

Dear R Community, Below is a problem with a simple ggplot2 graph. The code returns the error message below. Error: stat_count() must not be used with a y aesthetic. My code is below and the data is attached as a ?text? file. # Graph of the probabilities library(digest) library(DT) datatable(probability) str(probability) probability$Fertilizer <- as.factor(probability$Fertilizer)

Hi James, If you want to specify the y-values, you need to use stat="identity" as below: ggplot(probability, aes(x=Fertilizer, y=prob)) + geom_bar(stat="identity", aes(fill=Treatment)) best, huzefa On Tue, Apr 12, 2016 at 1:02 PM, James Henson <jfhenson1 at gmail.com> wrote: > Dear R Community, > > Below is a problem with a simple ggplot2 graph. The code

Hi all, I analysed my data with lme and after that I spent a lot of time for mean separation of treatments (post hoc). But still I couldn’t make through it. This is my data set and R scripts I tried. replication fertilizer variety plot height 1 level1 var1 1504 52 1 level1 var3 1506 59 1 level1 var4 1509 54 1 level1 var2 1510 48 2 level1 var1 2604 47 2 level1 var4 2606 51 2 level1 var3

Dear R folk, First off I want to thank those of you who responded with comments for my R quick and dirty stats tutorial. They've been quite helpful, and I'm in the process of revising them. When it comes to repeated measures MANOVA, I'm in a bit of a bind, however. I'm beginning to see that all of the documentation is written for psychologists, who have a slightly

Hi everyone! I tried to look up in the previous topics whether there was something similar to my question but I think there's nothing corresponding... Anyway! Here is my problem: I have a plot (qplot()) in which I want to insert axis labels with the following text: kg NH4+ N-equiv. * ha-1 yr-1 (well... 4 smaller and down and + as a power (the same for the -1) but I am sure you understood

r-help at stat.math.ethz.ch on Wednesday, October 26, 2005 at 6:00 AM -0500 wrote: Ronaldo, Try Harold's suggestion. The df still won't agree, because lmer (at least in its current version) just puts an upper bound on the df. But that should be OK, because all those t tests are approximations anyways, and you can get better confidence intervals (credible intervals, whatever) by using the

I would have expected the function getAnywhere to have behaved differently in the following: > search() [1] ".GlobalEnv" "file:C:/R/Rdata/miya/.Rdata" [3] "package:boot" "package:methods" [5] "package:ctest" "package:mva" [7] "package:modreg"

Hello everyone I am currently trying to conduct analysis of my graduate thesis data using a mixed effects model and I have reached an impass. When I try to conduct a multiple comparison, I get an error (See below): > fm3<- lme(abovegroundbiomass.m.2~medium*amelioration*fertilizer*treatment, random=~1|block/medium/amelioration/fertilizer) > tukeytest<-glht(fm3,

Hi there, I have a experimental design related question. I have done a experiment with 3 factors. The design matrix is similar to: > data.frame(Factor.1 = rep(rep(0:2, each = 3),3), Factor.2 = rep(c("U","S","N"), 9), Factor.3 = rep(0:2, each = 9)) Factor.1 Factor.2 Factor.3 1 0 U 0 2 0 S 0 3 0 N

Hello everyone. I am trying to plot data from a time series in R and have run across some problems. I was wondering if someone could help. I have data taken every 15 minutes of a list of rankings of 25 titles. I want to plot one of the article's titles's ranks as a time series. I've been trying if statements with no success and am pretty new to R. Can anyone help? Thank you. -- View

Hi everyone, I'm trying to plot some data across time. I have a list of articles, ranks, date/times, authors, etc. Someone suggested using zoo and someone suggested using ts. I'm pretty new at this and have been trying a simple if() plot() statement, but it doesn't seem to work. I keep getting an error that Title does not exist. Basically its saying that the headers do not exist as

Hi! How to calculate the correct SE of mean in a nested or spliplot anova? Nested example: --------------------- m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver)) > m Call: aov(formula = Glycogen ~ Treatment + Error(Treatment/Rat/Liver)) Grand Mean: 142.2222 Stratum 1: Treatment Terms: Treatment Sum of Squares 1557.556 Deg. of Freedom 2 Estimated

I'm not sure if this is a bug, or if I'm doing something wrong. =20 =46rom the worms dataframe, which is at in a file called worms.txt at =20 http://www.imperial.ac.uk/bio/research/crawley/therbook <http://www.imperial.ac.uk/bio/research/mjcraw/therbook/index.htm>=20 =20 the idea is to extract a subset of the rows, sorted in declining order of worm density, with only the maximum

Hello I am analysing aboveground biomass data from revegetation testplots which I constructed in a split-plot design using the function lme. For the experiment, the three factors are amelioration (2 levels), fertilizer (2 levels) and treatment (7 levels). Each testplot (block) has a singlereplicate of each treatment (total of 8 testplots). The blocks were constructed of topsoil. Each block was

Hello Everyone,   Does anyone have the R code for this book? I contacted the author about this and was told that it used to be on the book's website, but that the publisher had asked that it be taken down. I'm hoping that someone will have downloaded the code when it was available, and that they will be willing to share it with me.   The book is almost 1,000 pages long and I'd like to

Hi,) I am just getting started with R but have hit an early snag. I am working through Crawley (2008) The R Book and on page 6, 'Significance Stars', I am trying to enter the commands given. However, 'Gain.txt' does not seem to have been downloaded when I downloaded the R programme. I have searched for it in case it put it somewhere other than 'temp' but I cannot find it.

Windows 2000 . 5.00.2195 with Service Pack 1. R 1.5.1 Output from my split-split plot aov "alerted" me that I have done something wrong. I designed an experiment with all combinations of all levels of each treatment, but lost a little data (3 out of 192 plots). With the following data, I run the following model: > collim[c(1:6,187:192),c(1,3:6,9)] plot Litter Fert

I was surprised to find that just changing the base level of a factor variable changed the number of significant coefficients in the solution. I was surprised at this and want to know how I should choose the order of the factors, if the order affects the result. Here is the small example. It is taken from 'The R Book', Crawley p. 365. The data is at