Dear R,
The glmulti package no longer loads through the library() command, apparently
because of a problem with rJava.
I have today reinstalled R from scratch (updated to v2.14.0) and reinstalled
all packages from scratch and updated them all too. The problem is the same as
I found on v2.13.2.
If you google "barplot with error bars" you immediately find
http://svitsrv25.epfl.ch/R-doc/library/prada/html/barploterrbar.html .
Toby.
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Lathouri, Maria [m.lathour...@imperial
Hi RTSlider,
I suspect you rather need to use the lme command (or perhaps glmmPQL or lmer)
because you have a random predictor?
lme(fixed=LeafLength~AirTemp*SnowFreeDate,random=~1|Species)
See
http://socserv.mcmaster.ca/jfox/Books/Companion-1E/appendix-mixed-models.pdf
for a tutorial on lm
Hi Maria,
Instead of using errbar(), you can use the segments command which may be more
flexible. Following the example on p.80 (the errbar command) of
http://cran.ma.ic.ac.uk/web/packages/Hmisc/Hmisc.pdf
set.seed(1)
x=1:10
y=x+rnorm(10)
delta=runif(10)
barwidth=0.015
x11();plot(x,y)
segments
h lots of links for
further information.
Toby Marthews
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Jie TANG [totang...@gmail.com]
Sent: 04 February 2011 06:42
To: r-help@r-project.org
Subject: [R] how to learn more from
Dear All,
I've always used this code:
year=c(1948:1953,2000,2100,2200,2300)
numdays=ifelse((year%%4==0 & year%%100!=0) | year%%400==0,366,365)
> numdays
[1] 366 365 365 365 366 365 366 365 365 365
Toby
From: r-help-boun...@r-project.org [r-help-boun...
with glmmPQL()
Toby Marthews ouce.ox.ac.uk> writes:
>
> Dear R HELP,
>
> ABOUT glmmPQL and the anova command. Here is an example of a
> repeated-measures ANOVA focussing on the way
> starling masses vary according to (i) roost situation and
> (ii) time (two time points only)
ld, but I still recommend Fowler et
al. (1998) for those who are starting out with statistical testing: Fowler J,
Cohen L & Jarvis P (1998). Practical statistics for field biology (2nd ed.).
Wiley, Chichester, UK.
Toby Marthews
From: r-help-boun...@r-proje
didrog+state:results
):
class(lmeres2)="lme"
anova(object=lmeres2,test="Chisq")
numDF denDF F-value p-value
(Intercept)136 182.84356 <.0001
mnth 136 164.57288 <.0001
roostsitu 336 17.79263 <.0001
mnth:r
he same exp and ln fits using MS Excel and the numbers returned are
exactly as for the 'simpler ways', so nls is definitely in the wrong, but
the self-start is almost identical to the one for the logarithmic fit,
which works perfectly
I can't see my mistake at all. Can anyone
logical OR. The shorter form performs
elementwise comparisons in much the same way as arithmetic operators. The
longer form evaluates left to right examining only the first element of
each vector. Evaluation proceeds only until the result is determined."
Hope this helps,
Toby Marthews
L
orkspace loaded?
Toby Marthews
Le Mar 17 juin 2008 13:32, Michael Pearmain a écrit :
> Another newbie question.
>
> I've written a function and saved the file as Xtabs.R, in a central place
> on
> a network so others will be able ot use the function,
> My question is how do i
t;\n\tHello.\n")
Hello.
>
Is there any way to prevent strwrap doing this? I haven't found any
options for this.
Thanks very much (and thanks to Greg Snow for the tip about strwrap before).
Toby Marthews
==
Le Mer 11 juin 2008 18:00, Gre
a bit odd. However, R's benefits and advantages by far outweigh its little
oddities, as I presume you know since you are using it.
Effectively, I never use the round() command and always calculate using
the floor function.
Toby Marthews
Le Dim 15 juin 2008 11:26, Markus Didion a écrit :
Try
colscaledivs=100#colscaledivs=15 here is the R default
levelplot(z ~ x * y, g,xlab="x co-ordinate",ylab="y
co-ordinate",colorkey=TRUE,at=seq(from=-0.01,to=0.25,length=colscaledivs),col.regions=(col=gray((0:colscaledivs)/colscaledivs)))
Toby Marthews
Le
quot;First outer y axis label",side=2,outer=TRUE)
mtext("Second outer x axis label",side=3,outer=TRUE)
mtext("Second outer y axis label",side=4,outer=TRUE)
Cheers,
Toby Marthews
Le Jeu 12 juin 2008 15:32, Andrew Yee a écrit :
> Here's a naive question about axis()
>
I get word wrapping like this in Windows? Perhaps the attached
screenshots clarify this question.
THANKS FOR ANY HELP!
Toby Marthews
Previous relevant posts:
- The post from 2006 about Screen Wrapping
(http://tolstoy.newcastle.edu.au/R/help/06/05/26673.html) which Brian
Ripley answered was a
Dear R-help,
I am trying to create groupedData objects using the nlme library. I'm
missing something basic, I know:
Here is the first example in ch.1 of Pinheiro & Bates (2000):
library(nlme)
x2=Rail$travel;x1=Rail$Rail;eg1=data.frame(x1,x2);eg1gd=Rail
print(eg1gd)
x11();print(plot(eg1gd))
femod
18 matches
Mail list logo