Junho 2014 - R-br - DINF UFPR

EMV Gama
by ana paula coelho madeira 18 Jun '14

18 Jun '14

Bom dia! Estou verificando o ajuste da distribuição gama a um conjunto de dados. Tenho utilizado o programa Maple para calcular as estimativas dos parâmestros e feito o restante da análise no R. Gostaria de saber como faço calculo os estimadores de máxima verossimilhança da distribuição gama no R. Desde já agradeço. Att., Ana Paula

2 2

outputs de tabelas
by Daniela Recchia 17 Jun '14

17 Jun '14

Olá pessoal, pesquisei algumas funções para obter tabelas "mais bonitas" (como uma do excel, com linhas e tudo mais) porém só encontrei exemplos - vários - para casos de tabelas mais complexas. Tentei o paxote xtables e também o Hmisc mas não obtenho ainda o output que quero, poderia ser em pdf, html ou o que fosse - mas não dá certo pois tenho apenas uma tabela de frequencias e não um modelo. Por exemplo quero ter um output bonito para essa tabela: > table(groups_fit,groups_mbi) groups_mbi groups_fit [ 9.0,19.0) [19.0,26.4) [26.4,31.9) [31.9,38.0) [25.1,35.0) 31 15 1 0 [35.0,46.0) 50 51 28 4 [46.0,55.5) 11 34 35 33 [55.5,68.0) 2 6 27 39 [68.0,96.1] 0 0 2 9 groups_mbi groups_fit [38.0,53.0] [25.1,35.0) 0 [35.0,46.0) 0 [46.0,55.5) 11 [55.5,68.0) 28 [68.0,96.1] 50 tentei com o latex(tabular) mas ainda nao da certo, continua no mesmo formato. Alguém teria uma sugestão? Obrigada! Daniela -- Daniela Rodrigues Recchia M.Sc. in Statistics Tel. (Mobil): +49 01578 7588382 Germany “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

4 11

análise de reísduo
by Hélio Gallo Rocha 16 Jun '14

16 Jun '14

Boa noite listeiros. Uma dúvida de estatística básica. A questão é que ao analisar o CV dos resíduos entre dois modelos de estimativa de safra e o observado, o CV do modelo que melhor ajustou ficou muito, mas muito alto. Acredito que seja pelo fato de estar analisando resíduos, com valores negativos e positivos, com grande amplitude. Sendo que o modelo que melhor ajustou apresentou média baixa, tornando o CV alto. Pensando nisso, subtrai o menor valor de cada variável aos resíduos, tornando todos os valores positivos. Feito isto, o desvio padrão ficou igual, mas a média aumentou, diminuindo o CV, ficando no mesmo patamar dos estimadores. Este procedimento está correto? Se sim, quem poderia me indicar material que "validasse" este procedimento? Agradeço desde já as opiniões. -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

5 9

Alterar o tipo das variáveis de um data.frame
by Aline Lipsky 16 Jun '14

16 Jun '14

Bom dia prezados colegas, Eu tenho um data.frame com a seguinte estrutura: $ Bonito : int 0 0 0 0 0 0 0 0 0 0 ... $ Burriquete : chr "0" "0" "0" "0" ... $ Cabrinha : int 0 0 0 0 0 0 0 0 0 0 ... $ Cação : int 0 0 0 0 0 0 0 0 0 0 ... $ Cação.anjo : int 0 0 0 0 0 0 0 0 0 0 ... $ Camarão : Factor w/ 437 levels "0","1","10","100",..: 196 1 320 256 16 126 41 204 310 115 ... $ Camarão.barba.ruça: int 0 0 0 0 0 0 0 0 0 0 ... $ Camarão.rosa : int 0 0 0 0 0 0 0 0 0 0 ... $ Camarão.santana : int 0 0 0 0 0 0 0 0 0 0 ... Preciso que todas as variáveis sejam do tipo numérico. Como posso fazer essa alteração? Grata pela atenção, Aline

4 4

