[R] nested (hierarchical) anova

2010-04-15 Thread Héctor Villalobos
=FALSE) # P ## or from the second row in the two way anova: summary( aov(cChol ~ drug + dsource, data=zar) ) In order to avoid confussion in our students (and myself!) I wonder if we can obtain the desired results in one step? Thank you in advance for your kind help. Héctor -- Héctor

Re: [R] coping with a warning in if()

2010-01-20 Thread Héctor Villalobos
Thanks to all who responded. I've learn a lot from your comments. The best solution I can see is to use both 'missing()' and 'match.arg()' in my function. Héctor Hi, I'm sure this one is very easy I am trying to write a function where one of its arguments has two posible (strings)

[R] coping with a warning in if()

2010-01-19 Thread Héctor Villalobos
)) ){ stop(specify type of result 'simple' or 'complete') } else { res - rnorm(n) res } Thanks for your time Héctor -- Héctor Villalobos hvill...@ipn.mx CICIMAR - IPN A.P. 592. Col. Centro La Paz, Baja California Sur, MÉXICO. 23000 Tels. (+52 612) 122 53 44; 123 46 58; 123 47 34 ext

Re: [R] problems installing R packages

2009-12-02 Thread Héctor Villalobos
. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Héctor Villalobos hvill

Re: [R] package installation error

2009-09-01 Thread Héctor Villalobos
-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Héctor Villalobos hvill...@ipn.mx CICIMAR - IPN A.P. 592. Col. Centro La Paz, Baja California Sur, MÉXICO. 23000 Tels. (+52 612) 122 53 44; 123 46 58; 123 47 34 ext. 82425 Fax. (+52 612) 122

[R] manipulating text to generate different formulas to use in nls()

2009-08-10 Thread Héctor Villalobos
$period[3], f4=1/inip$period[4]) ) ## not run. These two will be used in res - nls(modl, data=data.frame(y, x), start=initP) Thank you in advance for any insigth. Héctor -- Héctor Villalobos hvill...@ipn.mx CICIMAR - IPN A.P. 592. Col. Centro La Paz, Baja California Sur, MÉXICO

Re: [R] biplot.princomp - changing score labels

2009-02-25 Thread Héctor Villalobos
the PCs. Anyway, thanks for the comment - and for designing biplot.princomp. Axel -- Héctor Villalobos hvill...@ipn.mx CICIMAR - IPN A.P. 592. Col. Centro La Paz, Baja California Sur, MÉXICO. 23000 Tels. (+52 612) 122 53 44; 123 46 58; 123 47 34 ext. 82425 Fax. (+52 612) 122 53 22

Re: [R] plotting arrows with different colors and varying head size

2009-01-20 Thread Héctor Villalobos
Thanks to Jim and Greg, Concerning the colors of the arrows, color.scale() function from Jim's solution seems more straigthforward to me. I'm trying now to include a proper legend with color.legend() function. Héctor On 20 Jan 2009 at 22:13, Jim Lemon wrote: Héctor Villalobos wrote: Dear

[R] plotting arrows with different colors and varying head size

2009-01-19 Thread Héctor Villalobos
=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods [7] base other attached packages: [1] car_1.2-9 hdf5_1.6.7 -- Héctor Villalobos hvill...@ipn.mx CICIMAR - IPN La Paz, Baja California Sur, MÉXICO [[alternative HTML version

Re: [R] as.numeric in data.frame, but only where it is possible

2008-11-27 Thread Héctor Villalobos
, reproducible code. -- Héctor Villalobos [EMAIL PROTECTED] CICIMAR - IPN La Paz, Baja California Sur, MÉXICO [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] nls() vs lm() estimates

2008-06-13 Thread Héctor Villalobos
Hi, I'm trying to understand why the coefficients a and b for the model: W = a*L^b estimated via nls() differs from those obtained for the log transformed model: log(W) = log(a) + b*log(L) estimated via lm(). Also, if I didn't make a mistake, R-squared suggests a better adjustment for the