- R-br - DINF UFPR

Criar um VETOR de um dataframe
by Manasses Nóbrega 12 Jun '18

12 Jun '18

Olá, pessoal! Este parece ser um problema fácil de resolver, mas não estou conseguindo!! Tenho um conjunto de dados (com quase 3000 linhas) que relaciona a resposta de um voluntário sob várias condições experimentais, cada uma dela com seus respectivos níveis. Exemplo: resp V1 V2 V3 0.35 1 1 1 0.38 1 2 1 0.65 1 3 1 0.49 1 4 1 -0.25 0 1 2 -1 0 2 2 -0.73 0 3 2 -0.47 0 4 2 Para realizar uma análise, por exemplo, vou precisar de *um vetor* com as respostas apenas para o nível 1 de V2 ou, por exemplo, de* um vetor* com as respostas apenas para o nível 2 de V3. Tentei a função aggregate mas não consegui. Vocês poderiam me ajudar? Desde já agradeço a todos. Um abraço, Manasses -- *I am a brain, Watson. The rest of me is a mere appendix. * *(Sherlock Holmes)*

3 3

Re: [R-br] como usar um layout gráfico do tipo 3x2 ? (Gerson Jr)
by Rafael Tieppo 12 Jun '18

12 Jun '18

a função grid.arrange() pode ajudar: library(gridExtra)grid.arrange(plot1, plot2) mais exemplos: https://www.r-graph-gallery.com/261-multiple-graphs-on-same-page/ Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: https://rafatieppo.github.io "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

1 0

como usar um layout gráfico do tipo 3x2 ?
by Cleber N.Borges 11 Jun '18

11 Jun '18

Olá a todos, Gostaria de fazer 5 gráficos em uma mesma janela... mas não entendi como usar as funções de layout... Gostaria de algo mais ou menos como coloco abaixo: ######## # primeira linha na janela gráfica => 3 gráficos plot( rnorm(10), col=1 ); plot( rnorm(10), col=2 ); plot( rnorm(10), col=3 ); # segunda linha na janela gráfica => 2 gráficos hist( rnorm(10), col=1 ); hist( rnorm(10), col=2 ); plot( rnorm(10), col=3 ); agradeço antecipadamente pela ajuda Cleber Borges --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

3 5

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