[R] glmulti fails because of rJava

2011-11-23 Thread Toby Marthews
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.

Re: [R] barplot with errorbars

2011-02-17 Thread Toby Marthews
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

Re: [R] Linear regressions: producing multiple outputs

2011-02-17 Thread Toby Marthews
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

Re: [R] add error bars in a plot

2011-02-09 Thread Toby Marthews
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

Re: [R] how to learn more from R

2011-02-04 Thread Toby Marthews
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

Re: [R] leap year and order function

2011-01-31 Thread Toby Marthews
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...

Re: [R] Using anova() with glmmPQL()

2011-01-20 Thread Toby Marthews
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)

Re: [R] Statistical formulas

2011-01-18 Thread Toby Marthews
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

[R] Using anova() with glmmPQL()

2011-01-17 Thread Toby Marthews
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

[R] nls() fails on a simple exponential fit, when lm() gets it right?

2008-08-29 Thread Toby Marthews
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

r-help@r-project.org

2008-06-18 Thread Toby Marthews
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

Re: [R] Calling functions

2008-06-17 Thread Toby Marthews
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

[R] How can I stop strwrap removing escape characters? (and solution to the screen wrapping question from 11/6)

2008-06-16 Thread Toby Marthews
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

Re: [R] round(1.5) = round(2.5) = 2?

2008-06-15 Thread Toby Marthews
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 :

Re: [R] Level Plot and Scale of Colorkey

2008-06-13 Thread Toby Marthews
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

Re: [R] controlling location of labels in axis()

2008-06-12 Thread Toby Marthews
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() >

[R] Word wrapping for character objects (WINDOWS R ONLY)

2008-06-11 Thread Toby Marthews
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

[R] Constructing groupedData objects in nlme - a little problem

2008-06-04 Thread Toby Marthews
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