Abril 2017 - R-br - DINF UFPR

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 17 Abr '17

17 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

Maximizar plot window no RStudio
by Elias Carvalho 13 Abr '17

13 Abr '17

Bom dia pessoal Alguém sabe como posso maximizar a janela plot do RSstudio via linha de comando ? Explicando: Estou fazendo uma rotina para verificação de outliers baseado no boxplot então ao examinar os dados eu gero os gráficos abaixo, no entanto, se a janela estiver pequena demais os gráficos não aparecem e ainda surge uma mensagem Error in plot.new() : figure margins too large in R <http://stackoverflow.com/questions/12766166/error-in-plot-new-figure-margin…> : -- *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"*

3 4

Combinar multiplos datasets - Inputação Múltipla com MICE
by Elias Carvalho 11 Abr '17

11 Abr '17

*Boa noite pessoal* *Eu preciso tratar dados missing do meu dataset, li bastante e o MICE me pareceu uma boa ferramente para isso, então p**rossegui da seguinte maneira:* mice(df.factor.variable, m=5, maxit = 50, method = 'polyreg', seed = 500) *e depois para juntar os 5 datasets usei:* complete.factor.data <- complete(factor.variable.imputed, 'long') *Como tenho 1999 registros no meu dataset e foram gerados 5 conjuntos de dados imputados, o resultado foi de 9995 observações (1999 * 5).* *Alguém sabe me dizer se o procedimento está correto ou há uma maneira de unir os 5 datasets imputados em um só gerando um dataset de 1999 registros sem NAs?* *ps: vou usar esse novo dataset sem missing para aprender uma rede Bayesiana* -- *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"*

1 0

código para construir tabela com dados pareados
by Mauricio Cardeal 09 Abr '17

09 Abr '17

Boa tarde! Por favor, como faço para construir uma tabela como a do modelo do final do script, onde o conteúdo das 4 caselas são contagens do número de pares que combinam melhora e não melhora de 2 tratamentos aplicados simultaneamente no mesmo individuo (perna direita e perna esquerda), por exemplo ? # CMR: par <- c(1,1,2,2,3,3,4,4,5,5,6,6,7,7) # refere-se ao número do individuo trata <- c(1,0,1,0,1,0,1,0,1,0,1,0,1,0) # tratamento 1 (A) ou 0 (B) melhora <- c(1,0,1,0,1,0,0,1,0,1,0,0,1,1) # resultado do tratamento 1 (melhora) 0 (não melhora) tratac <- factor(trata) levels(tratac)<-list("A"=1,"B"=0 ) melhorac <- factor(melhora) levels(melhorac)<-list("+"=1,"-"=0) data.frame(par,tratac,melhorac) par tratac melhorac 1 1 A + 2 1 B - 3 2 A + 4 2 B - 5 3 A + 6 3 B - 7 4 A - 8 4 B + 9 5 A - 10 5 B + 11 6 A - 12 6 B - 13 7 A + 14 7 B + # COMO FAÇO PARA GERAR ESSA TABELA A PARTIR DOS DADOS ACIMA? # tratamento B # melhora melhora # + - #tratamento A melhora + 1 3 # melhora - 2 1 Obrigado, Maurício Cardeal UFBA

2 3

Suavização do gráfico
by Mateus Dias Nunes 09 Abr '17

09 Abr '17

Bom dia, Eu estou com um problema para suavizar meu gráfico. Eu tenho dados horários e desejo plotar individualmente cada data (no período das 06h as 18h). Entretanto o gráfico tem ficado com muito ruído devido as observações serem horárias e existirem falhas (NA) nas observações. E por isso eu gostaria de suavizar o gráfico com as médias horárias, porém não estou conseguindo implementar isso. Espero que alguém possa me ajudar a resolver. dado utilizado: https://www.dropbox.com/s/b80zzzeaesnh8qp/dado.txt?dl=0 library(openair) mydata<-read.table("dado.txt", h = T, sep = ",") mydata$direta <- as.numeric(as.character(mydata$direta)) mydata$date <- as.POSIXct(strptime(mydata$date, format = "%d/%m/%Y %H:%M:%S",tz = "Etc/GMT-1")) means <- aggregate(mydata["direta"], format(mydata["date"],"%d-%m-%Y %H:%M:%S"), mean, na.rm = TRUE) means[ ,"date"] <- seq(min(mydata[ , "date"]), max(mydata[ , "date"]), length = nrow(means)) Figura1<-plot(mydata$date[2:721], mydata$direta[2:721], type = "l", col="red",xlab = "minutes",main = "01 de Janeiro de 2007", ylab = "Direct Solar Radiation (W/m²)") Figura1 dev.copy(device=jpeg,filename="Brasilia_01dejaneiro2007.jpg", width = 800, height = 500) dev.off () desde já, agradeço. _______________________________________________________________________________________ MATEUS DIAS NUNES Postgraduate Program in Meteorology - PGMET - CPTEC / INPE Center for Weather Forecasting and Climate Studies / National Institute for Space Research Cachoeira Paulista - SP - Brasil FONE: +55 (53) 98112-5154 SKYPE: omsndomateus(a)hotmail.com ____________________________________________ ___________________________________________ <https://mailtrack.io/> Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateus…> <#>

2 1

Duvida em estatística no R
by Tiago G. Santos 08 Abr '17

08 Abr '17

Boa tarde, Alguém sabe como faço para apresentar os resultados dos testes tukey ou dunnet como as letras que indicam significância estatística ..... -- Tiago Gonçalves Santos http://lattes.cnpq.br/7645251131511453 Doutorando em Fìsica pelo Instituto Tecnológico de Aeronáutica - ITA Laboratório de Caracterização de Dispositivos Semicondutores Instituto de Estudos Avançados - IEAv Trevo Coronel Aviador José Alberto Albano do Amarante, no. 1 12.228-001 São José dos Campos – SP – Brasil

2 1