- R-br - DINF UFPR

Problemas em cortar linhas iguais
by Amiko Bh 20 Fev '17

20 Fev '17

Pessoal: Dado o DF, com o seguinte nomes fictícios: *NOME* *FILIACAO* AENDER XXXXX AMORIM APARECIDA XXXXX AMORIM AGAPE XXXXX SILVA MARISA IRENEIA DE SOUZA SILVA AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA AGNALDO XXXXX SANTOS BORGES MIRTES XXXXX LUCIANO AGUIDA XXXXX LISBOA DIVA XXXXX LISBOA AGUINALDO XXXXX PAIXAO ROSANE XXXXX PAIXAO AILTON XXXXX DOS SANTOS PATRICIA XXXXX ARAUJO AILTON XXXXX ESTEVES MARIA XXXXX ESTEVES AILTON XXXXX JUNIOR FLAVIA XXXXX MARTINS Uso o seguinte para cortar os nomes repetidos: banco_sem_nome_repetido <- df[!duplicated(data.frame(df$NOME, df$FILIACAO)),] Meu problema: o script resulta o seguinte df: *NOME* *FILIACAO* AENDER XXXXX AMORIM APARECIDA XXXXX AMORIM AGAPE XXXXX SILVA MARISA IRENEIA DE SOUZA SILVA AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA AGNALDO XXXXX SANTOS BORGES MIRTES XXXXX LUCIANO AGUIDA XXXXX LISBOA DIVA XXXXX LISBOA AGUINALDO XXXXX PAIXAO ROSANE XXXXX PAIXAO AILTON XXXXX DOS SANTOS PATRICIA XXXXX ARAUJO AILTON XXXXX ESTEVES MARIA XXXXX ESTEVES AILTON XXXXX JUNIOR FLAVIA XXXXX MARTINS Observe-se que: AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA AGATIL XXXXX OLIVEIRA VANESSA XXXXX DA COSTA Não é cortada a linha porque após o OLIVEIRA tem um espaço. O banco é enorme e tem vários casos iguais ao citado. Alguma solução?

3 2

Re: [R-br] modelo misto aninhado
by Fernando Antonio de souza 19 Fev '17

19 Fev '17

Mande um CMR (código mínimo reproduzível) para que possamos ajudar. Att Em 19/02/2017 09:25, "Adriele Giaretta Biase via R-br" < r-br(a)listas.c3sl.ufpr.br> escreveu: Bom dia! como faço para estimar um modelo estatístico no R, que que tenho três fatores aninhados da seguinte forma: método, confinamento, gênero (sendo os animais a unidade experimental). Sendo os confinamentos aleatórios. E os seus desdobramentos dentro de cada nível dos fatores? -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

modelo misto aninhado
by Adriele Giaretta Biase 19 Fev '17

19 Fev '17

Bom dia! como faço para estimar um modelo estatístico no R, que que tenho três fatores aninhados da seguinte forma: método, confinamento, gênero (sendo os animais a unidade experimental). Sendo os confinamentos aleatórios. E os seus desdobramentos dentro de cada nível dos fatores? -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

1 0

anova em delineamento em dic desbalanceado
by Adriele Giaretta Biase 18 Fev '17

18 Fev '17

Pessoal, qual é a melhor função para testar um conjunto de dados desbalanceados, com fator qualitativo, ao acaso (delineamento inteiramente ao acaso), no R? Como usar o tukey para esse caso? --

2 1

pacote ribge
by Telmo Brugnara 16 Fev '17

16 Fev '17

Pessoal, Criei um pacote, ribge <https://github.com/tbrugz/ribge>, para consulta a algumas bases do IBGE. Funcionalidades/consultas atuais: - populacao dos municipios: 2000-2016 - pip dos municipíos - area dos municípios - indices de preços: INPC, IPCA, SINAPI, deflator PIB - consulta às series de <http://seriesestatisticas.ibge.gov.br/> - mapeamento dos códigos de municipios IBGE-TSE Espero que possa ser útil a alguns por aqui. Críticas e sugestões (e pull requests) são bem-vindos ;) Att, Telmo

