- R-br - DINF UFPR

Como escolher o melhor conjunto de dados imputador no MICE
by Elias Carvalho 21 Abr '17

21 Abr '17

Usando o pacote mice fiz o seguinte teste (imputação de dados): imp <- mice(nhanes) Isso gerou 5 conjuntos de dados imputados: imp$imp$bmi Eu gostaria de entender como escolher o melhor dataset imputado Por exemplo, para bmi (acima) qual das 5 colunas será a melhor escolha? Eu achei vários exemplos usando regressão logistica, mas eu não vou rodar uma regressão logistica, vou usar os dados para criar uma rede bayesiana. Alguém poderia me ajudar por favor? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

Re: [R-br] Plotar linhas em gráfico 3d no R
by sznelwar＠uol.com.br 21 Abr '17

21 Abr '17

1 0

Plotar linhas em gráfico 3d no R
by Laura Teixeira Rocha 20 Abr '17

20 Abr '17

Possuo um gráfico de dispersão 3d realizado com duas funções no R, a função 'plotly' e a ''scatterplot3d'. Os gráficos podem ser reproduzidos através do código: x <- c(0,50,100,150,200,250,300,350,400,450) y <- c(0,50,100,150,200,250,300,350,400,450) z <- c(1,2,1,1,2,1,2,1,2,1) plot_ly(x=x,y=y,z=z) scatterplot3d(x, y, z) Se possível, para ambas as funções, gostaria de saber como posso plotar, no plano deste gráfico, duas retas cortando o eixo x e duas retas cortando o eixo y, de modo que essas retas se encontrem no plano, formando uma região de intersecção entre as duas. Além disso, gostaria de saber se é possível preencher/ressaltar essa área de encontro com alguma cor. Muito obrigada! Laura.

1 0

Interface para usuário realizar chamada de função em R
by Paulo Nogueira Starzynski 20 Abr '17

20 Abr '17

Olá colegas de lista. Gostaria da opinião de vocês com relação a algo que envolve R, mas não apenas ele. Desenvolvi em R um código que automatiza alguns processos. Todo o trabalho foi encapsulado dentro de uma função principal (aqui chamada de *realiza.consulta*), que recebe um parâmetro com as configurações informando o que deve ser feito. Criei algo parecido para exemplificar: # # "Painel de Controle" do Usuário # configuracoes <- list( idUser <- 'usuario', idPassword <- 'senha123', consultas <- list(c("info1", "todos"), c("info2", "apenas alguns"), c("info3", NA), c("info4", "opcao1", "opcao2") ) ) # # Excucao do PROGRAMA PRINCIPAL source("/codigos.R") realiza.consulta(configuracoes) Como vocês podem ver, as configurações ficam armazenadas em um objeto do tipo lista, que é passado para a função principal. É de responsabilidade do usuário criar esse objeto e chamar a função realiza.consulta. A questão é que eu gostaria de desenvolver uma "máscara" para esse processo, onde os usuários informariam tais configurações e clicariam em um botão para iniciar o procedimento, sem ter de abrir o R ou RStudio. Vocês já fizeram algo semelhante para poder dividir o conhecimento? Abraços, Paulo

3 3

Experimento em faixa com medidas repetidas no tempo
by ROSIANA RODRIGUES ALVES 20 Abr '17

20 Abr '17

