- R-br - DINF UFPR

Re: [R-br] Distribuição exponencial mista
by Alisson Lucrécio 20 Set '16

20 Set '16

Cesar, Preciso solucionar essa equação f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)) para obter os paramentos da distribuição. Assim como solucionado pelo fitdistr do pacote MASS para as ditribuções gamma, webull e exponencial. Observe que fiz uma tentativa com o pacote Renext, mas sem sucesso. Muito obrigado. 2016-09-19 16:12 GMT-03:00 Alisson Lucrécio < alisson.lucrecio(a)ifgoiano.edu.br>: > Como solucionar essa equação f(x) = (alpha/beta1*exp(-x/beta1)) + > ((1-alpha)/beta2*exp(-x/beta2)) para obter os paramentos da distribuição. > Assim como é feito pelo fitdistr do pacote MASS para as ditribuções > gamma, webull e exponencial. Observe que fiz uma tentativa com o pacote > Renext, mas sem sucesso. > > Muito obrigado. > > > On 19 September 2016 at 16:05, Cesar Rabak <cesar.rabak(a)gmail.com> wrote: > >> . . . e sua dúvida seria? >> >> On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br < >> r-br(a)listas.c3sl.ufpr.br> wrote: >> >>> Caro Colegas, >>> >>> Bom dia. >>> >>> Algum colega da poderia me ajudar na solução de uma distribução >>> exponencial mista. >>> >>> # Script >>> >>> library(MASS) >>> >>> # Distribuição Gamma >>> >>> rg <- rgamma(100, shape = 1.5, rate = 1) >>> >>> fit_gamma <- fitdistr(rg,"gamma") >>> >>> par(bg="white", las=1, cex=1.1) >>> plot(density(rg, bw=0.5, cut=0), las=1, lwd=2, >>> xlim=c(0,5),col="steelblue") >>> >>> # Distribuição Weibull >>> >>> rw <- rweibull(100, scale=1, shape=1.5) >>> >>> par(bg="white", las=1, cex=1.1) >>> plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, >>> xlim=c(0,5),col="steelblue") >>> >>> fit_weibull <- fitdistr(rw,"weibull") >>> >>> # Distribuição Exponential >>> >>> ex <- rexp(100, rate = 1) >>> >>> par(bg="white", las=1, cex=1.1) >>> plot(density(ex, bw=0.5, cut=0), las=1, lwd=2, >>> xlim=c(0,5),col="steelblue") >>> >>> fit_exp <- fitdistr(ex, "exponential") >>> >>> plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, >>> xlim=c(0,5),col="steelblue") >>> lines(density(rg, bw=0.5, cut=0),col="green") >>> lines(density(ex, bw=0.5, cut=0),col="red") >>> >>> # Distribuição Exponencial mista >>> >>> library(Renext) >>> >>> # f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)) >>> >>> mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) >>> >>> fit_mixexp <- mom.mixexp2(mixexp) >>> >>> plot(density(mixexp, bw=0.5, cut=1), lwd=2, >>> xlim=c(0, 5),col="steelblue") >>> lines(density(ex, bw=0.5, cut=0),col="red") >>> >>> mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, >>> rate2), >>> start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = >>> nls.control(maxiter = 100)) >>> >>> Obrigado. >>> -- >>> Alisson Lucrecio da Costa >>> >>> _______________________________________________ >>> 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. >>> >> >> > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

2 1

Erro no RStudio
by ROSIANA RODRIGUES ALVES 20 Set '16

20 Set '16

Olá Boa Noite, Estou precisando de uma solução para um problema que descobri no RStudio. Todo vez que tento instalar um biblioteca , ou importar um arquivo de dados pela abas disponíveis no RStudio e o nome de usuário do computador ou o nome da pasta contém acentos aparece essa mensagem de erro "r code execution error". A solução que arrumei foi sempre instalar por linhas de comando, mas a maioria dos meus alunos preferem a forma "clique-aqui-clique-ali". 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 ____________________________________________________________________________ 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.

5 4

Usar R para Salvar Referencias
by Elias Carvalho 19 Set '16

19 Set '16

Boa tarde pessoal Alguém sabe se é possível fazer uma busca no google scholar e salvar os resultados automaticamente no mendeley, tudo via script ? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho*

3 2

Distribuição exponencial mista
by Alisson Lucrécio 19 Set '16

19 Set '16

Caro Colegas, Bom dia. Algum colega da poderia me ajudar na solução de uma distribução exponencial mista. # Script library(MASS) # Distribuição Gamma rg <- rgamma(100, shape = 1.5, rate = 1) fit_gamma <- fitdistr(rg,"gamma") par(bg="white", las=1, cex=1.1) plot(density(rg, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") # Distribuição Weibull rw <- rweibull(100, scale=1, shape=1.5) par(bg="white", las=1, cex=1.1) plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") fit_weibull <- fitdistr(rw,"weibull") # Distribuição Exponential ex <- rexp(100, rate = 1) par(bg="white", las=1, cex=1.1) plot(density(ex, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") fit_exp <- fitdistr(ex, "exponential") plot(density(rw, bw=0.5, cut=0), las=1, lwd=2, xlim=c(0,5),col="steelblue") lines(density(rg, bw=0.5, cut=0),col="green") lines(density(ex, bw=0.5, cut=0),col="red") # Distribuição Exponencial mista library(Renext) # f(x) = (alpha/beta1*exp(-x/beta1)) + ((1-alpha)/beta2*exp(-x/beta2)) mixexp <- rmixexp2(n = 100, prob1 = 0.5, rate1 = 1.0, rate2 = 3.0) fit_mixexp <- mom.mixexp2(mixexp) plot(density(mixexp, bw=0.5, cut=1), lwd=2, xlim=c(0, 5),col="steelblue") lines(density(ex, bw=0.5, cut=0),col="red") mod <- nls(mixexp ~ prob1 * dexp(x, rate1) + (1 - prob1) * dexp(x, rate2), start = list(prob1 = 1, rate1 = 1, rate2 = 3.0), control = nls.control(maxiter = 100)) Obrigado. -- Alisson Lucrecio da Costa

2 1

Divulgação de pacote para stemming
by Daniel Falbel 17 Set '16

17 Set '16

Boa tarde, Acabei de publicar no CRAN a versão 0.0.2 do pacote ptstem. https://github.com/dfalbel/ptstem Esse pacote implementa 3 algoritmos para fazer stemming na língua portuguesa. Espero que seja útil. Abraços, Daniel

1 0

Criação de um Fluxograma
by Abel Brasil Ramos da Silva 14 Set '16

14 Set '16

Olá a todos, Gostaria de saber como faço para produzir um fluxograma/diagrama de um data.frame, a ideia é produzir algo desse tipo, porém inserindo os totais de cada subamostra e a porcentagem em cada nó. P.S: A forma de saída pode ser qualquer uma. |Econômica /--S --- classe |Master |Executiva M---surv | |Econômica \--N --- classe |Master |Executiva sex--| |Econômica /--S --- classe |Master |Executiva F---surv | |Econômica \--N --- classe |Master |Executiva Exemplo seria um banco similar a esse: set.seed(1) sexo=sample(c(0,1),50,replace=T) sexo=factor(sexo,levels=c(0,1),label=c("Feminino","Masculino")) surv=sample(c(0,1),50,replace=T) surv=factor(surv,levels=c(0,1),label=c("Não","Sim")) classe=sample(c(0,1,2),50,replace=T) classe=factor(classe,levels=c(0,1,2),label=c("Econômica"," Master","Executiva")) dados=data.frame(sexo,surv,classe) s=table(dados$sexo);s sp=round(prop.table(s)*100,2);sp Msexo=subset(dados,sexo=="Masculino") cM=table(Msexo$classe);cM cpM=round(prop.table(c)*100,2);cpM Fsexo=subset(dados,sexo=="Feminino") cF=table(Fsexo$classe);cF cpF=round(prop.table(c)*100,2);cpF Agradeço pela ajuda desde já. Abs.

3 3

Relação de cidades IBGE x TSE
by Rodrigo Coster 12 Set '16

12 Set '16

Boa tarde, Alguém teria uma lista relacionando os códigos de cidade IBGE e TSE? Por exemplo, enquanto para o IBGE (e a maioria dos outros órgãos), Porto Alegre / RS tem o código 431490-2, enquanto para o TSE o código é 88013. Há um tempo atrás eu fiz essa relação através do nome das cidades, mas havia diferença na grafia de algumas cidades. []'s

3 3

Passando um objeto para "character"
by Pedro Emmanuel Alvarenga Americano do Brasil 09 Set '16

09 Set '16

Ola amigos de R, Estou tentando fazer um função para um plot mas sem obrigar que o usuário chame as variáveis de um banco entre aspas. Isso porque no conjunto de funções hierarquizadas, todas as outras não precisam que os argumentos sejam character, vai facilitar muito inclusive pra documentar. Eu achei a função Hmisc::Cs que faz o que eu gostaria, mas não gostaria de criar uma dependencia do Hmisc apenas por conta dessa função. Ao mesmo tempo, eu não consigo saber o que tem dentro da Cs() por que quando digito Cs no console, esse retorna apenas o namespace. Eu gostaria de clonar a Cs ou dar um outro jeito de não obrigar que usuario entrasse as variáveis como character. Eu dei uma olhada nas fuçoes parse, deparse e substitute no pacote base mas não consegui me desenrolar. Assim funciona mydata <- data.frame(Var1 = rnorm(10, 1 , 2), Var2 = rnorm(10, 2, 2)) myplot <- function(x, y, data){ plot(data[ , c(x, y)]) } myplot("Var1", "Var2", mydata) # Mas eu gostaria de algo assim... myplot <- function(x, y, data){ if ( is.character(x) & is.character(y)){ plot(data[ , c(x, y)]) } else { # essa parte não funciona plot(data[ , Hmisc::Cs(x, y)]) # essa é a parte que estou enrolado pq não quero depender de Hmisc } myplot(Var1, Var2, mydata) Eu pensei numa coisa do tipo ... with(mydata, plot(Var1 , Var2)) Mas dentro da função não funciona... myplot <- function(x, y, data){ if ( is.character(x) & is.character(y)){ plot(data[ , c(x, y)]) } else { with(data, plot(x, y)) } } myplot( Var1, Var2, mydata) Abraços fortes, Pedro Brasil

