- R-br - DINF UFPR

[Dúvida] Casas decimais
by Pedro Rafael 29 Jul '11

29 Jul '11

Pessoal qual é o objeto que eu posso acrescentar casas decimais à numeros inteiros, exemplo: se eu tenho 0 colocar 0.00, se eu tenho 2 colocar 2.00, etc. Obrigado.

3 3

Re: [R-br] Aplicações que necessitam de alto desempenho
by Junior Beleti 29 Jul '11

29 Jul '11

Olá Benilton, vou dar uma olhada em decomposicao espectral de matrizes. Em suas aplicações você já utiliza programação paralela então? Que pacotes você utiliza? Marcos.. Não é bem o que preciso, mas obrigado. E se tiver alguma aplicação que leve muito tempo pra processar, estou a disposição. Obrigado, abraço.. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

Letalidade para o Estado
by Pedro Rafael 29 Jul '11

29 Jul '11

Divirta-se. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

(sem assunto)
by Juraci Mendes Moreira 28 Jul '11

28 Jul '11

ola senhores, gostaria de me desligar da lista obridado.

2 1

Re: [R-br] Comparar informações de colunas e linhas
by Roney Fraga Souza 28 Jul '11

28 Jul '11

beniltoncarvalho, Muito obrigado pela ajuda. att Roney

1 0

Aplicações que necessitam de alto desempenho
by Junior Beleti 28 Jul '11

28 Jul '11

Boa tarde a todos. Como já mensionei em tópicos anteriores, tenho realizado pesquisa em programação paralela para o ambiente R. Já vi alguns tópicos sobre problemas de desempenho mas que foram solucionados com uma melhor programação e utilização de pacotes. Estou interessado em aplicações que realizem grande processamento e que talvez possam ter um melhor desempenho se executadas de forma paralela. Tenho realizado experimentos com uma aplicação, mas seria interessante ampliar tais experimentos para outras aplicações. Se alguém possuir tais aplicações e tiver interesse. Abraço. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

3 2

Cálculo da Previsão de Modelos ARIMA
by Paola Tame 28 Jul '11

28 Jul '11

Boa tarde Grupo, estou com problemas para interpretar as saídas da função arima. Por exemplo, usando a base de dados de exemplo do R "USAccDeaths". Ajustando o modelo: fit <- arima(USAccDeaths, order = c(0,1,1), seasonal = list(order=c(0,1,1))) Tenho a seguinte saída: Series: USAccDeaths ARIMA(0,1,1)(0,1,1)[12] Call: arima(x = USAccDeaths, order = c(0, 1, 1), seasonal = list(order = c(0, 1, 1))) Coefficients: ma1 sma1 -0.4303 -0.5528 s.e. 0.1228 0.1784 sigma^2 estimated as 99347: log likelihood = -425.44 AIC = 856.88 AICc = 857.32 BIC = 863.11 Utilizando a função *"predict"* para calcular os próximos 4 meses tenho o seguinte resultado: $pred Jan Feb Mar Apr 1979 8336.061 7531.829 8314.644 8616.869 Como vou precisar utilizar o modelo fora do R preciso entender como é feita essa previsão pois utilizando os coeficientes *ma1* e *sma1* estou obtendo resultados diferentes. Mesmo ajustando um modelo mais simples como um AR(1) para esses dados (apenas como exemplo, sem me preocupar com a qualidade do ajuste) os preditos que calculo não batem com os preditos retornados pela função "* predict*". Aparentemente a função de predição não é montada diretamente dos coeficientes estimados (p. ex. no caso de um modelo AR(1) X(t) = intercept + ar1*X(t-1)). Alguém poderia me ajudar? Obrigada, Paola Tame

2 1

Comparar informações de colunas e linhas
by Roney Fraga Souza 28 Jul '11

28 Jul '11