Boa noite, Tenho que analisar um experimento e estou com duvidas como analisar. Os experimento teve o objetivo de avaliar a produção e qualidade de matéria seca de uma pastagem ao longo de um ano, que foi conduzida em piquetes com duas alturas de pastejo e três lâminas de irrigação. Dois fatores (Altura de pastejo e Laminas de irrigação). Devido a dificuldade de casualização das laminas de irrigação o experimento foi conduzido em faixas. Já analisei e experimentos em faixa, e experimentos com medida repetida mas os dois ao mesmo tempo, nunca analisei, peço ajuda de alguém que já fez esse tipo de analise para me passar umas dicas. Abraços Rosiana Rodrig ues Alves Estat ística e Experi mentação Agropecuária Em brapa Pesca e Aquicu ltura Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Palmas/TO +55 (63) 3229-7805/7800 | Skype: cnpasa.rosiana http://cnpasa.sede.embrapa.br | twitter.com/embrapa Confira também: www.facebook.com/agrosustentavel ----- Mensagem original ----- De: r-br-request(a)listas.c3sl.ufpr.br Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quarta-feira, 19 de abril de 2017 12:00:17 Assunto: Digest R-br, volume 76, assunto 13 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=https%3a%… ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: Banco de Dados - Linhas e Colunas (Luis G. S. e Silva) 2. Como escolher o melhor conjunto de dados imputador no MICE (Elias Carvalho) ---------------------------------------------------------------------- Message: 1 Date: Tue, 18 Apr 2017 20:29:01 +0000 From: "Luis G. S. e Silva" <lgsilvaesilva(a)gmail.com> To: jose gomes de souza filho Gomes <zegomesufg(a)gmail.com>, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Banco de Dados - Linhas e Colunas Message-ID: <CALruoAEa75fUUwPbBtwmm+cxcZAk-YCXdU5m1uecviqUKLdsdw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" require(reshape2) dados <- read.csv2('dadosg.csv', encoding = 'latin1') dcast(cidade+parcela2+parcela3~item+ano, value.var = 'quantidade', data = dados) On Tue, Apr 18, 2017 at 4:27 PM jose gomes de souza filho Gomes via R-br < r-br(a)listas.c3sl.ufpr.br> wrote: > Boa tarde, > > Tem o seguinte banco de dados > > base de dados: > https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=https%3a%… > > *cidade* > > *parcela2* > > *parcela3* > > *ano* > > *item* > > *quantidade* > > *GOIÂNIA* > > *Ana* > > *serv* > > *2011* > > *baixados* > > *82331* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *2011* > > *baixados* > > *15733* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *2011* > > *baixados* > > *9989* > > *GOIÂNIA* > > *Ana* > > *serv* > > *2011* > > *Pendentes* > > *7282* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *2011* > > *Pendentes* > > *6819* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *2011* > > *Pendentes* > > *5420* > > *GOIÂNIA* > > *Ana* > > *serv* > > *2011* > > *novos* > > *4635* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *2011* > > *novos* > > *4270* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *2011* > > *novos* > > *4125* > > *GOIÂNIA* > > *Ana* > > *serv* > > *2012* > > *baixados* > > *82332* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *2012* > > *baixados* > > *15734* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *2012* > > *baixados* > > *9990* > > *GOIÂNIA* > > *Ana* > > *serv* > > *2012* > > *Pendentes* > > *7283* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *2012* > > *Pendentes* > > *6820* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *2012* > > *Pendentes* > > *5421* > > *GOIÂNIA* > > *Ana* > > *serv* > > *2012* > > *novos* > > *4636* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *2012* > > *novos* > > *4271* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *2012* > > *novos* > > *4126* > > > > > Gostaria que ficasse assim: > > > > *cidade* > > *parcela2* > > *parcela3* > > *baixados_2011* > > *baixados_2012* > > *Pendentes_2011* > > *Pendentes_2012* > > *novos_2011* > > *novos_2012* > > *GOIÂNIA* > > *Ana* > > *serv* > > *82331* > > *82332* > > *7282* > > *7283* > > *4635* > > *4636* > > *GOIÂNIA* > > *Paulo* > > *terc* > > *15733* > > *15734* > > *6819* > > *6820* > > *4636* > > *4271* > > *ANÁPOLIS* > > *Roberto* > > *serv* > > *9989* > > *9990* > > *5420* > > *5421* > > *4125* > > *4126* > > > > Obrigado. > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=https%3a%… > Leia o guia de postagem (https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=http%3a%2… e forneça > código mínimo reproduzível. -- Luís G. Silva e Silva -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=http%3a%2… ------------------------------ Message: 2 Date: Wed, 19 Apr 2017 00:34:50 -0300 From: Elias Carvalho <ecacarva(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Como escolher o melhor conjunto de dados imputador no MICE Message-ID: <CAH8E6UuWYW_xG_OKU0ft=UNWOFODNWX-c8e0kzxbi_qzV43xSw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Usando o pacote mice fiz o seguinte teste (imputação de dados): imp <- mice(nhanes) Isso gerou 5 conjuntos de dados imputados: imp$imp$bmi Eu gostaria de entender como escolher o melhor dataset imputado Por exemplo, para bmi (acima) qual das 5 colunas será a melhor escolha? Eu achei vários exemplos usando regressão logistica, mas eu não vou rodar uma regressão logistica, vou usar os dados para criar uma rede bayesiana. Alguém poderia me ajudar por favor? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"* -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=http%3a%2… -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: Captura de tela de 2017-04-19 00-30-32.png Tipo: image/png Tamanho: 14828 bytes Descrição: não disponível URL: <https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=http%3a%2… ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://imsva91-ctp.trendmicro.com:443/wis/clicktime/v1/query?url=https%3a%… ------------------------------ Fim da Digest R-br, volume 76, assunto 13 ***************************************** ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

2 1

Inclusão de estrutura de correlação função glmer do pacote lme4
by Cleber Iack 19 Abr '17

19 Abr '17

Prezados Eu estou utilizando a função glmer do pacote lme4 do R. A minha resposta é binomial. No Pacote bild, eu posso utilizar a opção Mc1 ou MC2 para estruturas de correlação, tendo em vista os meus dados serem longitudinais. Há alguma maneira de utilizando a funçao glmer, poder incluir esta estrutura? Obrigado

2 2

Banco de Dados - Linhas e Colunas
by jose gomes de souza filho Gomes 18 Abr '17

18 Abr '17

