- R-br - DINF UFPR

[Dúvida] Vantagens de usar Linux
by Pedro Rafael 17 Jul '11

17 Jul '11

Pessoal estou pensando seriamente sair do Windows para Linux. Sou estatísticos e trabalho com estatística, onde o que faço no trabalho corresponde a 50% usar o R, 25% usar o Excel e 25% escrever relatórios. Quase nunca, para não dizer nunca, necessitei de outros programas da estatística fora o R e o Gretl. A única coisa que tenho medo é de meus arquivos ou programas não estarem compatíveis para os outros usuários. Uma macro no Excel certamente não irá rodar no Calc e o contrário também é verdadeiro. O bom ganho que eu irei ter no Linux será a liberdade. Acredito que o R fica mais "disposto" trabalhando sobre um OS *UNIX o que também é um ganho legal. Alguém pode levantar prós e contras de ter o Linux para trabalhos estatísticos? Saudações, Pedro Rafael Diniz Marinho

3 2

sugestao aos que assinam a lista em modo 'email diario'
by Benilton Carvalho 15 Jul '11

15 Jul '11

Se assinantes da lista preferem assinar a lista em modo 'email diario' para evitar o acumulo de mensagem em suas Caixas de Entrada, segue uma sugestao que pode ser obvia para muitos: - Mude as configuracoes de sua assinatura para modo 'Email individual'; - No seu cliente de email (GMail, Thunderbird, etc), crie uma nova "pasta" (que recebera' os emails da lista); - Crie um filtro que para mover todas as mensagens enderecadas a lista para a pasta recem-criada (filtrar pelos destinatarios em "To:" e "Cc:" e' suficiente - "r-br(a)listas.c3sl.ufpr.br") Dessa maneira, voces podem ler as mensagens conforme segundo sua conveniencia e **importante** manter a organizacao por "threads/assunto". (notem que ao responder um 'email diario', todos os usuarios que fazem uso da organizacao por threads acabam prejudicados) b

3 6

Dúvida: Como responder correatmente emails do grupo
by Jose Claudio Faria 15 Jul '11

15 Jul '11

Pessoal, Uso o Gmail para responder emails do grupo, pois ainda não é possível via Nabble e não uso nenhum cliente de amail no meu computador. Antes usava o Thanderbird, mas ultimamente tenho preferido acesso direto na web via interface do Gmail. Contudo tenho uma dúvida, pois não tenho conseguido manter minha resposta agrupadas aos tópicos originai (thread). Por exemplo, postei recentemente um email com o seguinte assunto: Ajuda: tratamento de erros no R Gustavo Carvalho e Henrique Dallazuanna responderam e suas respostas ficaram agrupadas a pergunta original. OK! Gostaria de agradecer a ambos e mater minha resposta agrupada no tópico original, enviada por mim. Já testei as duas opções abaixo: Re: [R-br] Ajuda: tratamento de erros no R Re: Ajuda: tratamento de erros no R Mas as duas opções geram na lista minha resposta como um novo tópico, o que não é desejável. Alguma dica? Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 1

Eclipse + R, interface gráfica para R
by Roney Fraga Souza 15 Jul '11

15 Jul '11

Olá a todos. Essa não é uma dúvida, mas uma sugestão para aqueles que desejam trabalhar no R em um ambiente que permita aumentar a produtividade. Já analisando o resultado final, veja as seguintes imagens: http://i.imgur.com/twBWw.png http://i.imgur.com/CzuLC.png É necessário integrar o R a Integrated Development Environment (IDE) Eclipse, que se dá via o plugin StatET. O processo de instalação é chato mas o resultado final é satisfatório. O R fica com a mesma "cara" do Matlab. Segue dois links que mostram o passo a passo da instalação. http://navisan.com/Articles/EclipseRHTML.aspx http://lukemiller.org/index.php/2010/12/r-2-12-0-and-eclipse-with-statet-in… Vale notar que o processo é para windows, mas fiz no mac e está funcionando bem. Att Roney

1 0

Re: [R-br] Ajuda: tratamento de erros no R
by Jose Claudio Faria 15 Jul '11

15 Jul '11

Obrigado ao Gustavo Carvalho Henrique Dallazuanna pelas dicas. Henrique, muito bacana ter enviado um exemplo enxuto e didático, facilitou muito meu entendimento e na resolução dos problemas do pesquisador. Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

R em ambiente servidor
by Vinicius Brito Rocha 15 Jul '11

15 Jul '11

Senhores, preciso utilizar o R em um ambiente remoto, um servidor, com mais capacidade de processamento que minha maquina local. Alguem como posso fazer isso? Seria instalar o R e usar algum tipo aplicativo tipo tinn-r ou R-studio? Outra pergunta, qual q real diferença entre o R-gui o Rterm? Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

4 8

[Dúvida] Rodar código com hora marcada
by Pedro Rafael 15 Jul '11

15 Jul '11

