Outubro 2011 - R-br

Associação conjunta
by Edson Lira 18 Out '11

18 Out '11

Boa noite a todos, tenho um banco de dados com algumas variávei: AR(sim, não), AIN(sim, não, talvez), T_HIV(menos de 250, entre de 250 e 300, acima de 300), HISTO_A(normal, nia, nia I). Gostaria de avaliar a associação conjunta da variável AR com o as demais variáveis conjuntamente. Seria um modelo log-linear? Alguém poderia me dar alguma sugestão ou material para eu dar uma olhada? []'s Edson Lira Estatístico Manaus-Amazonas

3 3

Como obter a estrutura de um data.frame (linux)
by Jose Claudio Faria 18 Out '11

18 Out '11

Pessoal, Alguém poderia me informar uma função (no Linux) para exportar (por exemplo para a área de transferência) a estrutura de um data.frame. No Windows existe uma função (trExport) que acompanha o pacote TinnR que faz isso e não estou encontrando (ainda) nada similar no Linux. O objetivo é poder compartilhar dados reais (data.frame) pela lista de discussão sem a necessidade de conexões, read.table, etc. O que preciso é algo assim como o resultado para o data.frame BOD que fica no pacote datasets: structure(list(Time = c(1, 2, 3, 4, 5), demand = c(8,3, 10,3, 19, 16, 15.6, 19.8), .Names = c("Time", "demand"), row.names = c(NA, -6L), class = "data.frame") Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

Média e desvio de muitas variáveis de um fatorial!
by Ivan Bezerra Allaman 18 Out '11

18 Out '11

Bom dia senhores! Tenho um experimento fatorial 14x4, com 12 variáveis, e gostaria de saber se tem como colocar de uma só vez em um único quadro a média e o desvio padrão dos tratamentos e de todas as variáveis. Fiz algo como: medias <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,mean,na.rm=TRUE)) desvios <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,sd,na.rm=TRUE)) descrit <- list() for(i in 1:12){ descrit[[i]] <- cbind(medias[,i],desvios[,i]) } Mais sai como lista não sendo o ideal que procuro. Desde já agradeço pela atenção! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 3

Dois plots no mesmo frame
by Marcelo Claro de Souza 18 Out '11

18 Out '11

Gostaria de plotar duas regressões lineares no mesmo frame. Para tal plotei os pontos utilizando os comandos abaixo: x1=data$SLA1 y1=data$LLS 1 x2=data$SLA2 y2=data$LLS2 plot(x1, y1, type='p', xlim=range(x1,x2), ylim=range(y1, y2), xlab='x', ylab='y') points(x2, y2, type='p', col="red", xlab='x', ylab='y') reg1=lm(y1~x1) reg2=lm(y2~x2) Mas ao tentar aplicar as linhas de tendência utilizando abline(reg1) e abline(reg2) não obtive resultado algum. Como devo proceder para adicionar as linhas das regressões a essas figuras? Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Strings
by Rodrigo Sant'Ana 18 Out '11

18 Out '11

Caros Companheiros da R-BR, estou com um problema, que provavelmente, seja muito simples de resolver. Mas confesso, já quebrei a cabeça para resolver e num consegui chegar ao resultado que desejo. Abaixo, tento replicar meu problema para que vocês possam me ajudar! Muito obrigado pela atenção de todos! # Eis uma réplica simplificada da base de dados! dados <- data.frame(municipio = c("A", "B", "C", "D", "E"), pesc.emb = c("PESCADOR JOAO","EMBARCACAO ZE","PESCADOR FELIPE","EMBARCACAO AMARGURA II","EMBARCACAO STAR")) dados municipio pesc.emb 1 A PESCADOR JOAO 2 B EMBARCACAO ZE 3 C PESCADOR FELIPE 4 D EMBARCACAO AMARGURA II 5 E EMBACACAO STAR # Preciso criar uma coluna classificando a informação da linha em "PESCADOR" ou "EMBARCACAO", ou seja, preciso extrair a string "PESCADOR" ou "EMBARCACAO" da coluna pesc.emb. Já tentei utilizar todas as funções do pacote stringr, porém não obtive exito, pois preciso ambos na mesma coluna. library(stringr) dados$classify <- str_extract(dados$pesc.emb, "PESCADOR") dados$classify2 <- str_extract(dados$pesc.emb, "EMBARCACAO") dados municipio pesc.emb classify classify2 1 A PESCADOR JOAO PESCADOR <NA> 2 B EMBARCACAO ZE <NA> EMBARCACAO 3 C PESCADOR FELIPE PESCADOR <NA> 4 D EMBARCACAO AMARGURA II <NA> EMBARCACAO 5 E EMBARCACAO STAR <NA> EMBARCACAO sei que pode parecer simples, mas estou enfrentando dificuldades em solucionar este problema. Fico super agradecido pela atenção de todos, abs _______________ Rodrigo Sant'Ana Oceanógrafo

