Receita: Usando o RPy

Um exemplo simples de uso do módulo RPy para realizar testes estatísticos utilizando-se do pacote estatístico R. O conjunto de dados utilizado no exemplo é um arquivo tipo CSV comum, onde a primeira linha contém o nome dos dados.

Código

from rpy import *
###############Carrega os dados###################
#carrega o dataframe com seus dados
r('ed<-read.table("dados.csv", row.names=1,head=T,dec=".")')
#"disponibiliza" as colunas do dataframe para acesso direto
r('attach(ed)')
#################Testa normalidade################
#direciona a saida para o arquivo "qqnorm_rpy.png"
r('png("qqnorm_rpy.png")')
#testa a normalidade pelo qqnorm
r('qqnorm(altura)')
r('qqline(altura)')
#retorna a saida para o console
r('dev.off()')
################Gera um boxplot###################
#direciona a saida para o arquivo "boxplot_rpy.png"
r('png("boxplot_rpy.png")')
#gera o boxplot da altura para os diferentes sexos
r('boxplot(altura[sexo==0], altura[sexo==1], names=c("Feminino", "Masculino"))')
#retorna a saida para o console
r('dev.off()')
################Realiza um test T-student###########
#realiza o teste T-student para o modelo altura~sexo
print "\n\nTest T student \n\n"
test_t =r("t.test(altura~sexo)")
## Outra maneira de realizar um test T: r.test_t()
#imprime na tela os resultados do test T
for titulo, dados in test_t.iteritems():
    print titulo, dados

• Frederico G. Colombo Arnoldi