And... it works !! Thanks again !!
-Etienne
Gregory Warnes a écrit :
>
> Hi Etienne,
>
> The basic problem is that under the default conversion mode
> (BASIC_CONVERSION) all R objects are converted to roughly-equivalent
> python structures. As a consequence, the object 'av' isn't actually
> an R object, so r.summary(av) won't treat it as such. The simplest
> solution is to change the conversion mode to NO_CONVERSION adn then
> explicitly request conversion of an object when you need the python
> version. IE:
>
>> set_default_mode(NO_CONVERSION)
>> v = r.aov(r.formula("score~factor+Error(id_subject/factor)"),
>> data=Res)set_default_mode(BASIC_CONVERSION)
>> set_default_mode(BASIC_CONVERSION)
>> r.summary(av)
>
> -G
>
> On Dec 18, 2007, at 1:52PM , Etienne Gaudrain wrote:
>
>> Hi everyone,
>>
>> I'm new to RPy, and I came to this terrific module as I was used to
>> make some of my analyses in R, and I came to Python in replacement of
>> Matlab. Formerly, I manipulated data with Matlab, put it in a MySQL
>> database, and made my stats in R via ODBC. I'm now thinking about
>> jumping one step by calling R directly from Python with RPy.
>>
>> The analysis I almost always have to do is a repeated measure ANOVA.
>> The way I do this in R is :
>>
>>
>> # /after odbc connection and sql query, Res contains my data
>> /
>> library('stats')
>> av <- aov( score~factor+Error(id_subject/factor), data=Res)
>> summary(av)
>>
>>
>> Now I tried the same in RPy :
>>
>>
>> # /retrieve data from sql query, Res is a dictionnary
>> /
>> r.library('stats')
>> av = r.aov("score~factor+Error(id_subject/factor)", data=Res)
>>
>>
>> This fails saying that "Error" isn't defined in the dataframe...
>> After reading some R doc about GLM, I found that using the R function
>> formula() seemed to solve this problem:
>>
>> av = r.aov(r.formula("score~factor+Error(id_subject/factor)"), data=Res)
>> r.summary(av)
>>
>> However, a new problem rose in r.summary(). This function returns
>> something that isn't readable, and that does not contain the p
>> values, or anything similar. It seems that the r.summary_aov()
>> function might be adequat, but this function returns an Error saying
>> that there is a NaN somewhere...
>>
>> Does anybody have an advice on how to perform the repeated measure ANOVA?
>> Thanks!
>> -Etienne
>>
>>
>> PS : I use Windows XP, Python 2.5.1, Numpy 1.0.3.1 and RPy
>> 1.0.1-Numpy-py2.5 and R 2.6.1.
>>
>> -------------------------------------------------------------------------
>> SF.Net email is sponsored by:
>> Check out the new SourceForge.net Marketplace.
>> It's the best place to buy or sell services
>> for just about anything Open Source.
>> http://ad.doubleclick.net/clk;164216239;13503038;w?http://sf.net/marketplace<mime-attachment.txt>
>
> ------------------------------------------------------------------------
>
> -------------------------------------------------------------------------
> SF.Net email is sponsored by:
> Check out the new SourceForge.net Marketplace.
> It's the best place to buy or sell services
> for just about anything Open Source.
> http://ad.doubleclick.net/clk;164216239;13503038;w?http://sf.net/marketplace
> ------------------------------------------------------------------------
>
> _______________________________________________
> rpy-list mailing list
> rpy-list@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/rpy-list
>
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Etienne Gaudrain
Universite Claude Bernard LYON 1
CNRS - UMR5020, Neurosciences Sensorielles, Comportement, Cognition
50, avenue Tony Garnier
69366 LYON Cedex 07
FRANCE
Tél : 04 37 28 74 85
Fax : 04 37 28 76 01
Page web equipe : http://olfac.univ-lyon1.fr/unite/equipe-02/equipe-2-f.html
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-------------------------------------------------------------------------
SF.Net email is sponsored by:
Check out the new SourceForge.net Marketplace.
It's the best place to buy or sell services
for just about anything Open Source.
http://ad.doubleclick.net/clk;164216239;13503038;w?http://sf.net/marketplace
_______________________________________________
rpy-list mailing list
rpy-list@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rpy-list