Preenchimento e substituição de dados
by Leonardo Monteiro 16 Jun '14

16 Jun '14

Ola, pessoal, tenho um vetor de uma variavel 'x', que varia durante o ano, ou seja, tenho essa variavel de acordo com o numero do dia do ano (1,2,3,..., 366), num total de com 9862 linhas as quais estao presentes valores NA tb.. Aí, eu tenho um vetor com os valores médios dessa variavel por dia do ano ... entao esse vetor tem 366 linhas. Como eu poderia substituir os valores NA, considerando o valor medio de 'x', por dia do ano?? Obrigado, -- Leonardo Amaral Monteiro - Agronomist MSc in Agrometeorology PhD Student in Engineering Agricultural Systems Sao Paulo University - College of Agriculture "Luiz de Queiroz" - Brazil Landscape System Analysis - ZALF Institute, Müncheberg - Germany

3 2

Base de Dados do JSON para o R
by Lincoln Thiago Neves 16 Jun '14

16 Jun '14

Prezados, Gostaria de uma ajuda na manipulação de um arquivo JSON dentro do R. O arquivo tem as variáveis *_id*, *last_seen*, *section_views*, *urls *e *urls_seen *e segue a seguinte estrutura: { "_id" : "2minformatica(a)uol.com.br", "last_seen" : "2014-06-04T23:12:44.426Z", "section_views" : { "blogs．lancenet．com．br/palestra" : 1, "www．lancenet．com．br" : 3, "www．lancenet．com．br/flamengo" : 1, "www．lancenet．com．br/futebol" : 1, "www．lancenet．com．br/minuto" : 1, "www．lancenet．com．br/palmeiras" : 5 }, "urls" : [ "www.lancenet.com.br/minuto/Palmeiras-Alberto_Valentim_0_1148885234.html", " blogs.lancenet.com.br/palestra/2014/06/01/diogo-fica-mais-forte-no-ataque-do-palmeiras ", "www.lancenet.com.br/palmeiras", " www.lancenet.com.br/palmeiras/Uruguaios-citam-saida-Eguren-empresario_0_115… ", " www.lancenet.com.br/futebol/Contrato-Neilton-despede-Santos-Aprendi_0_11476… ", " www.lancenet.com.br/flamengo/Andre-Santos-retorna-Fenerbahce-proposta_0_115… ", " www.lancenet.com.br/palmeiras/Pratto-aguarda-Palmeiras-melhor-aceita_0_1150… " ], "urls_seen" : [ "www.lancenet.com.br/minuto/Palmeiras-Alberto_Valentim_0_1148885234.html", " blogs.lancenet.com.br/palestra/2014/06/01/diogo-fica-mais-forte-no-ataque-do-palmeiras ", "www.lancenet.com.br", "www.lancenet.com.br/palmeiras", " www.lancenet.com.br/palmeiras/Uruguaios-citam-saida-Eguren-empresario_0_115… ", " www.lancenet.com.br/futebol/Contrato-Neilton-despede-Santos-Aprendi_0_11476… ", " www.lancenet.com.br/flamengo/Andre-Santos-retorna-Fenerbahce-proposta_0_115… ", " www.lancenet.com.br/palmeiras/Pratto-aguarda-Palmeiras-melhor-aceita_0_1150… " ] } Eu preciso colocar este dados em uma 'tabela' da seguinte maneira: *_id* *last_seen* *section_views* *2minformatica(a)uol.com.br <2minformatica(a)uol.com.br>* 2014-06-04T23:12:44.426Z "blogs．lancenet．com．br/palestra" : 1, *2minformatica(a)uol.com.br <2minformatica(a)uol.com.br>* 2014-06-04T23:12:44.426Z "www．lancenet．com．br" : 3, *2minformatica(a)uol.com.br <2minformatica(a)uol.com.br>* 2014-06-04T23:12:44.426Z "www．lancenet．com．br/flamengo" : 1, *2minformatica(a)uol.com.br <2minformatica(a)uol.com.br>* 2014-06-04T23:12:44.426Z "www．lancenet．com．br/futebol" : 1, *2minformatica(a)uol.com.br <2minformatica(a)uol.com.br>* 2014-06-04T23:12:44.426Z "www．lancenet．com．br/minuto" : 1, *2minformatica(a)uol.com.br <2minformatica(a)uol.com.br>* 2014-06-04T23:12:44.426Z "www．lancenet．com．br/palmeiras" : 5 A princípio pretendo utilizar apenas as variáveis *_id*, *last_seen *e *section_view*. Desde já agradeço por qualquer ajuda. Abraços. -- Atenciosamente, Lincoln T Neves +55 21 96426-2949 lincolnthiago.neves(a)gmail.com www.lincolntneves.weebly.com