6 10

Diagrama de Vernn (Package eVenn)
by Íris Cardoso 17 Out '11

17 Out '11

Olá, pessoal Primeiramente, sou nova na lista e gostaria de dizer que é muito bom encontrar um grupo com uma troca tão grande e rica de informações. E, como toda mensagem tem um "interesse" por trás, gostaria de esclarecer uma dúvida: Apanhei do R um bom tempo para conseguir plotar um diagrama de Venn (package eVenn). Sempre acusava um erro do tipo ausência de "TRUE/FALSE" nos dados. Optei por alterar o único ponto diferente dos dados dos exemplos, os nomes das variáveis. Depois de alterar os nomes para ficarem iguais aos nomes da base do exemplo "Liste_1_.xxx." (onde xxx = soma de 1's na variável), consegui plotar o diagrama. Porém ele plotou a 1ª coluna (id das observações) como uma variável do diagrama. Excluí a referida coluna e obtive o Diagrama do jeito que queria. Gostaria de saber se alguém já conseguiu plotar usando o "nome verdadeiro" das variáveis ou se essa alteração é mesmo necessária. Achei estranho pois, se o nome da variável fosse mesmo necessário nesse formato, ele não deveria plotar a coluna de id. Desde já agradeço a atenção. P.S.: A coluna de id estava com o mesmo nome do exemplo "row.names" e o conteúdo estava configurado como character. Íris Cardoso Estatística e Metre em Logística

2 2

metamodelagem
by Edmar Alino 17 Out '11

17 Out '11

Srs, Estou tentando criar metamodelos para otimizacao e me recomendaram os seguintes pacotes: - Redes Neurais (rbf e mlp): AMORE - GaSP: Dicekriging. Como sou iniciante no R, qualquer sugestao/recomendacao/orientacao pertinente a esses pacotes ou outros pacotes mais apropriados serao muito bemvindos. Grato desde ja, Edmar

2 1

Ajuda
by Paloma Palmieri 17 Out '11

17 Out '11

Prezados, Boa tarde! Estou com o seguinte problema: Tenho um programa em MATLAB que gera um arquivo .prn que é utilizado por um arquivo compilado em C (que é chamado em um outro arquivo de extensão .bat). O que acontece é que eu tenho que rodar os programas individualmente e isso não é muito eficiente. Neste caso, gostaria de ver se alguém da lista sabe como eu chamo e rodo um executável (ou .bat) em C dentro do MATLAB. Desde já agradeço. Abraços, Paloma Palmieri Alves paloma.palmieri(a)gmail.com Cel.: (21) 8811-5230 (21) 7961-8718

2 1

ordenação de data para plot
by julio cesar araujo 15 Out '11

15 Out '11

Boa noite pessoal Tenho uma série temporal de retornos de ativos e gostaria de plotar um gráfico simples com data no eixo X e os valores no eixo Y.....como existem feriados no Brasil, essa série não é de segunda a sexta. Possuo +- 2000 datas, e quando vou plotar plot(data, valor) os dados do eixo X não ficam em ordem cronologica. O que fazer para que fiquem as datas no eixo X em ordem cronologica e os valores no eixo Y em linha? Atc

2 2

Re: [R-br] ordenação de data para plot
by Ivan Bezerra Allaman 15 Out '11

15 Out '11

É provável que seus dados não esteja formatado da maneira que o R entende, ou seja, YY/MM/DD, se eu não me engano. Caso queira formatar as datas da maneira que mais lhe convém, verifique as funções, ?as.Date Para uma ajuda mais refinada nos forneça um CRM. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0