3 4

Erro devtools no "build & reload" com Rstudio
by Pedro Emmanuel Alvarenga Americano do Brasil 09 Set '16

09 Set '16

Amigos de R, Uma luz por favor. Estou montando um pacote utilizando o Rstudio. Fiz algumas funções, estou usando o roxygen2 pra fazer a documentação. Estou tentando documentar as funções que já tenho e estava tudo indo bem até que o seguinte erro começou a impedir o meu progresso. O erro aparece na janela "build" quando pressiono "build & reload". Tentei desfazer as ultimas modificações que fiz para descobrir o que estava causando o erro mas sem sucesso. ==> devtools::document(roclets=c('rd', 'collate', 'namespace')) Updating DiagnosisMed documentation Loading DiagnosisMed Writing SS.Rd Error in dirname(path) : caminho muito longo Calls: suppressPackageStartupMessages ... roc_output.rd_roclet -> mapply -> <Anonymous> -> file.exists -> dirname Execu��o interrompida Exited with status 1. O que me deixa mais curioso é que eu não mudei a pasta do pacote e quando eu vou nas configurações do "build" para esse pacote, a pasta destino está correta. Alguma luz de como poderia acertar esse erro? Pedro Brasil

1 1

Re: [R-br] Parcela Subdividida
by Angélica Ricarte 08 Set '16

08 Set '16

Prezado Cesar, Não consegui compreender a mensagem de erro e nem o porquê que minha matriz é nula. Você pode ajudar? Agradeço e aguardo retorno. Obrigada. Att., Angélica Date: Mon, 5 Sep 2016 12:52:41 -0700 From: ml-node+s2285057n4666617h13(a)n4.nabble.com To: angelica.ricarte(a)hotmail.com Subject: Re: [R-br] Parcela Subdividida Angélica, Você entendeu a mensagem de erro? Não é possível analisar uma matriz nula.... 2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br <[hidden email]>: Prezados, Estou tentando realizar uma análise de variância para um experimentos em parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode me ajudar? Vejam como realizei as análises e os seus resultados: #ANOVA: > mod<-aov(PS[,4]~PS[,1]*PS[,2]+Error(PS[,1]:PS[,3]));summary(mod) Warning message: In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) : modelo Error() é singular Error: PS[, 1]:PS[, 3] Df Sum Sq Mean Sq F value Pr(>F) PS[, 1] 1 8.1 8.100 11.3 0.0152 * Residuals 6 4.3 0.717 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) PS[, 2] 4 108.9 27.23 24.94 3.02e-08 *** PS[, 1]:PS[, 2] 4 496.9 124.23 113.79 3.09e-15 *** Residuals 24 26.2 1.09 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1#Resíduos: > residuals(mod) NULL #Teste de Shapiro-Wilk: > res<-residuals(mod) > shapiro.test(res) Error: is.numeric(x) is not TRUE Obrigada. Att., Angélica. _______________________________________________ R-br mailing list [hidden email] 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. _______________________________________________ R-br mailing list [hidden email] 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. If you reply to this email, your message will be added to the discussion below: http://r-br.2285057.n4.nabble.com/R-br-Parcela-Subdividida-tp4666613p466661… To unsubscribe from R-br, click here. NAML

2 2

Parâmetros do plot
by Andre Oliveira 08 Set '16

08 Set '16

Bom dia,pesquisei e conversei com algumas pessoas e não resolvi o meu problema. Como tirar estes vetores desta imagem feita com o plot?require(bpca) fit=bpcaplot(bpca(iris[2:4],center =2),var.factor = 0.2) obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 2

(sem assunto)
by Elias Carvalho 07 Set '16

07 Set '16

Boa tarde a todos. Eu tenho uma matriz, chamada "matrix.results", como resultado de algumas análises: Consigo imprimir uma tabela no console com o comando: kable(head(matrix.results ), format = "rst", align = "r") Agora eu gostaria de salvar uma imagem desta tabela para inserir em um documento. No entanto eu gostaria de salvar como imagem uma tabela mais bem formatada, como esta: Alguém poderia ajudar ? *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho*

3 3

Parcela Subdividida
by Angélica Ricarte 05 Set '16

05 Set '16

Prezados, Estou tentando realizar uma análise de variância para um experimentos em parcelas subdivididas em DIC, mas, apesar de ser gerado um resultado, aparece uma mensagem de erro. Quando tento realizar o cálculo dos resíduos também dá erro e, consequentemente, não consigo fazer o teste de normalidade. Alguém pode me ajudar? Vejam como realizei as análises e os seus resultados: #ANOVA: > mod<-aov(PS[,4]~PS[,1]*PS[,2]+Error(PS[,1]:PS[,3]));summary(mod) Warning message: In aov(PS[, 4] ~ PS[, 1] * PS[, 2] + Error(PS[, 1]:PS[, 3])) : modelo Error() é singular Error: PS[, 1]:PS[, 3] Df Sum Sq Mean Sq F value Pr(>F) PS[, 1] 1 8.1 8.100 11.3 0.0152 * Residuals 6 4.3 0.717 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) PS[, 2] 4 108.9 27.23 24.94 3.02e-08 *** PS[, 1]:PS[, 2] 4 496.9 124.23 113.79 3.09e-15 *** Residuals 24 26.2 1.09 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1#Resíduos: > residuals(mod) NULL #Teste de Shapiro-Wilk: > res<-residuals(mod) > shapiro.test(res) Error: is.numeric(x) is not TRUE Obrigada. Att., Angélica.

2 1

Re: [R-br] Exclusão de casos por critério temporal
by sznelwar＠uol.com.br 04 Set '16

04 Set '16

1 0

Re: [R-br] Frequencia de linhas duplicadas
by sznelwar＠uol.com.br 04 Set '16

04 Set '16

1 0

Shapefile de bairros de Belém - PA
by Daniel Marcelino 04 Set '16

04 Set '16

Olá, Eu estou procurando por um shapefile com a malha de bairros de Belém do Pará, alguém já trabalhou em algum projeto que contenha essa malha? Att. Daniel

2 2

Unindo colunas de um data-frame
by Emerson Cotta Bodevan 03 Set '16

03 Set '16

Pessoal, boa noite. Tenho um conjunto de dados armazenados em um data.frame com três colunas, como segue: A B C 3 4 6 1 2 9 3 2 5 . . . Gostaria de gerar um vetor com os seguintes números: 346 129 325 ... Ou seja, os números do vetor possuem três algarismos. Eles são formados pela união dos trẽs números de cada linha do data frame. Obrigado por qualquer ajuda, *Emerson*

6 6

Leitura de range de aba de arquivo xlsx Excell no R
by Wilson Cohen 01 Set '16

01 Set '16

Leonard Assis Velho, instala e carrega o *"openxlsx"*. Depois abre o arquivo com *read.xlsx("arquivo.xlsx")* Espero ter ajudado! Um abraço! -- *Prof. MSc. Wilson Martins da Silva* *Mestre em Ecologia Aquática e Pesca - PPGEAP/ICB/UFPA* *Especialista em Biodiversidade Amazônica - UFPA/ATMDepto de Morfologia e Ciências Fisiológicas* *Centro de Ciências Biológicas e da Saúde* *Universidade do Estado do Pará* *Campus IX - Altamira - Pará* http://lattes.cnpq.br/2440617911389759 wilsoncelula(a)gmail.com, wilsoncelula(a)hotmail.com <wilsonms(a)ufpa.br>

1 0

Transformar distribuição de frequencia em dataset
by Ronaldo Câmara de Araújo 01 Set '16

01 Set '16

Bom dia listeiros, Gostaria da ajuda de vcs, procurei a solução aqui mas não vi post dessa natrureza. Utilizando o ftable ou xtabs executo o comando - ftable(cad2016$Mês.Ano ~ cad2016$SIGLADR + cad2016$SEXO) - o output é o que segue abaixo. Preciso disso em um data.frame, só que dessa forma, ou seja, as colunas sendo essas datas. Quando transformo isso em um data frame, ou seja, a <- data.frame(ftable(cad2016$Mês.Ano ~ cad2016$SIGLADR + cad2016$SEXO)) vem um data frame com essas colunas na linha. Se não tiver solução para a distribuição de frequencia abaixo, pelo menos no - ftable(cad2016$Mês.Ano ~ cad2016$SIGLADR) cad2016$Mês.Ano 2016-01-01 2016-02-01 2016-03-01 2016-04-01 2016-05-01 2016-06-01 2016-07-01cad2016$SIGLADR cad2016$SEXO AC F 1871 1871 1875 1887 1909 1913 1909 M 2232 2232 2244 2260 2266 2270 2268DR/ACR F 87 85 85 86 85 85 85 M 276 277 276 275 275 275 274DR/AL F 226 226 226 226 226 227 227 M 1007 1007 1007 1004 999 1000 996DR/AM F 200 200 198 198 198 198 198 Desde já, agradeço a quem possa me ajudar. Ronaldo Câmara de Araújo

2 1

Leitura de range de aba de arquivo xlsx Excell no R
by Maurício Cunha 31 Ago '16

31 Ago '16

Bom dia Pessoal, Estou tentando ler dados de uma planilha Excel com muitas abas, trabalhando com RStudio no Mac. Dei os seguintes comandos: > install.packages("XLConnect") trying URL ' https://cran.rstudio.com/bin/macosx/mavericks/contrib/3.3/XLConnect_0.2-12.… ' Content type 'application/x-gzip' length 5546911 bytes (5.3 MB) ================================================== downloaded 5.3 MB tar: Failed to set default locale The downloaded binary packages are in /var/folders/_4/f5_n3z0x24nc08f5whq81nl80000gn/T//RtmpxbUYbW/downloaded_packages > library(XLConnect) Loading required package: XLConnectJars JavaVM: requested Java version ((null)) not available. Using Java at "" instead. JavaVM: Failed to load JVM: /bundle/Libraries/libserver.dylib JavaVM FATAL: Failed to load the jvm library. Error : .onLoad failed in loadNamespace() for 'XLConnectJars', details: call: .jinit() error: JNI_GetCreatedJavaVMs returned -1 Error: package 'XLConnectJars' could not be loaded Ou seja o XLConnect instala mas a library dá problema. Alguém teria alguma orientação? Obrigado, Maurício