2 2

Encoding usando read.csv.sql()
by Leonardo Aguirre 16 Jun '14

16 Jun '14

Pessoal, Alguém sabe como definir o encoding do arquivo usando a função read.csv.sql? Tenho um banco em CSV com 19Gb e preciso ler já aplicando alguns filtros, e pensei em fazer isso usando a função read.csv.sql, do pacote sqldf. Entretanto, preciso definir o encoding="latin1" ao ler essa base. Obrigado, Leonardo.

5 8

Box-Cox para objetos da classe geotada
by Helder Gramacho 14 Jun '14

14 Jun '14

Boa tarde pessoal, Observando a ajuda da Transformação Box-Cox por meio do help(boxcox), vi que o objeto precisa ser da classe fórmula como no exemplo abaixo: boxcox(Volume ~ log(Height) + log(Girth), data = trees, lambda = seq(-0.25, 0.25, length = 10)) Entretanto, no tutorial disponível na página: http://leg.ufpr.br/geoR/tutorials/ctc.R vi que para objetos da classe geodata é possível realizar o teste Box-Cox, se é que não me enganei: ### <code r> ctc <- read.table("http://www.leg.ufpr.br/geoR/tutorials/datasets/ctc.dat", head=T) require(geoR) ctc <- as.geodata(ctc) require(MASS) boxcox(ctc) ### </code> É isso mesmo? Gostaria de entender como funciona o teste nos dois casos (fórmula e geodata). Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 2

Plotar categorias em clusplot
by Rodrigo Muniz 12 Jun '14

12 Jun '14

Caros membros da lista, gostaria de solicitar a ajuda de vocês mais uma vez, para plotar categorias em gráficos de cluster. Segue um exemplo. No exemplo abaixo temos um data frame com com as variáveis x,y e z, e tais variáveis estão dentro de categorias elencadas de A-F. A partir da função "kmeans" as variáveis foram agrupadas em 5 grupos, e plotados utilizando a função "clusplot". Eu gostaria de plotar no mesmo gráfico algum identificador das categorias, como cores, ou mesmo símbolos. Desculpem pelo exemplo simplista. Desde já agradeço. ex=data.frame(cat=rep(c("A", "B","C","D","E","F"), c(10, 20, 5, 15, 20,30)) , x=rnorm(100), y=rnorm(100), z=rnorm(100)) clus=scale(ex[,2:3]) fit=kmeans(clus, 5) clusplot(clus, fit$cluster, color=T, shade=T, label=2, lines=0) Att... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) E-mail - muniz.ra(a)usp.br Cel (19) 98300-4333 (Pessoal)

1 0

Análise Fatorial Exploratória de experimento em fatorial com parcela subdividida no tempo
by Alisson Lucrecio 12 Jun '14

12 Jun '14

Caro Colegas da r-br, Boa tarde. Eu posso fazer uma análise fatorial de dados de fertilidade de 2 amostragem (no mesmo local) de 3 solos que foram cultivado com cafeeiro com 4 doses lodo de esgoto com 3 blocos? Delineamento fatorial com parcela subdividida no tempo. 2 amostragem em data diferentes x 3 solos x 4 dose de lodo de esgoto com 3 blocos. Obrigado. -- Alisson Lucrecio da Costa

2 2