1 0

Fatorial incompleto
by Fernando Antonio de souza 14 Fev '17

14 Fev '17

Caros colegas Possuo um banco de dados de um experimento em que tenho os fatores experimentais a linhagem (LG) e Idade (faixa de idade). acontece que não possuo todas as combinações dos níveis dos fatores, como vocês podem ver na tabela abaixo. Isso seria um fatorial incompleto. Eu queria aproveitar os dados que tenho das linhagens na idade 1 (L03,L415,L65) e compará-las com as demais idades e linhagens. Como eu poderia proceder? DADOS<-read.csv("https://www.dropbox.com/s/yvd9cgz55462yeo/CAMILA.csv?raw=1 ",head=TRUE) DADOS<-transform(DADOS,LG=factor(LG),IDADE=factor(IDADE)) ℝ> table(CAMILA[,c(1,4)]) IDADE LG 1 2 3 L02 0 415 1258 L03 10 771 1052 L27 0 98 184 L415 31 181 1038 L62 0 163 570 L65 8 389 1226 -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

4 5

Fóruns e grupos
by Roberto Guevara 14 Fev '17

14 Fev '17

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

3 2

Benchmarking the unidades de saúde
by Pedro Emmanuel Alvarenga Americano do Brasil 14 Fev '17

14 Fev '17

Olá amigos de R, Uma colega e eu escrevemos um pacote que tem como missão principal fazer analises de benchmarking (comparações quanto a eficiência de previnir morte e reduzir tempo de internação) de unidades de saúde. As principais fontes são relacionadas a literatura de unidades de terapia intensiva. Estou escrevendo pra pedir aos amigos de R, baixarem, testarem e fazer algum comentário sobre possíveis bugs, aperfeiçoamentos ou novas funcionalidades, pois acho que o pacote já poderia ser submetido ao CRAN. *install.packages("ems", repos="http://R-Forge.R-project.org <http://R-Forge.R-project.org>")* As vezes a instalção do r-forge não funciona pq fazemos a sincronização de arquivos não finalizados, para copia de segurança e que impede r-forge de construir o pacote. A Version: *0.3.1* | Last change: *2017-02-10 21:32:52+01* | Rev.: *30. É a última que certamente funciona.* Nós estamos nos programando estudar outras coisas relacionadas com benchmarking (e.g. análise de envelopes), mas na área da saúde essa literatura é muito pobre. Então quem puder indicar textos ou livros de formas adicionais de fazer benchmarking eu agradeço. Abraço forte a todos, Pedro Brasil

1 0

Distribuição para regressão de resposta binária
by Marcos Bissoli 14 Fev '17

14 Fev '17