Olá pessoal como vão? Meu problema é o seguinte: Como rodar um script do R com hora marcada. Tenho um computador que quero deixar ligado e a medida que a base de dados for alimentada gostaria que o R se preocupasse em rodar novamente o código. Para mim seria suficiente estabelecer horas para que o R rodasse os comandos. Existe algum objeto já pronto que possa fazer isso de forma eficiente? Saudações, Pedro Rafael Diniz Marinho

3 2

problema rJava - xlsx
by Vinicius Brito Rocha 14 Jul '11

14 Jul '11

Srs. estou tentando instalar o pacote xlxs, mas recebo uma mensagem de erro no rJava. O que pode estar acontecendo? erro: Error : não é possível alocar vetor de tamanho 914.6 Mb Erro: carregamento do pacote/espaço de nomes falhou para 'rJava' o pacote tem apenas 650 kb -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

Re: [R-br] AJUDA - amova (package ade4)
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Olá Ricardo! Forneça-nos um CMR com uma pequena amostra dos seus dados para que possamos ajudá-lo. O erro me parece que está na forma de como vc está declarando no modelo. Pelo erro, é provável que não exista variação de 'samples' dentro de 'regions'. Veja melhor o seu modelo, e verifique se realmente há a necessidade de se colocar esta fonte de variação (se ela existe mesmo)! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Rubem, A função corr não é uma função específica para ser utilizada com o pacote boot. Pode ser utilizada em qualquer outra situação. Vide ?corr. A diferença da corr para cor está apenas num vetor de pesos que pode ser utilizado na função corr. Caso este vetor não seja declarado, o resultado da corr = cor. Quanto as processos de reamostragem, e lendo com mais atenção o artigo do DiCicci & Efron 1996, você tem toda a razão. O modo como estava sendo feito, estava equivocado, pois a distribuição bivariada ocorre para o par XiYi. Valeu Rubem, são essas chamadas de atenção que enriquecem nossos conhecimentos. Muito obrigado!! Abraço! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

AJUDA - amova (package ade4)
by Ricardo Fernandes Campos Junior 14 Jul '11

14 Jul '11

Olá, pessoal. Estou tentando usar a função "amova" do pacote ade4, mas não consigo. Comparei a estrutura dos meus objetos com a dos objetos do HELP da função, mas me parece tudo ok. A função roda tudo bem, mas recebo o seguinte resultado: amova(samples,sqrt(distance),structure) $call amova(samples = samp3, distances = sqrt(distance), structures = structurss) $results Df Sum Sq Mean Sq Between regions 2 0.1371102 0.06855508 Between samples Within regions 0 0.0000000 NaN Within samples 41 1.7349423 0.04231567 Total 43 1.8720525 0.04353610 $componentsofcovariance Sigma % Variations Between regions NaN NaN Variations Between samples Within regions NaN NaN Variations Within samples 0.04231567 NaN Total variations NaN NaN $statphi Phi Phi-samples-total NaN Phi-samples-regions NA Phi-regions-total NaN Abaixo o resultado de 'str' dos meus objetos. str(samples) 'data.frame': 19 obs. of 3 variables: $ AF : int 0 7 1 4 0 1 2 4 1 1 ... $ Caatinga: int 1 0 0 0 0 0 0 0 0 0 ... $ Cerrado : int 0 0 0 0 1 0 0 0 0 0 ... str(distance) Class 'dist' atomic [1:171] 0.0411 0.1007 0.0919 0.1007 0.0919 ... ..- attr(*, "Labels")= chr [1:19] "4" "12" "17" "19" ... ..- attr(*, "Size")= int 19 ..- attr(*, "call")= language as.dist.default(m = mat) ..- attr(*, "Diag")= logi FALSE ..- attr(*, "Upper")= logi FALSE str(structure) 'data.frame': 3 obs. of 1 variable: $ regions: Factor w/ 3 levels "AF","Caatinga",..: 1 2 3 -- Abraços, Ricardo Fernandes Campos Junior

1 0

conexao teradata
by Vinicius Brito Rocha 14 Jul '11

14 Jul '11

Srs. alguem sabe como funciona a conexao R ao banco de dados teradata? pelo que vi nao é via odbc. abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 3

Re: [R-br] Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Verdade!! Valeu. Agora uma coisa que estou quebrando a cabeça é a seguinte. Porque que o boot gerado pela função boot está extremamente diferente do boot gerado na 'unha'?. Vejamos um CMR. x <- c(0.75,2.55,1.40,0.80,2.60) y <- c(0.8249,0.797,0.7762,0.7797,0.7483) da <- data.frame(x,y) corr(da) mean(boot(da,corr,10000)$t) #Esta função eu fiz para evitar NAs, já que algumas amostragens são todas iguais, gerando um desvio padrão zero. sa <- function(x){ y <- sample(x,replace=T) if(y[1:1] == y[2:2] && y[1:1] == y[3:3] && y[1:1] == y[4:4] && y[1:1] == y[5:5]){ z <- sample(x,replace=T)}else{ y} } boot_da <- numeric(10000) n <- 10000) for(i in 1:n){ boot_da[[i]] <- corr(matrix(c(sa(da[,1]),sa(da[,2])),ncol=2)) } boot_da mean(boot_da) Se alguém puder ajudar de novo, agradeço. Abraços! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Ajuda: tratamento de erros no R
by Jose Claudio Faria 14 Jul '11