4 3

Avaliação da evolução do peso
by Edson Lira 31 Ago '16

31 Ago '16

Bom dia a todos da Lista! Abaixo o link de um arquivo teste que estou tentando fazer uma análise.https://www.datafilehost.com/d/UgWFEhooCR Os dados para análise são registros de comparecimento de pacientes em uma unidade de saúde. Cada vez que o paciente comparece na unidade uma das variáveis que é registrada é seu peso, portanto o mesmo paciente pode aparecer uma, duas ou mais vezes, cada uma aparição é um registro. Como faço para avaliar seu peso a cada vez que comparece? considerando com o peso inicial e com os pesos subsequentes registrados. Se um paciente comparece três vezes, seu peso será comparado três vezes (a segunda com a primeira e a terceira com a segunda e a terceira com a primeira). Obrigado a todos![ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 4

Shapefiles x latitude e longitude
by Jailson 31 Ago '16

31 Ago '16

Prezados, bom dia! Peço desculpas se minha dúvida já tenha sido abordada anteriormente, mas não encontrei ao pesquisar nas mensagens arquivadas. Tenho os shapefiles dos municípios do IBGE. Também possuo uma base de dados com cerca de 2000 registros, cada um com uma coordenada de latitude e longitude. Preciso verificar se a coordenada (latitude e longitude) de um registro qualquer pertence ao município a que ele deveria pertencer. Plotando as coordenadas no mapa eu consigo verificar visualmente, mas gostaria de fazer a verificação por programa, onde seriam listados apenas os registros cujas coordenadas não pertencessem ao respectivo município. Agradeço antecipadamente a quem puder me ajudar. Jailson

3 2

Análise de Sequências
by Paulo Nogueira Starzynski 30 Ago '16

30 Ago '16

Boa noite colegas de lista. Terei de abordar um problema de detecção de padrão na ocorrência de eventos ao longo do tempo, isto é, determinar por exemplo que a ocorrência do evento A está temporalmente relacionada ao evento B, e que este está ligado ao evento G. Assim, eu apontaria a sequência de eventos ABG como relevante. Minha base de dados contém uma coluna informando o id do evento, e outra a data/hora de sua ocorrência. São milhares de eventos diferentes e não há informação à priori de quais podem formar as sequências mais prováveis. Comecei a pesquisar alguns pacotes no R que abordam problemas desse tipo e cheguei no *arules* e *TraMineR*. Alguém já trabalhou questões como essa e pode sugerir algo mais? (em tempo, acredito que todas as opções de análise de sequência genética não podem me ajudar nesse problema) Abraços, Paulo Nogueira Starzynski

1 0

Visualizar grande data.frame
by Abel Brasil Ramos da Silva 30 Ago '16

30 Ago '16

Olá a Todos, Bom dia! Possuo um data.frame com 10000 observações e 1350 colunas. Gostaria de apenas visualizar esse dados no próprio R (ou R-Studio), que pacote e comando para executar essa tarefa. Com o comando View existe uma limitação da visualização. Agradeço a ajuda desde já. Abs.

4 9

R-Fiddle
by Raphael Saldanha 30 Ago '16

30 Ago '16

Olá pessoal, Não sei o pessoal conhece, mas não encontrei nenhuma postagem mencionando o R-Fiddle. http://www.r-fiddle.org Trata-se de uma página parecida com o JSFiddle, onde você pode escrever na linguagem R, rodar e ver os resultados, salvar e compartilhar o link. Muito útil na criação de código mínimo reproduzível. Segue um exemplo: http://www.r-fiddle.org/#/fiddle?id=UpgptFnA&version=2 Abs, Raphael Saldanha

1 0

Criando pacotes com o Rstudio
by Pedro Emmanuel Alvarenga Americano do Brasil 30 Ago '16

30 Ago '16

Amigos de R, Eu estou tentando organizar uma coleção de funções em um pacote utilizando o Rstudio. Ja fucei na documentação do Rstudio e alguns sitios por ai mas não consigui uma luz. A minha questão é que tenho tres funções em um mesmo arquivo *.R, que retornam a mesma análise mas por métodos diferentes. Mais ou menos assim... # bloco do roxygen2 #' Diferentes métodos para a mesma análise #' @aliases f1 f2 f3 #'... # bloco das funções f1 <- function(arg1, arg2, arg3) {...} f2 <- function(arg1, arg2, arg3) {...} f3 <- function(arg1, arg2, arg3) {...} O meu problema é que quando eu pressiono "Build & reload", o Rstudio não retorna qualquer aviso ou erro, mas so carrega no pacote a primeira função. Então, quando chamo a segunda ou terceira ele diz que a função não existe. Eu já pensei até em separar as funções em diferentes arquivos, mas gostaria muito de manter essa organização. Até o momento estou achando que preciso fazer alguma configuração do devtools ou do Rstudio pra poder fazer funcionar, mas não estou sabendo qual. Alguma luz de como resolver? Abraço forte, Pedro Brasil

2 2

Sub-amostragem
by Marcos Bissoli 30 Ago '16

30 Ago '16

Prezados, Peço desculpas de antemão se meu problema é de extrema trivialidade ou se o tema já fora aqui debatido. No entanto, pesquisei os arquivos do fórum e não encontrei, dessa vez, uma solução. Gerei um código simplificado de meu problema. Meu intuito é extrair uma sub-amostra de uma amostra de dados já coletados para fins de análise com melhor equilíbrio entre os grupos observados. Assim, elaborei algo semelhante com o abaixo. Id <- 1:20 X <- c("A","B","A","C","A","A","B","A","C","A","A","B","A","C","A","A","B","A","C","A") Dados <- data.frame(Id,X) IdDados <- sample(Dados$Id,5,replace = FALSE) DadosAm <- Dados[id==IdDados[1],] for (i in 2:5) DadosAm <- rbind(DadosAm,Dados[id==IdDados[i],]) DadosAm Este código funciona, e consigo extrair uma sub-amostra aleatória de n1=5 a partir de uma amostra inicial de n=20 em data.frame. No entanto, meu real problema é gerar uma sub-amostra de n1=88 em uma amostra inicial de n=1668. Mas, quando tento fazer com tais dimensões a sub-amostra gera uma série de NA's, que não existem na amostra original. Chequei o funcionamento na amostra real e maior, e percebi que o primeiro Id de IdDados não corresponde ao Id de Dados adicionado já no primeiro comando de criação de DadosAm. Creio que seja algum erro meu de implementação, mas cheguei a fazer testes com simulações de n e n1 maiores e deram certo. Mas quando vou para o meu banco real permanece o problema. Desde já, agradeço qualquer ajuda, e reitero minhas desculpas pela possível trivialidade da dúvida. Saudações acadêmicas, -- MARCOS BISSOLI Faculdade de Nutrição Universidade Federal de Alfenas Blog: bocademiamaldita.blogspot.com/ E-mail: mbissoli(a)gmail.com Twitter: #mbissoli Alfenas, Minas Gerais, Brasil *****Pense na Natureza antes de Imprimir***** Divulgue ON-LINE Eu apoio a ENEN "na luta por um Brasil sem fome" "por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan" (para cada povo sua própria língua, para todos os povos o Esperanto) E nunca votarei no PSDB/DEM!

5 6

r baixando dados inmet
by Rafael Tieppo 29 Ago '16

29 Ago '16

Algum tempo atrás o Alisson Lucrecio e o Éder Comunello estavam trabalhando em um tópico sobre obtenção de dados do INMET. Como precisei obter alguns dados hoje, resolvi testar o código, que por sinal funcionou perfeitamente.Como os dados são salvos em formato ```.html```, tentei montar um script organizar os dados em um ```data.frame```, assim como deixar apenas uma linha para cada dia de coletado, pois os dados do originais do INMET tem formato de duas linhas para cada dia de coleta. Segue o script com os comentários:obs: a função que ordena os dias é não tem uma boa performance, pois faz um bom tempo que fiz ela. #=================================================================# Rafael Tieppo# rafaelt(a)unemat.br# http://docente.unemat.br/rafaeltieppo/# 29-08-2016#================================================================= ##### Goal 1 ### Script to sign in - INMET### Historical Data### Original Source:### http://r-br.2285057.n4.nabble.com/R-br-r-baixando-dados-inmet-td4660459.htm… http://r-br.2285057.n4.nabble.com/R-br-RCurl-td4659610.html##### Goal 2### Ordering and filtering DATA from INMET### use oneline() function from### from https://github.com/rafatieppo/INMET_DATA_ORDER #************************************************************#************************************************************ # GOAL 1#************************************************************#************************************************************ #------------------------------------------------------------### Packageslibrary(RCurl)library(bitops)#------------------------------------------------------------ #------------------------------------------------------------### Logging INMET ### Login linkmyURL1 <- "http://www.inmet.gov.br/projetos/rede/pesquisa/inicio.php"### Data linkmyURL2 <- "http://www.inmet.gov.br/projetos/rede/pesquisa/gera_serie_txt.php?&mRelEsta…,"#------------------------------------------------------------ #------------------------------------------------------------### Access Data myParams=list( mCod="EMAIL", ### alterar! mSenha="PASSWORD", ### alterar! btnProcesso = " Acessar ")#------------------------------------------------------------ #------------------------------------------------------------### Getting Data myCurl <- getCurlHandle()curlSetOpt(cookiejar="cookies.txt", useragent="Chrome/10.0.648.133" , followlocation=TRUE, curl=myCurl)###"Mozilla/5.0" login <- postForm(myURL1, .params=myParams, curl=myCurl) DATA <- getURLContent(myURL2, curl=myCurl)#------------------------------------------------------------ #************************************************************#************************************************************ # GOAL 2#************************************************************#************************************************************ #------------------------------------------------------------### Tiding Data ### Using shell script to separate DATA from text#### grep "^83309" < ESTACAO_83309.html > ESTACAO_83303_DATA.csv### Shell script, get all lines that starts with "83309" (station number) SHELL_FUN = paste("grep", "^83309", "<", "ESTACAO_83309.html", ">", "ESTACAO_83303_DATA.csv", sep = ' ')### system(SHELL_FUN) DATA_83303 <- read.csv("ESTACAO_83303_DATA.csv", sep = ";", dec = ".") colnames(DATA_83303) <- c("Estacao", "Data", "Hora", "Precipitacao", "TempMaxima", "TempMinima", "Insolacao", "Evaporacao_Piche", "Temp_Comp_Media", "UR_Media", "Vel_Vent_Media") names(DATA_83303)#------------------------------------------------------------ #------------------------------------------------------------### Organizing DATA### Data from INMET has two lines for each### To put one line for each day, use oneline() function ### Calling ONE_LINE function### from https://github.com/rafatieppo/INMET_DATA_ORDER ### ATTENTION### To use oneline() is mandatory a data.frame with a specific cols data### order, as follow (names doen not matter, only sequence): ### Estacao; Data; Hora; Precipitacao; TempMaxima; TempMinima;### Insolacao; Umidade Relativa Media; Velocidade do Vento Media; DATA_83303 <- DATA_83303[c(1:7,10:11)] Estacao;Data;Hora;Precipitacao;TempMaxima;TempMinima;Insolacao;EvaporacaoPiche;Temp Comp Media;Umidade Relativa Media;Velocidade do Vento Media; ### Ordering with oneline() function DATA_83303_one_line <- one_line(DATA_83303) edit(DATA_83303_one_line) #------------------------------------------------------------ Rafael Tieppo State University of Mato Grosso - Department of Agricultural Engineering site: http://docente.unemat.br/rafaeltieppo/ blog: https://fuidebicicleta.wordpress.com"Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

