Displaying 20 results from an estimated 1000 matches similar to: "R-help Digest, Vol 32, Issue 26"
2005 Oct 24
3
lme and lmer syntax
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 -
2011 Jun 01
1
different results from lme() and lmer()
Hello R-help,
I'm studying an example in the R book.?
The data file is available from the link below.http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/fertilizer.txt
Could you explain Why the results from lme() and lmer() are different in the following case? In other examples, I can get the same results using the two functions, but not here...?
Thank you.Miya
library(lme4)library(nlme)#
2003 Oct 06
1
Re: Use of the Foreign package to import Stata files
Long ago (Sat, 2 Nov 2002), Bill Hart <w.hart at sbcglobal.net> wrote:
> An R newbie here. I am using R 1.6 currently and have
> (successfully, I think) installed the Foreign package.
> Tried to import a data file created with Stata 7.0
> SE. Had minor problems with syntax then R decided
> that my file was not really a Stata file. It rejected
> the file saying
2006 Jan 27
1
monochrome mosaic plot in vcd package
helpeRs,
I have a nice looking mosaic plot in an article to be published
soon. Sadly, the published version will be in black and white and so ruin
the advantage of the default shading scheme of tiles.
What would readers suggest as an alternative shading scheme? If I have a
black-and-white shading scheme graduated according to suitable cutoffs I
won't be able to tell positive from
2005 Jun 02
0
How to calculate the correct SE in a nested or spliplot anova?
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
2006 May 13
0
Maximum likelihood estimate ofbivariatevonmises-weibulldistribut
r-help at stat.math.ethz.ch on Saturday, May 13, 2006 at 6:00 AM -0500 wrote:
>>One of my friends recently wrote his PhD thesis from University of
>>Leeds under Kanti Mardia's direction.
I bet your friend was really angling for that.
--
Alan B. Cobo-Lewis, Ph.D. (207) 581-3840 tel
Department of Psychology (207) 581-6128 fax
University of Maine
Orono, ME 04469-5742 alanc
2005 Sep 22
3
anova on binomial LMER objects
Dear R users,
I have been having problems getting believable estimates from anova on a
model fit from lmer. I get the impression that F is being greatly
underestimated, as can be seen by running the example I have given below.
First an explanation of what I'm trying to do. I am trying to fit a glmm
with binomial errors to some data. The experiment involves 10
shadehouses, divided between
2009 Mar 23
1
Syntax in lmer...
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
2016 Apr 12
2
ggplot2
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)
2016 Apr 12
0
ggplot2
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
2011 Jan 19
3
lme-post hoc
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
2005 Nov 15
1
Repeates Measures MANOVA for Time*Treatment Interactions
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
2005 Apr 25
2
residuals in lmer
Does anyone know how to extract residuals in lmer?
Here's the error I get:
>
crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep:
pesticide)+(1|rep:pesticide:irrigation), crop.data)
> qqnorm(crop.lme)
Error in qqnorm.default(crop.lme) : y is empty or has only NAs
> resid(crop.lme)
NULL
Thanks!
--Jake
2011 Sep 22
1
How to do Multiple Comparisons for a Mixed Effects Model
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,
2006 Nov 22
1
differences between aov and lme
Hi,
we have a split-plot experiment in which we measured the yield of crop
fields. The factors we studied were:
B : 3 blocks
I : 2 main plots for presence of Irrigation
V : 2 plots for Varieties
N : 3 levels of Nitrogen
Each block contains two plots (irrigated or not) . Each plot is divided
into two secondary parcels for the two varieties.
Each of these parcels is divided into three subplots
2017 Jun 26
1
How to code the factor for ANOVA
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
2011 Jul 16
2
expression() and axis labels
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
2011 Jan 11
1
sorting question
Hi,
I have a data frame with variables a, b, c (character vars) and t (time var,
could be represented as POSIXct or character, depending on which is most useful.
The format is "yyyy-mm-dd hh:mm:ss CET"). Now, I want to sort the data frame in
ascending order by a, b, c and then in descending order by t.
Here's what I've got, but I'm not sure how to put the
2010 Jul 22
1
check menu button (tcltk)
Hi,
I am making a mock user interface in tcltk and I would like to add a 'check menu button' such as shown here: http://zetcode.com/tutorials/pygtktutorial/images/checkmenuitem.png
Does anybody know how to do this? I am quite new to R.
Cheers!!
Albert-Jan
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
All right, but apart from the sanitation, the medicine,
2010 Apr 29
1
UpdateLinks = FALSE
Hi,
I'm reading 100s of excel files and many of them contain links to external files (I hate that, but that aside). Every time such a file is opened, a menu pops up asking if I want to update the links. I never want to update the links. I used the macro recorder to see what code would be needed to suppress that message, but to no avail (I tried more variations, but one attempt is shown