14 Jul '11

Pessoal, Quem tiver experiência com tratamento de erros no R poderia fazer o favor de me colocar na veia boa: funções, pacotes, material didático, links, etc? Um pesquisador (aqui da UESC) está precisando proteger um trecho de seu script contra erros. Ele já detectou (e mapeou) os possíveis erros, mas não sabe como tratar. Ou seja, registrar o erro e continuar (quando possível) com suas simulações (bem demoradas relecionadas ao estudo de estabilidade de galáxias - Astrofísica). Minha experiência com tratamento de erros no R é bem simples (detectar erro, dar mensagem e parar) e insuficiente para sua necessidade. Nesse caso é necessário uma abordagem mais consistente e efetiva: 1. Proteger o bloco onde ocorrem os erros; 2. Uma vez que ocorrer verificar qual foi o erro; 3. Registrar sua ocorrência (por exemplo num arquivo_erros.txt) e prosseguir (quando/se ppossível) Toda boa informação é bem vinda! Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 2

Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Bom dia senhores!! Alguém que tem o livro do Davison, A.C. and Hinkley, D.V. (1997) poderia me disponibilizar o algorítimo para cálculo do intervalo de confianção BCa. A função boot.ci da library boot calcula estes intervalos com o parâmetro viesado, e eu não quero. Não consegui modificar a função também, pois ela está amarrada ao output da função boot. Desde já obrigado! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

4 3

Erro ao ler banco no mac os
by Edson Lira 14 Jul '11

14 Jul '11

Bom dia, lista. Alguém sabe como corrigir esse erro? Ou o que é esse erro? > ag<-read.xls("/Users/edsonlira/IBGE/Censo_2010/Agregado/Agre2010_sinopse_AM.xls",sheet=1) Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 2 > ag<-read.csv2("/Users/edsonlira/IBGE/Censo_2010/Agregado/Agre2010_sinopse_AM.csv") Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 2 Edson Lira Estatístico Manaus-Amazonas

1 1

Dúvida para montar um agrupamento
by Clayton Santos Delfino 14 Jul '11

14 Jul '11

Olá pessoal, boa noite. Fico feliz por ter encontrado uma lista de R em português pois estou começando a aprender e tenho várias dúvidas e espero contar com a preciosa ajuda de todos. Eu já consegui realizar algumas análises com R mas agora empaquei com o seguinte: tenho um data.frame com a seguinte estrutura: identificador municipio regiao principal 1 Sao Paulo 1 Sim 1 Sorocaba 1 Não 1 Curitiba 2 Não 1 Sao Paulo 1 Sim 1 Porto Alegre 3 Não 1 Cuiaba 4 Não 2 Sao Paulo 1 Não 2 Sorocaba 1 Não 2 Curitiba 2 Sim 2 Curitiba 2 Sim 3 Porto Alegre 3 Sim 3 Porto Alegre 3 Sim 3 Porto Alegre 3 Sim 3 Sorocaba 1 Não Eu o que eu estou tentando fazer e ainda não encontrei um meio seria gerar uma tabela com uma linha para cada identificador onde somente o municipio com maior frequencia, considerando-se "principal = Sim" fosse retornado, exemplo: identificador Municipio Principal (frequencia) 1 São Paulo 2 2 Curitiba 2 3 Porto Alegre 3 Se alguem puder me ajudar agredeço desde já. Um abraço a todos. Clayton

7 10

Problema com o pacote amélia no debian 5.0 amd64:
by GESER 13 Jul '11

13 Jul '11

Gostaria de fazer uma imputação no linux através do pacote Amélia II versão 1.2-18 e não estou conseguindo. O estranho é que o script roda perfeitamente no Windows. O script está abaixo: #### Início do script #### #Carrego o arquivo para ser imputado# umi99 <- read.table("data/umid99NA.txt", header=TRUE)[,-1] #Carrego o pacote Amélia# require(Amelia) nrep <- 3 system.time(a.out <- amelia(umi99, m=nrep)) #Há mais linhas de comando adiante, mas o problema acontece na linha acima# #### Fim do script #### O erro acontece sempre ao final de uma imputação e diz: Erro em dimnames(impdata$covMatrices)[[1]] <- prepped$theta.names : comprimento de 'dimnames' [1] deve ser o mesmo de 'dims' [3] Se necessário, posso enviar uma imagem do prompt anunciando o erro e o arquivo que carrego para os dados serem imputados. Obrigado, Fabio

4 5

Re: [R-br] Ajuda (possivelmente simples) para reestruturar um data.frame
by Jose Claudio Faria 13 Jul '11

13 Jul '11

Walmes e Benilton, Não deu tempo nem de eu ir tomar um café!!! ;) Muito elegantes as soluções: muito obrigado! Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 1

[Dúvida] Como plotar valores
by Pedro Rafael 13 Jul '11

13 Jul '11

Pessoal, como plotar o valor dos pontos em um gráfico de linha no R? Exemplo: x = 1:10 plot(x,type = "o",pch = 12) Isso vai me gerar um gráfico com 10 pequenos desenhos. Como colocar o seu valor próximo desses desenhos? Pedro