2 2

it's shocking
by Adevair Henrique da Fonseca 28 Ago '16

28 Ago '16

Hi! Have you seen these impressive pics? They are shocking! Just take a look <http://nkequatristi.elysburghauntedhouse.com/e4pwj> See you soon, Adevair Henrique da Fonseca

1 0

Fw: insomnia
by Adevair Henrique da Fonseca 28 Ago '16

28 Ago '16

Hey! When I saw that stuff I couldn't fall asleep and had insomnia for a couple of nights, so I decided to show it to you, just take a look <http://truhustune.tradefleets.com/e4sxuw> My best to you, Adevair Henrique da Fonseca

1 0

GOM.EM
by Anna Karoline R. da Cruz 26 Ago '16

26 Ago '16

Boa tarde, pessoal! Estou usando o "GoM.em" para um conjunto de dados. Minha dúvida é: como descrevo os perfis a partir do lambda que ele fornece? OBS: tenho algumas as variáveis com mais de duas categorias. > mod1$lambda [,1] [,2] [,3] [1,] 0.125000000 0.375000000 0.62500000 [2,] 0.999000000 0.758576148 0.94791392 [3,] 0.999000000 0.999000000 0.99900000 [4,] 0.999000000 0.192048304 0.99900000 [5,] 0.958999554 0.239871163 0.99900000 [6,] 0.999000000 0.582693490 0.60513053 [7,] 0.760374047 0.999000000 0.99900000 [8,] 0.001000000 0.046672897 0.90641452 [9,] 0.475789685 0.883268404 0.10884461 [10,] 0.355654569 0.232261520 0.05983177 [11,] 0.236404268 0.001038477 0.12928587 [12,] 0.009139838 0.364151602 0.28363626 [13,] 0.001000000 0.604041999 0.99900000 [14,] 0.268329215 0.339902463 0.99900000 [15,] 0.001000000 0.001000000 0.78822352 [16,] 0.031297098 0.338268504 0.99900000 [17,] 0.001000000 0.001000000 0.26163060 [18,] 0.999000000 0.999000000 0.99900000 [19,] 0.001000000 0.001000000 0.13074344 [20,] 0.999000000 0.999000000 0.99900000 [21,] 0.999000000 0.999000000 0.99900000 [22,] 0.999000000 0.999000000 0.99900000 [23,] 0.001000000 0.001000000 0.00100000 [24,] 0.768230734 0.999000000 0.47083018 [25,] 0.001000000 0.001000000 0.39344518 [26,] 0.532563903 0.999000000 0.00100000 [27,] 0.001000000 0.001000000 0.00100000 [28,] 0.289341873 0.999000000 0.21041057 [29,] 0.001000000 0.001000000 0.13074344 [30,] 0.001000000 0.790611117 0.00100000 [31,] 0.001000000 0.001000000 0.00100000 [32,] 0.531285647 0.999000000 0.00100000 [33,] 0.001000000 0.001000000 0.26163060 [34,] 0.283878725 0.999000000 0.00100000 Desde já agradeço a ajuda! Abs. Anna.

2 3

gorgeous
by david_ufvjm 26 Ago '16

26 Ago '16

Dear! That stuff is just gorgeous, you have to read about it more at <http://slaprislame.rldamodaran.com/e4xfwhrm> Hope this helps, david_ufvjm

1 0

Oportunidades PD/FAPESP
by Benilton Carvalho 24 Ago '16

24 Ago '16

Caros, No contexto do BRAINN (CEPID/FAPESP), temos duas oportunidades PD na UNICAMP, conforme descrições abaixo. Fico deveras agradecido caso possam também compartilhar com potenciais candidatos. Att, Benilton ## Oportunidade 1 O Instituto Brasileiro de Neurociência e Neurotecnologia (BRAINN), localizado no campus da UNICAMP em Campinas, SP. Procura candidatos a bolsa de Pós-Doutorado, financiada pela FAPESP, no âmbito de projeto CEPID. A bolsa tem duração Inicial de dois anos. O (a) candidato (a) deverá ter obrigatoriamente graduação em biologia ou ciências biológicas e experiência de pesquisa na área de análise de dados biológicos complexos. O (a) candidato (a) deve ter obrigatoriamente experiência com programação em R, modelagem matemática e estatística de sistemas biológicos complexos e análise de dados de proteômica. Candidatos (as) devem ter concluído o doutorado há não mais que 4 anos e demonstrar ser capaz de gerar produção científica, tendo publicado pelo menos dois artigos em periódicos científicos com revisão por pares. A inscrição será no período de 22/08/2016 a 05/09/2016 através do email: icendes(a)unicamp.br, sendo necessário encaminhar o Curriculum Vitae, Histórico Escolar da Pós-graduação e duas cartas de recomendação. ## Oportunidade 2 O Instituto Brasileiro de Neurociência e Neurotecnologia, localizado no campus da UNICAMP, busca candidatos altamente motivados e analiticamente capazes para oportunidade de pós-doutorado, vinculada a projeto CEPID financiado pela FAPESP. A bolsa tem duração inicial de dois anos. O candidato(a) deverá ter obtido o título de Doutor a menos de 4 anos em uma das seguintes áreas: Estatística, Bioestatística, Matemática, Biologia Computacional, Bioinformática ou áreas afins; e possuir interesse pelo trabalho na área de neurociências e, também, análise integrativa de dados -ômicos. Proficiência em programação (por exemplo, R, Python, Perl, C, C++ ou Java) e, pelo menos, três artigos em periódicos científicos com revisão por pares nos últimos três anos são requeridos. Excelente comunicação oral e escrita (Inglês, Português ou Espanhol) é indispensável. A inscrição será no período de 22/08/2016 a 05/09/2016 através do email: icendes(a)unicamp.br, sendo necessário encaminhar o Curriculum Vitae, Histórico Escolar da Pós-graduação e duas cartas de recomendação. http://www.fcm.unicamp.br/fcm/relacoes-publicas/saladeimprensa/cepid-brainn- oportunidades-de-pos-doutorado-em-biologia-computacional-e-bioestatistica https://support.bioconductor.org/p/86410/

1 0

Re: [R-br] O R não abre no linux mint
by Éder Comunello 24 Ago '16

24 Ago '16