Caros, Reescrevendo o problema da mensagem anterior, procurei ser mais preciso nos exemplos. Segue as dúvidas: Objetivo I: Considerando XX, XX <- data.frame("A"=c(8,11,14,11,11,8)) XX A 1 8 2 11 3 14 4 11 5 11 6 8 quero criar uma nova variável em XX, que retorne 1 quando um determinado valor em XX$A já foi repetido ou 0 se é a primeira vez que ele aparece. Nesse caso teria que ficar XX A nova_variavel1 1 8 1 2 11 1 3 14 0 4 11 1 5 11 1 6 8 1 A pergunta é, como criar essa variável (nova_variavel1)? Objetivo II: Considerando YY, YY <- data.frame("A"=c(8,11,14,11,11,8), "B"=c(14,16,33,14,16,14)) YY A B 1 8 14 2 11 16 3 14 33 4 11 14 5 11 16 6 8 14 quero criar agora uma nova variável que apresente valor 1 se existem valores repetidos em YY$A e YY$B simultaneamente, exemplo, na linha 1 de YY os valores (8 e 14) se repetem na linha 6, logo a nova variável terá valor 1 em ambas a linhas citadas, caso contrário assumirá valor 0. O resultado seria YY A B nova_variavel2 1 8 14 1 2 11 16 1 3 14 33 0 4 11 14 0 5 11 16 1 6 8 14 1 Como criar a variável (nova_variavel2)? Objetivo III: Considerando YY, quero criar uma nova variável que conte quantas vezes os valores de determinadas linhas se repetem. Como os valores das linhas 2 e 5 de YY$A e YY$B são iguais, então nas linhas 2 e 5 da nova variável aparecera o número 2, caso os valores aparecessem em 4 linhas distintas o valor da nova variável seria 4 nas linhas citadas. O resultado do exemplo YY seria YY A B nova_variavel3 1 8 14 2 2 11 16 2 3 14 33 1 4 11 14 1 5 11 16 2 6 8 14 2 Por fim, a última questão é, como criar a variável (nova_variavel3)? Atenciosamente Roney

2 1

Separar strings
by Edson Lira 28 Jul '11

28 Jul '11

Pessoal mandei esta mensagem para o Benilton e compartilho com vocês, tenho uma dúvida que ainda não trabalhei, vejam se vocês podem me ajudar. Vejam os dados abaixo: melhoria [1] água [2] água e saúde [3] água e saúde [4] educação e saúde [5] tudo [6] tudo [7] atendimento na prefeitura [8] educação e saúde [9] educação e saúde [10] tudo [11] tudo [12] emprego [13] alimentação, educação e saúde [14] saúde, transporte e água [15] água e educação [16] atendimento na prefeitura e água [17] segurança e saúde [18] tudo [19] tudo [20] educação e saúde [21] tudo [22] tudo [23] educação e segurança [24] água e saúde [25] tudo [26] transporte, água, educação e saúde [27] saúde e água [28] tudo [29] saúde [30] tudo Veja que as respostas para a variável melhoria são do tipo: (água e saúde), (transporte,água), etc. Gostaria de separar estas expressões em 2 ou três variáveis, por exemplo: var1 var2 var3 var4 água saúde transporte água educação saúde Tem alguma sugestão? Edson Lira Estatístico Manaus-Amazonas

4 4

Comparar informações de colunas e linhas
by Roney Fraga Souza 28 Jul '11

28 Jul '11

Pessoal, Gostaria de comparar informações de determinada tabela considerando suas linhas, mas não sei como fazer. Utilizando a notação abaixo para exemplificar, quero verificar se existe alguma linha da coluna A que apresenta valor repetido. Algo como, se a1 == an então retorna 1 em c1, caso contrário 0, do mesmo modo, se a2 == an retorna 1 em c2, caso contrário 0... Outra possibilidade útil é contar a quantidade de repetições de um "a" qualquer e retornar a soma da quantidade de repetições em c. Mas o que realmente seria útil era verificar se existe alguma linha onde a1==an e b1==bn, quando isso ocorrer retornar 1 em d1, caso contrário 0, do mesmo modo, se a2==an e b2==bn, retornar 1 em d2, caso contrário 0. A B C D 1 a1 b1 c1 d1 2 a2 b2 c2 d2 3 a3 b3 c3 d3 4 a4 b4 c4 d4 n an bn cn dn onde n vai de 1 a +infinito Atenciosamente Roney

3 2

mensagem de warning
by Edson Lira 27 Jul '11

27 Jul '11