Boa tarde, Tem o seguinte banco de dados base de dados: https://www.datafilehost.com/d/6c117eac *cidade* *parcela2* *parcela3* *ano* *item* *quantidade* *GOIÂNIA* *Ana* *serv* *2011* *baixados* *82331* *GOIÂNIA* *Paulo* *terc* *2011* *baixados* *15733* *ANÁPOLIS* *Roberto* *serv* *2011* *baixados* *9989* *GOIÂNIA* *Ana* *serv* *2011* *Pendentes* *7282* *GOIÂNIA* *Paulo* *terc* *2011* *Pendentes* *6819* *ANÁPOLIS* *Roberto* *serv* *2011* *Pendentes* *5420* *GOIÂNIA* *Ana* *serv* *2011* *novos* *4635* *GOIÂNIA* *Paulo* *terc* *2011* *novos* *4270* *ANÁPOLIS* *Roberto* *serv* *2011* *novos* *4125* *GOIÂNIA* *Ana* *serv* *2012* *baixados* *82332* *GOIÂNIA* *Paulo* *terc* *2012* *baixados* *15734* *ANÁPOLIS* *Roberto* *serv* *2012* *baixados* *9990* *GOIÂNIA* *Ana* *serv* *2012* *Pendentes* *7283* *GOIÂNIA* *Paulo* *terc* *2012* *Pendentes* *6820* *ANÁPOLIS* *Roberto* *serv* *2012* *Pendentes* *5421* *GOIÂNIA* *Ana* *serv* *2012* *novos* *4636* *GOIÂNIA* *Paulo* *terc* *2012* *novos* *4271* *ANÁPOLIS* *Roberto* *serv* *2012* *novos* *4126* Gostaria que ficasse assim: *cidade* *parcela2* *parcela3* *baixados_2011* *baixados_2012* *Pendentes_2011* *Pendentes_2012* *novos_2011* *novos_2012* *GOIÂNIA* *Ana* *serv* *82331* *82332* *7282* *7283* *4635* *4636* *GOIÂNIA* *Paulo* *terc* *15733* *15734* *6819* *6820* *4636* *4271* *ANÁPOLIS* *Roberto* *serv* *9989* *9990* *5420* *5421* *4125* *4126* Obrigado.

2 1

Re: [R-br] Suavização do gráfico
by sznelwar＠uol.com.br 16 Abr '17

16 Abr '17

1 0

Re: [R-br] Substituiçao de missing
by suzabio 16 Abr '17

16 Abr '17

Olá Leonardo, Obrigada! Deu certo, eu não estava conseguindo ver um jeito de substituir esses NAs, tudo que eu pensei deu errado.Um abraço e boa Páscoa! Suzana Enviado do meu dispositivo Samsung -------- Mensagem original -------- De : Leonardo Aguirre <leoaguirre32(a)gmail.com> Data: 16/04/2017 09:28 (GMT-03:00) Para: suzana de souza rodrigues <suzabio(a)yahoo.com.br>, "a lista Brasileira oficial de discussão do programa R." <r-br(a)listas.c3sl.ufpr.br> Assunto: Re: [R-br] Substituiçao de missing Como vc já calculou as medias, uma opção é jogar esse campo na tabela usando a função merge() e depois utilizar a função ifelse() para substituir os NA. Abs, Leonado Em 16 de abr de 2017 9:11 AM, "suzana de souza rodrigues via R-br" <r-br(a)listas.c3sl.ufpr.br> escreveu: Bom dia,Tenho um banco de dados: regiao semana notificados confirmados sexo local1 1 1 82 45 M A2 1 2 105 32 F A3 1 3 134 71 M A4 1 4 80 44 M A5 1 5 102 32 M A6 1 6 108 46 F A Formado por 50 regiões, a variável confirmados apresenta missings (NA).Preciso substituir os missings pela média de casos confirmados por região.Calculei a média por região por esse comando, sem considerar os missings: media <- bancos%>% group_by(regiao)%>% summarise(mean(na.exclude(confirmados)))as.data.frame(media) Com esse resultado: regiao mean(na.exclude(confirmados))1 1 50.151522 2 51.526043 3 69.907694 4 45.185575 5 61.658166 6 57.796887 7 47.305268 8 69.717959 9 68.4070410 10 62.2552111 11 58.00000 Só que não estou conseguindo trocar os missings por essas médias e preciso de ajuda nisso.Desde já agradeço,Suzana _______________________________________________ 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

Substituiçao de missing
by suzana de souza rodrigues 16 Abr '17

16 Abr '17

Bom dia,Tenho um banco de dados: regiao semana notificados confirmados sexo local1 1 1 82 45 M A2 1 2 105 32 F A3 1 3 134 71 M A4 1 4 80 44 M A5 1 5 102 32 M A6 1 6 108 46 F A Formado por 50 regiões, a variável confirmados apresenta missings (NA).Preciso substituir os missings pela média de casos confirmados por região.Calculei a média por região por esse comando, sem considerar os missings: media <- bancos%>% group_by(regiao)%>% summarise(mean(na.exclude(confirmados)))as.data.frame(media) Com esse resultado: regiao mean(na.exclude(confirmados))1 1 50.151522 2 51.526043 3 69.907694 4 45.185575 5 61.658166 6 57.796887 7 47.305268 8 69.717959 9 68.4070410 10 62.2552111 11 58.00000 Só que não estou conseguindo trocar os missings por essas médias e preciso de ajuda nisso.Desde já agradeço,Suzana

2 1