Sérgio, bom dia! Tentei reproduzir o erro, mas não tive sucesso. Usei o Linux Mint 13, que por padrão instala o R 2.14. Uma vez que essa versão do Mint é baseada no Ubuntu 12.04 (Precise Pangolin), segui as instruções do CRAN, atualizando para o R 3.2.5 com sucesso. <https://cran.r-project.org/bin/linux/ubuntu/README.html> Versões do R > 3.3.0 requerem Ubuntu > 12.04. Sendo o Mint baseado em versões iguais ou inferiores à 12.04, R > 3.2.5 não irá funcionar. ================================================ Éder Comunello Researcher at Brazilian Agricultural Research Corporation (Embrapa) DSc in Agricultural Systems Engineering (USP/Esalq) MSc in Environ. Sciences (UEM), Agronomist (UEM) --- Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00 > > > > > > Em 23 de agosto de 2016 22:18, Sergio Pinto via R-br < > r-br(a)listas.c3sl.ufpr.br> escreveu: > >> Boa noite Vinícius. >> Continuo na mesma. >> >> Att >> >> Sergio Mefiros Pinto >> Doutorando em Ciências Médicas >> Prof. de Educação Física >> >> Em 23/08/2016 21:04, "Vinícius Lionel Mateus" <vinynegrelli(a)gmail.com> >> escreveu: >> >>> Olá Sergio, >>> >>> Alguma novidade sobre seu problema? >>> >>> >>> Att, >>> Viny >>> >>> On 08/22/2016 07:43 PM, Sergio Pinto via R-br wrote: >>> >>> Prezado Éder e amigos listeiros. >>> Segui as instruções do link mas infelizmente o problema continua. >>> Se serve como ajuda o problema de fato começou comigo tentando atualizar >>> o R. >>> >>> Agradeço a ajuda >>> Att >>> Sergio Medeiros Pinto >>> Doutorando em Ciências Médicas >>> Prof. de Educação Física >>> >>> Em 22 de agosto de 2016 09:29, Éder Comunello <comunello.eder(a)gmail.com> >>> escreveu: >>> >>>> Sergio, bom dia! >>>> >>>> Encontrei a descrição de um erro similar no Ubuntu, causado pela >>>> existência de duas instalações do R, uma via source e outra via apt-get. Se >>>> lhe for útil, segue o link: >>>> http://grokbase.com/t/r/r-sig-debian/117nw5nch6/getting-conf >>>> used-with-two-versions-of-r%E2%80%8F%E2%80%8F >>>> >>>> >>>> ================================================ >>>> Éder Comunello >>>> Researcher at Brazilian Agricultural Research Corporation (Embrapa) >>>> DSc in Agricultural Systems Engineering (USP/Esalq) >>>> MSc in Environ. Sciences (UEM), Agronomist (UEM) >>>> --- >>>> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| >>>> ================================================ >>>> GEO, -22.2752, -54.8182, 408m >>>> UTC-04:00 / DST: UTC-03:00 >>>> >>>> >>>> >>>> >>>> Em 19 de agosto de 2016 16:04, Sergio Pinto via R-br < >>>> r-br(a)listas.c3sl.ufpr.br> escreveu: >>>> >>>>> Prezados >>>>> Essa é a saída que eu obtenho quando entro no R pelo terminal >>>>> >>>>> sergio@sergio-Infoway ~ $ R >>>>> /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file >>>>> or directory >>>>> /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file >>>>> or directory >>>>> >>>>> Acessei o /usr/local/bin/R , fui na linha 236 e aparece o seguinte: >>>>> . "${R_HOME}/etc${R_ARCH}/ldpaths" >>>>> >>>>> já na linha 259 aparece: >>>>> exec "${R_binary}" ${args} "${@}" >>>>> >>>>> Tem solução meu problema? >>>>> >>>>> Sergio Medeiros Pinto >>>>> Doutorando em Ciências Médicas UERJ >>>>> Prof. de Educação Física >>>>> >>>>> >>>>> Em 18 de agosto de 2016 19:35, Sergio Pinto <smp10bmc(a)gmail.com> >>>>> escreveu: >>>>> >>>>>> Éder boa noite! >>>>>> >>>>>> Apliquei os comandos que vocẽ me sugeriu ( um de cada vez). >>>>>> O problema persiste. >>>>>> Eu tenho um ícone do R no meu menu, mas ele não abre o programa. >>>>>> Quando digito R no terminal continua aparecendo: >>>>>> >>>>>> sergio@sergio-Infoway ~ $ R >>>>>> /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such >>>>>> file or directory >>>>>> /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file >>>>>> or directory >>>>>> >>>>>> Será que devo desinstalar o Linux Mint e instalar tudo do zero >>>>>> novamente? >>>>>> >>>>>> Sergio Medeiros Pinto >>>>>> Doutorando em Ciências Médicas UERJ >>>>>> Prof. de Educação Física >>>>>> >>>>>> Em 18 de agosto de 2016 17:20, Éder Comunello < >>>>>> comunello.eder(a)gmail.com> escreveu: >>>>>> >>>>>>> Sergio, boa tarde! >>>>>>> >>>>>>> Você instalou via apt-get? Já verificou se existem dependências? >>>>>>> >>>>>>> sudo apt-get update >>>>>>> sudo apt-get build-dep r-base >>>>>>> sudo apt-get install r-base r-recommended >>>>>>> >>>>>>> >>>>>>> ================================================ >>>>>>> Éder Comunello >>>>>>> Researcher at Brazilian Agricultural Research Corporation (Embrapa) >>>>>>> DSc in Agricultural Systems Engineering (USP/Esalq) >>>>>>> MSc in Environ. Sciences (UEM), Agronomist (UEM) >>>>>>> --- >>>>>>> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| >>>>>>> ================================================ >>>>>>> GEO, -22.2752, -54.8182, 408m >>>>>>> UTC-04:00 / DST: UTC-03:00 >>>>>>> >>>>>>> >>>>>>> >>>>>>> >>>>>>> Em 18 de agosto de 2016 00:50, Sergio Pinto via R-br < >>>>>>> r-br(a)listas.c3sl.ufpr.br> escreveu: >>>>>>> >>>>>>>> Prezado Marcelo >>>>>>>> Obrigado pela ajuda. >>>>>>>> Infelizmente o problema persiste. >>>>>>>> Continua aparecendo a mesma saída no terminal. >>>>>>>> Att >>>>>>>> Sergio Medeiros Pinto >>>>>>>> Doutorando em Ciências Médicas UERJ >>>>>>>> Prof. de Educação Fisica >>>>>>>> >>>>>>>> Em 17/08/2016 21:50, "Marcelo Laia" <marcelolaia(a)gmail.com> >>>>>>>> escreveu: >>>>>>>> >>>>>>>>> Sergio, bom dia! >>>>>>>>> >>>>>>>>> Encontrei a descrição de um erro similar no Ubuntu, causado pela >>>>>>>>> existência de duas instalações do R, uma via source e outra via apt-get. Se >>>>>>>>> lhe for útil, segue o link: >>>>>>>>> http://grokbase.com/t/r/r-sig-debian/117nw5nch6/getting-conf >>>>>>>>> used-with-two-versions-of-r%E2%80%8F%E2%80%8F >>>>>>>>> >>>>>>>>> >>>>>>>>> ================================================ >>>>>>>>> Éder Comunello >>>>>>>>> Researcher at Brazilian Agricultural Research Corporation (Embrapa) >>>>>>>>> DSc in Agricultural Systems Engineering (USP/Esalq) >>>>>>>>> MSc in Environ. Sciences (UEM), Agronomist (UEM) >>>>>>>>> --- >>>>>>>>> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| >>>>>>>>> ================================================ >>>>>>>>> GEO, -22.2752, -54.8182, 408m >>>>>>>>> UTC-04:00 / DST: UTC-03:00 >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> >>>>>>>>> Em 19 de agosto de 2016 16:04, Sergio Pinto via R-br < >>>>>>>>> r-br(a)listas.c3sl.ufpr.br> escreveu: >>>>>>>>> >>>>>>>>>> Prezados >>>>>>>>>> Essa é a saída que eu obtenho quando entro no R pelo terminal >>>>>>>>>> >>>>>>>>>> sergio@sergio-Infoway ~ $ R >>>>>>>>>> /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such >>>>>>>>>> file or directory >>>>>>>>>> /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such >>>>>>>>>> file or directory >>>>>>>>>> >>>>>>>>>> Acessei o /usr/local/bin/R , fui na linha 236 e aparece o >>>>>>>>>> seguinte: >>>>>>>>>> . "${R_HOME}/etc${R_ARCH}/ldpaths" >>>>>>>>>> >>>>>>>>>> já na linha 259 aparece: >>>>>>>>>> exec "${R_binary}" ${args} "${@}" >>>>>>>>>> >>>>>>>>>> Tem solução meu problema? >>>>>>>>>> >>>>>>>>>> Sergio Medeiros Pinto >>>>>>>>>> Doutorando em Ciências Médicas UERJ >>>>>>>>>> Prof. de Educação Física >>>>>>>>>> >>>>>>>>>> >>>>>>>>>> Em 18 de agosto de 2016 19:35, Sergio Pinto <smp10bmc(a)gmail.com> >>>>>>>>>> escreveu: >>>>>>>>>> >>>>>>>>>>> Éder boa noite! >>>>>>>>>>> >>>>>>>>>>> Apliquei os comandos que vocẽ me sugeriu ( um de cada vez). >>>>>>>>>>> O problema persiste. >>>>>>>>>>> Eu tenho um ícone do R no meu menu, mas ele não abre o programa. >>>>>>>>>>> Quando digito R no terminal continua aparecendo: >>>>>>>>>>> >>>>>>>>>>> sergio@sergio-Infoway ~ $ R >>>>>>>>>>> /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No >>>>>>>>>>> such file or directory >>>>>>>>>>> /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such >>>>>>>>>>> file or directory >>>>>>>>>>> >>>>>>>>>>> Será que devo desinstalar o Linux Mint e instalar tudo do zero >>>>>>>>>>> novamente? >>>>>>>>>>> >>>>>>>>>>> Sergio Medeiros Pinto >>>>>>>>>>> Doutorando em Ciências Médicas UERJ >>>>>>>>>>> Prof. de Educação Física >>>>>>>>>>> >>>>>>>>>>> Em 18 de agosto de 2016 17:20, Éder Comunello < >>>>>>>>>>> comunello.eder(a)gmail.com> escreveu: >>>>>>>>>>> >>>>>>>>>>>> Sergio, boa tarde! >>>>>>>>>>>> >>>>>>>>>>>> Você instalou via apt-get? Já verificou se existem dependências? >>>>>>>>>>>> >>>>>>>>>>>> sudo apt-get update >>>>>>>>>>>> sudo apt-get build-dep r-base >>>>>>>>>>>> sudo apt-get install r-base r-recommended >>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>>> ================================================ >>>>>>>>>>>> Éder Comunello >>>>>>>>>>>> Researcher at Brazilian Agricultural Research Corporation >>>>>>>>>>>> (Embrapa) >>>>>>>>>>>> DSc in Agricultural Systems Engineering (USP/Esalq) >>>>>>>>>>>> MSc in Environ. Sciences (UEM), Agronomist (UEM) >>>>>>>>>>>> --- >>>>>>>>>>>> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| >>>>>>>>>>>> ================================================ >>>>>>>>>>>> GEO, -22.2752, -54.8182, 408m >>>>>>>>>>>> UTC-04:00 / DST: UTC-03:00 >>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>>> Em 18 de agosto de 2016 00:50, Sergio Pinto via R-br < >>>>>>>>>>>> r-br(a)listas.c3sl.ufpr.br> escreveu: >>>>>>>>>>>> >>>>>>>>>>>>> Prezado Marcelo >>>>>>>>>>>>> Obrigado pela ajuda. >>>>>>>>>>>>> Infelizmente o problema persiste. >>>>>>>>>>>>> Continua aparecendo a mesma saída no terminal. >>>>>>>>>>>>> Att >>>>>>>>>>>>> Sergio Medeiros Pinto >>>>>>>>>>>>> Doutorando em Ciências Médicas UERJ >>>>>>>>>>>>> Prof. de Educação Fisica >>>>>>>>>>>>> >>>>>>>>>>>>> Em 17/08/2016 21:50, "Marcelo Laia" <marcelolaia(a)gmail.com> >>>>>>>>>>>>> escreveu: >>>>>>>>>>>>> >>>>>>>>>>>>> http://stackoverflow.com/questions/20967256/error-on-package >>>>>>>>>>>>> -and-r-paths-and-installations >>>>>>>>>>>>> >>>>>>>>>>>>> On 17/08/16 at 08:53pm, Sergio Pinto via R-br wrote: >>>>>>>>>>>>> > Prezados >>>>>>>>>>>>> > tenho um computador 32 bits com linux mint e instalei o >>>>>>>>>>>>> pacote r-base. >>>>>>>>>>>>> > O ícone do R foi criado mas não abre o programa. >>>>>>>>>>>>> > No terminal aparece a seguinte mensagem: >>>>>>>>>>>>> > >>>>>>>>>>>>> > sergio@sergio-Infoway ~ $ R >>>>>>>>>>>>> > /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No >>>>>>>>>>>>> such file or >>>>>>>>>>>>> > directory >>>>>>>>>>>>> > /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No >>>>>>>>>>>>> such file or >>>>>>>>>>>>> > directory >>>>>>>>>>>>> > >>>>>>>>>>>>> > >>>>>>>>>>>>> > O meu sistema é >>>>>>>>>>>>> > Linux sergio-Infoway 3.13.0-39-generic #66-Ubuntu SMP Tue >>>>>>>>>>>>> Oct 28 13:31:23 >>>>>>>>>>>>> > UTC 2014 i686 i686 i686 GNU/Linux >>>>>>>>>>>>> > >>>>>>>>>>>>> > Poderiam me ajudar? >>>>>>>>>>>>> > >>>>>>>>>>>>> > Desde já agradeço a atenção >>>>>>>>>>>>> > >>>>>>>>>>>>> > *Sergio Medeiros Pinto* >>>>>>>>>>>>> > Doutorando em Ciências Médicas UERJ >>>>>>>>>>>>> > Prof. de Educação Física >>>>>>>>>>>>> >>>>>>>>>>>>> > _______________________________________________ >>>>>>>>>>>>> > 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. >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> -- >>>>>>>>>>>>> Marcelo >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> >>>>>>>>>>>>> _______________________________________________ >>>>>>>>>>>>> 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. >>>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>>> >>>>>>>>>>> >>>>>>>>>> >>>>>>>>>> _______________________________________________ >>>>>>>>>> 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. >>>>>>>>>> >>>>>>>>> >>>>>>>>> >>> >>> _______________________________________________ >>> R-br mailing listR-br@listas.c3sl.ufpr.brhttps://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. >>> >>> >>> -- >>> Best regards, >>> >>> Vinícius Lionel Mateus, MSc (http://lattes.cnpq.br/6501001637020665) >>> Chemistry PhD Student >>> Atmospheric Chemistry Laboratory - Dep. Chemistry >>> PUC-Rio - Pontifical Catholic University of Rio de Janeiro >>> Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brazil CEP: 22451-900 >>> Phone: (+55) (21) 982 - 889 - 221 >>> Skype: vinicius.lionelhttp://www.qui.puc-rio.br/index.html >>> >>> >> _______________________________________________ >> 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

