Maio 2017 - R-br - DINF UFPR

R-Ladies Rio
by Celina Rebello 03 Mai '17

03 Mai '17

Pessoal, Aproveitando para divulgar: Qui, 18 de mai, 19:00 Faculdades Senac Rio Rua Santa Luzia, 735 - 4° andar - sala 402 - Próximo ao Metrô da Cinelândia., Rio de Janeiro Evento é gratuito, temos poucas vagas. Inscrições em: https://goo.gl/AOsCQp Quem estiver pelo Rio, vale muito ir.

1 0

Re: [R-br] Ajuda para elaboração de boxplot em painel
by Gilenio Borges Fernandes 03 Mai '17

03 Mai '17

Prezad@s Estou com dificuldade de elaborar um boxplot por grupo (ou painel), no formato do que está no arquivo pdf em anexo a esse E-mail. Segue o CMR do banco de dados (Documentação do SAS) O arquivo pdf pode ser acessado em https://www.dropbox.com/s/5q65geuq7e9m5ey/Dados%26grafico-DocSAS.pdf?dl=0 Desde já agradeço Parts=structure(list(Machine = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("A386", "A455", "C334"), class = "factor"), Sample = 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, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, 20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 23L, 23L, 23L, 23L, 24L, 24L, 24L, 24L), Diam = c(4.32, 4.55, 4.16, 4.44, 4.49, 4.3, 4.52, 4.61, 4.44, 4.32, 4.25, 4.5, 4.55, 4.15, 4.42, 4.49, 4.21, 4.3, 4.29, 4.63, 4.56, 4.61, 4.29, 4.56, 4.63, 4.3, 4.41, 4.58, 4.38, 4.65, 4.43, 4.44, 4.12, 4.49, 4.3, 4.36, 4.45, 4.56, 4.38, 4.51, 4.62, 4.67, 4.7, 4.58, 4.33, 4.23, 4.34, 4.58, 4.29, 4.38, 4.28, 4.41, 4.15, 4.35, 4.28, 4.23, 4.21, 4.3, 4.32, 4.38, 4.16, 4.28, 4.31, 4.59, 4.14, 4.18, 4.08, 4.21, 4.51, 4.2, 4.28, 4.19, 4.1, 4.33, 4.37, 4.47, 3.99, 4.09, 4.47, 4.25, 4.24, 4.54, 4.43, 4.38, 4.23, 4.48, 4.31, 4.57, 4.27, 4.4, 4.32, 4.56, 4.7, 4.65, 4.49, 4.38)), .Names = c("Machine", "Sample", "Diam"), class = "data.frame", row.names = c(NA, -96L)) Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860 Em 3 de maio de 2017 15:58, Gilenio Borges Fernandes <gilenio(a)ufba.br> escreveu: > Prezad@s > > Estou com dificuldade de elaborar um boxplot por grupo (ou painel), no > formato do que está no arquivo pdf em anexo a esse E-mail. Segue o CMR do > banco de dados (Documentação do SAS) > > Desde já agradeço > > > Parts=structure(list(Machine = structure(c(1L, 1L, 1L, 1L, 1L, 1L, > > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, > > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, > > 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, > > 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, > > 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("A386", "A455", > > "C334"), class = "factor"), Sample = 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, 9L, 9L, 9L, 9L, 10L, > > 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, > > 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, > > 16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, > > 20L, 20L, 20L, 20L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 23L, > > 23L, 23L, 23L, 24L, 24L, 24L, 24L), Diam = c(4.32, 4.55, 4.16, > > 4.44, 4.49, 4.3, 4.52, 4.61, 4.44, 4.32, 4.25, 4.5, 4.55, 4.15, > > 4.42, 4.49, 4.21, 4.3, 4.29, 4.63, 4.56, 4.61, 4.29, 4.56, 4.63, > > 4.3, 4.41, 4.58, 4.38, 4.65, 4.43, 4.44, 4.12, 4.49, 4.3, 4.36, > > 4.45, 4.56, 4.38, 4.51, 4.62, 4.67, 4.7, 4.58, 4.33, 4.23, 4.34, > > 4.58, 4.29, 4.38, 4.28, 4.41, 4.15, 4.35, 4.28, 4.23, 4.21, 4.3, > > 4.32, 4.38, 4.16, 4.28, 4.31, 4.59, 4.14, 4.18, 4.08, 4.21, 4.51, > > 4.2, 4.28, 4.19, 4.1, 4.33, 4.37, 4.47, 3.99, 4.09, 4.47, 4.25, > > 4.24, 4.54, 4.43, 4.38, 4.23, 4.48, 4.31, 4.57, 4.27, 4.4, 4.32, > > 4.56, 4.7, 4.65, 4.49, 4.38)), .Names = c("Machine", "Sample", > > "Diam"), class = "data.frame", row.names = c(NA, -96L)) > > > > Gilênio Borges Fernandes > Professor Associado IV (Aposentado) > Professor Adjunto A (Substituto) > Universidade Federal da Bahia > Instituto de Matemática e Estatística > Departamento de Estatística > Av. Adhemar de Barros, s/n – Ondina. > 40.170-110 - Salvador - BA, Brasil > Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 > Skype: gilenio.fernandes > Lattes: http://lattes.cnpq.br/6764860618464860 >