3 2

ajuda xlsReadWrite
by Vinicius Brito Rocha 13 Jul '11

13 Jul '11

Srs. alguem pode me mandar por email o pacote xlsReadWrite 1.54 do http://www.swissr.org/download estou no trabalho e nao consigo baixar. e ainda estou com o problema do rJava e não consigo usar o xlsx abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

1 0

Ajuda (possivelmente simples) para reestruturar um data.frame
by Jose Claudio Faria 13 Jul '11

13 Jul '11

Pessoal, Ajuda ai que estou travado! #Tenho isso (desconsiderar a numeração das linhas) tra r y 1 A 1 2.70 2 A 2 3.02 3 A 3 3.25 4 A 4 3.21 5 A 5 3.04 6 A 6 3.07 7 B 1 1.16 8 B 2 1.16 9 B 3 0.91 10 B 4 0.93 11 B 5 0.82 12 B 6 0.71 13 C 1 3.01 14 C 2 2.91 15 C 3 3.14 16 C 4 2.96 17 C 5 3.01 18 C 6 2.93 19 D 1 4.88 20 D 2 4.68 21 D 3 5.20 22 D 4 4.94 23 D 5 5.34 24 D 6 5.07 25 E 1 4.89 26 E 2 5.26 27 E 3 5.07 28 E 4 5.14 29 E 5 5.47 30 E 6 4.98 # Para economizar espaço em um relatório (que uso muito) preciso disso (idem): tra r y tra r y tra r y tra r y tra r y 1 A 1 2.70 B 1 1.16 C 1 3.01 D 1 4.88 E 1 4.89 2 A 2 3.02 B 2 1.16 C 2 2.91 D 2 4.68 E 2 5.26 3 A 3 3.25 B 3 0.91 C 3 3.14 D 3 5.20 E 3 5.07 4 A 4 3.21 B 4 0.93 C 4 2.96 D 4 4.94 E 4 5.14 5 A 5 3.04 B 5 0.82 C 5 3.01 D 5 5.34 E 5 5.47 6 A 6 3.07 B 6 0.71 C 6 2.93 D 6 5.07 E 6 4.98 # Fiz na "mão grande" usando "cbind" assim dados.tb <- cbind(dados.q1[1:6,], dados.q1[7:12,], dados.q1[13:18,], dados.q1[19:24,], dados.q1[25:30,]) Contudo, creio que deve ter uma solução mais racional e genérica para quando mudar o número de tratamentos e repetições! Tentei split, pacotes doBy e reshape, mas não fui bem sucedido. Para reprodução de possível solução pode usar 'df' abaixo: df <- data.frame(gl(n=5, k=6, lab=LETTERS[1:5]), rep(1:6, 5), 1:30) names(df) <- c('tra', 'r', 'y') Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 3

Interpolação
by Nelson Betzek 13 Jul '11

13 Jul '11

É possível fazer interpolação pelo inverso da distancia elevado a uma potencia, ou só por Krigagem?

3 2

Modelo linear generalizado misto
by Marcelo Cardoso mello 13 Jul '11

13 Jul '11

Pessoal da lista ... tenho um modelo logístico que segue glm(Doença~Raça+Fazenda+Medicamento) ou glmer(Doença~Raça+Medicamento+(1|Fazenda)) O que vcs acham de colocar fazenda (animais em uma amostra de fazndas) como um fator aleatório no modelo? OU então colocar os próprios indivíduos como um fator aleatório? É uma opção melhor? Posso colocar o intercepto como aleatório também? E a odds ratio como fica? ... posso estimar a odds dos fatores fixos do modelo e dai qual função estima odds com teste de significância a partir de um modelo glmer? Uso logistic.display da epicalc para glm e para glmer?. Na verdade não sei se a glmer é a melhor função para o meu caso. Alguem conhece algum material que possa me ajudar com exemplos no R?Help.Marcelo.

5 5

.net+R
by Vinicius Brito Rocha 12 Jul '11

12 Jul '11

Srs. existe uma maneira de eviar resultados de R diretamente para .net? att. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 2

Dúvida
by Edson Lira 12 Jul '11

12 Jul '11