Relações entre variáveis
by Elias Carvalho 24 Ago '16

24 Ago '16

Boa Dia Pessoal Gostaria de esclarecer uma dúvida: quando tenho dados contínuos eu devo fazer uma análise da relação linear entre as variáveis, quando são categóricos a relação a ser analisada é considerada polinomial (me corrijam se estiver errado). Como eu ainda estou procurando entender alguns pressupostos da estatística, pergunto: 1. Qual das duas relações produzem resultados mais fiéis / próximos da realidade ? 2. Embora ocorra uma perda de informação, quando tenho dados mistos (categóricos e contínuos) em um data set a melhor opção seria discretizar os dados contínuos em categóricos ? Eu pergunto isso porque estou trabalhando com redes Bayesianas e sempre me deparo com dados mistos que podem ser modelados por redes Bayesianas condicionais que modelam a probabilidade condicional de variáveis categóricas como percentuais por exemplo (30% sim e 70% não) e variáveis contínuas com regressão linear. Porém neste tipo de rede Bayesiana nós discretos não podem ter pais contínuos e isso pode não representar a realidade, neste caso é preciso optar pela discretização para modelar as relações não lineares. -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho*

1 0

Re: [R-br] O R não abre no linux mint
by Sergio Pinto 23 Ago '16

23 Ago '16

Prezados Essa é a saída que eu obtenho quando entro no R pelo terminal sergio@sergio-Infoway ~ $ R /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file or directory /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file or directory Acessei o /usr/local/bin/R , fui na linha 236 e aparece o seguinte: . "${R_HOME}/etc${R_ARCH}/ldpaths" já na linha 259 aparece: exec "${R_binary}" ${args} "${@}" Tem solução meu problema? Sergio Medeiros Pinto Doutorando em Ciências Médicas UERJ Prof. de Educação Física Em 18 de agosto de 2016 19:35, Sergio Pinto <smp10bmc(a)gmail.com> escreveu: > Éder boa noite! > > Apliquei os comandos que vocẽ me sugeriu ( um de cada vez). > O problema persiste. > Eu tenho um ícone do R no meu menu, mas ele não abre o programa. > Quando digito R no terminal continua aparecendo: > > sergio@sergio-Infoway ~ $ R > /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file or > directory > /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file or > directory > > Será que devo desinstalar o Linux Mint e instalar tudo do zero novamente? > > Sergio Medeiros Pinto > Doutorando em Ciências Médicas UERJ > Prof. de Educação Física > > Em 18 de agosto de 2016 17:20, Éder Comunello <comunello.eder(a)gmail.com> > escreveu: > >> Sergio, boa tarde! >> >> Você instalou via apt-get? Já verificou se existem dependências? >> >> sudo apt-get update >> sudo apt-get build-dep r-base >> sudo apt-get install r-base r-recommended >> >> >> ================================================ >> Éder Comunello >> Researcher at Brazilian Agricultural Research Corporation (Embrapa) >> DSc in Agricultural Systems Engineering (USP/Esalq) >> MSc in Environ. Sciences (UEM), Agronomist (UEM) >> --- >> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| >> ================================================ >> GEO, -22.2752, -54.8182, 408m >> UTC-04:00 / DST: UTC-03:00 >> >> >> >> >> Em 18 de agosto de 2016 00:50, Sergio Pinto via R-br < >> r-br(a)listas.c3sl.ufpr.br> escreveu: >> >>> Prezado Marcelo >>> Obrigado pela ajuda. >>> Infelizmente o problema persiste. >>> Continua aparecendo a mesma saída no terminal. >>> Att >>> Sergio Medeiros Pinto >>> Doutorando em Ciências Médicas UERJ >>> Prof. de Educação Fisica >>> >>> Em 17/08/2016 21:50, "Marcelo Laia" <marcelolaia(a)gmail.com> escreveu: >>> >>> http://stackoverflow.com/questions/20967256/error-on-package >>> -and-r-paths-and-installations >>> >>> On 17/08/16 at 08:53pm, Sergio Pinto via R-br wrote: >>> > Prezados >>> > tenho um computador 32 bits com linux mint e instalei o pacote >>> r-base. >>> > O ícone do R foi criado mas não abre o programa. >>> > No terminal aparece a seguinte mensagem: >>> > >>> > sergio@sergio-Infoway ~ $ R >>> > /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file >>> or >>> > directory >>> > /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file >>> or >>> > directory >>> > >>> > >>> > O meu sistema é >>> > Linux sergio-Infoway 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 >>> 13:31:23 >>> > UTC 2014 i686 i686 i686 GNU/Linux >>> > >>> > Poderiam me ajudar? >>> > >>> > Desde já agradeço a atenção >>> > >>> > *Sergio Medeiros Pinto* >>> > Doutorando em Ciências Médicas UERJ >>> > Prof. de Educação Física >>> >>> > _______________________________________________ >>> > 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. >>> >>> >>> >>> -- >>> Marcelo >>> >>> >>> >>> _______________________________________________ >>> 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. >>> >> >> >

3 4

Identificando ordenamento não-exato
by Abel Brasil Ramos da Silva 23 Ago '16

23 Ago '16

Olá a todos! Gostaria de saber se existe uma forma de avaliar se dois vetores numéricos tendem a possuir uma ordenação próxima. Por exemplo, Ex 1: x=c(1,2,3,4,5,6,7,8,9,10) y=c(1,2,4,3,5,6,7,8,9,10) #Troquei o 3 pelo 4 nesse exemplo 1 eu consideraria que existe certa similaridade (não exata!) na ordenação desse vetores. Ex 2: x=c(1,3,2,5,7,4,10,6,8,9) y=c(1,2,3,4,5,7,8,10,6,9) nesse exemplo 2 também existe certa similaridade (digamos da distância máxima de 2 casas). No vetor x o 7 ocupa a 5º casa, já no y ocupa na 6ª casa, ou seja distância de 1 casa. O 8 ocupa a 9ª casa no vetor x, já no y ocupa a 7ª casa, ou seja, distância de 2 casas. Agradeço pela ajuda! Att.

1 0

Superfície de Resposta
by Emerson Cotta Bodevan 23 Ago '16

23 Ago '16

Caros colegas, bom dia. Estou estudando sobre superfície de resposta. Alguém poderia me indicar pacotes sobre esse assunto no R? Também, sendo possível, pacotes para planejamento e análises de experimentos. Agradeço qualquer sugestão e indicação. Att., *Emerson*

3 3

Exclusão de casos por critério temporal
by Paulo Dick 23 Ago '16

23 Ago '16

