- R-br - DINF UFPR

Converter character em numeric
by Edmar Caldas 06 Jun '18

06 Jun '18

Boa Tarde! estou tentando converter uma variável character em numerica e não estou conseguindo, alguem pode me ajudar por favor? > typeof(mes_ano)[1] "character"> agregar<-as.numeric(as.character(mes_ano))> is.numeric(mes_ano)[1] FALSE Edmar

2 2

linhas de densidade com legenda
by Alice Peres 05 Jun '18

05 Jun '18

Bom dia, caros membros Estou tendo dificuldade, vocês podem me ajudar a a transformar o gráfico abaixo em 4 linhas de densidade com legenda? x=c(3.42,4.01,4.04,4.26,4.10,5.61,4.90,3.67,5.30,5.05,4.44,3.50,4.09,5.05,4.40,3.99,3.86, 4.99,2.86,2.18,3.61,4.45,5.01,3.11,3.53,4.70,4.03,3.23,2.87,3.38,3.04,4.59,1.45,3.91,5.46,5.70, 3.33,3.38,5.43,3.33,4.37,3.83,3.45,5.74,3.64,3.23,4.32,2.91,4.69,4.04,3.10,3.29,4.04,3.41,4.57, 2.97,5.46,4.02,5.62,4.94,3.23,4.39,3.67,3.99,2.60,3.05,2.81,2.61,3.06,4.21,3.19,2.30,2.84,2.60, 4.32,2.55,3.23,3.66,3.26,3.90,4.60,3.99,3.42,4.53,3.20,3.46,3.19) hist(rexp(1000,1/min(x)),freq=FALSE) hist(rgamma(1000,shape=mean(x)^2,rate=mean(x)),freq=FALSE,add=TRUE) hist(rexp(1000,1/mean(x)),freq=FALSE,add=TRUE) hist(rexp(1000,1/max(x)),freq=FALSE,add=TRUE)

2 1

Re: [R-br] ggplot2_apresentar equação no gráfico
by cayanportela＠hotmail.com 05 Jun '18

05 Jun '18

Você pode gerar um modelo utilizando a função lm e extrair os coeficientes. De acordo com https://stackoverflow.com/questions/30144785/how-to-extract-intercept-and-s… , no seu caso ficaria assim: db %>% group_by(type) %>% do({ mod = lm(y ~ exp(-x), data = .) data.frame(Intercept = coef(mod)[1], Slope = coef(mod)[2]) }) Cayan _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

3 3

Erro na análise de regressão
by Angélica Ricarte 05 Jun '18

05 Jun '18