Caros amigos e Benilton, tenho um banco csv que salvei como csv do mac. Ao rodar notei alguns problemas: 1.A acentuação; Outro 2. a mensagem abaixo: Mensagens de aviso perdidas: 1: it is not known that wchar_t is Unicode on this platform Usei o encoding='latin1' e fileEncoding='latin1', mas não consegui corrigir o problema Alguma idéia? Edson Lira Estatístico Manaus-Amazonas

2 1

Re: [R-br] Experimento fatorial com tratamento adicional sem normalidade
by Ivan Bezerra Allaman 27 Jul '11

27 Jul '11

Adriano, Esse é o caminho das pedras, http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/ Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

eliminar linha completa de amostra rejeitada
by hacker gostoso 27 Jul '11

27 Jul '11

Olá membros. Minha duvida é a seguinte: Tenho 2 arquivos, sendo que o arquivo A é uma matrix, e o outro arquivo um vetor B. O vetor B seria as amostras rejeitadas (qualidade). Como faço para eliminar as amostras B de A? Ex. A com 10 colunas e 500 linhas (i linhas). B (j linhas) (i in 1:length(A)) (j in 1:length(B)) Eliminar A[i,2] (2a coluna) SE A[i,2] = B[j] e gerando um terceiro arquivo "C" Desde já muito obrigado Ricardo Pereira

2 1

Re: [R-br] Experimento fatorial com tratamento adicional sem normalidade
by Ivan Bezerra Allaman 27 Jul '11

27 Jul '11

Adriano, Se a variável não tem distribuição normal e você não quer utilizar nenhuma transformação, terá que recorrer a metodologia dos Modelos Lineares Generalizados (GLM) e encontrar a família de distribuição que se encaixe aos seus dados. Eu o aconselho tentar uma transformação primeiro ser você não estiver familiarizado com a metodologia GLM. Já tentou boxcox? Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

3 3

big memory
by Vinicius Brito Rocha 27 Jul '11

27 Jul '11

Estou começando a utilizar o big.memory. pelo q ue vi ele só lê arquivos tipo matrix (big.matrix), via read.big.matrix. alguém sabe como importar em um formato data.frame, usando o big.memory? 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.

2 1

Experimento fatorial com tratamento adicional sem normalidade
by Adriano Carvalho Costa 27 Jul '11

27 Jul '11

Boa noite, solicitar a ajuda de voces para realizar uma análise de um experimento fatorial com um tratamento adicional sem normalidade. Utilizando a base de dados abaixo (Fonte: Walmes) e considerando que ela não tem distribuição normal como posso fazer a análise sem utilizar transformação? fa <- read.table("http://www.leg.ufpr.br/~walmes/cursoR/fat-adi.txt", header=TRUE) attach(fa) fa names(fa) concentração=concentraÃ.Ã.o fa <- transform(fa, concentração=factor(concentração), bloc=factor(bloc)) m0 <- aov(media~bloc+trat, fa) anova(m0) m1 <- aov(media~bloc+origem*concentração, data=subset(fa, trat!="TEST")) summary(m1) Obrigado!

1 0

teste: desconsiderar
by Jose Claudio Faria 26 Jul '11

26 Jul '11

test!

1 1

[Dúvida] Compara dois vetores.
by Pedro Rafael 26 Jul '11

26 Jul '11

como procurar valores entre vetores de tamalho diferente? Exemplo: x = c("a","L","k","w","t") y = c("t","a","3") x[x==y] = "Casa" Mensagens de aviso perdidas:In x == y : comprimento do objeto maior não é múltiplo do comprimento do objeto menor Gostaria de ver dentro de x e y e o que for igual subistituir pela palabra casa. É um exemplo bem simples mas só para que eu possa entender como resolver problemas dessa natureza. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 1

Re: [R-br] Erro ao gerar Operadores matemáticos na compilação sweave!
by Ivan Bezerra Allaman 26 Jul '11

26 Jul '11

Pela nossa senhora dos detalhes em Benilton. rsrsrsrsrs. Deu certo. Obrigado!! Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

Re: [R-br] Erro ao gerar Operadores matemáticos na compilação sweave!
by Ivan Bezerra Allaman 26 Jul '11

26 Jul '11

