Re: [R] Repeated Measures ANOVA and Missing Values in the data set

2015-06-19 Thread Stéphane Adamowicz
]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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. Stéphane

Re: [R] predict function in regression analysis

2015-05-06 Thread Stéphane Adamowicz
, self-contained, reproducible code. Stéphane Adamowicz Chercheur / Scientist stephane.adamow...@paca.inra.fr Centre PACA - UR 1115 PSH Tel. : +33 1 (0)4 32 72 24 35 Fax : +33 1 (0)4 32 72 24 32 228, route de l'Aérodrome 84 914 Avignon Cedex 9 France www.inra.fr

Re: [R] matrix manipulation question

2015-03-31 Thread Stéphane Adamowicz
Many thanks, Stéphane Le 30 mars 2015 à 10:42, peter dalgaard pda...@gmail.com a écrit : On 30 Mar 2015, at 09:59 , Stéphane Adamowicz stephane.adamow...@avignon.inra.fr wrote: However, in order to help me understand, would you be so kind as to give me a matrix or data.frame example

Re: [R] matrix manipulation question

2015-03-30 Thread Stéphane Adamowicz
Le 27 mars 2015 � 18:01, David Winsemius dwinsem...@comcast.net a �crit : On Mar 27, 2015, at 3:41 AM, St�phane Adamowicz wrote: Well, it seems to work with me. No one is doubting that it worked for you in this instance. What Peter D. was criticizing was the construction :

Re: [R] matrix manipulation question

2015-03-27 Thread Stéphane Adamowicz
Well, it seems to work with me. Y - as.matrix(airquality) head(Y, n=8) Ozone Solar.R Wind Temp Month Day [1,]41 190 7.4 67 5 1 [2,]36 118 8.0 72 5 2 [3,]12 149 12.6 74 5 3 [4,]18 313 11.5 62 5 4 [5,]NA NA 14.3 56

Re: [R] matrix manipulation question

2015-03-27 Thread Stéphane Adamowicz
Why not use complete.cases() ? data_no_NA - data[, complete.cases(t(data))==T] Le 27 mars 2015 à 06:13, Jatin Kala jatin.kala...@gmail.com a écrit : Hi, I've got a rather large matrix of about 800 rows and 60 columns. Each column is a time-series 800 long. Out of these 60 time

Re: [R] matrix manipulation question

2015-03-27 Thread Stéphane Adamowicz
Le 27 mars 2015 à 12:34, PIKAL Petr petr.pi...@precheza.cz a écrit : Very, very, very bad solution. as.matrix can change silently your data to unwanted format, complete.cases()==T is silly as Peter already pointed out. Perhaps, but it happens that in the original message, the

Re: [R] matrix manipulation question

2015-03-27 Thread Stéphane Adamowicz
example. Furthermore in my example no unwanted format occurred. You can Yes because data.frame was (luckily) numeric. Luck has nothing to do with this. I Chose this example on purpose … Stéphane __ R-help@r-project.org mailing list -- To

Re: [R] Removing objects except user-defined functions

2015-02-19 Thread Stéphane Adamowicz
There is a function keep() in package gdata for this purpose Le 19 févr. 2015 à 10:25, philippe massicotte pmassico...@hotmail.com a écrit : Dear R users. I would like to remove all object from my workspace except the function I have defined. However, is I use rm(list = ls()) everything

Re: [R] Using line (Tukey's robust line) on 3 observations?

2014-12-05 Thread Stéphane Adamowicz
Le 4 d�c. 2014 � 13:40, Tal Galili a �crit : By accident I came across the following example: x - 1:3 y - 1:3 line(x, y) # returns: Call: line(x, x) Coefficients: [1] -2 2 While when using 1:4, it will give the more reasonable 0,1 coefficients. I imagine this is in the

[R] Rép : ANOVA and permutation tests : beware of traps

2014-09-24 Thread Stéphane Adamowicz
packages that propose randomization (or bootstrap) methods ? Cheers, Stéphane ___ Stéphane Adamowicz stephane.adamow...@avignon.inra.fr ___ UR 1115 Plantes et Systèmes de Culture Horticoles (PSH

[R] ANOVA and permutation tests : beware of traps

2014-09-23 Thread Stéphane Adamowicz
Recently, I came across a strange and potentially troublesome behaviour of the lm and aov functions that ask questions about calculation accuracy. Let us consider the 2 following datasets dat1 dat2 : (dat1 - data.frame(Y=c(1:3, 10+1:3), F=c(rep(A,3), rep(B,3 Y F 1 1 A 2 2 A 3 3 A 4

[R] minor tick marks

2010-06-09 Thread Stéphane Adamowicz
major ticks obtained with the above command. The minor.tick function in library Hmisc gives an error when applied to log axes. Any solution ? Stéphane _ Stéphane Adamowicz INRA, unité PSH domaine St Paul, site agroparc 84914 Avignon, cedex 9 France