Prezados, boa noite. Estou com a seguinte dúvida, se puderem me ajudar, agradeço demais: Utilizando o pacote "ExpDes.pt", fiz a seguinte análise: >dados Oleo S 1 0.0 85.71 2 0.0 85.71 3 0.0 85.71 4 0.0 85.71 5 0.0 85.71 6 0.5 85.71 7 0.5 85.71 8 0.5 85.71 9 0.5 85.71 10 0.5 85.71 11 1.0 100.00 12 1.0 100.00 13 1.0 100.00 14 1.0 100.00 15 1.0 100.00 16 1.5 71.42 17 1.5 71.42 18 1.5 71.42 19 1.5 71.42 20 1.5 71.42 > dic(dados[,1],dados[,2],quali=F,sigF=0.05) ------------------------------------------------------------------------ Quadro da analise de variancia ------------------------------------------------------------------------ GL SQ QM Fc Pr>Fc Tratamento 3 2042 680.68 4.4071e+30 1.5357e-239 Residuo 16 0 0.00 Total 19 2042 ------------------------------------------------------------------------ CV = 0 % ------------------------------------------------------------------------ Teste de normalidade dos residuos Valor-p: 5.640644e-07 ATENCAO: a 5% de significancia, os residuos nao podem ser considerados normais! ------------------------------------------------------------------------ ------------------------------------------------------------------------ Teste de homogeneidade de variancia valor-p: NaN Error in if (pvalor.hvar < 0.05) { : valor ausente onde TRUE/FALSE necessário > dic(dados[,1],dados[,2],quali=F,sigF=0.05) Porque não consigo obter um resultado para a análise de regressão? Seria porque o CV é igual a zero e porque não há variância? Se sim, poderiam me dar uma explicação sobre a relação entre a variância e a regressão? Também tentei fazer o mesmo considerando os tratamentos qualitativos (quali=T) para obter um teste de médias, porém aparece a mesma mensagem de erro. OBS: para outras variáveis (que possuem variância e CV diferente de zero) deu certo a análise. Grata, aguardo retorno. Att., Angélica Batista.

2 3

IBGE - Distorção Malha de Setores Censitários Censo 2010
by Bruno Rezende 04 Jun '18

04 Jun '18

Prezados, 1. Estou trabalhando com a malha de setores censitários fornecida pelo IBGE (CENSO 2010). 2. Inicialmente, obtive a malha em formato .shp e carreguei no R. Notei que, ao projetar essa malha no Google Earth e no Leaflet, o desenho dos polígonos não coincidia com o desenho dos mapas. 3. Pensando que a distorção acima pudesse ser causada por algum erro meu, baixei o arquivo Google Earth .kmz dos setores censitários, fornecido diretamente pelo IBGE. 4. Contudo, os problemas são os mesmos: para minha cidade, algumas partes dos mapas se sobrepõem e algumas partes destoam do mapa "real". 5. Sabemos que um setor censitário é uma subdivisão do território com a finalidade de delimitar o trabalho do agente censeador a cerca de 250 a 350 domicílios. Uma divisa de setor censitário que cruze o meio de uma quadra não permite clareza na definição do setor para o agente que está em campo. A lógica infere que utilizar os próprios contornos do terreno, ou equipamentos públicos, facilita a divisão e clareza in loco. 6. Assim, deseja-se que, na maioria das vezes, as divisas do setor censitário e dos equipamentos públicos coincidam. 7. O problema é que o mapa dos setores não "casa" com o mapa "real". Assim, coletar uma coordenada em campo e projetar no shape file de setor censitário apresentará erro, inviabilizando análises. 8. Estou analisando a cidade de Passo Fundo - RS. Em Porto Alegre, capital do estado, essas distorções aparentam não ocorrer. 9. Tentei algumas transformações do CRS, mas sem sucesso. 10. Senhores, alguém de vocês já passo por problema semelhante? Saberiam indicar como corrigí-lo?

1 0

OpenBUGS no R
by Manasses Nóbrega 04 Jun '18

04 Jun '18

Prezados, Estou iniciando meus estudos em MCMC e uso do OpenBUGS no R. Ao tentar rodar um pequeno exemplo inicial (do livro do Kinas & Andrade, Introdução à Inferência Bayesiana com R), apareceu uma mensagem de erro que não consegui resolver. Abaixo segue o código que estou usando e as mensagens de erro que apareceram: lddados <- list(n=4,x=0) params <- list("teta") inits <- list(teta=0.5) sink("ru486_mcmc.txt") cat("model { teta ~ dunif(0,1) x ~ dbin(teta,n) }") sink() library(R2OpenBUGS) library(R2WinBUGS) library(BRugs) *Error :* .onLoad falhou em loadNamespace() para 'BRugs', detalhes: chamada: if (!is.na(ob$version)) msg <- paste(msg, "version", ob$version) else msg <- paste(msg, erro: argument is of length zero In addition: Warning messages: 1: package ‘BRugs’ was built under R version 3.3.3 2: In findOpenBUGS() : OpenBUGS 3.2.2 or greater must be installed (if so, this indicates missing registry keys of OpenBUGS). Setting the environment variable 'OpenBUGS_PATH' in advance of loading 'BRugs' overwrites the path. See ?loadOpenBUGS in order to load OpenBUGS manually. 3: In is.na(ob$version) : is.na() applied to non-(list or vector) of type 'NULL' Error: package or namespace load failed for ‘BRugs’ fit <- bugs(lddados,inits,params,model.file = "ru486_mcmc.txt", n.chains = 1,n.iter = 25000,n.burnin = 20000,n.thin = 5) *Error in* bugs(lddados, inits, params, model.file = "ru486_mcmc.txt", n.chains = 1, : Cannot find the OpenBUGS program Estou usando a versão 3.5 do R, e já tentei nas formas 32 e 64 bit. Também tentei o loadOpenBUGS e não tive êxito! Agradeço a ajuda de vocês na resolução desses erros. Um abraço a todos, Manasses -- *I am a brain, Watson. The rest of me is a mere appendix. * *(Sherlock Holmes)*

2 1

Re: [R-br] Digest R-br, volume 90, assunto 3
by Cayan Portela 04 Jun '18

04 Jun '18

Você pode extrair os coeficientes utilizando a função lm. De acordo com https://stackoverflow.com/questions/30144785/how-to-extract-intercept-and-s… ,no seu caso ficaria: db %>% group_by(type) %>% do({ mod = lm(y ~ exp(-x), data = .) data.frame(Intercept = coef(mod)[1], Slope = coef(mod)[2]) }) ________________________________ De: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> em nome de r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Enviado: segunda-feira, 4 de junho de 2018 15:00 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Digest R-br, volume 90, assunto 3 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://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flistas.in… 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. ggplot2_apresentar equação no gráfico (Henrique Boriolo Dias) ---------------------------------------------------------------------- Message: 1 Date: Mon, 4 Jun 2018 08:50:07 +0000 (UTC) From: Henrique Boriolo Dias <henrique.bdias(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] ggplot2_apresentar equação no gráfico Message-ID: <1597288900.8454811.1528102207591(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Olá pessoal, Sabem se é possível plotar a equação especifica (e talvez o R2) gerada pelo geom_smooth quando se usa a opção facet_wrap?Estou interessado no valor que multiplica o 'x' da equação.Alguém poderia me dar uma luz de qual o critério é usado para o ajuste dessa equação no geom_smooth? - Dados em anexo.- Script: db <- read.csv('data.csv', sep = ';', dec = '.', header = T) tiff(filename=paste0("graph.tiff"), width = 2000, height = 1000, res=300)ggplot(db, aes(y=y, x=x)) + geom_point(aes(shape=type))+ facet_wrap(~type) + geom_smooth(method = lm, formula = y ~ exp(-x), se = FALSE)dev.off() ObrigadoHenrique -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Flistas.inf…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: data.csv Tipo: application/vnd.ms-excel Tamanho: 5000 bytes Descrição: não disponível URL: <https://nam04.safelinks.protection.outlook.com/?url=http%3A%2F%2Flistas.inf…> ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flistas.in… ------------------------------ Fim da Digest R-br, volume 90, assunto 3 ****************************************

1 0

ggplot2_apresentar equação no gráfico
by Henrique Boriolo Dias 04 Jun '18

04 Jun '18

Olá pessoal, Sabem se é possível plotar a equação especifica (e talvez o R2) gerada pelo geom_smooth quando se usa a opção facet_wrap?Estou interessado no valor que multiplica o 'x' da equação.Alguém poderia me dar uma luz de qual o critério é usado para o ajuste dessa equação no geom_smooth? - Dados em anexo.- Script: db <- read.csv('data.csv', sep = ';', dec = '.', header = T) tiff(filename=paste0("graph.tiff"), width = 2000, height = 1000, res=300)ggplot(db, aes(y=y, x=x)) + geom_point(aes(shape=type))+ facet_wrap(~type) + geom_smooth(method = lm, formula = y ~ exp(-x), se = FALSE)dev.off() ObrigadoHenrique

1 0

Ajuda com amostra pequena
by Elias Carvalho 01 Jun '18

01 Jun '18

Bom dia Recebi um conjunto de dados (summary abaixo) para análise, porém vejo alguns problemas: 1. Amostra muito pequena, apenas 17 registros 2. Algumas variáveis com categorias ocorrendo apenas um ou duas vezes como sexo (SEX) por exemplo 3. Variável idade (IDA) e numero de alunos na sala (NAL) com muitas categoria 4. Variável gerenciamento (GBY) não possui variabilidade Esse é um banco de dados sobre bullying na escola, preciso analisar a capacidade do coordenador em definir bullying (DBY) reconhecer bullying (RBY) e gerenciar bullying (GBY) Alguem pode sugerir que tipo de ações eu poderia executar para conseguir fazer alguma análise sobre esses dados ou me dizer se é impossível fazer isso? -- *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 2

INDICADORES SOCIAIS OF TOPICS
by Edson Lira 01 Jun '18

01 Jun '18

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

4 4