Companheiros, bom dia. Gostaria que me ajudassem na solução do seguinte problema: Tenho uma base que possui mais de um registro por indivíduo, que preciso fazer uma limpeza a partir das variáveis mês e ano. Preciso retirar apenas os registros que possuem intervalo de tempo menor que 3 meses entre eles, priorizando os registros mais antigos. Ou seja, se tenho 4 registros para um mesmo indivíduo, nos meses 1, 3, 5 e 6 (mesmo ano), fico com as linhas dos meses 1 e 5, e excluo a linha dos meses 3 e 6. Segue sintaxe para gerar fragmento extraído desta base: structure(list(identificador = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L, 5L, 5L, 6L, 6L, 7L, 7L, 8L, 8L, 8L, 9L, 9L, 10L, 10L), .Label = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10"), class = "factor"), ano = c(2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L, 2016L), mes = c(5L, 6L, 2L, 7L, 3L, 7L, 2L, 6L, 2L, 5L, 5L, 6L, 6L, 7L, 1L, 3L, 6L, 3L, 6L, 2L, 7L)), .Names = c("identificador", "ano", "mes"), row.names = c(NA, 21L), class = "data.frame") Agradeço pela ajuda. Abraços *Paulo Dick* Estatístico / Epidemiologia em Saúde Pública Tel.: (55 21) 99591-2716

3 7

Frequencia de linhas duplicadas
by Wagner Tassinari 22 Ago '16

22 Ago '16

Saudações pessoal, Tenho um banco de dados do seguinte formato: ID CODOM 23 100 45 100 34 101 20 100 10 102 11 102 23 100 Gostaria de construir uma variável de contagem da variável CODOM, ou seja,no CODOM 100 existem 4 registros, no CODOM 101 apenas 1 registro e no CODOM 102 existem 2 registros. Saída deverá ser assim: ID CODOM FREQ 23 100 4 45 100 4 34 101 1 20 100 4 10 102 2 11 102 2 23 100 4 Obrigado, ----- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Cel: (21) 96488-5982 (WhatsApp) Skype: wagner.tassinari wtassinari(a)gmail.com ----------------------------------------------------------------------

2 2

Teste de Aderência
by ana paula coelho madeira 22 Ago '16

22 Ago '16

Boa tarde. Tenho um conjunto de dados nos quais estou trabalhando com diferentes fdp's. Gostaria de aplicar diferentes testes de aderência, entre eles o Anderson Darling. No entanto, consegui apenas o ks.test. No pacote Nortest encontrei o ad.test para testar a normalidade, que não é o caso. Alguém tem algum material para me indicar? Desde já agradeço. Ana Paula

2 1

Teste T-pareado
by Simone D. Sartorio 22 Ago '16

22 Ago '16

Pessoal, boa noite! Surgiu uma dúvida conceitual, que não estou conseguindo resolver. Estava preparando minhas aulas e ao conferir um exercício no R que inicialmente tinha resolvido com esses comandos (Solução 1) observei que pelo p-value aceitaria a H0: não há diferença entre as médias. Contudo, pela minha resolução manual, deveria rejeitar H0. Logo, propus a solução 2, exatamente como realizei manualmente e dai, os resultados foram iguais. Gostaria de saber, qual é a diferença da solução 1 para a solução 2? Para mim, a sintaxe de ambas soluções deveria retornar a mesma conclusão. Problema:Uma distribuidora de combustíveis deseja verificar se um novo tipo de gasolina é eficaz na revitalização de motores velhos.Selecionou-se 12 automóveis de um mesmo modelo com mais de 8 anos de uso e, após regulagem dos motores, verifica-se o consumo de combustível. Em seguida, o carro é abastecido com o novo tipo de combustível durante 15 semanas e uma nova aferição é feita (km/litro). Como o desempenho dos automóveis foi medido antes e depois das 15 semanas, é razoável assumir que exista alguma dependência entre as variáveis. Gostaria de testar de realizar o teste de hipótese ao nível de 5%:H0: mD = 0 (O novo combustível não aumenta o rendimento) Ha: mD > 0 (O novo combustível aumenta o rendimento), onde mD é a verdade média da diferença do rendimento apos e antes. Aqueles que puderem me ajudar a entender, agradeço desde já!!!att.Simone ### Solução 1 consumo<- c(8.1, 7.9, 6.8, 7.8, 7.6, 7.9, 5.7, 8.4, 8.0, 9.5, 8.0, 6.8, 11.6, 8.8, 9.9, 9.5, 11.6, 9.1, 10.6, 10.8, 13.4, 10.6, 10.5, 11.4)momento<- factor(c(rep(‘Antes’,12), rep(‘Apos’,12))); momento t.test(consumo ~ momento, paired = T, conf.level = 0.95, alternative=‘greater’) # Paired t-test##data: consumo by momento#t = -6.5396, df = 11, p-value = 1#alternative hypothesis: true difference in means is greater than 0#95 percent confidence interval:# -3.7495 Inf#sample estimates:#mean of the differences # -2.941667 ### Solução 2X <- c(8.1, 7.9, 6.8, 7.8, 7.6, 7.9, 5.7, 8.4, 8.0, 9.5, 8.0, 6.8)# antesY <- c(11.6, 8.8, 9.9, 9.5, 11.6, 9.1, 10.6, 10.8, 13.4, 10.6, 10.5, 11.4)# aposD <- Y-X # Diferença > t.test(D, paired = F, conf.level = 0.95, alternative='greater')## One Sample t-test##data: D#t = 6.5396, df = 11, p-value = 2.097e-05#alternative hypothesis: true mean is greater than 0#95 percent confidence interval:# 2.133833 Inf#sample estimates:#mean of x # 2.941667 ***---***---***---***---***---***---***---***---***---***---***---***---**** Simone Daniela Sartorio* Professora Adjunta II no CCA/UFSCar, Araras/SP.* Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP.* Licenciada em Matemática - UNESP/Rio Claro.* Contatos: * tel: (19)3543-2664* cel: (19)98182-0586 - Tim***---***---***---***---***---***---***---***---***---***---***---***---*** Tenha um bom dia! ;)

4 3

Re: [R-br] Teste T-pareado
by Ivan Bezerra Allaman 22 Ago '16

22 Ago '16

A resposta é simples! Você está testando a hipótese erra na solução 1. Perceba no output, que a diferença está negativa, ou seja, o teste está fazendo ANTES - DEPOIS. Logo, se nós esperamos que DEPOIS melhore o desempenho, fazendo mais km/litro, então a DIFERENÇA é negativa e a hipótese a ser testada é: H0: md >= 0 Ha: md < 0 That is it!

2 1

erro no loop
by Edimeire Alexandra Pinto 22 Ago '16

22 Ago '16

Oi Pessoal. Como não achei um pacote que trabalhasse com previsões dinâmicas,estou eu mesma tendo de montar isso. Já dei alguns passos, agora resta entender“for” de quando coloco for(j in 1:20), tendo em vista que o R diz que os índicesestão fora do limite, sendo que na realidade dava para realizar as previsões.Além disso, o tamanho das previsões, em predict, era para ser do tamanho doloop de 20+j+1 até 84 (total de observações) e isso não ocorre... Suponha que eu tenha a base de dados Canada (é só copiar ecolar os comandos), meu objetivo rodarvários modelos dynlm para cada amostra que começa em 1 e vai de 20+j. Asprevisões de cada modelo rodado devem ser de 20+j+1 até 76, pois a janela ésempre fixa de tamanho 8.Gente, qualquer ajuda é muito bem vinda. require(vars)require(dynlm)data(Canada)colnames(Canada)summary(Canada)for(j in 1:20){ deflator_prev<-predict(dynlm(e-prod~ I(L(e,1)-rw)+ I(U-rw) +L(U,1),data=ts(Canada[(1:20+j),],start=c(1980,01),freq=4)), newdata=ts(Canada[(20+j+1:20+j+8),],start=c(1985,02),freq=4 )) } deflator_prev

1 0

Gerar series sintéticas de precipitação no R
by Alisson Lucrécio 22 Ago '16

22 Ago '16

Caro Colegas, Bom dia. Algum colega do tem um script para geração de séries sintéticas de precipitação no R. Baseado nessa metodologia. Segue o link. https://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source= web&cd=1&cad=rja&uact=8&ved=0ahUKEwjTy-XKgtXOAhWIgZAKHWO3BA8QFggcMAA& url=ftp%3A%2F%2Fftp.ufv.br%2Fdea%2FGPRH%2Fteses%2Fds_ oliveira%2FVicente%2520de%2520Paulo%2520Santos%2520de% 2520Oliveira%2520-%2520Tese%2520DS.doc&usg=AFQjCNEWA- 2QE2NG05QQWr0sUUbG5Csw4g&bvm=bv.129759880,d.Y2I Muito obrigado. -- Alisson Lucrecio da Costa

1 0

Geração de series sintéticas de precipitação no R
by Alisson Lucrécio 22 Ago '16

22 Ago '16

Caro Colegas, Bom dia. Algum colega do tem um script para geração de séries sintéticas de precipitação no R. Baseado nessa metodologia. Segue o link. https://www.google.com.br/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&cad=rja&… Muito obrigado. -- Alisson Lucrecio da Costa

1 0

previsão dinâmica
by Edimeire Alexandra Pinto 19 Ago '16

19 Ago '16

OI gente. Preciso de uma ajudinha. Sei que parece bobo, mas nãoconsigo estender a janela de previsões sequenciais com largura de , com largura fixa de tamanho 8.Criar um modelo com sample de 1:70 e prever 71 até 78,depois reestimar a equação com sample de 1:72 e previsão de 73:80, reestimar omodelo de 1:73 e previsão de 74:83 eassim por adiante até fechar o range em 84, última observação.Vejamos um data frame já existente no R como nome de Canadaque são séries que vão do primeiro trimestre de 1980 ao último trimestre de2004:require("dynlm") library("vars") data("Canada") summary(Canada) class(Canada)colnames(Canada) suponha que eu definaque minha amostra para estimar o modelo com sequencia em j que vai até 10,assim teríamos um modelo estimado de 01/1980 a 3/1996, depois outros modelos estimado de 01/1980a 4/1996, depoismais outro modelo estimado de 01/1980 a 1/1997, assim por adiante até 4/1998, ou seja, paraem 4/1998. No R seria mais ou menos assim:eq<-list() for(j in1:10){ eq(i)<-dynlm(e~L(e,1)+prod+rw+u,data=window(Canada,start=c(1980,1),end=c(1996,2),freq=4,extend=j)) o problema é que a previsão tem de ser sequencial tanto nostart quanto no end e com tamanho de fixo de 8, e também com uma janela que vaiindo sequencialmente em j de 1 a 18. Por exemplo, a previsão teria de serassim: a primeira de 04/1996 a 4/1998 (observe que o primeiro modelo termina em3/1996, conforme falei antes), depois 01/1997 a 2/1999 (observe que o primeiromodelo termina em 1/1997, conforme falei antes), depois faria outra previsão de2/1997 a 4/1999 e assim por adiante até chegar na última observação que é04/2000. Note que aumentamos sempre de 8+j em que j vai de 1 a 18. No R, não consigo desenvolver, pois seria algo mais ou menosassim: ## sei que está errado, mas é só para ter ideiafor( i in1:18){ predict(eq(i), newdata=window(Canada,start=c(1996,4+i),end=c(1998,4+i+8),freq=4,extend=i) ## 8 é o tamanho fixo da janela de previsão Se tiver ficado confuso, desculpe-me, mas quis apenasmostrar que quero um previsão dinâmica e não estática, tendo e, vista que aequação que sempre estimada à medida que o tamanho da amostra vai variando.Qualquer dica é válida!Obrigada, gente.

3 3

Re: [R-br] Digest R-br, volume 68, assunto 17
by ana paula coelho madeira 19 Ago '16

19 Ago '16

Prezado Alexandro, Os comandos em questão tentam apenas a normalidade. Tenho interesse em ver ficar a aderência de uma distribuição qualquer. Agradeço > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 68, assunto 17 > To: r-br(a)listas.c3sl.ufpr.br > Date: Wed, 17 Aug 2016 12:00:02 -0300 > > 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://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > 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. previsão dinâmica (Edimeire Alexandra Pinto) > 2. Re: Anderson Darling (Alexandro (Yahoo)) > 3. OFF TOPIC - Banco de dados (Andre Oliveira) > 4. Re: OFF TOPIC - Banco de dados (luc.souza(a)usp.br) > 5. RES: OFF TOPIC - Banco de dados (Leonard Mendonça de Assis) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 16 Aug 2016 21:53:38 +0000 (UTC) > From: Edimeire Alexandra Pinto <economatistica(a)yahoo.com.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] previsão dinâmica > Message-ID: > <228175548.15741202.1471384418129.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > OI gente. > > Preciso de uma ajudinha. Sei que parece bobo, mas nãoconsigo estender a janela de previsões sequenciais com largura de , com largura fixa de tamanho 8.Criar um modelo com sample de 1:70 e prever 71 até 78,depois reestimar a equação com sample de 1:72 e previsão de 73:80, reestimar omodelo de 1:73 e previsão de 74:83 eassim por adiante até fechar o range em 84, última observação.Vejamos um data frame já existente no R como nome de Canadaque são séries que vão do primeiro trimestre de 1980 ao último trimestre de2004:require("dynlm") > > library("vars") > > data("Canada") > > summary(Canada) > > class(Canada)colnames(Canada) > > > > suponha que eu definaque minha amostra para estimar o modelo com sequencia em j que vai até 10,assim teríamos um modelo estimado de 01/1980 a 3/1996, depois outros modelos estimado de 01/1980a 4/1996, depoismais outro modelo estimado de 01/1980 a 1/1997, assim por adiante até 4/1998, ou seja, paraem 4/1998. No R seria mais ou menos assim:eq<-list() > > for(j in1:10){ > > eq(i)<-dynlm(e~L(e,1)+prod+rw+u,data=window(Canada,start=c(1980,1),end=c(1996,2),freq=4,extend=j)) > > o problema é que a previsão tem de ser sequencial tanto nostart quanto no end e com tamanho de fixo de 8, e também com uma janela que vaiindo sequencialmente em j de 1 a 18. Por exemplo, a previsão teria de serassim: a primeira de 04/1996 a 4/1998 (observe que o primeiro modelo termina em3/1996, conforme falei antes), depois 01/1997 a 2/1999 (observe que o primeiromodelo termina em 1/1997, conforme falei antes), depois faria outra previsão de2/1997 a 4/1999 e assim por adiante até chegar na última observação que é04/2000. Note que aumentamos sempre de 8+j em que j vai de 1 a 18. No R, não consigo desenvolver, pois seria algo mais ou menosassim: > > ## sei que está errado, mas é só para ter ideiafor( i in1:18){ > > predict(eq(i), newdata=window(Canada,start=c(1996,4+i),end=c(1998,4+i+8),freq=4,extend=i) ## 8 é o tamanho fixo da janela de previsão > > Se tiver ficado confuso, desculpe-me, mas quis apenasmostrar que quero um previsão dinâmica e não estática, tendo e, vista que aequação que sempre estimada à medida que o tamanho da amostra vai variando.Qualquer dica é válida!Obrigada, gente. > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160816/39b43423/atta…> > > ------------------------------ > > Message: 2 > Date: Tue, 16 Aug 2016 22:04:57 +0000 (UTC) > From: "Alexandro (Yahoo)" <vl.alexandro(a)yahoo.com.br> > To: ana paula coelho madeira <apcmadeira(a)hotmail.com>, a lista > Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Anderson Darling > Message-ID: > <1102947158.17317609.1471385097903.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > #Segue um exemplo!require(ADGofTest)#?ad.test > set.seed( 123 )x <- runif( 100 ) > ad.test( x )$p.valuead.test( x, pnorm, 0, 1 )$p.value #Espero ter ajudado! > Alexandro Vieira Lopes > > > > De: ana paula coelho madeira via R-br <r-br(a)listas.c3sl.ufpr.br> > Para: lista R <r-br(a)listas.c3sl.ufpr.br> > Enviadas: Segunda-feira, 15 de Agosto de 2016 15:28 > Assunto: [R-br] Anderson Darling > > Boa tarde. > Tenho um conjunto de dados nos quais estou trabalhando com diferentes fdp's. Gostaria de aplicar diferentes testes de aderência, entre eles o Anderson Darling. No entanto, consegui apenas o ks.test. No pacote Nortest encontrei o ad.test para testar a normalidade, que não é o caso. Alguém tem algum material para me indicar? > Desde já agradeço. > Ana Paula > _______________________________________________ > 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� c�igo m�imo reproduz�el. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160816/268f052c/atta…> > > ------------------------------ > > Message: 3 > Date: Wed, 17 Aug 2016 03:21:29 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] OFF TOPIC - Banco de dados > Message-ID: > <811409596.15775329.1471404089455.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Boa noite,tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > obrigado > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/1a282242/atta…> > > ------------------------------ > > Message: 4 > Date: Wed, 17 Aug 2016 06:47:50 -0300 (BRT) > From: luc.souza(a)usp.br > To: Andre Oliveira <andreolsouza(a)yahoo.com.br>, "a lista Brasileira > oficial de discussão do programa R." <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] OFF TOPIC - Banco de dados > Message-ID: <3E69A4AD-23EB-49AC-B77D-2914AD7D210C(a)usp.br> > Content-Type: text/plain; charset="utf-8" > > Na conab tem uma série histórica de vários produtos e no ibge Tb. > > Enviado do meu iPhone > > > Em 17 de ago de 2016, às 00:21, Andre Oliveira via R-br <r-br(a)listas.c3sl.ufpr.br> escreveu: > > > > Boa noite, > > tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > > obrigado > > > > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > > _______________________________________________ > > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/6a558ec2/atta…> > > ------------------------------ > > Message: 5 > Date: Wed, 17 Aug 2016 09:37:57 -0300 > From: Leonard Mendonça de Assis <assis.leonard(a)gmail.com> > To: 'a lista Brasileira oficial de discussão do programa R.' > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] RES: OFF TOPIC - Banco de dados > Message-ID: <00ad01d1f884$30927090$91b751b0$(a)gmail.com> > Content-Type: text/plain; charset="utf-8" > > André, eu sei que o CEPEA/USP disponibiliza séries temporais com a cotação de diversos produtos agrícolas. Área plantada eu desconheço algum órgão que possa ter isso. Uma tentativa seria ministério da Agricultura ou censo agropecuário do ibge > > > > De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Andre Oliveira via R-br > Enviada em: quarta-feira, 17 de agosto de 2016 00:21 > Para: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Assunto: [R-br] OFF TOPIC - Banco de dados > > > > Boa noite, > > tenho um projeto com aluno de ensino médio sobre produzir interface em shiny para visualizar dados. É uma feira onde o tema é alimentação e pensei em explorar banco de dados de produção de milho, soja ou outro produtos por hectare ao longo de anos no Brasil, via shiny. Alguém saberia dizer onde posso encontrar banco de dados disponíveis com varias culturas e vários anos? Este aluno é do técnico em informática e vem trabalhando com estatística onde há uma integração entre as áreas. > > > > obrigado > > > > > > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160817/bcb2784f/atta…> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 68, assunto 17 > *****************************************

2 1

Re: [R-br] O R não abre no linux mint
by Sergio Pinto 18 Ago '16

18 Ago '16

Prezado Marcelo Obrigado pela ajuda. Infelizmente o problema persiste. Continua aparecendo a mesma saída no terminal. Att Sergio Medeiros Pinto Doutorando em Ciências Médicas UERJ Prof. de Educação Fisica Em 17/08/2016 21:50, "Marcelo Laia" <marcelolaia(a)gmail.com> escreveu: http://stackoverflow.com/questions/20967256/error-on- package-and-r-paths-and-installations On 17/08/16 at 08:53pm, Sergio Pinto via R-br wrote: > Prezados > tenho um computador 32 bits com linux mint e instalei o pacote r-base. > O ícone do R foi criado mas não abre o programa. > No terminal aparece a seguinte mensagem: > > sergio@sergio-Infoway ~ $ R > /usr/local/bin/R: line 236: /usr/local/lib/R/etc/ldpaths: No such file or > directory > /usr/local/bin/R: line 259: /usr/local/lib/R/bin/exec/R: No such file or > directory > > > O meu sistema é > Linux sergio-Infoway 3.13.0-39-generic #66-Ubuntu SMP Tue Oct 28 13:31:23 > UTC 2014 i686 i686 i686 GNU/Linux > > Poderiam me ajudar? > > Desde já agradeço a atenção > > *Sergio Medeiros Pinto* > Doutorando em Ciências Médicas UERJ > Prof. de Educação Física > _______________________________________________ > 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. -- Marcelo

2 1

Girar rótulo da lattice
by Andre Oliveira 18 Ago '16

18 Ago '16

Olá,como posso girar os rótulos da lattice? par(las=1) não deu certo! require(lattice) par(las=1) xyplot(GP~factor(Ano)|S,pch=19, ylab="Ganho de peso (g)",xlab="", main="Ganho de peso por sexo") André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

1 0