Prezados, De antemão peço desculpas se desvio o tópico da lista. Mas creio que o tema da mensagem é minimamente transversal aos aqui tratados. Tenho uma variável resposta binária. Como a frequência da resposta é alta (38,11%), teóricos da Estatística aplicada à Epidemiologia sugerem que não seja usada uma regressão logística. Neste caso (de alta prevalência do desfecho), a primeira opção deveria ser uma log-binomial. Mas (e isso não é raro de ocorrer), minha log-binomial não apresentou convergência. Quando não há convergência, os teóricos sugerem uma regressão de Poisson com variância robusta. Entretanto, como meus dados sugerem subdispersão, optei por um modelo de quasi-poisson. Isso já deu certo em outras análises que fiz para terceiros. Inclusive, tenho conseguido adaptar a variância robusta ao modelo de quasi-poisson. Mas justamente agora, com os dados de minha tese... O diagnóstico visual está, ao meu ver, péssimo, para ajuste. A imagem anexa é do modelo de quasi-poisson. Mas experimentei todos os acima citados (logística e Poisson) e o gráfico não diferiu muito. [image: Imagem inline 1] A dúvida é... Há alguma outra alternativa de técnica de regressão que eu poderia tentar? Minhas variáveis explicativas são diversas, em quantidade e tipo (há contínuas, ordinais e binárias). Ou será (embora eu ache pouco provável) que este gráfico não significa um grande incômodo? Fiz o teste de qui-quadrado da deviance residual e estranhamente o valor p está resultando em 1, tanto para Poisson quanto para quasi-Poisson. Um outro fato estranho é o pseudo R² de Nagelkerke ter acusado 20%: todas as outras minhas variáveis resposta não passaram de 12%. Não sei se é correto (consultei bibliografia que sugeria isso para a regressão logística), mas apliquei um teste de Hosmer e Lemeshow e ele acusou um bom ajuste do modelo, também (p = 0,2718). Até uma curva de ROC eu fiz e a área está grande no gráfico (mais uma técnica que não sei se deve ser aplicada além da regressão logística,). Seguem alguns resultados, caso possa ajudar em algo. Desde já agradeço qualquer comentário. E reforço minhas desculpas caso eu tenha desviado do tópico além do esperado, e desde já acato qualquer negativa em prosseguir o debate. Nesse caso, se possível, aceitaria sugestões de boas listas para debates nesse nível onde eu pudesse me inscrever. Há braços, Marcos Bissoli Faculdade de Nutrição Unifal-MG > Mod1 <- glm(Tabagismo~.,data = TabModelagem,family = quasipoisson)> summary(Mod1) Call: glm(formula = Tabagismo ~ ., family = quasipoisson, data = TabModelagem) Deviance Residuals: Min 1Q Median 3Q Max -1.4867 -0.7821 -0.5889 0.5349 1.6624 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -1.245e+00 8.738e-01 -1.424 0.154644 factor.SexoDic.1 5.800e-01 8.273e-02 7.011 4.11e-12 *** factor.Branca.1 -8.332e-01 7.836e-01 -1.063 0.287863 factor.Negra.1 -8.210e-01 7.987e-01 -1.028 0.304185 factor.Parda.1 -9.009e-01 7.863e-01 -1.146 0.252163 factor.Amarela.1 -1.089e+00 8.481e-01 -1.284 0.199466 factor.SemReligiao.1 -9.670e-02 1.888e-01 -0.512 0.608566 factor.Catolica.1 -4.813e-01 1.862e-01 -2.585 0.009863 ** factor.Espirita.1 -1.235e-01 2.181e-01 -0.566 0.571230 factor.Evangelica.1 -9.177e-01 2.429e-01 -3.779 0.000166 *** factor.AfroBrasileira.1 6.068e-01 4.303e-01 1.410 0.158794 factor.Turno.1 1.534e-03 1.034e-01 0.015 0.988169 factor.Aposentado.1 -4.516e-02 1.055e-01 -0.428 0.668597 factor.OcupaEstDiApenasDesemp.1 7.249e-02 1.411e-01 0.514 0.607474 factor.ComFamilia.1 -4.323e-01 2.128e-01 -2.031 0.042444 * factor.ComOutParentes.1 -5.029e-01 3.517e-01 -1.430 0.153011 factor.Republica.1 8.985e-03 1.959e-01 0.046 0.963429 factor.Sozinho.1 -2.475e-01 2.236e-01 -1.107 0.268673 factor.Pensao.1 -8.439e-01 4.000e-01 -2.110 0.035106 * factor.OutroMoradia.1 -5.262e-01 3.353e-01 -1.569 0.116880 factor.RU.1 -1.937e-01 1.059e-01 -1.830 0.067589 . factor.praec4.1 -1.583e-01 2.666e-01 -0.594 0.552951 IdadeA 3.787e-02 9.381e-03 4.037 5.79e-05 *** escola 8.576e-02 3.441e-02 2.492 0.012836 * RendaPC 4.045e-05 1.313e-05 3.080 0.002119 ** Dist 2.605e-05 1.296e-04 0.201 0.840689 PraecSoma 2.419e-02 3.086e-02 0.784 0.433427 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for quasipoisson family taken to be 0.6036898) Null deviance: 834.67 on 1135 degrees of freedom Residual deviance: 706.16 on 1109 degrees of freedom AIC: NA Number of Fisher Scoring iterations: 5 -- MARCOS BISSOLI Faculdade de Nutrição Universidade Federal de Alfenas Blog: bocademiamaldita.blogspot.com/ E-mail: mbissoli(a)gmail.com Twitter: #mbissoli Alfenas, Minas Gerais, Brasil *****Pense na Natureza antes de Imprimir***** Divulgue ON-LINE Eu apoio a ENEN "na luta por um Brasil sem fome" "por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan" (para cada povo sua própria língua, para todos os povos o Esperanto) E nunca votarei no PSDB/DEM!