Gera este erro Benilton, Writing to file teste.tex Processing code chunks with options ... 1 : term verbatim pdf Erro: chunk 1 Error in teste.Rnw:15:26: símbolo inesperado 14: plot(1:10,1:10) 15: text(3,8,expression('se' x ^ Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

Re: [R-br] Erro ao gerar Operadores matemáticos na compilação sweave!
by Benilton Carvalho 26 Jul '11

26 Jul '11

da forma q eu mandei, ele gera a saida como voce espera (ie., com o simbolo matematico). qual eh o erro q vc tem ao executar conforme minha sugestao? b 2011/7/26 Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br>: > Não sei te responder Fernando!! > Benilton, > não roda da forma como você colocou!! dá um erro ao compilar o arquivo. > (arquivo.Rnw). Rodaria se eu colocasse, > text(4,8,expression('se ~x <= 3')) > > Mais da forma como sai (literalmente), eu não quero. Gostaria que saísse > como um operador matemático mesmo. Se alguém souber como, eu agradeço. > Allaman > (S,f,P) > > Ivan Bezerra Allaman > Doutor em Zootecnia/UFLA > email e msn - ivanalaman(a)yahoo.com.br > > ________________________________ > De: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br; Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Enviadas: Terça-feira, 26 de Julho de 2011 10:36 > Assunto: Re: [R-br] Erro ao gerar Operadores matemáticos na compilação > sweave! > > use as ferramentas da forma que elas foram desenhadas para serem usadas. > > text(4,8,expression('se '~x <= 3)) > > b > > 2011/7/26 Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br>: >> Bom dia senhores!! >> Estou tendo um problema ao gerar um símbolo no arquivo compilado com o >> sweave. Quando eu rodo normalmente o script, o símbolo sai direitinho no >> gráfico. No entanto, quando eu compilo o arquivo, o símbolo não sai. >> Alguém >> teria uma solução para isso? Segue um CRM. >> Primeiro apenas no R. >> plot(1:10,1:10) >> text(3,8,expression('se x \u2264 3')) >> #Tudo certo aqui!! >> Agora com o código sweave! >> \documentclass[11pt]{article} >> \usepackage[brazil]{babel} >> \usepackage[T1]{fontenc} >> \usepackage[utf8]{inputenc} >> \usepackage[top=4.5cm,left=4.5cm,right=4.5cm,bottom=4.5cm]{geometry} >> \usepackage{graphicx} >> \usepackage{amsthm} >> \usepackage{amsfonts,amssymb,amsxtra,empheq} >> \usepackage[mathscr]{eucal} >> \begin{document} >> <<fig=T,echo=F>>= >> plot(1:10,1:10) >> text(3,8,expression('se x \u2264 3')) >> @ >> \end{document} >> Percebam no arquivo que o operador matemático não sai. Se alguém puder >> ajudar, agradeço. >> Allaman >> (S,f,P) >> >> >> Ivan Bezerra Allaman >> Doutor em Zootecnia/UFLA >> email e msn - ivanalaman(a)yahoo.com.br >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código >> mínimo reproduzível. >> > > > > -- > Successful people ask better questions, and as a result, they get > better answers. (Tony Robbins) > > > -- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins)

1 0

Erro ao gerar Operadores matemáticos na compilação sweave!
by Ivan Bezerra Allaman 26 Jul '11

26 Jul '11

Bom dia senhores!! Estou tendo um problema ao gerar um símbolo no arquivo compilado com o sweave. Quando eu rodo normalmente o script, o símbolo sai direitinho no gráfico. No entanto, quando eu compilo o arquivo, o símbolo não sai. Alguém teria uma solução para isso? Segue um CRM. Primeiro apenas no R. plot(1:10,1:10) text(3,8,expression('se x \u2264 3')) #Tudo certo aqui!! Agora com o código sweave! \documentclass[11pt]{article} \usepackage[brazil]{babel} \usepackage[T1]{fontenc} \usepackage[utf8]{inputenc} \usepackage[top=4.5cm,left=4.5cm,right=4.5cm,bottom=4.5cm]{geometry} \usepackage{graphicx} \usepackage{amsthm} \usepackage{amsfonts,amssymb,amsxtra,empheq} \usepackage[mathscr]{eucal} \begin{document} <<fig=T,echo=F>>= plot(1:10,1:10) text(3,8,expression('se x \u2264 3')) @ \end{document} Percebam no arquivo que o operador matemático não sai. Se alguém puder ajudar, agradeço. Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

3 2

[Dúvida] Sequência de String no R
by Pedro Rafael 26 Jul '11

26 Jul '11

Pessoal como gerar uma sequencia do tipo A00, A01, A02, ..., A99 no R? Existe algum objeto que faça isso? -- Saudações, Pedro

4 3

ARFIMA
by Ana Paula Sousa 26 Jul '11

26 Jul '11

Alguém já trabalhou com o modelo ARFIMA para séries temporais? Atenciosamente, Ana Sousa

3 3

erro na leitura do banco
by Edson Lira 25 Jul '11

25 Jul '11

Alguém saberia me indicar como corrigir este erro? Erro em type.convert(data[[i]], as.is = as.is[i], dec = dec, na.strings = character(0L)) : string multibyte inválida em '<a0>gua' Estou usando um banco csv convertido de um outro computador. Edson Lira Estatístico Manaus-Amazonas

2 1

[Dúvida] Diferenças entre .RData e .data
by Pedro Rafael 25 Jul '11

25 Jul '11

Pessoal qual a diferença de um arquivo .RData para o .data. Em geral tenho tabelas no trabalho em csv ou outro banco de dados e quando passo para o R salvo sempre em .RData. Qual o melhor formato para salvar os objetos seja eles matrizes, vetores, data.frames no R, .RData ou .data? Pedro

2 1

Filtro em massa (subset)
by Vinicius Brito Rocha 25 Jul '11

25 Jul '11

Srs, tenho uma data.frame onde preciso realizar um filtro do mesmo, a partir de uma das variáveis categóricas. Porém tenho uma lista de restrições dessa variável à filtrar. Escrever uma a uma cada restrição dará muito trabalho. Alguém sabe como posso fazer isso de forma mais inteligente? 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 6

ANOVA - dados desbalanceados
by Vinicius Brito Rocha 25 Jul '11

25 Jul '11

Caros amigos de R, como faço uma ANOVA onde minhas classes são desbalanceadas? 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.

3 4

mjca
by Lincoln Thiago Neves 25 Jul '11

25 Jul '11

Amigos, Gostaria de uma ajuda com relação ao "plot" no output da função mjca de Análise de Correspondência. Estou construindo um mapa usando o plot na saída dessa função, mas queria separar os pontos (pois eles pertencem a grupos diferentes) por símbolos diferentes e gostaria que o valor da contribuição não apareça. Quanto aos símbolos já compreendi um pouco a utilização do "pch" para colocar os pontos da maneira que quero, mas ele sempre muda todos os pontos. E em relação a contribuição eu também tentei utilizar o contrib ="none", mas não tive sucesso. Alguem poderia me auxiliar ou me indicar um bom mateiral que fale com relação a esse tipo de gráfico? Muito obrigado! -- Att; Lincoln T. Neves MSN: lincolnthiago(a)hotmail.com E-mail: lincolnthiago.neves(a)gmail.com Web: https://sites.google.com/site/consultoriaeaulasdeestatistica/home Cel: 8883-3632

2 1

Re: [R-br] Grupo de discussão Latex
by Ivan Bezerra Allaman 24 Jul '11

24 Jul '11

Já cadastrei!! Valeu Laia! Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

1 0

Re: [R-br] Grupo de discussão Latex
by Ivan Bezerra Allaman 24 Jul '11

24 Jul '11

Se os nobres colegas não se importarem, poderíamos discutir latex aqui mesmo, colocando no corpo da mensagem o seguinte: '' [OFF-TOPIC:latex] mensagem ''. Caso não concordarem, dou todo apoio ao Laia. Allaman (S,f,P) Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

6 7

artigo sobre o pacote ScottKnott enviado ao Rjournal
by Enio Jelihovschi 24 Jul '11

24 Jul '11

José Claudio Faria e eu ( Enio Jelihovschi) enviamos um artigo ao Rjournal sobre o Pacote ScottKnott que nós escrevemos. A primeira resposta que recebemos, os referees e o próprio editor, Heather Turner, fizeram sugestões que nós compreendemos e fizemos as mudanças sugeridas. Após o reenvio do artigo recebemos a resposta do editor, que não nem sequer enviou aos referees, recusando o artigo. No nosso entender esta resposta final do editor não faz sentido, e, no nosso entender nem é compreensível do ponto de vista científico. Nós gostaríamos de discutir com a comunidade brasileira do R, escutar as opiniões, para depois enviar ao R-help e ao corpo editorial completo do Rjournal. Estamos enviando em anexo aos interessados cópias dos e-mails do editor assim como nossa resposta final e as cópias da primeira e da segunda versão do artigo. Obrigado Enio Jelihovschi José Claudio Faria

6 7

Re: [R-br] PLOT SERIE TEMPORAL
by Rafael Felice Fan Chen 23 Jul '11

23 Jul '11

Boa tarde, Como eu faço a sazonalidade mensal de uma serie de preços?

2 1

Cluster! Padronização das variáveis...
by Thiago De paula protásio 23 Jul '11

23 Jul '11

Prezados membros da R-Br; Novamente peço a ajuda de vcs...Preciso padronizar as variáveis para posterior análise de agrupamento (cluster). Contudo, gostaria de saber se existe alguma função implementada no R que possa me auxiliar. Muito obrigado pela ajuda! 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

3 2

Re: [R-br] Cluster! Padronização das variáveis...
by Thiago De paula protásio 23 Jul '11

23 Jul '11

Muito obrigado Fernando; Era exatamente isso que eu estava procurando...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

Grupo de discussão Latex
by Alexandre Santos 23 Jul '11

23 Jul '11

Bom Dia Pessoal, Bom minha dúvida não tem haver com R, mas como boa parte dos usuários R usam também Tex, gostaria de saber se alguém saberia me indicar um bom e ativo grupo de discussão de Latex como o r-br, Obrigado, Alexandre

2 1

Re: [R-br] ANOVA - dados desbalanceados
by Ivan Bezerra Allaman 22 Jul '11

22 Jul '11

Bom dia Walmes!! Dá uma olhadinha novamente no seu exemplo de dados extremamente desbalanceados. Acho que há um equívoco na utilização da drop1, pois como o próprio help da função informa, ela testa a hipótese do tipo II, e no caso do seu exemplo, como nós já sabemos, a hipótese mais correta a ser testada é a do tipo IV. Ainda, no help da função tem a seguinte mensagem: ''Only use these functions with data containing missing values with great care''. Pode ser que eu esteja equivocado, por que não? Mais como eu uso diariamente os seus exemplos, e tenho certeza que muitas pessoas também recorrem-nos, acho que vale a pena você dar uma olhadinha lá. 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

4 4

[Dúvida] Rodando Script R no terminal sem entrar no prompt.
by Pedro Rafael 22 Jul '11

22 Jul '11

Pessoal existe alguma forma de rodar o script do R no terminal Linux sem ser necessário "entar no R"? Tipo assim, referenciar o código do R no diretório que ele se encontra e mandar ele rodar sem ser necessário entrar no prompt de comando do R. Obrigado, Pedro Rafael

4 3

Problema em regressão segmentada!!
by Ivan Bezerra Allaman 22 Jul '11

22 Jul '11

Boa tarde senhores!! Estou tendo dificuldades ao ajustar uma regressão segmentada com a library ''segmented'. Não entendo qual a dificuldade que a função está encontrando para ajustar uma regressão aos meus dados. Creio que os breaks pontos chutados não estão tão fora da realidade. Segue abaixo um CMR real. Y <- c(18.335, 19.456, 19.383, 19.542, 20.718, 21.531) X <- c(0,3,6,9,12,15) plot(X,Y,type='b') mod <- lm(Y ~ X) library(segmented) mod.seg <- segmented(mod, seg.Z=~X, psi=list(X=c(3,8))) Se alguém puder ajudar, agradeço. (S,f,P) Allaman Ivan Bezerra Allaman Doutor em Zootecnia/UFLA email e msn - ivanalaman(a)yahoo.com.br

2 1

alterar cores apenas na area grafico
by jose coelho 22 Jul '11

22 Jul '11

Vivam a todos Gostaria de saber com alterar a cor apenas da área do gráfico feita com plot (x,y). já tentei com plot(x,y,bg="gray) não dá erro mas também não funciona!! ao usar par(bg = "gray") fica toda a janela a cinzento e não apenas a área do gráfico!! Como posso fazer isso? Muito obrigado -- jmmcoelho

2 1

Re: [R-br] xtabs ou equivalente
by Jose Claudio Faria 22 Jul '11

22 Jul '11

Test! -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

PLOT SERIE TEMPORAL
by Enéias Façanha 22 Jul '11

22 Jul '11

PessoALL, Utilizando os dados abaixo, como eu faço para plotar essas duas séries, "MILLE_META" e "MILE_RECEITA", no mesmo gráfico, inserindo os meses no eixo x? MILLE_META MILLE_RECEITA JAN 339.155,42 275.649,07 FEV 320.427,98 267.601,66 MAR 333.703,99 279.127,36 ABR 318.293,90 184.393,12 MAI 342.795,95 262.414,92 JUN 328.074,05 268.703,78 Obrigado, -- Enéias Façanha

2 1

Minimização Condicionada
by Luciano Ramos Gonçalves 22 Jul '11

22 Jul '11

Gostaria de saber qual é a melhor função para realizar uma minimização com restrições: nlm(), optim(), optimize() ou constrOptim(), ou outra? Abaixo vai a função objeto, e as restrições descritas. Qualquer dica já ajuda. O meu problema maior foi aplicar as restrições aos métodos de otimização. Obs.: Salvo engano, tudo abaixo é plenamente reproduzível. # Tenho a função abaixo. LogLik <- function(parametros,dados){ alpha <- parametros[1] delta <- parametros[2] epsilon.c <- parametros[3] epsilon.v <- parametros[4] mu <- parametros[5] B <- dados[,2] S <- dados[,3] logl <- sum(-epsilon.c-epsilon.v+((max(c(B,S))+min(c(B,S)))/2)* (log((epsilon.v/(mu+epsilon.v)))+ log((epsilon.c/(mu+epsilon.c))))+ B*log(mu+epsilon.c)+S*log(mu+epsilon.v))+ sum( log( (alpha)*(1-delta)*exp(-mu)* ((((epsilon.v/(mu+epsilon.v))))^ (S-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.c/(mu+epsilon.c))^ (-((max(c(B,S))+min(c(B,S)))/2)))+ (alpha)*(delta)*exp(-mu)* ((((epsilon.c/(mu+epsilon.c))))^ (B-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.v/(mu+epsilon.v))^ (-((max(c(B,S))+min(c(B,S)))/2)))+ (1-alpha)* ((epsilon.c/(mu+epsilon.c))^ (B-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.v/(mu+epsilon.v))^ (S-((max(c(B,S))+min(c(B,S)))/2))) )) return(logl) } # Supondo os valores BS <- cbind(1,500,501) theta <- c(0.4,0.5,40,40,50) # Quero minimizar a função: LogLik(theta,BS) # Com as seguintes restrições, que eu não consegui implementar (as restrições são # aos valores de theta, ou "parametros" na função, em ordem): # 0 <= alpha <= 1 ("alpha" é um valor entre 0 e 1) # 0 <= delta <= 1 ("delta" é um valor entre 0 e 1) # epsilon.c > 20 ("epsilon.c" é maior que 20) # epsilon.v > 20 ("epsilon.v" é maior que 20) # mu > 30 ("mu" é maior que 30) Grato, Luciano.

1 0

[Dúvida] Salvando arquivo CSV sem aspas nas strings
by Pedro Rafael 21 Jul '11

21 Jul '11

Pessoal estou salvando um data.frame com o objeto write.table(objeto, file = "nome_do_arquivo.csv", sep = ";"). Contudo quero que as strings não contenham aspas. Como fazer para que o csv não contenham nas strings? Saudações, Pedro

3 2

xtabs ou equivalente
by Vinicius Brito Rocha 21 Jul '11

21 Jul '11

Senhores, estou usando o xtabs, no intuito de fazer alto parecido com o que se faz em uma tabela dinamica tenho um data.frame com dois fatores e uma variavel resposta continua estou fazendo algo do tipo: xtabs(y~A+B,data=base) tenho como resposta uma tabela com os níveis de A e B e uma estatística dentro de cada célula. Minha pergunta é: Como faço para mudar a estatística do xtabs? média, quantidade, mediana,.. não tem um parâmetro FUN dentro do xtabs. tem algum pacote com uma função equivalente que faça isso? 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

R is Hot!
by Alejandro C. Frery 21 Jul '11

21 Jul '11

Boa noite, Recomendo a leitura da matéria "R is Hot!" http://t.co/MFs3Bz9 Alejandro -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

Incluir formula em um função
by Leo Bastos 20 Jul '11

20 Jul '11

Olá pessoal, Tenho algumas funções que eu gostaria de incluir uma formula como parametro de entrada. Da mesma forma que as funcoes lm e glm. Alguem poderia me indicar um exemplo, ou indicar um artigo/tutorial/link/whatever para que eu possa incluir uma formula com entrada de uma funcao? Por exemplo, suponha que eu tenho um data frame com um vetor de resposta y e algumas covariáveis. E suponha q eu quero uma função q retorne uma estimativa dos coeficientes de uma regressão linear. (sim, eu sei q o lm faz isso, mas eu gostaria de um exemplo de função onde eu possa usar uma formula como parametro de entrada). Uma função com um vetor y e uma matrix X é trivial, mas como fazer usando uma formula como parametro de entrada? Agradeço antecipadamente qq ajuda. Abraço, Leo Bastos

2 1

Fwd: [Dúvida] Ambiente para programação em R
by Pedro Rafael 19 Jul '11

19 Jul '11

Eu já conhecia o RKward e é bem interessante. Recentemente conheci o RStudio e acredito que ele está sendo mantido com mais regularidade ou estou enganado? Pedro ---------- Mensagem encaminhada ---------- De: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> Data: 18 de julho de 2011 08:51 Assunto: Re: [R-br] [Dúvida] Ambiente para programação em R Para: "jcfaria [via R-br]" <ml-node+3675065-198716963-223914(a)n4.nabble.com> Eu já conhecia o RKward e é bem interessante. Recentemente conheci o RStudio e acredito que ele está sendo mantido com mais regularidade ou estou enganado? Pedro Em 18 de julho de 2011 08:44, jcfaria [via R-br] < ml-node+3675065-198716963-223914(a)n4.nabble.com> escreveu: > Eu > > uso o vim. > > > Eu geralmente > > programo com uma janela para o vim e duas rodando o R: uma para as > > páginas de ajuda e outra para o código. > > Gustavo, > > Um opção muito interessante (eu prefiro assim) é habilitar o > vim/vim-r-plugin para abrir os arquivos de help do R dentro do próprio > vim (<leader>rh) em uma nova tab (para alternar entre tabs: gt (next > Tab) gT (prior Tab)) . > > Vantagens: > 1. Fica bem mais fácil para ler o help (sintaxe); > 2. Fica bem mais fácil de navegar no help (marcas, atalhos do vim, etc); > 3. Dá para interagir com o R (com todos os recursos) a partir do vim, > por exemplo, executar linha por linha os exemplos para entender a > função que está no help; > 4. As tabs vão se abrindo sequencialmente (de uma vc pode abrir outra, etc) > > > Teste colocando a linha abaixo no seu .vimrc: > > let vimrplugin_vimpager = "tabnew" " Open help always in new tab > > Abs, > -- > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > Jose Claudio Faria > Estatistica - prof. Pleno > UESC/DCET/Brasil > [hidden email] <http://user/SendEmail.jtp?type=node&node=3675065&i=0> > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > _______________________________________________ > R-br mailing list > [hidden email] <http://user/SendEmail.jtp?type=node&node=3675065&i=1> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Duvida-Ambiente-para-programacao-em-… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

4 3

[Dúvida] Ambiente para programação em R
by Pedro Rafael 18 Jul '11

18 Jul '11

Pessoal qual um bom ambiente para gerenciar projetos bem como para programação em R no Linux? Pedro

4 4

Produzindo um Codigo Minimo Reproduzivel (CMR)
by Benilton Carvalho 18 Jul '11

18 Jul '11

Na esperanca de que melhores CMRs sejam empregados na lista, adaptei um material que descreve como apropriadamente criar um CMR. Para detalhes, visitem: http://www.leg.ufpr.br/doku.php/software:rbr e sigam ate' o final da pagina. Ou, para uma versao "resumida": http://gist.github.com/1088208 b

2 1