2 1

geoR > Krige.bayes com trend e como obter DIC
by Wagner Wolff 03 Mai '17

03 Mai '17

Olá pessoal da lista! Estou com dificuldade de implementar na função krige.bayes() o atributo trend. Primeiramente, por uma abordagem frequentista escolhi um modelo de trend, quando tento usar o mesmo na função krige.bayes() ocorre alguns erros com questões matriciais. Alguem tem alguma dica para corrigir isso? Outra questão cé omo obtendo o DIC para os modelos que testarei na função krige.bayes(). Por exemplo, modelo com trend e componentes espaciais, modelo só com trend e modelo só com os componentes espacias. A figura abaixo mostra um exemplo de um artigo que fez essas comparações utilizando binom.krige.bayes(). Fonte: Mapping, Bayesian Geostatistical Analysis and Spatial Prediction of Lymphatic Filariasis Prevalence in Africa -- *Wagner Wolff, **PhD* "*Luiz de Queiroz**" College of Agriculture,* University of São Paulo Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil Phone: +55 19 982385582 <+55%2019%2098238-5582> http://orcid.org/0000-0003-3426-308X https://github.com/wwolff7 http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1

1 0

Comparações Múltiplas
by Isabela M. 03 Mai '17

03 Mai '17

Boa tarde, Meu nome é Isabela, faço iniciação científica e trabalho com análise multivariada, mais especificamente com a MANOVA. Alguém poderia me ajudar com script do R, para realizar a comparação múltipla multivariada, no caso para 5 variáveis? Desde já agradeço!

5 7

Banco de Dados
by jose gomes de souza filho Gomes 02 Mai '17

02 Mai '17

Boa tarde, Tem o seguinte banco de dados *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.

3 2

Extrair P-value do resultado de uma imputação com MICE
by Elias Carvalho 01 Mai '17

01 Mai '17

Boa noite pessoal Estou testando o pacote mice para fazer imputação múltipla de dados missing, conforme exemplo: library(mice) imp <- mice(nhanes2) fit <- with(imp, lm(chl~age+bmi)) summary(fit) O resultado se resume em: ## summary of imputation 1 : Call: lm(formula = chl ~ age + bmi) Residuals: Min 1Q Median 3Q Max -46.37 -19.23 -11.24 12.66 83.69 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 32.134 53.202 0.604 0.55232 age2 49.058 18.153 2.702 0.01334 * age3 62.837 19.886 3.160 0.00472 ** bmi 4.907 1.795 2.733 0.01246 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 34.89 on 21 degrees of freedom Multiple R-squared: 0.3835, Adjusted R-squared: 0.2954 F-statistic: 4.354 on 3 and 21 DF, p-value: 0.01557 Eu gostaria de extrair o p-value: 0.01557 dessa estrutura, mas não consegui com os métodos convencionais, pois a estrutura parece diferente da estrutura gerada por um lm ou glm. Alguém poderia me ajudar -- *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 5

Avaliação geral do modelo de regressão linear e logística
by Elias Carvalho 01 Mai '17

01 Mai '17

Ola pessoal Eu estou iniciando estudos sobre regressão linear e regressão logística vi diversos materiais e videos sobre os dois métodos, mas vi pouca coisa sobre como avaliar qual o melhor modelo (talvez esteja pesquisando errado). Estou interessado nisso porque estou usando o pacote MICE para imputar dados missing, que por sua vez gera n conjuntos de dados para cada variável e usa regressão linear e logística para nos ajudar a inspecionar cada dado imputado. Sobre regressão linear, pelo que entendi posso avaliar a qualidade geral do modelo, vendo se o p-value é significativo, Erro médio é menor e R-Squared maior ou achar um meio termo entre essas duas medidas, me corrijam por favor se eu estiver errado. Com relação a regressão logística, encontrei materiais diversos, cada um fazendo uma recomendação diferente, o que me deixou confuso de qual a melhor forma de avaliar a qualidade deste modelo. Alguém poderia me dar uma dica e fornecer links com referências atuais e confiáveis sobre o assunto ? -- *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