Boa tarde R-br lista, estou com um pequeno problema. Tenho uma lista de exames que associo ao laboratório respectivo, vejam um exemplo abaixo: Descrição lab Quantidade 1 17 ALFA HIDROXI PROGESTERONA <NA> 1 2 ACIDO URICO Bioquímica 1 3 ACIDO URICO Bioquímica 1 4 ACIDO URICO Bioquímica 1 5 ACIDO URICO Bioquímica 1 6 ACIDO URICO Bioquímica 4 7 ACIDO URICO Bioquímica 9 8 ACIDO URICO Bioquímica 1 9 ACIDO URICO Bioquímica 1480 10 ACIDO URICO Bioquímica 37 11 ACIDO VALPROICO <NA> 1 12 ALANINA AMINOTRANSFERASE (ALT/GPT) IFCC <NA> 1 13 ALBUMINA Bioquímica 205 14 ALBUMINA Bioquímica 17 15 ALFA-FETO-PROTEINA HORMÔNIOS 3 16 ALFA-FETO-PROTEINA HORMÔNIOS 2 17 AMILASE Bioquímica 77 18 AMILASE Bioquímica 4 19 ANDROSTENEDIONA <NA> 1 20 ANDROSTENEDIONA <NA> 1 Para gerar os dados de contagem dos exames por laboratório. Estou implementando uma rotina para otimizar o trabalho de consolidação dos dados. O problema é que alguns exames que já estão associados ao respectivo laboratório quando não são realizados no período, eles não têm registros, mas no script estão relacionados, vejam a rotina abaixo: library(memisc) acj$lab<-recode(acj$Lab, "Bioquímica" <-c("ACIDO URICO","ALBUMINA","BILIRRUBINA DIRETA", "BILIRRUBINA INDIRETA","BILIRRUBINA TOTAL E FRAÇÕES","BILIRRUBINAS TOTAL", "COLESTEROL", "COLESTEROL HDL", "COLESTEROL LDL", "COLESTEROL VLDL", "CREATININA","DHL - DESIDROGENASE LÁCTICA", "FERRITINA", "FERRO SÉRICO", "FERRO SÉRICO LABTEST", "FOSFATASE ALCALINA", "FOSFATASE ÁCIDA FRAÇÃO PROSTÁTICA", "FOSFATASE ÁCIDA PROSTÁTICA (RIE)", "FOSFORO", "GAMA GT- GAMA-GLUTAMIL-TRANSFERASE", "GLICOSE", "GLICOSE COM FLUORETO", "GLICOSE POS-PRANDIAL", "GLICOSE-6-FOSFATO DEHIDROGENASE (GGFD)", "GLOBULINA", "HEMOGLOBINA GLICOSILADA", "LIPIDIOS TOTAIS", "LIPIDOGRAMA", "MAGNESIO", "POTÁSSIO", "PROTEINAS TOTAIS", "SATURAÇÃO DE TRANSFERRINA","SÓDIO","TRANSAMINASE OXALACÉTICA - AST", "TRANSAMINASE PIRÚVICA - ALT", "TRANSAMINASES", "TRANSAMINASES", "TRANSFERRINA", "ALFA-1 GLICOPROTEINA ÁCIDA", "ALT-TRANSAMINASE PIRÚVICA", "AMILASE", "AST-TRANSSAMINASE OXALACÉTICA", "CÁLCIO","CAPACIDADE DE OXIDAÇÃO DO FERRO","COLESTEROL TOTAL","COLESTEROL TOTAL E FRAÇÕES", "ESTUDO DO FERRO", "FOSFATASE ÁCIDA TOTAL", "TESTE DE ROG", "TRIGLICERÍDEOS", "URÉIA", "LÍTIO"), "Citoquímica" <-c("FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM")) Quando rodo esta rotina o R mostra a mensagem: Erro em factor(y, levels = 1:max(y), labels = newcodes) : invalid labels; length 13 should be 1 or 12 Além disso: Mensagens de aviso perdidas: In recode(acj$Lab, "Bioquímica" <- c("ACIDO URICO", "ALBUMINA", : recoding "Citoquímica" <- c("FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM") has no consequences O problema é que não existe nenhum exame associado ao laboratório de Citoquímica (FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM). ou seja não foi realizado nenhum exame, portanto, não aparecem no banco, como resolver este problema? []'s Edson Lira Estatístico Manaus-Amazonas

1 0

Equação do plano
by Helber Freitas 12 Jul '11

12 Jul '11

Olá, Pessoal. Fico feliz em poder contar com o grupo. Sou novo por aqui e aproveito para agradecer e parabenizar a iniciativa. Preciso de uma função que forneça a equação do plano dados 3 pontos distintos (x1,y1,z1), (x2,y2,z2) e (x3,y3,z3). Grato, Helber

3 2

RJAVA JRI (problema ao plotar gráficos)
by Emmanuel Arnhold 12 Jul '11

12 Jul '11

Estou tentando fazer o R "conversar" com o JAVA. Estou usando RJAVA e JRI. Tudo certo quando rodo algumas análises. O problema é que os gráficos não estão aparecendo na janela gráfica. Alguma dica?? O que pode estar errado?? Agradeço desde já. Emmanuel

2 2

Re: [R-br] Derivadas Parciais e Ponto Crítico!
by Thiago De paula protásio 11 Jul '11

11 Jul '11

Obrigado Walmes pela dica do pacote... Resolveu o meu problema....a função rms retorna os pontos críticos e os autovalores da equação característica: IΒ −λI = 0 E a análise dos sinais dos autovalores são fundamentais.... Novamente agradeço-lhe!Abraço! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Derivadas Parciais e Ponto Crítico!
by Thiago De paula protásio 11 Jul '11

11 Jul '11

Pessoal; Recorro novamente a ajuda de vcs...Estou trabalhando com a seguinte superfície de resposta.... z = 31.13884 -2.23631*x +0.04441*x^2 + 1.25980*y -0.01134*y^2 Contudo, gostaria de determinar no R os pontos críticos e verificar se a natureza do ponto: sela, máximo ou mínimo. A ideia consiste nas derivadas parciais....Dessa forma, gostaria de saber se existe alguma função no R que me auxilie nesse problema. Muito obrigado!!! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

ajuste com nls
by Gilson Sanchez 11 Jul '11

11 Jul '11

Pessoal mais uma vez procuro ajuda de vc's, consegui ajustar um modelo com nls, mas agora o que quero saber como eu faço para estimar o coeficiente (R2) e ver se a equação ajustou corretamente. obrigado -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

3 2

[Notícias-Dúvidas] Novidades de interfaces gráficas no R.
by Pedro Rafael 11 Jul '11

11 Jul '11

Pesquisando mais sobre algumas novidades do R vi que existem vários projetos de construção de interfaces gráficas para o R fora o RStudio, R Comander, Tinn-R e RKward que sãos os mais conhecidos. Existe o Red-R que tenta facilitar a vida do usuário através de uma programação visual, uma espécie de SAS GUIDE que é uma funcionalidade do SAS bastante elogiada pelos usuários SAS. Observe o vídeo do Red-R em: http://www.youtube.com/watch?v=A4WwhOnuC10&feature=player_embedded. Existe um outro projeto também livre chamado R Analytic Flow, veja: http://www.ef-prime.com/products/ranalyticflow_en/. Amigos vocês conhecem mais alguns projetos como esse? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico (Secretaria de Saúde de Estado da Paraíba).

2 1

Re: [R-br] Identidade de modelos!!!
by Thiago De paula protásio 10 Jul '11

10 Jul '11

Walmes; No meu caso aqui ajustei modelos não lineares para 11 tratamentos e gostaria de saber se posso utilizar um único modelo para todos eles.... Muito obrigado pela informação.... Abraço! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Identidade de modelos!!!
by Thiago De paula protásio 10 Jul '11

10 Jul '11

Prezados membros da R-br; Venho novamente pedir a ajuda de vcs...Gostaria de saber como posso proceder ao teste de identidade de modelos no R?Ajustei algumas equações aqui para 6 tratamentos e visualmente percebi certa similaridade entre elas....dessa forma, gostaria de proceder a algum teste para testar a hipótese de que apenas uma equação seria necessária para todos os tratamentos.... Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Extrair células de uma matrix
by Diego Pujoni 10 Jul '11

10 Jul '11

Olá Pessoal, Tenho uma matriz onde cada coluna contém valores de probabilidade variando de 0 a 1. Como posso excluir os valores 0's e 1's de cada coluna, ficando apenas com as probabilidades intermediárias? Toda coluna possui apenas um valor 0 e apenas um valor 1, mas eles estão em linhas diferentes em cada coluna. Fiz um exemplo abaixo > A = matrix (c(1,.5,.8,.9,0.,.8,.5,0,1,.4,.8,0,.7,.3,1),5,3) > A [,1] [,2] [,3] [1,] 1.0 0.8 0.8 [2,] 0.5 0.5 0.0 [3,] 0.8 0.0 0.7 [4,] 0.9 1.0 0.3 [5,] 0.0 0.4 1.0 Quero obter essa matriz > A [,1] [,2] [,3] [1,] 0.5 0.8 0.8 [2,] 0.8 0.5 0.7 [3,] 0.9 0.4 0.3 Obrigado Diego PJ

2 1

regressão linear pelo origem
by Gilson Sanchez 08 Jul '11

08 Jul '11

Pessoal mais uma vez procuroa juda de vc's, bom a verdade é pra um colega ele é novo com o R, se alguem puder ajudar. Estou precisando da ajuda de você, sou novo no grupo e estou tentando fazer um grafico de regressao linear simples ajustando o intercepto da reta deve passar pela origem dos eixos (0,0) mas não estou conseguindo. O CMR que estou trabalhando é o seguinte. # REGRESSÃO EI30 VS DS require(XLConnect) reg<-loadWorkbook("Estadistica Final.xls", create= FALSE) reg1<-readWorksheet(reg, sheet="Plan1") reg1 attach(reg1) cor.test(EI30, DS, alternative = c("two.sided"), method = c("pearson"), conf.level = 0.95) reglin1<-lm(DS~EI30, data=reg1) anova(reglin1) reglin1 summary(reglin1) # GRAFICOS par(mfrow = c(3,1), mar=c(4.5,5,2,2)) plot(EI30,DS, type="p", pch=15, lwd=1.5, ylim=c(0,100), xlim=c(0,50), cex.main=1.2, cex.lab=0.9, cex.axis=0.9, xlab="", ylab="",cex.lab=2, main="Descoberto", font.main=2) abline(reglin1, lwd=2.0, col="black", lty=1) ## pode-se usar lty=9 legend("bottomright", bty="n", lwd=c(1.5), pch=c(15), col=c("black"), legend=c(expression(y=='1,516x'~+~'9,241'~','~R^2=='0,58'^'*')), cex=1) plot(EI30, LS, pch=16, col="black", ylim=c(0,3), xlim=c(0,50), cex.main=1.2, cex.lab=0.9, cex.axis=0.9, xlab="", ylab=expression(bold(Perda~do~solo~(g~L^{-1}))), font.lab=2, main="Liteira", font.main=2) abline(reglin2, lwd=2.0, col="black", lty=2) legend("bottomright", bty="n", lwd=c(1.5), pch=c(16), col=c("black"), legend=c(expression(y=='0,018x'~+~'0,323'~','~R^2=='0,16'^'*')), cex=1) Aguardo suas consideraçoes. Att. Omar Cubas Encinas Instituto Nacional de Pesquisas da Amazônia / INPA Aluno de Mestrado em Agricultura do Tropico Úmido / ATU Manaus - AM - Brasil (092) 9102-2542 -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

3 3

Multivariada - Clustering
by Emerson Cotta Bodevan 08 Jul '11

08 Jul '11

Prezados, qual pacote do R é mais adequado para análise multivariada? Especificamente, clustering (gráficos de dendograma,etc). Abs, -- *Emerson*

4 4

Kappa
by Edson Lira 08 Jul '11

08 Jul '11

Boa tarde R-Br lista!!! Para as tabelas abaixo, como faço para calcular o índice de concordância de Kappa e sua signficância. TABELA 1 OURO negative positive OB1 negative 449 3 positive 6 7 TABAELA2 OURO negative positive OB2 positive 1 8 TABELA 3 OURO negative positive OB3 negative 0 1 positive 1 7 Estou usando o pacote Psych, mas não vi nenhuma medida de signficância do Kappa. Alguem tem alguma sugestão? Edson Lira Estatístico Manaus-Amazonas

4 3

Re: [R-br] funcao tableStack do epicalc
by Pedro Emmanuel Alvarenga Americano do Brasil 07 Jul '11

07 Jul '11

Paula, Algumas dicas pra poder rodar essa funcao. Em primeiro lugar ela foi feita pra rodar tudo dE uma vez,nao faz sentido chamar pra grupos de variageis diferentes. Nessa funcao as coisas dao pau por conta do formato de uma ou mais variaveis. Nao use character, somente factor. Confira se as factor possui algum nivel que nao tem dado. Faca a funcao progressivamente inserindo as variaveis uma a uma. De uma olhada nos exemplos pra ver te da uma - paula - luz. :-) Pedro Brasil via android (:)= Em 07/07/2011 17:49, "Paula Mendes Luz" <luzpaulamendes(a)gmail.com>escreveu: oi benilton com certeza. por acaso não estou com o código em mãos agora e por isso tentei reproduzir em palavras... mas entendo q é mais complicado assim. então depois envio outro email. p. 2011/7/7 Benilton Carvalho <beniltoncarvalho(a)gmail.com>: > Paula, > > se vc puder seguir as recomendacoes do guia de postagem (vide rodape') > e fornecer um ...

1 0

funcao tableStack do epicalc
by Paula Mendes Luz 07 Jul '11

07 Jul '11

oi pessoal estou tentando usar a funcao tableStack da library epicalc mas tá complicado... os erros que recebo nao sao mto informativos... alguem tem alguma dica? tem algum tipo de variavel que a funcao nao aceita? acho que o problema tá nisso (apesar de q os erros sao diversos)... lendo o help acho que entendi que a funcao só aceita um tipo de var de cada vez... entao se tenho vars numericas e character, devo separar as vars... por ex, tenho idade (numeric), gender (character: 1/2), education (character: 1/2/3/4/9), altura (numeric), centro (character: 1/2) faco o tableStack de idade e altura usando o by = centro depois o tableStack de gender e education, by = centro, com vars.to.factor de gender e education (preciso transformar o centro tb?) ex de erros recebidos incluem: object 'table0' not found e object 'test.method' not found qq ajuda é bem vinda! valeu! paula.

2 2

Re: [R-br] R utilizando apenas 1 processador para rodar
by Junior Beleti 06 Jul '11

06 Jul '11

Olá Paulo Nogueira, como já foi mencionado, o pacote multicore possui funções de "alto nivel" para utilizar o paralelismo. A sua aplicação necessita de alto desempenho? É que estou trabalhando com processamento paralelo e talvez possa ajudar. Abraço -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

7 13

Grau de concordância
by Edson Lira 06 Jul '11

06 Jul '11

Bom dia R-Br lista, tenhos as observações abaixo: g_ob1 g_ob2 wbc 1 15 3 0.00 2 35 40 0.09 3 200 200 0.10 4 310 300 0.30 5 60 60 0.09 6 860 800 1.60 7 506 320 0.60 8 125 138 0.24 9 4 0 0.01 Desejo checar a concocrdância entre ( g_ob1, g_ob2); (g_ob1,wbc), (g_ob2,wbc) Nos pacotes DiagnosisMed e psych só trabalham com positivo e negativo,certo, errado. É isso ou estou enganado? Alguém poderia me indicar uma solução para isso? []'s. Edson Lira Estatístico Manaus-Amazonas

5 6

Plot matriz em 3D
by Tiago Vieira 06 Jul '11

06 Jul '11

Boa noite a todos, eu tenho uma matriz 10 X 50 e quero plotar um gráfico em que x é o termos da linha, y os termos das colunas e z o conteúdo da matriz. Vocês podem me dar uma dica,? Desde já, agradeço! Tiago Vieira.

2 1

Enc: Diagnosis rotina
by Edson Lira 06 Jul '11

06 Jul '11

Boa tarde R-Br lista, estou encaminhando essa mensagem que havia mando para o endereço antigo da lista. Acho que foi a pressa. []'s Edson Lira Estatístico Manaus-Amazonas - ---- Mensagem encaminhada ----- De: Edson Lira <edinhoestat(a)yahoo.com.br> Para: Grupo R <r_stat(a)yahoogrupos.com.br> Enviadas: Quarta-feira, 6 de Julho de 2011 10:09 Assunto: Diagnosis rotina Pessoal tenho essa tabela, ob2 negative positive ob1 negative 455 0 positive 2 8 estou tentando aplicar a rotina abaixo: require(DiagnosisMed) #diagnosis(a,b,c,d, CL = 0.95, print = TRUE, plot = FALSE) a A number representing True Positives from a 2x2 table. Also, TP could either a 2x2 table (see below) or a column in a dataset representing the reference standard. b A number representing False Negatives from a 2x2 table. If TP is a column in a dataset FN should also be a columns in a dataset, however representing the index test. c A number representing False Positives from a 2x2 table d A number representing True Negatives from a 2x2 table diagnosis(455,0,1,8, CL = 0.95, print = TRUE, plot = FALSE) que me está me mostrando o erro abaixo: Erro em uniroot(function(or) { : f() values at end points not of opposite sign Se substituir o 0 pelo 1, roda, como faço para corrigir isso? []'s Edson Lira Estatístico Manaus-Amazonas

1 0

Re: [R-br] R-studio - hot key
by Ivan Bezerra Allaman 06 Jul '11

06 Jul '11

Vinicius, Talvez lhe seja útil, http://rstudio.org/docs/using/keyboard_shortcuts Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Fazendo busca por mensagens na lista R-Br
by Paulo Nogueira 06 Jul '11

06 Jul '11

Aproveitando o gancho, como entrei faz pouco tempo na R-Br não recebi as primeiras mensagens no meu email. Existe um local na web para procurar as mensagens antigas? obrigado Paulo 2011/7/6 Gustavo Carvalho <gustavo.bio(a)gmail.com> > (...) Tem um post do Benilton em que ele mostra as situações > mais comuns e compara o tempo de execução delas. Dê uma procurada. > >

3 2

R-studio - hot key
by Vinicius Brito Rocha 06 Jul '11

06 Jul '11

Tem como criar uma tecla de atalho para rodar o código selecionado no R-Studio? sempre usei o F10 (por causa do S-plus) abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 1

R utilizando apenas 1 processador para rodar
by Paulo Nogueira 06 Jul '11

06 Jul '11

Bom dia R-Br. Eu participava da antiga R-Stat e demorei um pouco para migrar para o novo ambiente mas estou por aqui faz alguns dias. Tenho uma questão com relação à priorização que o R faz relacionada aos processadores disponíveis no computador. Tenho um notebook com 4 núcleos de processamento mas percebo que quando rodo algo pesado no R apenas um deles é utilizado. O escolhido varia de vez em vez. Alguém sabe como fazer ele utilizar melhor o processamento disponível ou que caminho devo seguir para descobrir? Acho que o ideal seria ter os processadores trabalhando juntos. Utilizo o R através do RStudio (http://rstudio.org/) rodando em Ubuntu 10.10. Obrigado. Paulo Nogueira Starzynski

6 8

FEDORA ou UBUNTU?
by Leonard Assis 06 Jul '11

06 Jul '11

Qual dos dois é melhor pra rodar o R em um i5 notebook com 6GB de RAM e hd de 640 gb (Samsung sd1) Leonard

10 12

desenhos de blocos incompletos com restrição
by Bernardo Rangel Tura 06 Jul '11

06 Jul '11

Bom dia Estou precisando de uma ajuda para gerar uma tabela de blocos incompletos no R. Meu problema é o seguinte tenho 9 tratamentos: A1,A2,A3 B1,B2,B3 C1,C2,C3 Preciso gerar blocos com 6 tratamentos cada mas com três limitações 1- Cada bloco tem que ter pelo menos um representante de cada letra (ou seja não pode haver o bloco: A1,A2,A3,B1,B2,B3) 2- Cada bloco tem que ter pelo menos um representante de cada numero (ou seja não pode haver o bloco A1,A2,B2,B1,C2,C1) 3- Não pode haver tratamentos repetidos (ou seja não pode haver o bloco: A1,A2,A1,B1,C2,B3) O total de observações deve ser em torno de 300. Se não houvesse as restrições seria fácil require(crossdres) find.BIB(9,50,6) Porém com esta restrições (1,2) não estou conseguindo alguém pode me ajudar? -- []s Tura

2 3