[R] please Help me on a repeated measures anova

David Winsemius dwinsemius at comcast.net
Thu Dec 30 16:22:26 CET 2010


On Dec 30, 2010, at 10:11 AM, David Winsemius wrote:

>
> On Dec 30, 2010, at 5:38 AM, soileil wrote:
>
>>
>> So far I've done this:
>>
>>
>>> tableau = read.table("Classeur1.txt", h=T)
>>> tableau
>>   Temps.s. Modalite  AchE
>> 1         15       ER   0,1
>>
> snipped data
>> 103      120   SED+ED  0,19
>> 104      120   SED+ED 0,172
>>> attach (tableau)
>>
>>
>>> tableau <- lm (AchE ~ Temps.s. * Modalite, random = ~ 1 | individu)
>> Erreur dans storage.mode(y) <- "double" :
>> la modification du mode de stockage
>
> Even without an English translation, it seems very likely that the  
> AchE variable is not of numeric class,  but rather a factor. You may  
> want to redo your data input step with a setting for the decimal  
> separator set properly for the conventions of your locale settings.  
> See help(read.table)

And the other _big_  problem that I did not notice at first is that  
you are overwriting your data with the model output. This isn't  
exactly the problem addressed by fortune("dog") but an even more  
severe case of identity conflation. This is even a further example of  
why it is a bad idea to use the attach function.

>
> You may also want to look at the settings for locale if you will  
> typically be using "," as a decimal separator.
>
> ?options
> ?locales
>
>
>> d'un objet 'factor' n'est pas
>> autorisée
>> De plus : Messages d'avis :
>> 1: In model.response(mf, "numeric") :
>> l'utilisation de type="numeric" avec une réponse de type facteur sera
>> ignorée
>> 2: In lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
>> les arguments surnuméraires random sont ignorés.
>>
>> What is my mistake?
>> -- 
>> View this message in context: http://r.789695.n4.nabble.com/please-Help-me-on-a-repeated-measures-anova-tp3159868p3168134.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
> -- 
>
> David Winsemius, MD
> West Hartford, CT
>
> ______________________________________________
> R-help at 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.

David Winsemius, MD
West Hartford, CT



More information about the R-help mailing list