[R-br] Assunto: Re: JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES
Muito obrigado. Enviado do Yahoo Mail no Android Em qua, 30 30e jan 30e 2019 às 22:56, Marcus Nunes por (R-br) escreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES
Cesar, vc tem o link? Não consegui encontrar. [ ]'sProf. Edson Lira, Me Estatístico Manaus-Amazonas Em terça-feira, 29 de janeiro de 2019 10:21:39 AMT, Cesar Rabak por (R-br) escreveu: Edson, Veja uma trédi recente nesta lista sobre o uso da função merge ela parece servir exatamente para o que você deseja. Na trédi você poderá ver que podes escolher entre a merge da biblioteca base ou a do pacote data.table. HTH--Cesar Rabak On Fri, Jan 25, 2019 at 6:40 PM Edson Lira por (R-br) wrote: Boa tarde a todos!!! Tenho dois arquivos: um com os CEP's de pacientes de um hospital (57000 registros) sem as coordenadas, o outro arquivo que consegui no link abixo http://results.openaddresses.io/?runs=all#runs No link acima tem POSTCODES (CEP"S) de todo o mundo. Baixei o do Amazonas que tem o CEP DE MANAUS (SEPAREI ESSE ARQUIVO com 554000 registros) Quero fazer uma rotina aonde eu possa juntar o arquivo de consultas sem as coordenadas com o arquivo de CEP'S com as coordendas. (LAT LON), ou seja, quero fazer um arquivo mantendo a dimensão do arquivo das consultas. Espero ter sido claro. Abaixo o link dos arquivos (um com CEP e o outro com CEP e Coordenadas) https://www.datafilehost.com/d/b93da0f3 https://www.datafilehost.com/d/b8f65984 [ ]'s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES
Boa tarde a todos!!! Tenho dois arquivos: um com os CEP's de pacientes de um hospital (57000 registros) sem as coordenadas, o outro arquivo que consegui no link abixo http://results.openaddresses.io/?runs=all#runs No link acima tem POSTCODES (CEP"S) de todo o mundo. Baixei o do Amazonas que tem o CEP DE MANAUS (SEPAREI ESSE ARQUIVO com 554000 registros) Quero fazer uma rotina aonde eu possa juntar o arquivo de consultas sem as coordenadas com o arquivo de CEP'S com as coordendas. (LAT LON), ou seja, quero fazer um arquivo mantendo a dimensão do arquivo das consultas. Espero ter sido claro. Abaixo o link dos arquivos (um com CEP e o outro com CEP e Coordenadas) https://www.datafilehost.com/d/b93da0f3 https://www.datafilehost.com/d/b8f65984 [ ]'s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] RES: COMPARAR TRATAMENTOS
Obrigado Antonio Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em quinta-feira, 29 de novembro de 2018 08:56:06 AMT, Antonio Silva escreveu: #yiv8513670777 #yiv8513670777 -- _filtered #yiv8513670777 {panose-1:2 4 5 3 5 4 6 3 2 4;} _filtered #yiv8513670777 {font-family:Calibri;panose-1:2 15 5 2 2 2 4 3 2 4;} _filtered #yiv8513670777 {font-family:Garamond;panose-1:2 2 4 4 3 3 1 1 8 3;}#yiv8513670777 #yiv8513670777 p.yiv8513670777MsoNormal, #yiv8513670777 li.yiv8513670777MsoNormal, #yiv8513670777 div.yiv8513670777MsoNormal {margin:0cm;margin-bottom:.0001pt;font-size:11.0pt;font-family:sans-serif;}#yiv8513670777 a:link, #yiv8513670777 span.yiv8513670777MsoHyperlink {color:blue;text-decoration:underline;}#yiv8513670777 a:visited, #yiv8513670777 span.yiv8513670777MsoHyperlinkFollowed {color:#954F72;text-decoration:underline;}#yiv8513670777 .yiv8513670777MsoChpDefault {} _filtered #yiv8513670777 {margin:70.85pt 3.0cm 70.85pt 3.0cm;}#yiv8513670777 div.yiv8513670777WordSection1 {}#yiv8513670777 Enviado do Email para Windows 10 De: Edson Lira por (R-br) Enviado:quarta-feira, 28 de novembro de 2018 23:46 Para: a lista Brasileira oficial de discussão do programa R. Cc:Edson Lira Assunto: [R-br] COMPARAR TRATAMENTOS Boa noite, os dados abaixo referem-se a 10 tratamentos e ao estagio de uma doença. O grau 2 é o estágio maior. doenca Trat 0 1 2 Total GM1 3 2 1 6 GM2 0 2 4 6 GM3 1 3 2 6 GM4 0 3 2 5 GM5 0 4 2 6 GM6 0 1 5 6 GM7 1 1 4 6 GM8 0 0 6 6 GM9 0 4 1 5 GM10 0 0 6 6 Fiz um teste de associação (Exato de Fisher) e verifiquei que ha uma associação entre o tratamento e o estagio da doença. Agora gostaria de verificar qual dos tratamentos é o mais eficiente. Posso testar também: G1xG2 G3xG4 G5xG6 G7xG8 G9XG10 Alguém tem alguma sugestão? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Caro Edson Eu montaria a matriz desta forma, a seguir, faria uma análise não paramétrica e faria a escala de notas começando do 1, eliminando o dígito 0. Trat rep nota GM0 d1 0 GM1 d1 3 GM2 d1 0 GM3 d1 1 GM4 d1 0 GM5 d1 0 GM6 d1 0 GM7 d1 1 GM8 d1 0 GM9 d1 0 GM10 d1 0 GM0 d2 1 GM1 d2 2 GM2 d2 2 GM3 d2 3 GM4 d2 3 GM5 d2 4 GM6 d2 1 GM7 d2 1 GM8 d2 0 GM9 d2 4 GM10 d2 0 GM0 d3 2 GM1 d3 1 GM2 d3 4 GM3 d3 2 GM4 d3 2 GM5 d3 2 GM6 d3 5 GM7 d3 4 GM8 d3 6 GM9 d3 1 GM10 d3 6 Espero ter ajudado. Brito. ___ R-br mailing list R-br@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.
[R-br] COMPARAR TRATAMENTOS
Boa noite, os dados abaixo referem-se a 10 tratamentos e ao estagio de uma doença. O grau 2 é o estágio maior. doencaTrat 0 1 2 Total GM1 3 2 1 6GM2 0 2 4 6GM3 1 3 2 6GM4 0 3 2 5GM5 0 4 2 6GM6 0 1 5 6GM7 1 1 4 6GM8 0 0 6 6GM9 0 4 1 5GM10 0 0 6 6 Fiz um teste de associação (Exato de Fisher) e verifiquei que ha uma associação entre o tratamento e o estagio da doença. Agora gostaria de verificar qual dos tratamentos é o mais eficiente. Posso testar também: G1xG2G3xG4G5xG6G7xG8G9XG10 Alguém tem alguma sugestão? [ ]'s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] Eliminar NA no gráfico
Boa noite a todos, usando a rotina abaixo: library(ggpubr)p <- ggboxplot(hg2, x = "Grau", y = "glicose", color = "Grau", palette = "jco", add = "jitter",na.rm=T)p + stat_compare_means(aes(group = Grau), na.rm=T) Tenho a seguinte mensagem de erro: Don't know how to automatically pick scale for object of type quosure/formula. Defaulting to continuous.Don't know how to automatically pick scale for object of type quosure/formula. Defaulting to continuous.Error in validDetails.text(x) : objeto do tipo 'pairlist' não é coercível para tipo 'double'Além disso: Warning message:Removed 2 rows containing missing values (geom_point). Quando ploto o gráfico abaixo, aparece NA, como faço para eliminar o NA e qual o problema com a minha primeira rotina? O banco está nesse link:https://www.datafilehost.com/d/15ca7130 [ ]'sProf. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Assunto: Classificar um conjunto de dados
###CONTAGEM DE DOADORESCOM ESSA ROTINA ABAIXO, VC ORDENA DO MAIOR PARA O MAIOR AS DUAS VARIÁVEISd_cad0_10<- tip_doa[order(tip_doa$NUMEROGENTE,tip_doa$DATAHORADOACAO, decreasing=TRUE) ,] NA ROTINA ABAIXO VC MANTEM AS DUPLICADAS E MANTEM SOMENTE A DE MAIOR FREQUÊNCIAdc_0_10<- d_cad0_10[!duplicated(d_cad0_10$NUMEROGENTE),] USEI ESSA ROTINA EM UM BANCO DE DOADORES AQUI DE MANAUS PARA MANTER SOMENTE SOMENTE ELE UMA VEZ, E COMO ORDENEI ELE MANTINHA O PRIMEIRO REGISTRO DOS DUPLICADOS E O PRIMEIRO ERA O DE MAIOR FREQUÊNCIA. ACHO QUE ERA ISSO. TESTE ESSA. [ ]'sProf. Edson Lira, Me Estatístico Manaus-Amazonas Em quinta-feira, 8 de novembro de 2018 18:29:43 AMT, Fernando Souza por (R-br) escreveu: Cara criei esse comando que aproxima bem do que você deseja,Porém a última situaçao (Local=D) não estou conseguindo programar, veja se tem alguma idéia que possa complementar. att dados<-data.frame(Local=c('A','A','A','B','C','C','C','C','D','D'),Frequencia=c(2,2,1,3,5,5,5,5,3,4),Data=c("13/05/2015","13/05/2015","13/05/2015","13/05/2015","13/05/2015","10/05/2015","13/05/2015","13/05/2015","13/05/2015","13/05/2015"),qtde=c(500,700,700,400,300,400,400,900,550,500))dados$Data <-as.Date(dados$Data,"%d/%m/%Y")str(dados)dados # install.packages("plyr")library(plyr) arranjo<-function(x){ valido <- c() for(i in 1:nrow(x)){ if(x[i,2]==max(x[,2])){ valido[i]<-1 if( x[i,3]>= max(x[,3]) & x[i,2]==max(x[,2])){ valido[i] <-1 }else{ valido[i] <- 0 } if(x[i,4]>= max(x[,4])){ valido[i] <- 1 }else{ valido[i] <- 0 } }else{ valido[i]=0 } } return(mutate(x,valido=valido))} ddply(dados,.(Local, Data),arranjo) On Nov 7 2018, at 11:29 am, Yury Duarte por (R-br) wrote: Olá Daniel! Utilizei a abordagem que você sugeriu e criei dois laços para fazer a checagem.Ainda não está exatamente como gostaria mas o filtro já está sendo feito da maneira que preciso.Segue a estrutura utilizada. Obrigado! Abraços dados:area_id dia mes ano frequencia valora 15 12 2012 10 5a 22 12 2012 15 10a 15 12 2012 12 15a 24 11 2013 15 20b 8 12 2013 10 5b 14 10 2012 12 10c 7 11 2012 15 5c 19 11 2013 12 10c 19 10 2014 12 15c 19 10 2013 15 20d 5 10 2014 10 5d 23 10 2014 10 10d 5 12 2012 10 15d 5 12 2012 15 20d 5 12 2013 12 25d 6 11 2014 12 30 estrutura:#Remover Objetos do R#rm(list = ls())set.seed(29031989) raiz = 'C:\\Users\\Desktop\\'a = read.table(paste0(raiz, 'teste_classificacao_final.txt'), header = T, sep = '\t')colunas = colnames(a)locais = data.frame(unique(a$area_id)) for (i in 1:length(locais[,1])){ #nome = paste0("area_", locais[i,1]) #area = assign(nome, a[which(a$area_id==locais[i,1]),]) area = a[which(a$area_id==locais[i,1]),] anos = data.frame(unique(area$ano)) for (x in 1:length(anos[,1])) { area_ano = area[which(area$ano==anos[x,1]),] frequencia_max = max(area_ano$frequencia) valor_max = max(area_ano$valor) area_ano_freq = area_ano[which(area_ano$frequencia==frequencia_max),] area_ano_freq_val = area_ano_freq[which(area_ano_freq$frequencia==frequencia_max),] write.table(area_ano_freq_val, paste0(raiz, 'teste', '.txt'), sep = "\t", append = T, col.names = F, row.names = F) }} b = read.table(paste0(raiz, 'teste.txt'), header = F, sep = '\t')colnames(b) = colunas Yury DuarteEngenheiro Agrônomo - ESALQ/USP Em qua, 7 de nov de 2018 às 09:26, Daniel Guimarães Tiezzi por (R-br) escreveu: Uma opção eh usar Uma estrutura de repetição para checar as variaveis e retornar o valor q deseja Daniel On Wed, 7 Nov 2018, 09:22 Yury Duarte por (R-br) escreveu: Dá uma olhada no aggregate Enviado do Yahoo Mail no Android Em ter, 6 6e nov 6e 2018 às 15:25, Yury Duarte por (R-br) escreveu:Boa tarde colegas listeiros, como vão? Tenho um problema e estou com certa dificuldade em resolver de forma prática. Meu conjunto de dados é o seguinte: | local | frequencia | data | quantidade | valido | | A | 2 | 13/05/2015 | 500,00 | 0 | | A | 2 | 13/05/2015 | 700,00 | 1 | | A | 1 | 13/05/2015 | 700,00 | 0 | | B | 3 | 13/05/2015 | 400,00 | 1 | | C | 5 | 13/05/2015 | 300,00 | 0 | | C | 5 | 10/05/2015 | 400,00 | 0 | | C | 5 | 13/05/2016 | 400,00 | 1 | | C | 5 | 13/05/2015 | 900,00 | 1 | | D | 3 | 13/05/2015 | 550,00 | 0 | | D | 4 | 13/05/2015 | 500,00 | 1 | Para esse conjunto, gostaria de gerar a coluna 'valida', que me classifica os dados da seguinte maneira: para cada local e para cada ano, retornar apenas o local de maior frequência, mais recente e de maior quantidade (nessa ordem).Pen
[R-br] Assunto: Classificar um conjunto de dados
Dá uma olhada no aggregate Enviado do Yahoo Mail no Android Em ter, 6 6e nov 6e 2018 às 15:25, Yury Duarte por (R-br) escreveu: Boa tarde colegas listeiros, como vão? Tenho um problema e estou com certa dificuldade em resolver de forma prática. Meu conjunto de dados é o seguinte: | local | frequencia | data | quantidade | valido | | A | 2 | 13/05/2015 | 500,00 | 0 | | A | 2 | 13/05/2015 | 700,00 | 1 | | A | 1 | 13/05/2015 | 700,00 | 0 | | B | 3 | 13/05/2015 | 400,00 | 1 | | C | 5 | 13/05/2015 | 300,00 | 0 | | C | 5 | 10/05/2015 | 400,00 | 0 | | C | 5 | 13/05/2016 | 400,00 | 1 | | C | 5 | 13/05/2015 | 900,00 | 1 | | D | 3 | 13/05/2015 | 550,00 | 0 | | D | 4 | 13/05/2015 | 500,00 | 1 | Para esse conjunto, gostaria de gerar a coluna 'valida', que me classifica os dados da seguinte maneira: para cada local e para cada ano, retornar apenas o local de maior frequência, mais recente e de maior quantidade (nessa ordem).Pensei em gerar um objeto para cada local a partir do which(), mas como meu conjunto de dados tem mais de 150 mil linhas e mais de 2000 locais, achei que essa abordagem não seria muito prática nem rápida de ser processada. Desde já, agradeço pela colaboração de todos! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
[R-br] Comparar semanas (experimento longitudinal)
Boa noite listeiros! Estou trabalhando com um bd (o pesquisador não conversou comigo antes) aonde preciso comparar a evolução de uma medida ao longo do tempo. São 10 grupos com 5 tratamentos, em determinados grupos tenho 5 ou 6 sujeitos. As semanas são: 1 e em alguns grupos vai a 7,8, 15 ou 16 semanas os grupos que têm 15, 16 semanas, tem um tratamento até a 9 semana, e outro tratamento até a 15 ou 16 semana. O objetivo do pesquisador é comparar a evolução e o efeito dos tratamentos. Alguém me sugere uma rotina para fazer a evolução entre a primeira e última semana? O banco está no link abaixo. https://www.datafilehost.com/d/f906462d [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] SOMA DE LINHAS
Bom dia a todos! Tenho um banco aonde da coluna 19 a coluna 29 tenho respostas do tipo: Concordância, Indiferente e Discordância, com um banco de 29 linhas. O meu objetivo é somar por linha (por cada indivíduo, a quantidade e o percentual de Concordância), ou seja acrescentar mais duas colunas por indivíduo. Abaixo as respostas de um indivíduo (uma linha). Concordância Concordância Concordância Concordância Concordância Concordância Concordância Discordância Concordância Concordância Concordância Discordância Indiferente Indiferente Indiferente Concordância Concordância Discordância Indiferente Indiferente Indiferente [ ]'sProf. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] Dados ausaentes (NA)
Boa tarde, a todos, estou com um problema com essa tabela, dentre outras.Vejam que há 38 observações em branco (seriam NAs). $DispSecNao Qtde Perc 38 79.17DispAcionamento, DispParadaEmergencia, DispInterTrav 1 2.08DispParadaEmergencia 1 2.08DispParadaEmergencia, DispInterTrav 7 14.58DispPartida, DispParadaEmergencia, DispInterTrav 1 2.08Total 48 100.00 Usei a rotina abaixo para eliminar os NAs. b<-read.csv2("C:\\Users\\Edson Lira\\Google Drive\\CONSULTA\\CLAUDIO_SILVIO\\banco1.csv")b[b==""]<-NA Para outras tabelas que existem registro em branco, o NA foi eliminado das frequências, exceto na tabela acima. Alguém tem alguma sugestão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Planejamento e Análise de Experimentos
Faça a busca no Google e você verá uma lista de possibilidades R BOOK FREE DOWNLOAD Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em sexta-feira, 15 de junho de 2018 19:22:13 AMT, Mauro Sznelwar via R-br escreveu: Tente de novo, já pedi para compartilhar com o seu e-mail e do grupo. Boa tarde! Eu não consegui acessar o livro, o link não funcionou. Obrigada! Gislaine Em 14/06/18 22:28, Mauro Sznelwar via R-br escreveu: Veja este livro lhe é útil, tem um capítulo grande sobre Análise de Variância, que é planejamento de experimentos.https://www.dropbox.com/home?preview=The+R+Book.pdf Prezados, boa tarde. Estou estudando o assunto e gostaria de fazê-lo utilizando o R. Alguém poderia indicar um bom livro/material sobre planejamento e análise de experimentos (preferencialmente utilizando o R) voltado para área de Farmácia? Agradeço qualquer indicação. Att., Emerson___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] INDICADORES SOCIAIS OF TOPICS
Caros amigos, boa tarde, um amigo meu Economista me fez uma pergunta que resolvi compartilhar com vocês. "Existe alguma técnica em Estatística, específica para tratamento de Indicadores Sociais?" Eu diria que não, mas gostaria de saber se vocês conhecem algo do tipo. [ ['s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] Assunto: Re: TRABALHAR COM O CEP PARA CRIAR LATITUTDE E LONGITUDE EM UM BANCO
Obrigado Nani e Elias. Elias, as linhas são duplicadas em relação aos bairros, pois cada linha é um prontuário. Enviado do Yahoo Mail no Android Em sáb, 19 19e mai 19e 2018 às 8:02, nani via R-brescreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
[R-br] TRABALHAR COM O CEP PARA CRIAR LATITUTDE E LONGITUDE EM UM BANCO
Boa tarde Listeiros, estou estudando a elaboração de mapas e me deparei com uma situação que ainda não consegui resolver. Se eu tenho um banco de dados com CEPs, como faço para a partir do CEP criar as variáveis lat lomg. No google mapas tem esse serviço, quando eu tenho um banco pequeno é fácil de resolver, mas quando o banco é grande, só para assinantes. Suponha um banco com 16000 linhas, ai é dificultoso. Eu pesquisei e encontrei um pacote chamado cepR, mas não consegui com o material que está disponível resolver o problema. Alguém poderia sugerir algo? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] Fw: MAPAS
Caros amigos, reenviando. Prof. Edson Lira, Me Estatístico Manaus-Amazonas - Mensagem encaminhada - De: Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>Para: Cleber N.Borges via R-br <r-br@listas.c3sl.ufpr.br>Enviado: terça-feira, 17 de abril de 2018 11:24:20 GMT-4Assunto: [R-br] MAPAS Bom dia caros amigos, nesses links tem o script e o banco que estou usando para confeccionar mapas de sarampo no município de Manaus. Eu consegui plotar o mapa por bairro. Plotei os mapas de cada região (6) separadas. Meu problema: quero montar um mapa único das 6 regiões , no script vcs verão que são 6 zonas geográficas separadas. Só não consegui montar u mapa único das 6 regiões, sem plotar os nomes dos bairros aonde o meu poligono seria as 6 regiões. SCRIPT Download mapa mao sar.R | | | | Download mapa mao sar.R | | | BANCO Download manaus sar 63 2016.xlsx | | | | Download manaus sar 63 2016.xlsx | | | [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] MAPAS
Bom dia caros amigos, nesses links tem o script e o banco que estou usando para confeccionar mapas de sarampo no município de Manaus. Eu consegui plotar o mapa por bairro. Plotei os mapas de cada região (6) separadas. Meu problema: quero montar um mapa único das 6 regiões , no script vcs verão que são 6 zonas geográficas separadas. Só não consegui montar u mapa único das 6 regiões, sem plotar os nomes dos bairros aonde o meu poligono seria as 6 regiões. SCRIPT Download mapa mao sar.R | | | | Download mapa mao sar.R | | | BANCO Download manaus sar 63 2016.xlsx | | | | Download manaus sar 63 2016.xlsx | | | [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] Assunto: Re: dificuldade em carregar pacote já instalado
Use dep=T Outra forma, é vc baixar o pacote para uma pasta no seu PC e instalar pacotes ZIP, se não me engano Enviado do Yahoo Mail no Android Em sáb, 31 31e mar 31e 2018 às 17:47, Jobenil Luiz Magalhães Júnior via R-brr-br@listas.c3sl.ufpr.br> escreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Pacote ExpDes.pt
Como você resolveu? Até para sabermos.qual sugestão foi implementada por vc. Só uma observação. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em quinta-feira, 15 de março de 2018 16:10:02 GMT-4, Lucas Stempkowski <lucas_stempkow...@hotmail.com> escreveu: #yiv4488972980 #yiv4488972980 -- P {margin-top:0;margin-bottom:0;}#yiv4488972980 Resolvido, obrigado a todos. -- - Lucas Antonio Stempkowski Eng. Agronômo - UFFS Mestrando em Produção Vegetal CAV/UDESC +55 54 99153-0124/ 98448-4633 -- De: Edson Lira <edinhoes...@yahoo.com.br> Enviado: segunda-feira, 12 de março de 2018 10:23:17 Para: Rafael Pertille; Lucas Stempkowski; a lista Brasileira oficial de discussão do programa R.; a lista Brasileira oficial de discussão do programa R. Assunto: Re: [R-br] Pacote ExpDes.pt Parece que o pacote não está instalado, ou você não usou;library(ExpDes.pt) [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Domingo, 11 de Março de 2018 19:53, Rafael Pertille via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Confere se está instalado o pacote, pois o nome está certo... Em 11 de março de 2018 20:13, Lucas Stempkowski via R-br<r-br@listas.c3sl.ufpr.br> escreveu: Olá, Estou tentando rodar o pacote ExpDes.pt , porém como resposta tenho: Error in library(ExpDes.pt) : there is no package called ‘ExpDes.pt’ Alguém pode me ajudar?Obrigado ---Lucas Antonio StempkowskiEng. Agronômo - UFFSMestrando em Produção Vegetal CAV/UDESC+55 54 99153-0124/ 98448-4633-- __ _ R-br mailing list R-br@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. -- Rafael Henrique PertilleTécnico em AgropecuáriaAcadêmico de Agronomia - UTFPR Campus Pato Branco (46) 99770049___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Instalação do pacote brms
As vezes eu baixo o pacote que necessito zipado e instalo de um pasta local. Bom dia a todos!!! Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em quarta-feira, 14 de março de 2018 17:44:48 GMT-4, Jônatan via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Olá André, no seu e-mail só apareceram avisos durante a instalação do pacote, e as suas mensagens foram as mesmas que apareceram na minha sessão.A minha instalação foi demorada, você aguardou até o final? ela acusou algum erro, além daqueles warnings de sua mensagem? Instalei com sucesso o pacote rstan. Coloquei as saídas (muito extensa) nesse gist: https://gist.github.com/lhmet/9b715830ea226545aefa85d845620a40 Testei o exemplo da vinheta e os reproduzi com sucesso. Sugiro que verifique também se tens todos pré-requisitos reportados em https://github.com/stan-dev/rstan/wiki/Installing-RStan-on-Mac-or-Linux Att. sessionInfo() R version 3.4.3 (2017-11-30) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 14.04.5 LTS Matrix products: default BLAS: /usr/lib/openblas-base/libblas.so.3 LAPACK: /usr/lib/lapack/liblapack.so.3.0 locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=pt_BR.UTF-8 LC_NAME=pt_BR.UTF-8 [9] LC_ADDRESS=pt_BR.UTF-8 LC_TELEPHONE=pt_BR.UTF-8 [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=pt_BR.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] rstan_2.17.3 StanHeaders_2.17.2 inmetr_0.2.6.9000 ggmap_2.7.900 [5] forcats_0.2.0 stringr_1.3.0 dplyr_0.7.4 purrr_0.2.4 [9] readr_1.2.0 tidyr_0.8.0 tibble_1.4.2 ggplot2_2.2.1.9000 [13] tidyverse_1.2.1 loaded via a namespace (and not attached): [1] tagcloud_0.6 TH.data_1.0-8 colorspace_1.3-3 rjson_0.2.15 [5] gdalUtils_2.0.1.13 easypackages_0.1.0 modeltools_0.2-21 . . . (omitido) ver gist. Att. 2018-03-14 7:53 GMT-03:00 Andre Oliveira via R-br <r-br@listas.c3sl.ufpr.br>: boa dia,estou a procura de um pacote para estimar componentes de variâncias via Bayesiana. Já testes alguns outros recomendados por professores, mas, me falaram da brms. Acontece que são tantos erros que não consegui instalar no Ubuntu. A rstan é uma dependência que não instala por nada. Alguém tem algumas sugestões? R version 3.4.3 (2017-11-30) -- "Kite-Eating Tree" Copyright (C) 2017 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R é um software livre e vem sem GARANTIA ALGUMA. Você pode redistribuí-lo sob certas circunstâncias. Digite 'license()' ou 'licence()' para detalhes de distribuição. R é um projeto colaborativo com muitos contribuidores. Digite 'contributors()' para obter mais informações e 'citation()' para saber como citar o R ou pacotes do R em publicações. Digite 'demo()' para demonstrações, 'help()' para o sistema on-line de ajuda, ou 'help.start()' para abrir o sistema de ajuda em HTML no seu navegador. Digite 'q()' para sair do R. [Workspace loaded from ~/.RData] > install.packages("rstan") Installing package into ‘/home/andre/R/x86_64-pc- linux-gnu-library/3.4’ (as ‘lib’ is unspecified) tentando a URL 'https://cran-r.c3sl.ufpr.br/ src/contrib/rstan_2.17.3.tar. gz' Content type 'application/x-gzip' length 747202 bytes (729 KB) == downloaded 729 KB * installing *source* package ‘rstan’ ... ** package ‘rstan’ successfully unpacked and MD5 sums checked ** libs g++ -I/usr/share/R/include -DNDEBUG -I"../inst/include" -I"../inst/include/boost_not_ in_BH" -I"`"/usr/lib/R/bin/Rscript" --vanilla -e "cat(system.file('include', 'src', package = 'StanHeaders'))"`" -DBOOST_RESULT_OF_USE_TR1 -DBOOST_NO_DECLTYPE -DBOOST_DISABLE_ASSERTS -DFUSION_MAX_VECTOR_SIZE=12 -DBOOST_PHOENIX_NO_VARIADIC_ EXPRESSION -I"/home/andre/R/x86_64-pc- linux-gnu-library/3.4/Rcpp/ include" -I"/home/andre/R/x86_64-pc- linux-gnu-library/3.4/ RcppEigen/include" -I"/home/andre/R/x86_64-pc- linux-gnu-library/3.4/BH/ include" -I"/home/andre/R/x86_64-pc- linux-gnu-library/3.4/ StanHeaders/include" -fpic -g -O2 -fstack-protector --param=ssp-buffer-size=4 -Wformat -Werror=format-security -D_FORTIFY_SOURCE=2 -g -c chains.cpp -o chains.o .In file included from /home/andre/R/x86_64-pc-linux- gnu-library/3.4/BH/include/ boost/spirit/include/phoenix_ limits.hpp:11:0, from /home/andre/R/x86_64-pc-linux- gnu-library/3.4/BH/include/ boost/spirit/home/support/ meta_compiler.hpp:16, from /home/andre/R/x86_64-pc-linux- gnu-library/3.4/BH/include/ boost/spirit/home/qi/meta_ compiler.hpp:14,
Re: [R-br] Pacote ExpDes.pt
Parece que o pacote não está instalado, ou você não usou; library(ExpDes.pt) [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Domingo, 11 de Março de 2018 19:53, Rafael Pertille via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Confere se está instalado o pacote, pois o nome está certo... Em 11 de março de 2018 20:13, Lucas Stempkowski via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Olá, Estou tentando rodar o pacote ExpDes.pt , porém como resposta tenho: Error in library(ExpDes.pt) : there is no package called ‘ExpDes.pt’ Alguém pode me ajudar?Obrigado ---Lucas Antonio StempkowskiEng. Agronômo - UFFSMestrando em Produção Vegetal CAV/UDESC+55 54 99153-0124/ 98448-4633-- __ _ R-br mailing list R-br@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. -- Rafael Henrique PertilleTécnico em AgropecuáriaAcadêmico de Agronomia - UTFPR Campus Pato Branco (46) 99770049___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Assunto OFF
Isso César, gostaria de calibrar o processo. e a partir daí monitora-lo. [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Segunda-feira, 5 de Fevereiro de 2018 16:06, Cesar Rabak <cesar.ra...@gmail.com> escreveu: OK. Você quer criar uma curva de "calibração" do processo. Até aqui estou usando a palavra "processo" no sentido industrial do termo. O quê você precisa fazer Edson, é modelá-lo em um processo estatístico. A menos que o produto (no sentido da saída do processo "industrial") seja uma novação inaudita, provavelmente com pesquisa sobro o domínio do problema você encontrará alguma indicação de qual seria a modelagem mais apropriada e a partir disso você pode tentar criar uma estimação das perdas esperadas para esse processo ("industrial"), admitindo, obviamente que ele esteja OK. HTH 2018-02-05 18:18 GMT-02:00 Edson Lira <edinhoes...@yahoo.com.br>: Caros, eu só tenho disponível os dados referentes a janeiro de 2016 a abril de 2017. Gostaria de construir um gráfico (fase 1) para estabelecer um padrão de perdas no processo. O período de agosto a janeiro de 2017 apresenta valores muito acima da média dos outros dois períodos (janeiro a julho de 2016 e fevereiro a abril de 2017). Tentei ajustar uma média ponderada para os períodos, estimei um valor médio que substitui nos valores que se apresentam muito maiores do que o normal. Porém., ao calcular as MR incorro em diferenças zeradas e esse fato influencia nos cálculos, alguém tem alguma sugestão ? Nesse período, janeiro de 2016 a abril de 2017, pretendo estimar um valor padrão de perdas. Muito obrigado a todos!!! [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Segunda-feira, 5 de Fevereiro de 2018 14:10, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Prez@do Edson, SE você se refere a um "gráfico de controle" no sentido usado em Qualidade, gráfico de média com uma medida de dispersão e no caso geral, assume-se que o processo siga a "Lei Normal" (gaussiana). O problema estaria no seu processo que "visualmente" estaria "fora de controle". A sua dúvida advém do resultado ou pelo fato de o seu processo deveria estar correto e precisa gerar um gráfico de controle apropriado para indicar essa condição? --Cesar Rabak 2018-02-03 23:30 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: Caros amigos, boa noite, Estou com um problema na elaboração de um gráfico de controle. Estou analisando perdas em um período de 16 meses. Abaixo os dados. jan/16 fev/16 mar/16 abr/16 mai/16 jun/16 jul/16 2.096,23 3.532,09 1.995,22 2.294,34 6.471,51 6.161,56 3.687,96 ago/16 set/16 out/16 nov/16 dez/16 jan/17 30.469,19 52.931,60 39.768,36 21.855,01 16.399,77 17.521,40 fev/17 mar/17 abr/17 7.673,57 1.125,67 1.652,13 Vejam que de janeiro a julho de 2016, tenho uma grandeza, de agosto de 2016 a janeiro de 2017 a grandeza é bem maior e de fev de 2017 a abril a grandeza está na proporção do primeiro período. Tentei ajustar usando médias ponderadas pelas grandezas, porém ao usar a média encontrada no período divergente, ao calcular as diferenças para a média da amplitude tenho valores iguais e as diferenças zeram no período que uso a média ponderada e então fica difícil estabelecer os limites. Os limites são calculados da seguinte forma LSC=xbar +2*(MR/d2 ) d2=1,128 LSC=xbar - 2*(MR/d2) MR= soma(xi-xi-1)/m m= número de diferenças Alguém teria uma idéia de como ajustar os dados para fazer um gráfico de controle adequando? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Assunto OFF
Caros, eu só tenho disponível os dados referentes a janeiro de 2016 a abril de 2017. Gostaria de construir um gráfico (fase 1) para estabelecer um padrão de perdas no processo. O período de agosto a janeiro de 2017 apresenta valores muito acima da média dos outros dois períodos (janeiro a julho de 2016 e fevereiro a abril de 2017). Tentei ajustar uma média ponderada para os períodos, estimei um valor médio que substitui nos valores que se apresentam muito maiores do que o normal. Porém., ao calcular as MR incorro em diferenças zeradas e esse fato influencia nos cálculos, alguém tem alguma sugestão ? Nesse período, janeiro de 2016 a abril de 2017, pretendo estimar um valor padrão de perdas. Muito obrigado a todos!!! [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Segunda-feira, 5 de Fevereiro de 2018 14:10, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Prez@do Edson, SE você se refere a um "gráfico de controle" no sentido usado em Qualidade, gráfico de média com uma medida de dispersão e no caso geral, assume-se que o processo siga a "Lei Normal" (gaussiana). O problema estaria no seu processo que "visualmente" estaria "fora de controle". A sua dúvida advém do resultado ou pelo fato de o seu processo deveria estar correto e precisa gerar um gráfico de controle apropriado para indicar essa condição? --Cesar Rabak 2018-02-03 23:30 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: Caros amigos, boa noite, Estou com um problema na elaboração de um gráfico de controle. Estou analisando perdas em um período de 16 meses. Abaixo os dados. jan/16 fev/16 mar/16 abr/16 mai/16 jun/16 jul/16 2.096,23 3.532,09 1.995,22 2.294,34 6.471,51 6.161,56 3.687,96 ago/16 set/16 out/16 nov/16 dez/16 jan/17 30.469,19 52.931,60 39.768,36 21.855,01 16.399,77 17.521,40 fev/17 mar/17 abr/17 7.673,57 1.125,67 1.652,13 Vejam que de janeiro a julho de 2016, tenho uma grandeza, de agosto de 2016 a janeiro de 2017 a grandeza é bem maior e de fev de 2017 a abril a grandeza está na proporção do primeiro período. Tentei ajustar usando médias ponderadas pelas grandezas, porém ao usar a média encontrada no período divergente, ao calcular as diferenças para a média da amplitude tenho valores iguais e as diferenças zeram no período que uso a média ponderada e então fica difícil estabelecer os limites. Os limites são calculados da seguinte forma LSC=xbar +2*(MR/d2 ) d2=1,128 LSC=xbar - 2*(MR/d2) MR= soma(xi-xi-1)/m m= número de diferenças Alguém teria uma idéia de como ajustar os dados para fazer um gráfico de controle adequando? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Assunto OFF
Caros amigos, boa noite, Estou com um problema na elaboração de um gráfico de controle. Estou analisando perdas em um período de 16 meses. Abaixo os dados. jan/16 fev/16 mar/16 abr/16 mai/16 jun/16 jul/16 2.096,23 3.532,09 1.995,22 2.294,34 6.471,51 6.161,56 3.687,96 ago/16 set/16 out/16 nov/16 dez/16 jan/17 30.469,19 52.931,60 39.768,36 21.855,01 16.399,77 17.521,40 fev/17 mar/17 abr/17 7.673,57 1.125,67 1.652,13 Vejam que de janeiro a julho de 2016, tenho uma grandeza, de agosto de 2016 a janeiro de 2017 a grandeza é bem maior e de fev de 2017 a abril a grandeza está na proporção do primeiro período. Tentei ajustar usando médias ponderadas pelas grandezas, porém ao usar a média encontrada no período divergente, ao calcular as diferenças para a média da amplitude tenho valores iguais e as diferenças zeram no período que uso a média ponderada e então fica difícil estabelecer os limites. Os limites são calculados da seguinte forma LSC=xbar +2*(MR/d2 ) d2=1,128 LSC=xbar - 2*(MR/d2) MR= soma(xi-xi-1)/m m= número de diferenças Alguém teria uma idéia de como ajustar os dados para fazer um gráfico de controle adequando? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Recode rotina
Valeu Fernando, outra boa sugestão a ser testada [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Quarta-feira, 31 de Janeiro de 2018 15:44, Fernando Souza <nandodeso...@gmail.com> escreveu: A função revalue do pacote plyr é ótima para isso dados <- expand.grid(Resposta=c("otimo","bom","regular","ruim","NaoUtilizou","pessimo"),Questao=1:10) head(dados) library(plyr) dados$novaclass <- revalue(dados$Resposta, c("otimo"="positiva","bom" = "positiva","regular"="","NaoUtilizou" = "","ruim" ="negativa", "pessimo" = "negativa")) Em 31 de janeiro de 2018 17:54, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Obrigado Thiago e Leonard, vou testar a sugestão de vcs. [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em quarta-feira, 31 de janeiro de 2018 13:12:19 GMT-4, Leonard Assis via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Eu, particularmente, usaria o recode, mas jogaria em uma outra variável Leonard Em 31 de jan de 2018 1:48 PM, "Tiago Fragoso via R-br" <r-br@listas.c3sl.ufpr.br> escreveu: Olá Você pode usar a função 'mapvalues' do pacote plyr também. No seu caso, ficaria ads <- plyr::mapvalues(ads, from = c('péssimo','ruim','otimo',' bom','regular'), to = c('Negativa','Negativa',' Positiva','Positiva','')) 2018-01-31 13:05 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: > Já resolvi, com a rotina abaixo. > > ads <- apply(ads, 2, function(x) {x <- recode(x,"'pessimo'='Negativa' ; > 'ruim'='Negativa'; > 'otimo'='Positiva';'bom'=' Positiva';'nao > utilizou o servico'='';'regular'=''"); x}) > > Se alguém tiver uma solução mais elegante, eu gostaria de saber. > > Bom dia a todos! > > [ ]'s. > Prof. Edson Lira, Me > Estatístico > Manaus-Amazonas > > > Em Quarta-feira, 31 de Janeiro de 2018 9:45, Edson Lira via R-br > <r-br@listas.c3sl.ufpr.br> escreveu: > > > Bom dia caros amigos, tenho um pequeno problema que ainda não consegui > encontrar a solução. > > É o seguinte: Estou aplicando um questionário para avaliar a satisfação do > cliente em um hospital público. > > Na avaliação dos clientes estou usando uma escala de likert com as respostas > para algumas questões. > > ótimo, bom, regular, ruim, péssimo e não utilizou o serviço > > Gostaria de fazer uma mudança nessas respostas usando uma rotina do R. > > Gostaria de mudar da seguinte forma para todo o banco de uma única vez. > > otimo e bom receberia o label "positiva" regular e não utilizou o serviço o > label " " > e ruim e péssimo receberia o label "negativa". > > Alguém tem alguma sugestão? > > Bom diae muito obrigado!!! > Prof. Edson Lira, Me > Estatístico > Manaus-Amazonas > __ _ > R-br mailing list > R-br@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. > > > __ _ > R-br mailing list > R-br@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. __ _ R-br mailing list R-br@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 fornea cdigo mnimo reproduzvel. __ _ R-br mailing list R-br@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. __ _ R-br mailing list R-br@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. -- = Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodeso...@gmail.com Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ == ___ R-br mailing list R-br@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.
Re: [R-br] Recode rotina
Obrigado Thiago e Leonard, vou testar a sugestão de vcs. [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em quarta-feira, 31 de janeiro de 2018 13:12:19 GMT-4, Leonard Assis via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Eu, particularmente, usaria o recode, mas jogaria em uma outra variável Leonard Em 31 de jan de 2018 1:48 PM, "Tiago Fragoso via R-br" <r-br@listas.c3sl.ufpr.br> escreveu: Olá Você pode usar a função 'mapvalues' do pacote plyr também. No seu caso, ficaria ads <- plyr::mapvalues(ads, from = c('péssimo','ruim','otimo',' bom','regular'), to = c('Negativa','Negativa',' Positiva','Positiva','')) 2018-01-31 13:05 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: > Já resolvi, com a rotina abaixo. > > ads <- apply(ads, 2, function(x) {x <- recode(x,"'pessimo'='Negativa' ; > 'ruim'='Negativa'; > 'otimo'='Positiva';'bom'=' Positiva';'nao > utilizou o servico'='';'regular'=''"); x}) > > Se alguém tiver uma solução mais elegante, eu gostaria de saber. > > Bom dia a todos! > > [ ]'s. > Prof. Edson Lira, Me > Estatístico > Manaus-Amazonas > > > Em Quarta-feira, 31 de Janeiro de 2018 9:45, Edson Lira via R-br > <r-br@listas.c3sl.ufpr.br> escreveu: > > > Bom dia caros amigos, tenho um pequeno problema que ainda não consegui > encontrar a solução. > > É o seguinte: Estou aplicando um questionário para avaliar a satisfação do > cliente em um hospital público. > > Na avaliação dos clientes estou usando uma escala de likert com as respostas > para algumas questões. > > ótimo, bom, regular, ruim, péssimo e não utilizou o serviço > > Gostaria de fazer uma mudança nessas respostas usando uma rotina do R. > > Gostaria de mudar da seguinte forma para todo o banco de uma única vez. > > otimo e bom receberia o label "positiva" regular e não utilizou o serviço o > label " " > e ruim e péssimo receberia o label "negativa". > > Alguém tem alguma sugestão? > > Bom diae muito obrigado!!! > Prof. Edson Lira, Me > Estatístico > Manaus-Amazonas > __ _ > R-br mailing list > R-br@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. > > > __ _ > R-br mailing list > R-br@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. __ _ R-br mailing list R-br@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 fornea cdigo mnimo reproduzvel. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Recode rotina
Já resolvi, com a rotina abaixo. ads <- apply(ads, 2, function(x) {x <- recode(x,"'pessimo'='Negativa'; 'ruim'='Negativa'; 'otimo'='Positiva';'bom'='Positiva';'nao utilizou o servico'='';'regular'=''"); x}) Se alguém tiver uma solução mais elegante, eu gostaria de saber. Bom dia a todos! [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Quarta-feira, 31 de Janeiro de 2018 9:45, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Bom dia caros amigos, tenho um pequeno problema que ainda não consegui encontrar a solução. É o seguinte: Estou aplicando um questionário para avaliar a satisfação do cliente em um hospital público. Na avaliação dos clientes estou usando uma escala de likert com as respostas para algumas questões. ótimo, bom, regular, ruim, péssimo e não utilizou o serviço Gostaria de fazer uma mudança nessas respostas usando uma rotina do R. Gostaria de mudar da seguinte forma para todo o banco de uma única vez. otimo e bom receberia o label "positiva" regular e não utilizou o serviço o label " "e ruim e péssimo receberia o label "negativa". Alguém tem alguma sugestão? Bom diae muito obrigado!!! Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Recode rotina
Bom dia caros amigos, tenho um pequeno problema que ainda não consegui encontrar a solução. É o seguinte: Estou aplicando um questionário para avaliar a satisfação do cliente em um hospital público. Na avaliação dos clientes estou usando uma escala de likert com as respostas para algumas questões. ótimo, bom, regular, ruim, péssimo e não utilizou o serviço Gostaria de fazer uma mudança nessas respostas usando uma rotina do R. Gostaria de mudar da seguinte forma para todo o banco de uma única vez. otimo e bom receberia o label "positiva" regular e não utilizou o serviço o label " "e ruim e péssimo receberia o label "negativa". Alguém tem alguma sugestão? Bom diae muito obrigado!!! Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Problema ao gerar um mapa com dados estatísticos
Cesar, esta ai o problema, no arquivo que estou usando o merge a variável comum está com os nomes todos acentuados corretameente, por exemplo no arquivo da shape ele aparede assim ALVAR├âES, no aruqivo que estou usando aparede Alvarâes. Usei a rotina library(descr) am$NM_MUNICIP <- toUTF8(am$NM_MUNICIP, "IBM850") para tentar corrigir a acentuação quando uso o summary para ver depois da rotina aparece como segue abaixo. Por isso o merge dá problema. ALVAR├âES AMATUR├ü ANAM├â ANORI 1 1 1 1 APU├ì ATALAIA DO NORTE AUTAZES BARCELOS 1 1 1 1 BARREIRINHA BENJAMIN CONSTANT BERURI BOA VISTA DO RAMOS 1 1 1 1 BOCA DO ACRE BORBA CAAPIRANGA CANUTAMA 1 1 1 1 CARAUARI CAREIRO CAREIRO DA V├üRZEA COARI 1 1 1 1 CODAJ├üS EIRUNEP├ë ENVIRA FONTE BOA 1 1 1 1 GUAJAR├ü HUMAIT├ü IPIXUNA IRANDUBA 1 1 1 1 Corrigi meu banco de dados de acordo com o arquivo da shape e resolvi o problema. Tenho outras dúvidas, mas vamos descansar no final de semana. Bom final de semana a todos e muito obrigado pela dica Cesar, foi a partir dai que entendi o que estava acontecendo. [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Sábado, 13 de Janeiro de 2018 17:04, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Cesar, retorna NA, esse problema aparece depois do merge, estou desconfiando que é a relação dos municípios, que diverge na hora do merge. [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Sexta-feira, 12 de Janeiro de 2018 19:03, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Edson, O que a chamada: > summary(am$idh2000) retorna? 2018-01-12 13:25 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: Caros amigos, abaixo os arquivos (links) que estou usando para gerar um mapa do estado do Amazonas com o IDH Municipal de 2000 (está 2005, mas é 2000). Quando rodo a rotina para fazer a categorização do IDH (cut) me dá a mensagem: > am$idhc2000 <- cut(am$idh2000, 4) Error in seq.int(rx[1L], rx[2L], length.out = nb) : 'from' must be a finite number Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf > Não consegui descobrir o que está errado. Bom dia!!! Download amazonas idh 2005.R | | | Download amazonas idh 2005.R | | | Download amazonas calha.csv | | | Download amazonas calha.csv | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.cpg | | | Download 13MUE250GC SIR.cpg | | | https://www.datafilehost.com/ d/27e20ca4 # SHAPE [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Problema ao gerar um mapa com dados estatísticos
Cesar, retorna NA, esse problema aparece depois do merge, estou desconfiando que é a relação dos municípios, que diverge na hora do merge. [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Sexta-feira, 12 de Janeiro de 2018 19:03, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Edson, O que a chamada: > summary(am$idh2000) retorna? 2018-01-12 13:25 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: Caros amigos, abaixo os arquivos (links) que estou usando para gerar um mapa do estado do Amazonas com o IDH Municipal de 2000 (está 2005, mas é 2000). Quando rodo a rotina para fazer a categorização do IDH (cut) me dá a mensagem: > am$idhc2000 <- cut(am$idh2000, 4) Error in seq.int(rx[1L], rx[2L], length.out = nb) : 'from' must be a finite number Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf > Não consegui descobrir o que está errado. Bom dia!!! Download amazonas idh 2005.R | | | Download amazonas idh 2005.R | | | Download amazonas calha.csv | | | Download amazonas calha.csv | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.cpg | | | Download 13MUE250GC SIR.cpg | | | https://www.datafilehost.com/ d/27e20ca4 # SHAPE [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Problema ao gerar um mapa com dados estatísticos
Caros amigos, abaixo os arquivos (links) que estou usando para gerar um mapa do estado do Amazonas com o IDH Municipal de 2000 (está 2005, mas é 2000). Quando rodo a rotina para fazer a categorização do IDH (cut) me dá a mensagem: > am$idhc2000 <- cut(am$idh2000, 4) Error in seq.int(rx[1L], rx[2L], length.out = nb) : 'from' must be a finite number Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf > Não consegui descobrir o que está errado. Bom dia!!! Download amazonas idh 2005.R | | | Download amazonas idh 2005.R | | | Download amazonas calha.csv | | | Download amazonas calha.csv | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.shx | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.dbf | | | Download 13MUE250GC SIR.cpg | | | Download 13MUE250GC SIR.cpg | | | https://www.datafilehost.com/d/27e20ca4 # SHAPE [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Reshape
Valeu Delcio, muito obrigado!!! [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Quinta-feira, 16 de Novembro de 2017 12:13, Delcio Rudinei Bortolanza via R-br <r-br@listas.c3sl.ufpr.br> escreveu: #yiv6713012464 #yiv6713012464 -- P {margin-top:0;margin-bottom:0;}#yiv6713012464 dados<-read.table(header = F, stringsAsFactors = F, text = "n id UF x2002 x2003 x2004 x2005 x2006 x2007 x2008 x2009 x2011 x2012 x2013 x20141 1 Acre 192 189 311 311 333 327 338 354 357 381 337 3732 2 Alagoas 1251 1269 1292 1357 1413 1410 1355 1401 1344 1391 1429 14173 3 Amapá 204 198 239 262 252 294 317 299 302 311 332 3554 4 Amazonas 958 1050 1433 1528 1502 1488 1531 1673 1689 1720 1771 19085 5 Bahia 6885 6901 7167 7336 7292 7338 7665 7792 7602 7416 7530 78796 6 Ceará 3705 3871 3944 4120 4143 4138 4314 4380 4121 4184 4215 4319 ")dadosdados1<-t(dados)dados1 dados2<-dados1[c(-1,-2),]dados2 dados3<-as.data.frame(dados2)rownames(dados3)<-NULLcolnames(dados3)<-NULLdados3 De: R-br <r-br-boun...@listas.c3sl.ufpr.br> em nome de Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> Enviado: quinta-feira, 16 de novembro de 2017 11:16 Para: a lista Brasileira oficial de discussão do programa R. Assunto: [R-br] Reshape Bom dia caros amigos, estou trabalhando com a base de dados abaixo(somente 6 linhas): id UF x2002 x2003 x2004 x2005 x2006 x2007 x2008 x2009 x2011 x2012 x2013 x2014 1 1 Acre 192 189 311 311 333 327 338 354 357 381 337 373 2 2 Alagoas 1251 1269 1292 1357 1413 1410 1355 1401 1344 1391 1429 1417 3 3 Amapá 204 198 239 262 252 294 317 299 302 311 332 355 4 4 Amazonas 958 1050 1433 1528 1502 1488 1531 1673 1689 1720 1771 1908 5 5 Bahia 6885 6901 7167 7336 7292 7338 7665 7792 7602 7416 7530 7879 6 6 Ceará 3705 3871 3944 4120 4143 4138 4314 4380 4121 4184 4215 4319 Estou usando o reshape para tentar transformar as linhas em colunas, ou seja, cada uma UF seria uma coluna, e assim com as demaisComo gostaria que ficasse: Acre Alagoas Amapá Amazonas Bahia Ceará ano 192 1251 3705 2002 189 1269 3871 2003 311 1293 3944 2004 ... Estou usando a rotina: pea<-reshape(pea1, varying=c("x2002","x2003","x2004","x2005","x2006","x2007", "x2008","x2009","x2011","x2012","x2013","x2014"), v.names="Medida", timevar="UF", times=c("Acre","Alagoas","Amapá","Amazonas","Bahia", "Ceará","Distrito Federal","Espírito Santo", "Goiás","Maranhão","Mato Grosso","Mato Grosso do Sul", "Minas Gerais","Pará","Paraíba","Paraná"," Pernambuco", "Piauí","Rio de Janeiro","Rio Grande do Norte", "Rio Grande do Sul","Rondônia","Roraima", "Santa Catarina","São Paulo","Sergipe","Tocantins"), new.row.names=1:27, direction="wide") que está me dando o erro: Erro em varying[, i] : número incorreto de dimensões Alguém tem alguma sugesstão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Reshape
Bom dia caros amigos, estou trabalhando com a base de dados abaixo(somente 6 linhas): id UF x2002 x2003 x2004 x2005 x2006 x2007 x2008 x2009 x2011 x2012 x2013 x2014 1 1 Acre 192 189 311 311 333 327 338 354 357 381 337 373 2 2 Alagoas 1251 1269 1292 1357 1413 1410 1355 1401 1344 1391 1429 1417 3 3 Amapá 204 198 239 262 252 294 317 299 302 311 332 355 4 4 Amazonas 958 1050 1433 1528 1502 1488 1531 1673 1689 1720 1771 1908 5 5 Bahia 6885 6901 7167 7336 7292 7338 7665 7792 7602 7416 7530 7879 6 6 Ceará 3705 3871 3944 4120 4143 4138 4314 4380 4121 4184 4215 4319 Estou usando o reshape para tentar transformar as linhas em colunas, ou seja, cada uma UF seria uma coluna, e assim com as demaisComo gostaria que ficasse: Acre Alagoas Amapá Amazonas Bahia Ceará ano 192 1251 3705 2002 189 1269 3871 2003 311 1293 3944 2004 ... Estou usando a rotina: pea<-reshape(pea1, varying=c("x2002","x2003","x2004","x2005","x2006","x2007", "x2008","x2009","x2011","x2012","x2013","x2014"), v.names="Medida", timevar="UF", times=c("Acre","Alagoas","Amapá","Amazonas","Bahia", "Ceará","Distrito Federal","Espírito Santo", "Goiás","Maranhão","Mato Grosso","Mato Grosso do Sul", "Minas Gerais","Pará","Paraíba","Paraná"," Pernambuco", "Piauí","Rio de Janeiro","Rio Grande do Norte", "Rio Grande do Sul","Rondônia","Roraima", "Santa Catarina","São Paulo","Sergipe","Tocantins"), new.row.names=1:27, direction="wide") que está me dando o erro: Erro em varying[, i] : número incorreto de dimensões Alguém tem alguma sugesstão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Corrigir ausência de dados
Valeu Cesar e Elias Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Terça-feira, 14 de Novembro de 2017 18:27, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Edson, Considerando a última msg, sobre o problema em geral: Não existe uma única fórmula para resolver esse problema. No jargão atual a tentativa de resolvê-lo, além do óbvio de deletar os casos com as variáveis faltantes ou considerar uma a uma, há, a assim denominada, imputação. A respeito disso há desde tutoriais simples até obras (com custo até salgado) inteiras dedicadas a esse assunto. Como sempre, existem pacotes do R para isso, que consideram o usuário já sapiente sobre as técnicas que automatiza. HTH--Cesar Rabak 2017-11-14 12:30 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: Caros amigos, tenho uma série de dados de pea e pnea que têm uma ausência de dados no ano de 2010. Como faço para corrigir essa ausência? TrêsDuas dúvidas: 1. Como organizar a planilha visto que devo corrigir a ausência por estado para corrigir a ausência de pea e pnea para o ano de 2010 por estado?2. Qual a técnica que me aconselham a usar?3. E como fazê-lo com uma única rotina? Abaixo o link da planilha Download pea pnea 02 14.csv | | | Download pea pnea 02 14.csv | | | [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Corrigir ausência de dados
Obrigado Elias, mas supondo que houvesse esse problema, como resolvê-lo? Pensei em um modelo de regressão, como terei que fazer por estrado, como fazer em única rotina? Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Terça-feira, 14 de Novembro de 2017 16:29, Elias T Krainski via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Acho que você precisa procurar por dados do censo pois 2010 é ano de censo e não está incluido em séries de PNADs. Aliás, essas estatísticas são ainda mais precisas para esse ano. Elias On 14/11/17 23:14, Edson Lira via R-br wrote: Caros amigos, tenho uma série de dados de pea e pnea que têm uma ausência de dados no ano de 2010. Como faço para corrigir essa ausência? ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Corrigir ausência de dados
Ninguém tem alguma idéia? [ ['s Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Terça-feira, 14 de Novembro de 2017 10:31, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Caros amigos, tenho uma série de dados de pea e pnea que têm uma ausência de dados no ano de 2010. Como faço para corrigir essa ausência? TrêsDuas dúvidas: 1. Como organizar a planilha visto que devo corrigir a ausência por estado para corrigir a ausência de pea e pnea para o ano de 2010 por estado?2. Qual a técnica que me aconselham a usar?3. E como fazê-lo com uma única rotina? Abaixo o link da planilha Download pea pnea 02 14.csv | | | Download pea pnea 02 14.csv | | | [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Corrigir ausência de dados
Caros amigos, tenho uma série de dados de pea e pnea que têm uma ausência de dados no ano de 2010. Como faço para corrigir essa ausência? TrêsDuas dúvidas: 1. Como organizar a planilha visto que devo corrigir a ausência por estado para corrigir a ausência de pea e pnea para o ano de 2010 por estado?2. Qual a técnica que me aconselham a usar?3. E como fazê-lo com uma única rotina? Abaixo o link da planilha Download pea pnea 02 14.csv | | | Download pea pnea 02 14.csv | | | [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] Assunto: Re: ajuda
Antes vc deve verificar qual a estrutura da variável dt_nascstr(dt_nasc) Quando o banco está em Excell, geralmente ela não vem como character, vc deve converter a variável em character, depois transforma-la em dar, aí fazer sua operação como foi sugerido anteriormente. Enviado do Yahoo Mail no Android Em ter, 24 24e out 24e 2017 às 19:32, Fernando Souza via R-brr-br@listas.c3sl.ufpr.br> escreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Anova para amostras de tamanhos diferentes
Dá uma olha nesse material Posts sobre desbalanceado em Ridículas | | | Posts sobre desbalanceado em Ridículas Posts sobre desbalanceado escritos por Walmes Zeviani | | | Bom dia a todos! Prof. Edson Lira, Me Estatístico Manaus-Amazonas Em Terça-feira, 4 de Julho de 2017 6:52, Fernando Souza via R-br <r-br@listas.c3sl.ufpr.br> escreveu: install.packages(car)library(car)Anova(modelo, type="lll") Em 03/07/2017 23:20, "Fabiana Goncalves Barbosa via R-br" <r-br@listas.c3sl.ufpr.br> escreveu: Olá Pessoal Como faço uma ANOVA para amostras de tamanhos diferentes? __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Assunto: Criar novo dataframe a partir de uma estrutura
Não seria o caso de subset? Enviado do Yahoo Mail no Android Em seg, 19 19e jun 19e 2017 às 18:21, Elias Carvalho via R-brr-br@listas.c3sl.ufpr.br> escreveu: Boa tarde pessoal Tenho um dataframe com a seguinte esturutra: data.frame':1999 obs. of 4 variables: $ RAC: Factor w/ 5 levels "amarela","branca",..: 1 1 1 1 1 1 2 2 2 2 ... $ ECV: Factor w/ 7 levels "casada","desquitada",..: 1 1 6 6 7 1 6 6 6 6 ... $ DCF: Factor w/ 3 levels "não","não sabe",..: NA NA 3 NA NA NA NA NA 3 3 ... $ ERG: Factor w/ 3 levels "ALTO RISCO","BAIXO RISCO",..: 2 1 3 1 1 1 2 2 2 2 ... Gostaria de criar um novo dataframe com o mesmo numero de linhas(1999) mesma quantidade de variáveis(4) e mesmo número de niveis dessas variáveis tipo fator. Achei o comando para criar a estrutura: df2 <- data.frame(matrix("", ncol = ncol(df1), nrow = nrow(df1))) mas não achei como criar com o mesmo tipo de variável e mesma quantidade de níveis, no caso factor -- In Jesu et Maria ObrigadoProf. Elias Carvalho "Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Assunto: Assunto: Re: Ajuda: Append
Ler no mesmo CRM Enviado do Yahoo Mail no Android Em ter, 4 4e abr 4e 2017 às 18:01, Edson Lira via R-brr-br@listas.c3sl.ufpr.br> escreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Assunto: Re: Ajuda: Append
Você pode também ler os dois bancos no mesmo vem e criar dois novos data.frame (df1,df2) com as variáveis de interesse e fazer um rbind. df1=data.frame(lazer=bv1$lazer, trabalho=bv1$trab) df2=data.frame(lazer=bv2$lazer, trabalho=bv2$trab) dfn=rbind(df1,df2) Acho que assim vai resolver, faça o teste. Caso contrário, deve usar o merge, que é mais apropriado. Enviado do Yahoo Mail no Android Em ter, 4 4e abr 4e 2017 às 17:01, FHRB Toledo via R-brr-br@listas.c3sl.ufpr.br> escreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Assunto: Assunto: [R-Br] Alguém sabe sobre o suporte a vetores longos (> 2 GB)?
ff Enviado do Yahoo Mail no Android Em dom, 2 2e abr 2e 2017 às 21:26, Edson Lira via R-brr-br@listas.c3sl.ufpr.br> escreveu: ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Assunto: [R-Br] Alguém sabe sobre o suporte a vetores longos (> 2 GB)?
Baixa o pacote FF, vê se ele resolve. Enviado do Yahoo Mail no Android Em sáb, 1 1e abr 1e 2017 às 17:39, Maurício Camargo via R-brr-br@listas.c3sl.ufpr.br> escreveu: Olá a todos. Preciso alocar um vetor com mais de 2 GB. Alguém sabe como anda o suporte a LongVectors do R? Obrigado. Maurício Camargo___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Fóruns e grupos
Procure no Google"Foruns de Estatística usando o software R" Boa noite. Edson Lira Estatístico Manaus-Amazonas Em Sábado, 11 de Fevereiro de 2017 12:10, Roberto Guevara via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Colegas usuários de R, estou em busca de fóruns na internet ou grupos de whatsapp sobre Matemática e/ou Estatística (incluindo sobre R). Poderiam me indicar algum? Obrigado -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600 ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Trabalhando com datas
Muito obrigado Cesar e Ivan. A todos os listeiros um bom dia!!! Edson Lira Estatístico Manaus-Amazonas Em Segunda-feira, 6 de Fevereiro de 2017 19:02, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Edson, Não sabendo exatamente para quê você necessita a diferença, talvez lhe convenha ter este aditamento: a diferença após a conversão das cadeias de caracteres com strptime() pode ser feita diretamente:> hht[2]- hht[1]Time difference of 30 mins se você quiser apenas o valor numérico, pode-se obtê-lo por:> as.double(hht[2]- hht[1]) 30se as unidades forem importantes:> units(hht[2]- hht[1])[1] "mins" Ademais, se você quiser "normalizar" as unidades, para evitar o seguinte: hh <- c("10:00","11:30")hht <- strptime(hh, format = '%H:%M', 'GMT')as.double(hht[2] - hht[1])1.5units((hht[2] - hht[1])"hours" você pode fazeras.numeric(hht[2] - hht[1], units="mins")90 HTH--Cesar Rabak 2017-02-06 15:44 GMT-02:00 Luis Iván Ortiz Valencia via R-br <r-br@listas.c3sl.ufpr.br>: hh <- c("09:20","09:50")hht <- strptime(hh, format = '%H:%M', 'GMT')difftime(hht[2],hht[1])Time difference of 30 mins 2017-02-06 12:00 GMT-02:00 Edson Lira via R-br <r-br@listas.c3sl.ufpr.br>: Bom dia listeiros, tenho esse banco: | dt_sol | dt_atend | hora_ent | hora_sai | | 01/02/2016 | 02/02/2016 | 09:20 | 09:50 | | 01/02/2016 | 03/02/2016 | 09:00 | 10:44 | | 02/02/2016 | 02/02/2016 | 09:30 | 11:00 | | 02/02/2016 | 05/02/2016 | 10:00 | 12:00 | | 02/02/2016 | 12/02/2016 | 10:30 | 12:00 | | 02/02/2016 | 16/02/2016 | 10:00 | 11:30 | | 02/02/2016 | 18/02/2016 | 09:30 | 10:16 | As variáveis hora_ent e hora_saida como fator. Gostaria de calcular a diferença entre o tempo da entrada para a saída Edson Lira Estatístico Manaus-Amazonas __ _ R-br mailing list R-br@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-g uia) e forneça código mínimo reproduzível. -- Luis Iván Ortiz Valencia Doutor Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc. __ _ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] boxplot com outliers, como colocar os valores?
algumas formas pie(frota, main="Frota 2009 - Niterói_RJ", init.angle=180)text(locator(length(names(frota))),rotulos) text(450,100,"reta de regressão") Veja esse link abaixohow to add text in my boxplot? | | | | || | | | | | how to add text in my boxplot? I have plotted box plot for my data using given command and trying to write text (p-value ) in it. Using : b... | | | | Edson Lira Estatístico Manaus-Amazonas Em Segunda-feira, 6 de Fevereiro de 2017 10:15, Christina Grupioni via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Olá listeiros! Estou precisando incluir o valor dos outliers no boxplot que estou plotando junto com o histograma. E também necessito incluir o valor dos quartis. É possível? O código que estou usando é o seguinte. mat <- matrix(c(1,2,0,0), 2)matlayout(mat, c(3.5,1), c(1,3))par(mar=c(0.5, 4.5, 0.5, 0.5))boxplot(campo$ef, horizontal=TRUE, axes=FALSE, col=(cores))par(mar=c(4.5, 4.5, 0.5, 0.5))par(cex.axis=0.8,cex.main=1)hist(campo$ef,col=(cores),main="Histograma da Eficiência de colheita",xlab=("Eficiência (%)"),ylab="Frequencia",adj=0,breaks=12,border="white", labels=qcap)axis(1,at=seq(50,100,5))axis(2,at=seq(0,22,2))density.default(campo$ef)d<-density(campo$ef)par(new=TRUE)plot(d,ann=FALSE, axes=FALSE)-- Abaixo aqui meus dados. Não sei se está correto, sou nova mesmo no uso e me indicaram usar o comando dput para incluir os dados no post: dput(campo)structure(list(trat = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L), cacho = c(34L, 35L, 15L, 33L, 24L, 28L, 25L, 29L, 14L, 8L, 5L, 13L, 20L, 31L, 32L, 27L, 9L, 10L, 11L, 12L, 26L, 23L, 22L, 21L, 6L, 4L, 7L, 2L, 19L, 18L, 17L, 16L), carret = structure(c(1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L), .Label = c("nao", "sim"), class = "factor"), local = structure(c(1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L), .Label = c("MG", "MS"), class = "factor"), est = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("maduro", "vez"), class = "factor"), A = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("continuo", "repique" ), class = "factor"), H = c(4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L), tempo = c(87L, 127L, 35L, 100L, 59L, 27L, 30L, 64L, 44L, 31L, 52L, 75L, 40L, 63L, 62L, 32L, 43L, 35L, 21L, 57L, 19L, 112L, 113L, 73L, 27L, 24L, 29L, 51L, 41L, 50L, 39L, 88L), frutos.bons = c(114L, 138L, 170L, 222L, 83L, 81L, 125L, 134L, 76L, 155L, 240L, 189L, 55L, 295L, 281L, 369L, 74L, 67L, 126L, 214L, 62L, 10L, 32L, 185L, 83L, 83L, 77L, 208L, 42L, 54L, 41L, 50L), Kg.frutos.bons = c(5.9, 7.6, 4.7, 10.8, 4.1, 5.2, 7.8, 5.6, 5.5, 8, 9.1, 4.8, 2.9, 11, 13.2, 15.5, 4, 2.1, 5.1, 5.1, 2.7, 2.4, 2.5, 8.4, 3.2, 4, 4.7, 9.9, 1.9, 2.5, 4.4, 5.7), colhidos = c(146L, 187L, 189L, 288L, 110L, 113L, 128L, 180L, 148L, 162L, 249L, 246L, 80L, 349L, 397L, 422L, 82L, 104L, 131L, 240L, 66L, 88L, 96L, 235L, 84L, 91L, 97L, 258L, 50L, 109L, 204L, 250L), dleves = c(6L, 30L, 12L, 37L, 14L, 24L, 3L, 24L, 51L, 5L, 5L, 50L, 22L, 41L, 86L, 43L, 2L, 29L, 5L, 13L, 4L, 44L, 23L, 40L, 1L, 6L, 18L, 20L, 6L, 39L, 101L, 48L), dgraves = c(26L, 19L, 7L, 29L, 13L, 8L, 0L, 22L, 21L, 2L, 4L, 7L, 3L, 13L, 30L, 10L, 6L, 8L, 0L, 13L, 0L, 34L, 41L, 10L, 0L, 2L, 2L, 30L, 2L, 16L, 62L, 152L), naocolhidos = c(7L, 4L, 17L, 8L, 2L, 9L, 0L, 16L, 10L, 3L, 8L, 23L, 2L, 8L, 4L, 1L, 5L, 4L, 1L, 9L, 0L, 19L, 87L, 6L, 1L, 1L, 2L, 21L, 6L, 13L, 10L, 19L), carga = c(153L, 191L, 206L, 296L, 112L, 122L, 128L, 196L, 158L, 165L, 257L, 269L, 82L, 357L, 401L, 423L, 87L, 108L, 132L, 249L, 66L, 107L, 183L, 241L, 85L, 92L, 99L, 279L, 56L, 122L, 214L, 269L ), tamanho = c(0.76, 0.84, 0.77, 0.97, 0.64, 0.71, 1.04, 0.89, 0.94, 0.88, 0.86, 0.82, 0.81, 1, 0.95, 1.09, 0.83, 0.74, 0.84, 0.63, 1.02, 0.54, 0.87, 0.86, 0.87, 0.79, 0.82, 0.78, 0.74, 0.69, 0.57, 0.66), altura = c(2.4, 2.6, 3, 2.67, 2.26, 2.36, 2.45, 2.98, 2.97, 2.39, 3.13, 2.28, 1.76, 2.66, 2.51, 3.8, 3.52, 2.31, 2.86, 2.97, 2.39, 3.34, 3.42, 2.38, 3.33, 3.02, 2.47, 3.46, 1.58, 3.56, 3.4, 3.71), cap = c(244.137931, 215.4330709, 483.4285714, 388.8, 250.1694
[R-br] Trabalhando com datas
Bom dia listeiros, tenho esse banco: | dt_sol | dt_atend | hora_ent | hora_sai | | 01/02/2016 | 02/02/2016 | 09:20 | 09:50 | | 01/02/2016 | 03/02/2016 | 09:00 | 10:44 | | 02/02/2016 | 02/02/2016 | 09:30 | 11:00 | | 02/02/2016 | 05/02/2016 | 10:00 | 12:00 | | 02/02/2016 | 12/02/2016 | 10:30 | 12:00 | | 02/02/2016 | 16/02/2016 | 10:00 | 11:30 | | 02/02/2016 | 18/02/2016 | 09:30 | 10:16 | As variáveis hora_ent e hora_saida como fator. Gostaria de calcular a diferença entre o tempo da entrada para a saída Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Script Odds Ratio para o multinom ou vglm
Veja se ajuda | | | CodigosR | | | CodigosR Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 13 de Janeiro de 2017 9:07, Walmes Zeviani via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Prezado, Ajude àqueles dispostos a te ajudar. Forneça CMR. À disposição. Walmes. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Avaliação da evolução do peso
Eu gostaria de fazer que a estrutura do banco ficasse assim: id sexo dt_nasc dt1_acom dt2acom dt3_acom peso1 peso2 peso3... [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 10:49, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Sim, o problema é que ja me trouxeram o banco com essa estrutura, a idéia e avaliar o paciente em todos os seus comparecimentos, evidentemente em datas diferentes. Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 10:23, Elias T. Krainski via R-br <r-br@listas.c3sl.ufpr.br> escreveu: On 08/30/2016 10:40 AM, Edson Lira via R-br wrote: Como faço para avaliar seu peso a cada vez que comparece? considerando com o peso inicial e com os pesos subsequentes registrados. Se um paciente comparece três vezes, seu peso será comparado três vezes (a segunda com a primeira e a terceira com a segunda e a terceira com a primeira). Análise longitudinal? Elias ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Avaliação da evolução do peso
Sim, o problema é que ja me trouxeram o banco com essa estrutura, a idéia e avaliar o paciente em todos os seus comparecimentos, evidentemente em datas diferentes. Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 10:23, Elias T. Krainski via R-br <r-br@listas.c3sl.ufpr.br> escreveu: On 08/30/2016 10:40 AM, Edson Lira via R-br wrote: Como faço para avaliar seu peso a cada vez que comparece? considerando com o peso inicial e com os pesos subsequentes registrados. Se um paciente comparece três vezes, seu peso será comparado três vezes (a segunda com a primeira e a terceira com a segunda e a terceira com a primeira). Análise longitudinal? Elias ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Visualizar grande data.frame
Se vc selecionar as variáveis de interesse, fica mais fácil de usar. Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 10:00, Abel Brasil Ramos da Silva <abelbrasi...@gmail.com> escreveu: Olá, obrigado pela ajuda! Conheço os comandos head e tail. A intenção era ver todo o banco de dados mesmo, por exemplo, numa janela fora do R (se for possível!). Valeu. Em 30 de agosto de 2016 10:56, Edson Lira <edinhoes...@yahoo.com.br> escreveu: Ou use: bd[1:100, c(1,2,3,7)] Vc vai visualizar as 100 primeiras linhas com as variáveis 1,2,3 e 7 Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 9:41, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Baixa o pacote ff e tentar com ele. Transforma antes para csv. Geralmente faço isso. ´[ ]'1s Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 9:33, Abel Brasil Ramos da Silva via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Olá a Todos, Bom dia! Possuo um data.frame com 1 observações e 1350 colunas. Gostaria de apenas visualizar esse dados no próprio R (ou R-Studio), que pacote e comando para executar essa tarefa. Com o comando View existe uma limitação da visualização. Agradeço a ajuda desde já. Abs. __ _ R-br mailing list R-br@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� c�igo m�imo reproduz�el. __ _ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Visualizar grande data.frame
Ou use: bd[1:100, c(1,2,3,7)] Vc vai visualizar as 100 primeiras linhas com as variáveis 1,2,3 e 7 Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 9:41, Edson Lira via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Baixa o pacote ff e tentar com ele. Transforma antes para csv. Geralmente faço isso. ´[ ]'1s Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 30 de Agosto de 2016 9:33, Abel Brasil Ramos da Silva via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Olá a Todos, Bom dia! Possuo um data.frame com 1 observações e 1350 colunas. Gostaria de apenas visualizar esse dados no próprio R (ou R-Studio), que pacote e comando para executar essa tarefa. Com o comando View existe uma limitação da visualização. Agradeço a ajuda desde já. Abs. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
[R-br] Avaliação da evolução do peso
Bom dia a todos da Lista! Abaixo o link de um arquivo teste que estou tentando fazer uma análise.https://www.datafilehost.com/d/UgWFEhooCR Os dados para análise são registros de comparecimento de pacientes em uma unidade de saúde. Cada vez que o paciente comparece na unidade uma das variáveis que é registrada é seu peso, portanto o mesmo paciente pode aparecer uma, duas ou mais vezes, cada uma aparição é um registro. Como faço para avaliar seu peso a cada vez que comparece? considerando com o peso inicial e com os pesos subsequentes registrados. Se um paciente comparece três vezes, seu peso será comparado três vezes (a segunda com a primeira e a terceira com a segunda e a terceira com a primeira). Obrigado a todos![ ]'s. Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Superfície de Resposta
Veja se ajuda https://www.r-project.org/conferences/useR-2011/TalkSlides/Invited/Gromping-Design_of_Experiments.pdfhttps://cran.r-project.org/doc/contrib/Vikneswaran-ED_companion.pdfhttps://cran.r-project.org/web/views/ExperimentalDesign.html [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 23 de Agosto de 2016 8:39, Emerson Cotta Bodevan via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Caros colegas, bom dia. Estou estudando sobre superfície de resposta. Alguém poderia me indicar pacotes sobre esse assunto no R? Também, sendo possível, pacotes para planejamento e análises de experimentos. Agradeço qualquer sugestão e indicação. Att., Emerson ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Concordância
Bom dia a todos! Muito obrigado Leonardo, garimpando na net, achei algo relacionado a Bland e Altmann. Valeu!!! Edson Lira Estatístico Manaus-Amazonas Em Segunda-feira, 1 de Agosto de 2016 17:02, Leonardo Fontenelle via R-br <r-br@listas.c3sl.ufpr.br> escreveu: Se você quer estudar a concordância entre dois aparelhos quanto a uma medida (variável numérica), deve estar interessado no trabalho de Bland e Altmann: https://www-users.york.ac.uk/~mb55/meas/ba.htm Leonardo Ferreira Fontenelle Em Sex 22 jul. 2016, às 13:21, Edson Lira via R-br escreveu: Caros amigos, tenho um aparelho que o aparelho que sabemos estar fazendo as medidas corretas. Temos outros aparelhos que foram feitas 3 medidas para cada.Como testar a concordância entre as medidas do aparelho ouro com os outros aparelhos? Busquei o epiR, mas não consegui implementar um teste para verificar a sensibilidade, especificidades dos aparelhos testados com o ouro. Não sei se usar Kappa resolve. faixa eqouro eqtest A 56 42 43 43 50-60 B 54 41 42 42 C 53 39 41 40 D 87 66 67 70 70-95 E 81 65 63 66 F 90 67 68 67 G 210 148 145 147 115-400 H 117 92 88 90 I 141 120 117 114 J 200 157 156 155 Alguém tem alguma sugestão? Um abraço. Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
[R-br] Concordância
Caros amigos, tenho um aparelho que o aparelho que sabemos estar fazendo as medidas corretas. Temos outros aparelhos que foram feitas 3 medidas para cada.Como testar a concordância entre as medidas do aparelho ouro com os outros aparelhos? Busquei o epiR, mas não consegui implementar um teste para verificar a sensibilidade, especificidades dos aparelhos testados com o ouro. Não sei se usar Kappa resolve. faixa eqouro eqtest | | A | 56 | 42 | 43 | 43 | | 50-60 | B | 54 | 41 | 42 | 42 | | | C | 53 | 39 | 41 | 40 | | | | | | | | | | D | 87 | 66 | 67 | 70 | | 70-95 | E | 81 | 65 | 63 | 66 | | | F | 90 | 67 | 68 | 67 | | | | | | | | | | G | 210 | 148 | 145 | 147 | | 115-400 | H | 117 | 92 | 88 | 90 | | | I | 141 | 120 | 117 | 114 | | | J | 200 | 157 | 156 | 155 | Alguém tem alguma sugestão? Um abraço. Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Assunto: Re: Trabalhar com planilha de dados muito grande
Experimente o pacote: ff Olha esse material abaixo. http://ff.r-forge.r-project.org/bit_WU_Vienna2010.pdf O nome é esse mesmo. Edson Lira Estatístico Manaus-Amazonas Em Sábado, 4 de Junho de 2016 18:24, Leonardo Ferreira Fontenelle <leonar...@leonardof.med.br> escreveu: Não tenho experiência com o banco de dados do Censo, apenas com o da PNAD. Em geral, dá para você converter os dados das instruções de leitura do SAS para parâmetros a serem fornecidos para a leitura do arquivo TXT, que está em formato de largura fixa (FWF). A partir daí você tem dois caminhos possíveis. Um é jogar os dados num banco de dados como já foi discutido. Isso lhe possibilita trabalhar com mais dados do que a abordagem seguinte, que é a que eu conheço melhor. Outro caminho, para o caso de você querer trabalhar com poucas variáveis, é você converter, com a função fwf2csv() do pacote descr, o arquivo TXT num arquivo CSV. Esse arquivo CSV, por sua vez, pode ser lido pela função fread() do pacote data.table. A função fread permite ler apenas as variáveis desejadas, de forma a fazer caber na memória RAM todas as observações desejadas. Se você ainda não conhece data.table, leia a introdução oficial. Caso você opte pelo segundo caminho, pode adaptar este script para seus próprios fins: fwf2csv.SAScii <- function(fwffile, csvfile, input, verbose = getOption("verbose")) { require("SAScii") require("descr") # Usuário precisa conferir antes se a codificação de caracteres está correta... # https://github.com/ajdamico/SAScii/issues/5 dic <- within(parse.SAScii(input), { end <- cumsum(width) begin <- c(1, 1 + cumsum(width)[-length(width)]) }) fwf2csv(fwffile = fwffile, csvfile = csvfile, verbose = verbose, names = dic$varname, begin = dic$begin, end = dic$end) } Onde "input" é o arquivo ".sas" com as instruções de leitura para o SAS. Saudações, Leonardo Ferreira Fontenelle Em Sáb 4 jun. 2016, às 16:46, Edimeire Alexandra Pinto escreveu: Rodnei, no e-mail que mandei antes, so resumindo, os dados do censo estão em txt em meu computador e preciso selecionar algumas variáveis e exportar p o r, usando o pacote sqlite ou outro que souber Gente, se alguém mais souber eu agradeço a ajuda Enviado do Yahoo Mail no Android Em 12:38 Sáb, 4 de jun de PM, Rodrigo Coster <rcos...@gmail.com> escreveu: André, Na documentação do ultimo censo escolar tem instruções de leitura em R, usando um pacote especial para grandes bancos de dados. Escrevi uma postagem sobre leitura desses bancos usando SQlite (http://rcoster.blogspot.com.br/2014/02/lendo-grandes-bancos-de-dados.html). Dependendo do ano que tu quiser, eu tenho os comandos []'s 2016-06-04 8:18 GMT-03:00 Andre Oliveira <andreolso...@yahoo.com.br>: bom dia, estou precisando filtrar informação de arquivos de microdados do INEP, tem arquivos csv e txt com até 6gb a título basicamente de estudos. Sou leigo neste tema e gostaria de saber qual melhor forma de proceder. Me recomendaram usar um banco de dados neste caso (SQLite, MySQL, PostgreSQL ou Oracle), mas não sei por onde começar. Alguém tem material que julgue importante e que possa compartilhar comigo? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Moving Average, Modelos mistos
Time Series Analysis - With Applications in R | Jonathan D. Cryer | Springer | | | | | | | | | Time Series Analysis - With Applications in R | Jonathan D. Cryer | SpringerTime Series Analysis With Applications in R, Second Edition, presents an accessible approach to understanding time series models and their applications | | | | Visualizar em www.springer.com | Visualizado por Yahoo | | | | | Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 25 de Maio de 2016 16:58, Cesar Rabak <cesar.ra...@gmail.com> escreveu: Fernando, A melhor resposta à questão " Faz sentido uma ordem de X em ARMA?" é responder às seguintes questões: A ARMA é a melhor abordagem para a análise do fenômeno/experimento que estou realizando? O modelo subjacente de uma ordem X indica que processo de geração desses dados, e esse processo faz sentido em relação ao processo físico pelo qual os dados são gerados? HTH--Cesar Rabak 2016-05-25 13:32 GMT-03:00 Fernando Antonio de souza <nandodeso...@gmail.com>: Boa tarde, colegas Estou trabalhando com modelos mistos, com dados de medidas repetidas no tempo.Ao definir a estrutura da matriz de covariância, testei a estrutura autoregressiva de ordem 1 AR1(), Moving average ARMA(p=0,q=?), onde testei para ordem de 1 a 6 (acima de 6 o modelo não convergiu). Acontece que com a ordem 6 obtive o melhor modelo (segundo anova()), embora nem todas as correlações foram resolvidas, conforme o gráfico (anexo) Ainda não me sinto confortal em trabalhar com estas estruturas. Faz sentido uma ordem de 6 em ARMA? Vocês teriam algum material que me explique o uso destas estruturas? Atenciosamente === Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim) e-mail:nandodeso...@gmail.com Lattes: http://lattes.cnpq.br/6519538815038307 blog: https://producaoanimalcomr.wordpress.com/ ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Multiplots usando lattice
Não sei se vai funcionar no lattice.Use: par(mfrow=c(2,2)) Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 13 de Maio de 2016 9:33, Éder Comunello <comunello.e...@gmail.com> escreveu: Marcos, Tentou adicionar o argumento more=T nas instruções print de de g2 e g3? Éder ComunelloAgronomist (UEM), MSc in Environ. Sciences (UEM) DSc in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa)Dourados, MS, Brazil || GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 Em 12 de maio de 2016 23:14, marcos paulo <marcospaulo_agron...@hotmail.com> escreveu: Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica,usando o pacote lattice pra confecção desses plots. Porém não estou conseguindo colocar mais que dois na mesma janela, usando a função print. Segue exemplo: g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="NDVI", main='Vitrine 2014/2016', strip=strip.custom(bg="gray90")) g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="NDVI", strip=strip.custom(bg="gray90")) g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="Clorofila", strip=strip.custom(bg="gray90")) g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="Clorofila", strip=strip.custom(bg="gray90")) print(g1, split=c(1,1,1,4), more=TRUE)print(g2, split=c(1,2,1,4))print(g3, split=c(1,3,1,4))print(g4, split=c(1,4,1,4)) Cordialmente, Marcos Paulo. Técnico Agrícola em Zootecnia Bacharel em Agronomia Mestrando em Produção Vegetal - Universidade Federal de GoiásEstagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO Contato: (62) 85075783 http://lattes.cnpq.br/4322347592884852 | | Livre de vírus. www.avast.com. | ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] CRIAR NOVO DATA FRAME COM LINHAS SELECIONADAS
Acho que essa rotina está faltando uma pequena correção teste<-subset(banco, SE_TEM_MBA %in% NÃO) Faça assim: nao <-subset(banco, banco$SE_TEM_MBA %in% c("NÃO")) sim<-subset(banco, banco$SE_TEM_MBA %in% c("SIM")) [ ]'s.Edson Lira Estatístico Manaus-Amazonas Em Domingo, 17 de Abril de 2016 9:17, Daniel Tiezziescreveu: Amiko Tente isso: teste <- banco[banco$SE_TEM_MBA == “NÃO”,] daniel Daniel Tiezzi, MD, PhDProfessor AssociadoDepartamento de Ginecologia e ObstetríciaSetor de Mastologia e Oncologia GinecológicaFaculdade de Medicina de Ribeirão Preto - USPTel.: 16 3602-2488e-mail: dtie...@fmrp.usp.br On Apr 17, 2016, at 1:39 PM, Jobenil Luiz Magalhães Júnior wrote: Veja os exemplos aqui.https://stat.ethz.ch/R-manual/R-devel/library/base/html/subset.htmlAt.Jobenil Junior Enviado do meu iPad Em 16 de abr de 2016, às 21:58, Amiko Bh escreveu: Prezados, Tenho o seguinte banco (com aproximadamente 15000 linhas): | SE_TEM_MBA | MOTIVO_DO_MBA | ATO_INFRACIONAL_MBA | | NÃO | NSA | NSA | | NÃO | NSA | NSA | | SIM | NSA | NSA | | NÃO | NSA | NSA | | NÃO | NSA | NSA | | SIM | NSA | NSA | | NÃO | NSA | NSA | | NÃO | NSA | NSA | | NÃO | NSA | NSA | Gostaria de criar um novo data frame selecionando apenas a linhas com o dado NÃO na coluna SE_TEM_MBA. Usei os dois comandos abaixo e vejam o que acontece: 1) > teste<-subset(banco, SE_TEM_MBA %in% NÃO) Erro em match(x, table, nomatch = 0L) : objeto 'NÃO' não encontrado 2) teste<-subset(banco, SE_TEM_MBA %in% "SIM") Quando uso este retorna um data frame com os títulos da coluna, mas sem observação alguma. Muito obrigado. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Referenciar o R
> citation() To cite R in publications use: R Core Team (2014). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL http://www.R-project.org/. A BibTeX entry for LaTeX users is @Manual{, title = {R: A Language and Environment for Statistical Computing}, author = {{R Core Team}}, organization = {R Foundation for Statistical Computing}, address = {Vienna, Austria}, year = {2014}, url = {http://www.R-project.org/}, } We have invested a lot of time and effort in creating R, please cite it when using it for data analysis. See also‘citation("pkgname")’ for citing R packages. Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 14 de Abril de 2016 9:16, Rodrigo Campos <rod...@gmail.com> escreveu: Olá Marcelo, To cite R in publications use (para citar o R em publicações use): R Core Team (2015). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/. abr. 2016-04-14 10:03 GMT-03:00 Marcelo Menezes <biologou...@gmail.com>: Gente bom dia! eu gostaria de saber como posso referenciar o R em meu trabalho. Alguémpode me ajudar por favor? Att Marcelo Menezes AlvesGraduando em Ciências Biológicas pela UniversidadeFederal de Campina GrandeTécnico em Suporte e Manutenção em Informática peloInstituto Federal da Paraíba.Lattes: http://lattes.cnpq.br/1992059133849245 ___ R-br mailing list R-br@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. -- Rodrigo Campos www.pleuston.com ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Dados Missing
Elias, acho que vc precisa especificar melhor o que quer fazer. Olhe essa rotina. bd[bd==""]<-NA Elimina os missings. É só isso? [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Segunda-feira, 11 de Abril de 2016 10:12, Elias Carvalho <ecaca...@gmail.com> escreveu: Alguém da lista poderia me da uma dica de como tratar dados missing ? Multiple Imputation seria a técnica mais moderna ? Obrigado Elias ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
Re: [R-br] Como remover um pacote ???
remove.packages {utils} Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 18 de Março de 2016 16:12, Rodrigo Coster <rcos...@gmail.com> escreveu: A função é remove.packages() , com s no final: remove.packages("Rgooglemaps") 2016-03-18 17:06 GMT-03:00 Edmar Caldas <edmar_cal...@yahoo.com.br>: Pessoal não estou conseguindo remover um pacote. remove.package (Rgooglemaps) Erro: não foi possível encontrar a função "remove.package" Edmar ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
[R-br] Dúvidas em condicional composto
Caros listeiros, tenho essa configuração em um banco eem dpm brpm 1 Alto Alto Médio 2 Alto Alto Alto 3 Alto Alto Médio 4 Alto Alto Alto 5 Alto Alto Médio 6 Alto Alto Baixo 7 Alto Alto Baixo Se duas dessas variáveis apresentarem o nível alto, digo que há um estresse no indivíduo. Por exemplo se eem ou dpm for alto, digo que há um estresse, basta uma cominação qualquer de duas ou três variáveis apresentando alto, digo que há estresse. eem=alto e dpm= alto ->estresseeem=alto e brpm= alto ->estressedpm=alto e brpm= alto ->estressedpm=medio e brpm= alto ->sem estresse Não sei como implementar no R. Alguém poderia me sugerir alguma rotina? [ ]'s Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Dúvidas em condicional composto
É isso mesmo Rodrigo, é que se das 3, se pelo menos duas apresentarem o nível alto, afirma-se que tem estresse, caso contrário não. [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 26 de Janeiro de 2016 11:19, Rodrigo Coster <rcos...@gmail.com> escreveu: Edson, O teu exemplo não bate com o perguntado (no exemplo só aparecem 2 variáveis, nos dados/texto tem 3). Caso for: sempre que pelo menos 2 for alto, ter estresse, o comando abaixo deve servir: dados <- read.table(text="eem dpm brpmAlto Alto MédioAlto Alto AltoAlto Alto MédioAlto Alto AltoAlto Alto MédioAlto Alto BaixoAlto Alto BaixoAlto Medio Medio", sep = ' ', header = TRUE) ifelse(rowSums(dados == 'Alto') > 1, 'Estresse', 'Sem Estresse') 2016-01-26 13:02 GMT-02:00 Edson Lira <edinhoes...@yahoo.com.br>: Caros listeiros, tenho essa configuração em um banco eem dpm brpm 1 Alto Alto Médio 2 Alto Alto Alto 3 Alto Alto Médio 4 Alto Alto Alto 5 Alto Alto Médio 6 Alto Alto Baixo 7 Alto Alto Baixo Se duas dessas variáveis apresentarem o nível alto, digo que há um estresse no indivíduo. Por exemplo se eem ou dpm for alto, digo que há um estresse, basta uma cominação qualquer de duas ou três variáveis apresentando alto, digo que há estresse. eem=alto e dpm= alto ->estresseeem=alto e brpm= alto ->estressedpm=alto e brpm= alto ->estressedpm=medio e brpm= alto ->sem estresse Não sei como implementar no R. Alguém poderia me sugerir alguma rotina? [ ]'s Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Comandos
No google busque : "Teste de médias com o software R""Gráficos com o R" Boa sorte!!! Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 22 de Outubro de 2015 1:44, sergio simao joao <sergiosimaoj...@live.com> escreveu: Bom dia. Acabo de ser novo usuário do programa R, e estou tendo problema de encontrar comandos que permitem comparar graficamente a significância de dois ou mais grupos de medias.Pedir desculpa desde já porque ser uma questão tao básica. Continuação de optimo dia ___ R-br mailing list R-br@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� c�igo m�imo reproduz�el. ___ R-br mailing list R-br@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.
[R-br] Qualidade de vida
Estou tentando analisar um banco que mede qualidade de vida. Abaixo um script usado no SPSS (http://www.ufrgs.br/psiquiatria/psiq/whoqol86.html)Instrumentos de Avaliação de Qualidade de vida - OMS - versão em português | | | | | | | | | | | Instrumentos de Avaliação de Qualidade de vida - OMS - v...Versão em português dos Instrumentos de Avaliação de Qualidade de Vida - WHOQOL-100 e WHOQOL-bref - da Organização Mundial de Saúde. Projeto de... | | | | Visualizar em www.ufrgs.br | Visualizado por Yahoo | | | | | STEPS FOR CHECKING AND CLEANING DATA AND COMPUTING FACET AND DOMAIN SCORES Recode q3 q4 q26 (1=5) (2=4) (3=3) (4=2) (5=1). (scores based on a 4-20 scale) compute dom 1= (mean.6 (q3, q4, q10, q15, q16, q17, q18)) * 4. compute dom 2= (mean.5 (q5, q6, q7, q11, q19, q26)) * 4 compute dom 3= (mean.2 (q20,q21,q22)) * 4 compute dom 4= (mean.6 (q8, q9, q12, q13, q14, q23, q24, q25)) * 4 compute overall= (mean.2 (q1, q2)) * 4 (scores transformed to a 0-100 scale) compute dom1b= (dom1 - 4) * (100/16). compute dom2b= (dom2 - 4) * (100/16). compute dom3b= (dom3 - 4) * (100/16). compute dom4b= (dom4 - 4) * (100/16). compute q1b= (q1 - 1) * (100/16). compute q2b= (q2 - 1) * (100/16). Gostaria de replicar este script, alguém poderia me ajudar? O link do banco: Download qv.csv | | | | | | | | | | | Download qv.csvDownload qv.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehos... | Visualizado por Yahoo | | | | | [ ]'s Edson LiraEstatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Fatorial 2^4 com ponto central
Procura quick R no Google Enviado do Yahoo Mail no Android De:ana paula coelho madeira apcmade...@hotmail.com Data:13:03 seg, 22 de jun de PM Assunto:[R-br] Fatorial 2^4 com ponto central Prezados colegas, boa tarde. Tenho um planejamento fatorial 2^4 com um ponto central ( 4 repetições), totalizando 20 experimentos. Alguém tem sugestão de como fazer a análise no R (obs: consegui fazer no Statistica, mas queria apresentar no R)? Desde já agradeço Att., Ana Paula ___ R-br mailing list R-br@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.
Re: [R-br] Teste de Kruskal Wallis
Prezados, bom dia. Preciso fazer uma comparação entre4 tratamentos via teste de Kruskal Wallis. Utilizei o comando: kruskal.test(Resposta~trat, dados) Kruskal-Wallis rank sum test data: Resposta by trat Kruskal-Wallis chi-squared = 14.602, df = 3, p-value = 0.00219 indicando haver diferenças significativas. Como fazer comparações múltiplas entre os grupos? Desde já agradeço. Att., Ana Paula ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Pivot tables / Tabelas dinâmicas.
Meu caro, para gerar tabelas simples uso a rotina abaixo: Seja o banco: dfdf[df==]-NAtabelas-function(x){ Qtde-table(x) Perc-100*prop.table(table(x)) tab-cbind(Qtde,Perc) Total-colSums(tab) rbind(tab,Total) } apply(df,2,tabelas) Você terá todas tabelas com Qtde e perc de seu banco. Para gerar tabelas com duas ou mais entradas, uso:v12-100*prop.table(table(df$var1,df$var2))ftable(addmargins(v12, c(1,2), FUN=list(Total = sum))) [ ]'sEdson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 3 de Junho de 2015 23:00, Paulo Abreu abre...@me.com escreveu: As funções melt e dcast fazem a mágica. Em 03/06/2015, à(s) 10:47, Diego Miro d.miro1...@gmail.com escreveu: Alexandre, dependendo do que você quer fazer há algumas opções. Procure pelas funções xtable() e aggregate() No pacote reshape2 tem a função cast() também. Um abraço. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Dúvida ao juntar duas bases de dados
Diego, imaginando que é da forma abaixo seus bancos, eu faria o seguinte dados2-data.frame(ident=dados$id, varquero=var$dados)dados3-data.frame(ident=dados$id, varquero=var$dados1) var$dados=var$dados1 dfq-rbind(dados2,dados3) Veja que seu banco terá a soma da dimensão dos dois bancos. [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Sábado, 16 de Maio de 2015 13:05, Victor Eduardo victorduc...@gmail.com escreveu: Oi Diego, eu tentei isso, e ele simplesmente aumenta o número de linhas. Sobe para mais de 20 mil linhas na hora do merge Em 16 de maio de 2015 13:32, Diego Miro d.miro1...@gmail.com escreveu: Victor, tente usar o argumento (all.x = TRUE) dentro do merge. Dessa forma você estará fazendo um Left Join. Ele vai retornar a tabela dados inteira, com todos os seus registros e apenas os registros de dados1 que tenham um id pertencente a coluna id de dados. Caso não encontre o id em dados, retorna NA em dados1. merge(x = dados, y = dados1, by = 'Id', all.x = TRUE) Em 16 de maio de 2015 12:00, Fernando Antonio de souza nandodeso...@gmail.com escreveu: Veja a função merge()Em 16/05/2015 11:35, Victor Eduardo victorduc...@gmail.com escreveu: Bom dia, pessoal! Estou com um problema na hora de juntar a base de dados. Estou com duas bases com as seguintes dimensões: dim(dados)[1] 2056 11 dim(dados1)[1] 1782 2 Ambas as bases possuem uma variável chamada id. O que queria é formar um subconjunto a partir de dados1 baseado nesses ids, ou seja, aproveitaria as variáveis da base dados com as variáveis em dados1. Eu tentei com o merge e não está dando certo. Existe alguma função que faça isso? Abraços e ótima semana! ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Merge com mais de 2 bancos
Como o Walmes escreveu anteriormente, se os bancos que você quer juntar tem em comum somente o id, ou outra vriável em comum seu problema está resolvido. Se você tem dois ou mais bancos com a mesma estrutura (mesmas variáveis), basta usar o rbind. Exemplodf1 tem as variáveis id, sexo idade altura com 200 observações df2 tem as variáveis id, sexo idade altura com 100 observaçõesdf3 tem as variáveis id, sexo idade altura com 50 observações ndf-rbind(df1,df2,df3) O banco ndf terá sexo idade altura com 350 observações. [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 17 de Março de 2015 18:06, walmes . walmeszevi...@gmail.com escreveu: Eu acredito que já deve existir funções para essa tarefa em algum pacote. Bem, de qualquer forma, o código disponível no link abaixo pode ser útil para você. http://www.leg.ufpr.br/doku.php/ridiculas#como_fazer_a_justaposicao_de_varios_dataframes À disposição. Walmes. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] MCA gráfico
Bom dia caros amigos. Novamente retomo um problema que ainda não consegui resolver. Download MCA COMUNICACAO.docx # nesse link tem o script o gráfico e as dúvidas | | | | | | | | | | | Download MCA COMUNICACAO.docxDownload MCA COMUNICACAO.docx free. File hosted free by DataFileHost. | | | | Visualizar em www.datafileho... | Visualizado por Yahoo | | | | | Download mj1.csv # os dados | | | | | | | | | | | Download mj1.csvDownload mj1.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafileho... | Visualizado por Yahoo | | | | | No geral, necessito melhorar a aparência e a escala do gráfico e ainda não consegui. [ ]'s. Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
[R-br] GGPLOT
Meus caros amigos, estou tentando usar um exemplo de aplicação de MCA, script abaixo: #Usando MCA# load packagesrequire(FactoMineR)require(ggplot2) # load data teadata(tea)teanames(tea)# select these columnsnewtea = tea[, c(Tea, How, how, sugar, where, always)] # take a peekhead(newtea)names(newtea)str(newtea)cats = apply(newtea, 2, function(x) nlevels(as.factor(x)))cats # apply MCAmca1 = MCA(newtea, graph = FALSE) # list of resultsmca1 # table of eigenvaluesmca1$eig #We can use the package ggplot2() to get a nice plot:# data frame with variable coordinatesmca1_vars_df = data.frame(mca1$var$coord, Variable = rep(names(cats), cats)) # data frame with observation coordinatesmca1_obs_df = data.frame(mca1$ind$coord) par(cex.axis=0.5)# plot of variable categoriesggplot(data=mca1_vars_df,cex.axis=0.5, cex.lab=0.5, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df))) + geom_hline(yintercept = 0, colour = gray70) + geom_vline(xintercept = 0, colour = gray70) + geom_text(aes(colour=Variable)) + ggtitle(MCA plot of variables using R package FactoMineR) # MCA plot of observations and categoriesggplot(data = mca1_obs_df, aes(x = Dim.1, y = Dim.2)) + geom_hline(yintercept = 0, colour = gray70) + geom_vline(xintercept = 0, colour = gray70) + geom_point(colour = gray50, alpha = 0.7) + geom_density2d(colour = gray80) + geom_text(data = mca1_vars_df, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df), colour = Variable)) + ggtitle(MCA plot of variables using R package FactoMineR) + scale_colour_discrete(name = Variable) Como faço para mudar o tamanho da fonte do título, dos eixos, legenda, etc.? Tentei os comandos cex, cex.axis, cex.lab, mas não consegui. Alguém tem alguma sugestão? [ ]'s.Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] GGPLOT
Ok vou tentar. [ ]'s.Edson Lira Estatístico Manaus-Amazonas Em Segunda-feira, 23 de Fevereiro de 2015 17:07, Tropidurus Torquatus t.torqua...@gmail.com escreveu: Edson, Para construir meus gráficos com ggplot2 normalmente eu uso e abuso desses 2 sites: Acho que esse vai atender a suas dúvidas: http://www.cookbook-r.com/Graphs/ Esse abaixo é completo… mas pode gerar mais dúvidas: http://ggplot2.org Se não me falha a memória, podem ser usados duas formas… alteração direta na inserção dos dados mas aparente não seria a forma que você está utilizando (xlab/ylab), outra forma alterando o “theme”. Em 23/02/2015, à(s) 17:48, Edson Lira edinhoes...@yahoo.com.br escreveu: Meus caros amigos, estou tentando usar um exemplo de aplicação de MCA, script abaixo: #Usando MCA# load packagesrequire(FactoMineR)require(ggplot2) # load data teadata(tea)teanames(tea)# select these columnsnewtea = tea[, c(Tea, How, how, sugar, where, always)] # take a peekhead(newtea)names(newtea)str(newtea)cats = apply(newtea, 2, function(x) nlevels(as.factor(x)))cats # apply MCAmca1 = MCA(newtea, graph = FALSE) # list of resultsmca1 # table of eigenvaluesmca1$eig #We can use the package ggplot2() to get a nice plot:# data frame with variable coordinatesmca1_vars_df = data.frame(mca1$var$coord, Variable = rep(names(cats), cats)) # data frame with observation coordinatesmca1_obs_df = data.frame(mca1$ind$coord) par(cex.axis=0.5)# plot of variable categoriesggplot(data=mca1_vars_df,cex.axis=0.5, cex.lab=0.5, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df))) + geom_hline(yintercept = 0, colour = gray70) + geom_vline(xintercept = 0, colour = gray70) + geom_text(aes(colour=Variable)) + ggtitle(MCA plot of variables using R package FactoMineR) # MCA plot of observations and categoriesggplot(data = mca1_obs_df, aes(x = Dim.1, y = Dim.2)) + geom_hline(yintercept = 0, colour = gray70) + geom_vline(xintercept = 0, colour = gray70) + geom_point(colour = gray50, alpha = 0.7) + geom_density2d(colour = gray80) + geom_text(data = mca1_vars_df, aes(x = Dim.1, y = Dim.2, label = rownames(mca1_vars_df), colour = Variable)) + ggtitle(MCA plot of variables using R package FactoMineR) + scale_colour_discrete(name = Variable) Como faço para mudar o tamanho da fonte do título, dos eixos, legenda, etc.? Tentei os comandos cex, cex.axis, cex.lab, mas não consegui. Alguém tem alguma sugestão? [ ]'s.Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Transformar banco
Valeu Felipe e Listeiros, eu consegui também usando matriz. [ ]'s.Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 3 de Fevereiro de 2015 10:16, Felipe felipe.e.barle...@gmail.com escreveu: Talvez isso resolva: Transformando variáveis = y-sapply(dados,function(x) all(x%in%c('1','2','3','4','5'))) dados[,y]-sapply(dados[,y], as.numeric) ___ R-br mailing list R-br@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.
Re: [R-br] Correlação
Veja se é isso o que você precisa. #Examples ## EXAMPLE 1: ## Kidney samples from 291 salmon were split with one half of the ## samples sent to each of two laboratories where an IFAT test ## was run on each sample. The following results were obtained: ## Lab 1 positive, lab 2 positive: 19 ## Lab 1 positive, lab 2 negative: 10 ## Lab 1 negative, lab 2 positive: 6 ## Lab 1 negative, lab 2 negative: 256 dat - as.table(matrix(c(19,10,6,256), nrow = 2, byrow = TRUE)) colnames(dat) - c(L1-pos,L1-neg) rownames(dat) - c(L2-pos,L2-neg) epi.kappa(dat, method = fleiss, alternative = greater, conf.level = 0.95) ## The z test statistic is 11.53 (P 0.01). We accept the alternative ## hypothesis that the kappa statistic is greater than zero. ## The proportion of agreements after chance has been excluded is ## 0.67 (95% CI 0.56 to 0.79). We conclude that, on the basis of ## this sample, that there is substantial agreement between the two ## laboratories. ## EXAMPLE 2 (from Watson and Petrie 2010, page 1170): ## Silva et al. (2007) compared an early pregnancy enzyme-linked immunosorbent ## assay test for pregnancy associated glycoprotein on blood samples collected ## from lactating dairy cows at day 27 after artificial insemination with ## transrectal ultrasound (US) diagnosis of pregnancy at the same stage. ## The results were as follows: ## ELISA positive, US positive: 596 ## ELISA positive, US negative: 61 ## ELISA negative, US positive: 29 ## ELISA negative, Ul negative: 987 dat - as.table(matrix(c(596,61,29,987), nrow = 2, byrow = TRUE)) colnames(dat) - c(US-pos,US-neg) rownames(dat) - c(ELISA-pos,ELISA-neg) epi.kappa(dat, method = watson, alternative = greater, conf.level = 0.95) ## The proportion of agreements after chance has been excluded is ## 0.89 (95% CI 0.86 to 0.91). We conclude that that there is substantial ## agreement between the two pregnancy diagnostic methods. [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 6 de Janeiro de 2015 6:20, Leonard Assis assis.leon...@gmail.com escreveu: Acredito que para seu problema, correlação ajude pouco. Ideal seria 'concordância'. Procure sobre agreement coefficient no r searchEm 05/01/2015 16:57, Gustavo Faria gustavo...@gmail.com escreveu: Boa tarde, estou com dificuldade (iniciantes ¬¬ hehe) em fazer uma correlação com dados qualitativos, meus dados são 3 metodologias de identificação de organismos, e eu gostaria de saber o quão parecido e/ou divergente eles foram nas identificações realizadas, como poderia montar essa matriz? ja que tenho somente as identificações realizadas pelas 3 metodologias -- --- Gustavo de Faria Silva Biólogo herpetólogo Crbio 87421/04-P Pós-Graduado em Consultoria e Licenciamento Ambiental Mestrando em Zoologia de Vertebrados Telefone: (31) 9579-8473 (Vivo) ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Correlação
Não esqueça de instalar o epiR install.packages(epiR,dep=T) Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 6 de Janeiro de 2015 7:22, Edson Lira edinhoes...@yahoo.com.br escreveu: Veja se é isso o que você precisa. #Examples ## EXAMPLE 1: ## Kidney samples from 291 salmon were split with one half of the ## samples sent to each of two laboratories where an IFAT test ## was run on each sample. The following results were obtained: ## Lab 1 positive, lab 2 positive: 19 ## Lab 1 positive, lab 2 negative: 10 ## Lab 1 negative, lab 2 positive: 6 ## Lab 1 negative, lab 2 negative: 256 dat - as.table(matrix(c(19,10,6,256), nrow = 2, byrow = TRUE)) colnames(dat) - c(L1-pos,L1-neg) rownames(dat) - c(L2-pos,L2-neg) epi.kappa(dat, method = fleiss, alternative = greater, conf.level = 0.95) ## The z test statistic is 11.53 (P 0.01). We accept the alternative ## hypothesis that the kappa statistic is greater than zero. ## The proportion of agreements after chance has been excluded is ## 0.67 (95% CI 0.56 to 0.79). We conclude that, on the basis of ## this sample, that there is substantial agreement between the two ## laboratories. ## EXAMPLE 2 (from Watson and Petrie 2010, page 1170): ## Silva et al. (2007) compared an early pregnancy enzyme-linked immunosorbent ## assay test for pregnancy associated glycoprotein on blood samples collected ## from lactating dairy cows at day 27 after artificial insemination with ## transrectal ultrasound (US) diagnosis of pregnancy at the same stage. ## The results were as follows: ## ELISA positive, US positive: 596 ## ELISA positive, US negative: 61 ## ELISA negative, US positive: 29 ## ELISA negative, Ul negative: 987 dat - as.table(matrix(c(596,61,29,987), nrow = 2, byrow = TRUE)) colnames(dat) - c(US-pos,US-neg) rownames(dat) - c(ELISA-pos,ELISA-neg) epi.kappa(dat, method = watson, alternative = greater, conf.level = 0.95) ## The proportion of agreements after chance has been excluded is ## 0.89 (95% CI 0.86 to 0.91). We conclude that that there is substantial ## agreement between the two pregnancy diagnostic methods. [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 6 de Janeiro de 2015 6:20, Leonard Assis assis.leon...@gmail.com escreveu: Acredito que para seu problema, correlação ajude pouco. Ideal seria 'concordância'. Procure sobre agreement coefficient no r searchEm 05/01/2015 16:57, Gustavo Faria gustavo...@gmail.com escreveu: Boa tarde, estou com dificuldade (iniciantes ¬¬ hehe) em fazer uma correlação com dados qualitativos, meus dados são 3 metodologias de identificação de organismos, e eu gostaria de saber o quão parecido e/ou divergente eles foram nas identificações realizadas, como poderia montar essa matriz? ja que tenho somente as identificações realizadas pelas 3 metodologias -- --- Gustavo de Faria Silva Biólogo herpetólogo Crbio 87421/04-P Pós-Graduado em Consultoria e Licenciamento Ambiental Mestrando em Zoologia de Vertebrados Telefone: (31) 9579-8473 (Vivo) ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] Problema com Dates no Macbook pro
Valeu Marcus, funcionou, não sei o que estava errado. Valeu!!![ ]'sEdson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 2 de Janeiro de 2015 20:32, Marcus Nunes marcus.nu...@gmail.com escreveu: Edson, Se entendi corretamente o teu pedido, o que tu deseja fazer está no código abaixo. Rodei o código no Mac OS 10.10.1 e R 3.1.1 Abs, ## rm(list = ls()) require(ff)cb - read.csv2(~/Desktop/TESTE_DATA.csv) cb[cb==] - NA #cb$chega-format(as.Date(cb$ARRIVAL, format=%d/%b/%Y))cb$chega - as.Date(cb$ARRIVAL,%d/%m/%Y) cb$ent_ano - format(cb$chega, %Y) # anocb$ent_mes - format(cb$chega, %m) # anocb$ent_dia - format(cb$chega, %d) # ano 2015-01-02 21:53 GMT-02:00 Edson Lira edinhoes...@yahoo.com.br: Estou com problemas na conversão de datas. Quero separar dia, mês e ano. Abaixo o link de uma amostra dos dados com as variáveis de datas Download TESTE DATA.csv | | | | | | | | | | | Download TESTE DATA.csvDownload TESTE DATA.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Link do Script Download CB.R | | | | | | | | | | | Download CB.RDownload CB.R free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Feliz 2015 para os amigos!!! Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list R-br@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. -- Marcus Nunes http://marcusnunes.me/ ___ R-br mailing list R-br@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.
[R-br] Problema com Dates no Macbook pro
Estou com problemas na conversão de datas. Quero separar dia, mês e ano. Abaixo o link de uma amostra dos dados com as variáveis de datas Download TESTE DATA.csv | | | | | | | | | | | Download TESTE DATA.csvDownload TESTE DATA.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Link do Script Download CB.R | | | | | | | | | | | Download CB.RDownload CB.R free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Feliz 2015 para os amigos!!! Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Homogeneidade de variancias
Veja se te ajuda | bartlett.test {stats} | R Documentation | Bartlett Test of Homogeneity of Variances Description Performs Bartlett's test of the null that the variances in each of the groups (samples) are the same. Usage bartlett.test(x, ...) ## Default S3 method: bartlett.test(x, g, ...) ## S3 method for class 'formula' bartlett.test(formula, data, subset, na.action, ...) Arguments | x | a numeric vector of data values, or a list of numeric data vectors representing the respective samples, or fitted linear model objects (inheriting from class lm). | | g | a vector or factor object giving the group for the corresponding elements of x. Ignored if x is a list. | | formula | a formula of the form lhs ~ rhs where lhs gives the data values and rhs the corresponding groups. | | data | an optional matrix or data frame (or similar: see model.frame) containing the variables in the formulaformula. By default the variables are taken from environment(formula). | | subset | an optional vector specifying a subset of observations to be used. | | na.action | a function which indicates what should happen when the data contain NAs. Defaults togetOption(na.action). | | ... | further arguments to be passed to or from methods. | Details If x is a list, its elements are taken as the samples or fitted linear models to be compared for homogeneity of variances. In this case, the elements must either all be numeric data vectors or fitted linear model objects, g is ignored, and one can simply use bartlett.test(x) to perform the test. If the samples are not yet contained in a list, usebartlett.test(list(x, ...)).Otherwise, x must be a numeric data vector, and g must be a vector or factor object of the same length as x giving the group for the corresponding elements of x. Value A list of class htest containing the following components: | statistic | Bartlett's K-squared test statistic. | | parameter | the degrees of freedom of the approximate chi-squared distribution of the test statistic. | | p.value | the p-value of the test. | | method | the character string Bartlett test of homogeneity of variances. | | data.name | a character string giving the names of the data. | References Bartlett, M. S. (1937). Properties of sufficiency and statistical tests. Proceedings of the Royal Society of London Series A160, 268–282. See Also var.test for the special case of comparing variances in two samples from normal distributions; fligner.test for a rank-based (nonparametric) k-sample test for homogeneity of variances; ansari.test and mood.test for two rank based two-sample tests for difference in scale. Examples require(graphics) plot(count ~ spray, data = InsectSprays) bartlett.test(InsectSprays$count, InsectSprays$spray) bartlett.test(count ~ spray, data = InsectSprays) Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 2 de Janeiro de 2015 10:44, Andre Oliveira andreolso...@yahoo.com.br escreveu: Qual seria o procedimento para para subdividida? Obrigado pessoal! André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES Em Quinta-feira, 1 de Janeiro de 2015 23:23, Fernando Antonio de souza nandodeso...@gmail.com escreveu: VC pode criar uma variável interação e fazer o Bartlett. Teste sobre ela Dados$interação-interaction(fator1,fator2) bartlett.test(resid(modelo2), dados$interaction)Em 01/01/2015 18:57, Andre Oliveira andreolso...@yahoo.com.br escreveu: Mas neste caso como aplicar o bartlett.test? bartlett.test(modelo2$residuals,fat1) # Homogeneidade de variâncias bartlett.test(modelo2$residuals,fat2) # Homogeneidade de variâncias Só assim? Por fatores separadamente? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES Em Quinta-feira, 1 de Janeiro de 2015 18:40, henrique jose de paula alves paula alves jpahenri...@gmail.com escreveu: O teste de Bartley deve funcionar bem. Em 1 de janeiro de 2015 17:52, Andre Oliveira andreolso...@yahoo.com.br escreveu: Pessoal boa tarde, dado meu modelo em fatorial. modelo=aov(respsota~factor1*factor2) Como avaliar homogeneidade de variâncias? O que fiz aqui foi ..par(mfrow=c(2,2)) boxplot(resp~fat1) boxplot(resp~fat2) boxplot(resp~fat1*fat2) interaction.plot(fat1,fat2,resp,ylab=Médias,main=Interação Fat1 vs Fat2) Mas está muito subjetivo. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES ___ R-br mailing list R-br@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. -- Henrique José de Paula AlvesPós-graduando em
Re: [R-br] RES: Off
text/html; charset=iso-8859-1: Unrecognized ___ R-br mailing list R-br@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.
[R-br] Off
Caros amigos, vê se alguém poderia me ajudar. Tenho o seguinte problema. Uma empresa recebeu um lote de 20 tvs, a empresa rejeita o lote se o mesmo apresenta pelo menos 4 tvs com defeito, a probabilidade de uma TV apresentar defeito é 1%. Qual a probabilidade de rejeitar o lote todo?Estou sem calculadora do CEL e sem PC. Enviado do Yahoo Mail para iPhone ___ R-br mailing list R-br@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.
Re: [R-br] RES: [OFF-topic] Processar uma base de dados grande - ENEM 2012 - 6 milhões de registros - só 40 segundos
Experimenta o pacote ff Edson Lira Estatístico Manaus-Amazonas Em Domingo, 26 de Outubro de 2014 16:22, Mauro Sznelwar sznel...@uol.com.br escreveu: !--#yiv1316119439 _filtered #yiv1316119439 {font-family:Calibri;panose-1:2 15 5 2 2 2 4 3 2 4;}#yiv1316119439 #yiv1316119439 p.yiv1316119439MsoNormal, #yiv1316119439 li.yiv1316119439MsoNormal, #yiv1316119439 div.yiv1316119439MsoNormal {margin:0cm;margin-bottom:.0001pt;font-size:12.0pt;font-family:Times New Roman, serif;}#yiv1316119439 a:link, #yiv1316119439 span.yiv1316119439MsoHyperlink {color:blue;text-decoration:underline;}#yiv1316119439 a:visited, #yiv1316119439 span.yiv1316119439MsoHyperlinkFollowed {color:purple;text-decoration:underline;}#yiv1316119439 span.yiv1316119439EstiloDeEmail17 {font-family:Calibri, sans-serif;color:#1F497D;}#yiv1316119439 .yiv1316119439MsoChpDefault {font-size:10.0pt;} _filtered #yiv1316119439 {margin:70.85pt 3.0cm 70.85pt 3.0cm;}#yiv1316119439 div.yiv1316119439Section1 {}--Existeeste pacote para o Windows? Tentei baixar e disse que não está disponível paraa versão R 3.1 Pessoal, eu fiz um teste com o MonetDB importando um arquivocsv do censo demográfico 2010 do estado de São Paulo, que tem 2.2 gb, 224colunas e 4038218 linhas. Importou em 109.1149 segundos. Meu computador (MacBook Pro) tem processador i7 da primeirageração, 16 gb de ram e ssd. ##-- library(MonetDB.R) #criando um banco de dados startscript - monetdb.server.setup( database.directory = paste0( getwd() , /MonetDB ) , , dbname =db1, dbport =50002 ) pid - monetdb.server.start(startscript) conn - dbConnect(MonetDB.R(), monetdb://localhost:50002/db1) # 224 colunas # 4038218 linhas # arquivo com 2.2 GB # 109.149 segundos system.time({ monetdb.read.csv(conn, 'censo00sp1.csv',censo00sp, 4038218) }) ##-- Att Roney | | Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa. | ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
Re: [R-br] tabela com frequencia relativa e absoluta
tabelas-function(x){ Qtde-table(x) #absoluta Perc-100*prop.table(table(x)) #relativa Tab-cbind(Qtde,Perc) Total-colSums(Tab) rbind(Tab,Total) } options(digits=3) apply(bd,2,tabelas) Espero ter ajudado[ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 22 de Outubro de 2014 21:55, lutipilo...@yahoo.com.br lutipilo...@yahoo.com.br escreveu: !--#yiv2873321972 p.yiv2873321972MsoListParagraph, #yiv2873321972 li.yiv2873321972MsoListParagraph, #yiv2873321972 div.yiv2873321972MsoListParagraph {margin-top:0in;margin-right:0in;margin-bottom:0in;margin-left:.5in;margin-bottom:.0001pt;}#yiv2873321972 p.yiv2873321972MsoNormal, #yiv2873321972 li.yiv2873321972MsoNormal, #yiv2873321972 div.yiv2873321972MsoNormal {margin:0in;margin-bottom:.0001pt;}#yiv2873321972 p.yiv2873321972MsoListParagraphCxSpFirst, #yiv2873321972 li.yiv2873321972MsoListParagraphCxSpFirst, #yiv2873321972 div.yiv2873321972MsoListParagraphCxSpFirst, #yiv2873321972 p.yiv2873321972MsoListParagraphCxSpMiddle, #yiv2873321972 li.yiv2873321972MsoListParagraphCxSpMiddle, #yiv2873321972 div.yiv2873321972MsoListParagraphCxSpMiddle, #yiv2873321972 p.yiv2873321972MsoListParagraphCxSpLast, #yiv2873321972 li.yiv2873321972MsoListParagraphCxSpLast, #yiv2873321972 div.yiv2873321972MsoListParagraphCxSpLast {margin-top:0in;margin-right:0in;margin-bottom:0in;margin-left:.5in;margin-bottom:.0001pt;line-height:115%;}--Olá pessoal, tem algum comando para fazer uma tabela única direto no R com frequencia relativa e absoluta? E tb valores de chi² e p.Estou usando os comandos abaixo, mas o resultado aparece em tabelas individuais. table(q11,q13) rowPercents(table(q11,q13)) f - chiSquare(q13 ~ q11) f Atenciosamente,Luciane Pilotto Enviado do Email do Windows De: Hélio Gallo Rocha Enviado: quarta-feira, 22 de outubro de 2014 09:07 Para: r-br Mauro, bom dia. Não sei qual seu objetivo no uso desta função, assim vou te passar o histórico e se quiser dando idéias, será bem vindoEstou estudando geoestatística e é interessante saber que regiões tem seus atributos parecidos.A idéia inicial era utilizar a função skater do pacote spdep. Como tenho pontos teria de criar poligonos, que foi realizado com a função deldir, e listar os poligonos. A listagem era um problemão, pois teria de crira poligono por poligono, como no CRM inicial.Problema resolvido com a preciosa ajuda do Rafael. que sugeriu usar a função poly2nb do pacote sddep Com a função deldir apenas gera poligonos assim foi necessário inserir dados aos mesmo e torná-los um SpatialPolygonsDataFrame Dai pra frente o código está em estudo, veja:require(deldir) ;require(spdep)## 1. Coordenadasn - 40 set.seed(47)x - runif(n)y - runif(n)plot(x,y)## 2. Dados set.seed(1389)dados - data.frame(rnorm(n,10,sd=5), rnorm(n,15,sd=5),rnorm(n,20,sd=5))colnames(dados)- paste(var,1:3,sep=)row.names(dados) - paste(reg, 1:n, sep = )dadosclass(dados) ## 3. gerando coordenadas dos poligonos com função deldird - deldir(x,y) ttt - triang.list(d)plot(ttt)plot(ttt,border=red,showrect=TRUE)sss - tile.list(d)plot(sss)plot(ttt,add=TRUE,border=blue) ## 4. Gerando poligonos### Rafael Tieppo G R A N D E A J U D A ! ! ! ! FUN_VOR_MAT - function (A,TT) #TT é o arquivo, N o número do POL { AA - 1 BB - 2 NROW - length (A[[TT]]$x) + 1 #add primeiro ponto VORONOI_MATRIX - matrix(0,NROW,2) VORONOI_MATRIX[1:NROW,AA] - A[[TT]]$x[1:NROW] VORONOI_MATRIX[1:NROW,BB] - A[[TT]]$y[1:NROW] VORONOI_MATRIX[NROW,AA] - A[[TT]]$x[1] #+ 1 add primeiro ponto VORONOI_MATRIX[NROW,BB] - A[[TT]]$y[1] #+ 1 add primeiro ponto if(TT == NROW){ print(17/10/2014, [hidden email] )} return(VORONOI_MATRIX) } ## Agora vc vai usar a funcao:AREAS- lapply(1:length(sss), function (x) ( FUN_VOR_MAT(sss,x)) ) summary(AREAS)## 5. Formando PolygonsPOLS- lapply(1:length(AREAS), function(x) Polygons(list(Polygon(AREAS[[x]])), ID = paste(reg, x, sep = ))) summary(POLS) class(POLS) SPol - SpatialPolygons(POLS) class(SPol) plot(SPol) bbox(SPol) ## 6. Inserindo dados aos poligonos SPolDF - SpatialPolygonsDataFrame(SPol, dados) plot(SPolDF) ## 7. Código para visualizar áreas parecidas ## em estudo(coordenadas=cbind(x,y))class(coordenadas) ## rn - sapply(slot(SPolDF, polygons), function(x) slot(x, ID)) knn1 - knearneigh(coordenadas, 1) knn2 - knearneigh(coordenadas, 2) nb1 - knn2nb(knn1, row.names=rn) nb2 - knn2nb(knn2, row.names=rn) ##1diffs - diffnb(nb2, nb1)plot(SPolDF, border=grey)plot(nb1, coordenadas, add=TRUE)plot(diffs, coordenadas, add=TRUE, col=red)title(main=paste(Differences (red) in Columbus GAL weights (black), and polygon generated queen weights, sep=\n)) ##2 xxx - poly2nb(SPolDF, queen=FALSE)dxxx - diffnb(xxx, nb1)plot(SPolDF, border = grey)plot(nb1, coordenadas, add = TRUE)plot(dxxx, coordenadas, add = TRUE, col = red)title(main=paste
Re: [R-br] RES: Problema
Ao converter para csv, para leitura use o pacote ff que foi criado, se não me engano, para ler arquivo muito grande. install.package(ff) Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 16 de Outubro de 2014 23:30, Mauro Sznelwar sznel...@uol.com.br escreveu: Eu sugeriria que convertesse o seu Excell xls ou xl em csv e carregasse no R com read.csv (‘nome do arquivo.csv’) Boa noite, pessoal! Bem, sou iniciante no uso do R e estou com um problema. Estou trabalhando com muitos dados que estou exportando do excel e estou utilizando o pacote XLConnect e quando coloco o programa para executar ele inicia normalmente, quanto chega em determinado trecho aparece o seguinte erro: Error: cannot allocate vector size Gb Tenho 8 de memória ram, sistema de 64 btis e o windows 7 O que devo fazer??? Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Ordenamento entre colunas
Uma alternativa é um novo banco, novo-data.frame(Velho,A1, A2, A3, A4, A5, B1, B2, B3) [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 17 de Outubro de 2014 8:13, Giselle Davi giselle_d...@yahoo.com.br escreveu: Prezados, Como eu consigo ordenar as seguintes colunas de uma matriz no R : A5, A3, A1, A4, A2, B2, B1, B3 ? O objetivo é que elas fiquem assim: A1, A2, A3, A4, A5, B1, B2, B3 Utilizei o comando order, mas o ordenamento ocorre somente dentro das colunas e não entre como eu preciso. Desde já agradeço a atenção. Att, Giselle ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
[R-br] boxplots
Caros amigos da lista, tenho um banco com as variáveis: certos, fase(antes/depois) e ocupa(com 5 categorias). A minha pergunta é: Tem como fazer boxplots de acertos para as fases antes e depois por ocupação. A idéia é colocar a caixa antes e depois lado a lado. O banco para análise está no link abaixo: Download novo.csv Download novo.csv Download novo.csv free. File hosted free by DataFileHost. Visualizar em www.datafileho... Visualizado por Yahoo Um grande abraço. Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] mapa: rosa dos ventos
Aline, essa função é do pacote sp. Você já tentou reinstalar esse pacote? Ou então tenta uma busca na internet compassRose e no R ??compassRose Acho que não ajudei muito mas seriam as primeiras coisas que tentaria para depois buscar outra saída. [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 9 de Outubro de 2014 12:36, Aline Lipsky line_l...@hotmail.com escreveu: Boa tarde caros colegas, Estou plotando um mapa com o seguinte script: library(maptools) library(maps) library(mapdata) library(mapproj) library(sp) SHP22 - readShapeSpatial(hidrografia_poligonos_utm22.shp) #Plotar o mapa x11() par(mar=c(2,2,2,4.2)) # ajustando as margens (bottom, left, top, right) plot(SHP22, col='light blue', axes=F, ylim=c(6441000,653),xlim=c(38,41)) #adicionando eixo no mapa axis(4, las=1, cex.axis=0.9) #adicionando latitude do lado direito do mapa axis(1, las=1, cex.axis=0.9) #adicionando longitude abaixo do mapa Mas quando tento plotar a rosa dos ventos com a seguinte função: compassRose(43,6525000, rot=0, cex=0.65) aparece o erro: Error: could not find function compassRose. Alguém saberia me dizer o que posso fazer? Grata desde já pela atenção, Att., Aline ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Kruskal Wallis - identificar o tratamento diferente
Veja se ajuda. library(agricolae) data(corn) attach(corn) str(corn) comparison-kruskal(observation,method,group=TRUE, main=corn) comparison-kruskal(observation,method,p.adj=bon,group=FALSE, main=corn) detach(corn) Edson Lira Estatístico Manaus-Amazonas Em Sexta-feira, 26 de Setembro de 2014 11:09, Daniela Recchia daniela_recc...@yahoo.com.br escreveu: Boa tarde! Já pesquisei um pouco sobre isso mas não encontrei nada, tenho 3 variáveis X,Y,Z (não normais) e estou usando o kruskal-wallis para compara-las. Recebo então que elas são estatisticamente diferentes mas não consigo saber qual delas, e a X diferente da Y, é a X diferente da Z, é a Y diferente da Z? Alguém sabe se existe alguma nova função que fornece essa informação extra na anova não paramétrica? Abraços, 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 ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Codificação
Gisele, na lista foi adotado a rotina de não anexar arquivos aos emails e sim postar os arquivos, por exemplo, no site http://www.datafilehost.com/ Busque no Google Avaliar dados de genética no software R Olha o que eu achei: http://geneticapopulacional.blogspot.com.br/ fastStructure fastStructure é um algoritmo para inferir a estrutura da população a partir de grandes volumes de dados genótipo SNP. Baseia-se em um quadro Bayesian variacional para posterior inferência e é escrito em python2.x. Aqui, resumimos como configurar este pacote de software, compilar os scripts de C e Cython e executar o algoritmo em um teste simulado dataset genótipo. Boa sorte! Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 3 de Setembro de 2014 10:45, Giselle Davi giselle_d...@yahoo.com.br escreveu: Prezados Colaboradores, Tenho um arquivo com uma matriz formada por letras ( SNP: bases nitrogenadas: A, C, T, G) e preciso codificá-la para números, onde as letras A e T recebam 1 e as letras C e G recebam 0, formando assim uma matriz de genótipos Gn x m = {0;1;2} onde n é o número de indivíduos e m o número de SNPs. Anexado a este e-mail encontra-se parte dessa matriz com dimensão 6 x 30 ( Facilitando, portanto, o envio) que necessita dessa codificação. Matematicamente, a sequência de SNP é atribuída ( A qual eu não precisei construir, pois a sequência foi originada por genotipagem ) como: Dado um valor para PA (probabilidade de ocorrer a base A) e para PC (probabilidade de C), de modo que PA= PT, PC = PG e PA+PC = 0; 5, então as bases são geradas do seguinte modo: Se y existe (E virado para direita) [0;PA), a base é A; Se y existe [PA;0:5), a base é C; Se y existe [0:5; (0:5+PC)), a base é G; Se y existe [(0:5+PC);1], a base é T. Cada indivíduo possui um par de cromossomos homólogos, que se mantêm juntos por ligações entre as bases A e T e ligações entre as bases C e G. A codificação das sequências ( SNPs) deve ser feita, portanto, do seguinte modo: Contando o alelo de interesse (no caso, A) e seu correspondente no pareamento (T) em cada fita numa posição do genótipo, e atribuir 1 quando um dos dois (A ou T) ocorre e 0 se não ocorrem. Esta codificação resulta em: 0 - CC, CG, GC, GG)P(0) = 4*PC elevado a 2 1 - AC, AG, TC, TG, CA, GA, CT, GT)P(1) = 8*PA*PC 2 - AA, AT, TA, TT )P(2) = 4*PA elevado a 2 Sendo assim, preciso utilizar estas informações matemáticas no R de modo a proceder a codificação. Alguém poderia me ajudar ? Desde já agradeço. Att, Giselle ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] RES: RES: Arquivos com extensão json.
Muito obrigado!!! Vou tentar sua sugestão. [ ]'s Edson Lira Estatístico Manaus-Amazonas Em Domingo, 31 de Agosto de 2014 17:51, Leonard de Assis assis.leon...@gmail.com escreveu: Lira, O site do opendatakit explica como é a estrutura da ferramenta, inclusive do banco de dados (Sim, eles guardam tudo em um banco de dados quando você sobe o xml do smartphone pro servidor). Em 2012 minha empresa desenvolveu internamente uma solução baseada em opendatakit e pelo que eu lembro, o servidor de agregação era mysql (O nosso está hospedado na amazon). Eu lembro que tivemos problemas com encoding, mas apenas referente aos questionários em xml. Resolvemos isso colocando na 1ª linha do xml o encoding em utf-8 e resolveu Pode ser útil também este arquivo: http://opendatakit.org/downloads/download-info/odk-csvconverter-v1-4-3-production-jar/ De:R-br [mailto:r-br-boun...@listas.c3sl.ufpr.br] Em nome de Edson Lira Enviada em: sábado, 30 de agosto de 2014 21:49 Para: r-br@listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: Arquivos com extensão json. Assis, como vejo isso no odk? Tem algum comando? [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Sábado, 30 de Agosto de 2014 19:35, Leonard de Assis assis.leon...@gmail.com escreveu: Lira, O ODK usa um bd Mysql para armazenar tanto os questionários quanto as respostas. Verifique no local onde você criou o servidor de agregação, qual é o ‘encoding’ do bd mysql que está gerenciando o processo De:R-br [mailto:r-br-boun...@listas.c3sl.ufpr.br] Em nome de Benilton Carvalho Enviada em: quinta-feira, 28 de agosto de 2014 09:12 Para: r-br; Edson Lira Assunto: Re: [R-br] Arquivos com extensão json. pacote rjson Em 28 de agosto de 2014 01:52, Edson Lira edinhoes...@yahoo.com.br escreveu: Caros amigos, uso o aplicativo ODKAggregate para gerenciar dados coletados via smartfones, ele me permite exportar os dados coletados no formato csv, porém, sua conversão tem problemas com a acentuação, já tentei usar a opcão fileEncod e não consegui corrigir o problema. Quando exporto no formato json, o problema da acentuação se resolve, consigo abrir no notepad, e ver uma estrutura tipo de programação. No link abaixo, tem o conjunto de dados exportado pelo ODKAggregate. http://www.datafilehost.com/d/d60c05e4 Alguém poderia dar uma olhada e verificar se há uma saída analisar esses dados com o uso do R? [ ]'s. Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] RES: Arquivos com extensão json.
Assis, como vejo isso no odk? Tem algum comando? [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Sábado, 30 de Agosto de 2014 19:35, Leonard de Assis assis.leon...@gmail.com escreveu: Lira, O ODK usa um bd Mysql para armazenar tanto os questionários quanto as respostas. Verifique no local onde você criou o servidor de agregação, qual é o ‘encoding’ do bd mysql que está gerenciando o processo De:R-br [mailto:r-br-boun...@listas.c3sl.ufpr.br] Em nome de Benilton Carvalho Enviada em: quinta-feira, 28 de agosto de 2014 09:12 Para: r-br; Edson Lira Assunto: Re: [R-br] Arquivos com extensão json. pacote rjson Em 28 de agosto de 2014 01:52, Edson Lira edinhoes...@yahoo.com.br escreveu: Caros amigos, uso o aplicativo ODKAggregate para gerenciar dados coletados via smartfones, ele me permite exportar os dados coletados no formato csv, porém, sua conversão tem problemas com a acentuação, já tentei usar a opcão fileEncod e não consegui corrigir o problema. Quando exporto no formato json, o problema da acentuação se resolve, consigo abrir no notepad, e ver uma estrutura tipo de programação. No link abaixo, tem o conjunto de dados exportado pelo ODKAggregate. http://www.datafilehost.com/d/d60c05e4 Alguém poderia dar uma olhada e verificar se há uma saída analisar esses dados com o uso do R? [ ]'s. Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Read .spt e export como .csv
?write.csv ## Not run: ## To write a CSV file for input to Excel one might use x - data.frame(a = I(a \ quote), b = pi) write.table(x, file = foo.csv, sep = ,, col.names = NA, qmethod = double) ## and to read this file back into R one needs read.table(foo.csv, header = TRUE, sep = ,, row.names = 1) ## NB: you do need to specify a separator if qmethod = double. ### Alternatively write.csv(x, file = foo.csv) read.csv(foo.csv, row.names = 1) ## or without row names write.csv(x, file = foo.csv, row.names = FALSE) read.csv(foo.csv) ## To write a file in MacRoman for simple use in Mac Excel 2004/8 write.csv(x, file = foo.csv, fileEncoding = macroman) ## or for Windows Excel 2007/10 write.csv(x, file = foo.csv, fileEncoding = UTF-16LE) Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 27 de Agosto de 2014 9:08, Alisson Lucrécio alisson.lucre...@ifgoiano.edu.br escreveu: Caros Colegas, Bom dia. Tenho um arquivo .spt (.SPT) e que quero transformar para .csv (.CSV), algum colega saberia como fazer isso no R? Obrigado. -- Alisson Lucrecio da Costa ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
Re: [R-br] Arquivos com extensão json.
Obrigado Benilton Edson Lira Estatístico Manaus-Amazonas Em Quinta-feira, 28 de Agosto de 2014 8:12, Benilton Carvalho beniltoncarva...@gmail.com escreveu: pacote rjson Em 28 de agosto de 2014 01:52, Edson Lira edinhoes...@yahoo.com.br escreveu: Caros amigos, uso o aplicativo ODKAggregate para gerenciar dados coletados via smartfones, ele me permite exportar os dados coletados no formato csv, porém, sua conversão tem problemas com a acentuação, já tentei usar a opcão fileEncod e não consegui corrigir o problema. Quando exporto no formato json, o problema da acentuação se resolve, consigo abrir no notepad, e ver uma estrutura tipo de programação. No link abaixo, tem o conjunto de dados exportado pelo ODKAggregate. http://www.datafilehost.com/d/d60c05e4 Alguém poderia dar uma olhada e verificar se há uma saída analisar esses dados com o uso do R? [ ]'s. Edson Lira Estatístico Manaus-Amazonas ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.
[R-br] Arquivos com extensão json.
Caros amigos, uso o aplicativo ODKAggregate para gerenciar dados coletados via smartfones, ele me permite exportar os dados coletados no formato csv, porém, sua conversão tem problemas com a acentuação, já tentei usar a opcão fileEncod e não consegui corrigir o problema. Quando exporto no formato json, o problema da acentuação se resolve, consigo abrir no notepad, e ver uma estrutura tipo de programação. No link abaixo, tem o conjunto de dados exportado pelo ODKAggregate. http://www.datafilehost.com/d/d60c05e4 Alguém poderia dar uma olhada e verificar se há uma saída analisar esses dados com o uso do R? [ ]'s. Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Erro em q$limits
Walmes, peço desculpas pela falta de detalhes, estou usando o pacote qcc (http://cran.r-project.org/web/packages/qcc/index.html) Link dos dados: http://www.datafilehost.com/d/6655a016 A rotina completa abaixo: require(qcc) old - qcc.options() # save defaults qcc.options(cex.stats) # get a single parameter qcc.options(cex.stats=0.8) # change parameters qcc.options(cex.axis=0.7) # change parameters qcc.options(bg.margin=white) par(cex=0.7,cex.lab=0.8,cex.axis=0.7) qcc(M04, type=xbar,center=60,title=,nsigmas=2,cex=0.7,cex.lab=0.8,cex.axis=0.7, xlab=Dia, ylab=Tempo (min),cex=0.7,label.limits = c(, ), restore.par = FALSE) axis(side = 4, at = q$limits[1], labels = 52, tick = T, cex.axis = 0.7, las = 1, line = -0.5) axis(side = 4, at = q$limits[2], labels = 68, tick = T, cex.axis = 0.7, las = 1, line = -0. [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Terça-feira, 19 de Agosto de 2014 7:44, Rodrigo Coster rcos...@gmail.com escreveu: Edson, Acredito que o erro seja que tu não possui a variável q definida, e por isso não consegue desenhar o axis (q, por padrão, é uma função, por isso a mensagem de erro) 2014-08-18 22:09 GMT-03:00 walmes . walmeszevi...@gmail.com: Prezado Edson, O seu código é não reproduzível. Você não forneceu meios para reproduzir os dados nem mencionou a que pacote pertence a função qqc(). Alguém bem intencionado em ajudá-lo deve proceder como? Gerar dados que você poderia fornecer, procurar pela origem da função que você usa e achar a fonte do erro? Será que isso não é ser muito otimista? Seja atencioso com seus potenciais colaboradores, forneça meios de reproduzir o seu problema, forneça um CMR. À disposição. Walmes. ___ R-br mailing list R-br@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. ___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.
[R-br] Erro em q$limits
Tenho um bd com a seguinte estrutura: head(M04) A B C D 1 60.9460 62.2680 58.5355 60.3665 2 57.3075 57.6275 85.6005 65.5345 3 73.4505 48.5370 71.5110 65.6680 4 58.2090 67.5285 64.3855 71.5910 5 53.2815 61.4225 51.4050 65.2255 7 60.7725 72.0105 71.2575 54.3345 Alguém poderia me indicar qual o erro que estou cometendo na rotina abaixo: par(cex=0.7,cex.lab=0.8,cex.axis=0.7) qcc(M04, type=xbar,center=60,title=,nsigmas=2,cex=0.7,cex.lab=0.8,cex.axis=0.7, xlab=Dia, ylab=Tempo (min),cex=0.7,label.limits = c(, ), restore.par = FALSE) axis(side = 4, at = q$limits[1], labels = 52, tick = T, cex.axis = 0.7, las = 1, line = -0.5) axis(side = 4, at = q$limits[2], labels = 68, tick = T, cex.axis = 0.7, las = 1, line = -0.5) O erro na janela do Log é: axis(side = 4, at = q$limits[1], labels = 52, tick = T, cex.axis = 0.7, las = 1, line = -0.5) Erro em q$limits : objeto de tipo 'closure' não possível dividir em subconjuntos Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.
Re: [R-br] Gráfico de Ishikawa
Muito obrigado!!! [ ]'s. Edson Lira Estatístico Manaus-Amazonas Em Quarta-feira, 13 de Agosto de 2014 3:47, Tropidurus Torquatus t.torqua...@gmail.com escreveu: Edson, Eu coloquei os nomes entre “ “ e com espaço: cause.and.effect(cause = list(Método de foto cópia = c('Colocação do Original', 'Tempo de secagem'), Liquido = c('Qualidade', 'Qualidade quando novo', 'Sujeira'), Papel de cópia = c('Qualidade do papel', 'Foto sensibilidade'), Condições ambientais = c('Mesa suja', 'Mãos sujas'), Papel original = c('Dobras', 'Nitidez', 'Transparência'), Copiadora = c('Velocidade', 'Condições do enrolamento', 'Potência da lâmpada')), title=' Gestão da Qualidade \n Diagrama de Causas e Efeito', effect = 'Cópia de má qualidade', cex = c(0.7, 0.5, 1)) O plot saiu do jeito que você deseja. Atenciosamente, T.T. Em 13/08/2014, à(s) 07:27, Edson Lira edinhoes...@yahoo.com.br escreveu: Caros liteiros, estou tentando usar a rotina abaixo: library(qcc) # Ativar o pacote qcc cause.and.effect(cause= list(Método.de.foto.cópia= c('Colocação do Original', 'Tempo de secagem'), Liquido= c('Qualidade', 'Qualidade quando novo', 'Sujeira'), Papel.de.cópia= c('Qualidade do papel', 'Foto sensibilidade'), Condições.ambientais= c('Mesa suja', 'Mãos sujas'), Papel.original= c('Dobras', 'Nitidez', 'Transparência'), Copiadora= c('Velocidade', 'Condições do enrolamento', 'Potência da lâmpada')), title=' Gestão da Qualidade \n Diagrama de Causas e Efeito', effect= 'Cópia de má qualidade', cex= c(0.7, 0.5, 1)) Como faço para excluir os pontos que separaram as palavras no gráfico, por exemplo como faço para aparecer Método de foto cópia e não Método.de.foto.cópia. A rotina exige o ponto entre as palavras das causas principais. Bom dia a todos. [ ]'s. Edson Lira Estatístico Manaus-Amazonas___ R-br mailing list R-br@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.___ R-br mailing list R-br@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.