4 7

Acrescentar nova coluna com base em palavras de outra coluna
by Amikobh 11 Fev '17

11 Fev '17

Pessoal: Trabalho um banco de dados com a seguinte variável, dentre muitas outras: *NOME_DA_ESCOLA* CASA ESPERENCA III COLEGIO METODO COLEGIO MUNICIPAL MARCONI COLEGIO MUNICIPAL MARCONI COLEGIO MUNICIPAL MARCONI COLEGIO MUNICIPAL MARCONI COLEGIO SALESIANO COLEGIO SALESIANO COLEGIO SALESIANO ESCOLA ESTADUAL PROFESSORA AMELIA DE CASTRO MONTEIRO ESCOLA ESTADUAL PROFESSORA AMELIA DE CASTRO MONTEIRO ESCOLA MUNICIPAL BELO HORIZONTE ESCOLA MUNICIPAL BELO HORIZONTE ESCOLA ESTADUAL AFONSO PENA ESCOLA ESTADUAL AFONSO PENA ESCOLA ESTADUAL ALISSON PEREIRA GUIMARAES ESCOLA ESTADUAL ALISSON PEREIRA GUIMARAES ESCOLA ESTADUAL ALISSON PEREIRA GUIMARAES 1) Pretendo acrescentar mais 1 variável com base nas palavras "MUNICIPAL" e "ESTADUAL" deste banco para que fique assim: *NOME_DA_ESCOLA* *TIPO_DE_ESCOLA* CASA ESPERENCA III PARTICULAR COLEGIO METODO PARTICULAR COLEGIO MUNICIPAL MARCONI MUNICIPAL COLEGIO MUNICIPAL MARCONI MUNICIPAL COLEGIO MUNICIPAL MARCONI MUNICIPAL COLEGIO MUNICIPAL MARCONI MUNICIPAL COLEGIO SALESIANO PARTICULAR COLEGIO SALESIANO PARTICULAR COLEGIO SALESIANO PARTICULAR ESCOLA ESTADUAL PROFESSORA AMELIA DE CASTRO MONTEIRO ESTADUAL ESCOLA ESTADUAL PROFESSORA AMELIA DE CASTRO MONTEIRO ESTADUAL ESCOLA MUNICIPAL BELO HORIZONTE MUNICIPAL ESCOLA MUNICIPAL BELO HORIZONTE MUNICIPAL ESCOLA ESTADUAL AFONSO PENA ESTADUAL ESCOLA ESTADUAL AFONSO PENA ESTADUAL ESCOLA ESTADUAL ALISSON PEREIRA GUIMARAES ESTADUAL ESCOLA ESTADUAL ALISSON PEREIRA GUIMARAES ESTADUAL ESCOLA ESTADUAL ALISSON PEREIRA GUIMARAES ESTADUAL 2) Há como fazê-lo no R. Grande Abraço.

2 1