Maio 2015 - R-br - DINF UFPR

Instalação - Ubuntu 64 bits
by Andre Oliveira 10 Mai '15

10 Mai '15

Pessoal apesar de não usar sempre instalo o Rcmdr para apresentar a quem tenha resistência ao uso do R. O fato de ser 64 meu linux e R pode haver problemas na instalação do Rcmdr? 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

Re: [R-br] Fatorial tiplo com dados adicionais
by Tiago Souza Marçal 10 Mai '15

10 Mai '15

Te indico o pacote ExpDes se você quer uma rotina simples. Att. Tiago. Alisson Lucrécio <alisson.lucrecio(a)ifgoiano.edu.br> escreveu: Caro Colegas da lista r-br, Preciso de uma ajudar para montar uma análise de variância para um fatorial tripo (Solo x Solução de Extração x Concentração) e um fator adicional água. Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 níveis de solo e 1 nível de solução (Água). Segue um exemplo abaixo. fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, 0.50))) Rep <- apply(data.frame(Solo = rep(1:5, each = 18), Sol = rep(rep(c(0.5, 1), each = 9), 5), Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), 1, sum) fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), Sol = rep("Água", 15)) Rep <- rep(1:5, each = 3) adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) Muito 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.

1 0

Fatorial tiplo com dados adicionais
by Alisson Lucrécio 09 Mai '15

09 Mai '15

Caro Colegas da lista r-br, Preciso de uma ajudar para montar uma análise de variância para um fatorial tripo (Solo x Solução de Extração x Concentração) e um fator adicional água. Os dados fatoriais são compostos por 5 níveis de solos, 2 níveis de solução e 3 níveis de concentração. Os dados adicionais são compostos por 5 níveis de solo e 1 nível de solução (Água). Segue um exemplo abaixo. fat_data <- data.frame(Solo = gl(5, 18, la = LETTERS[1:5]), Sol = gl(2, 9, length = 90, la = c("PO4", "SO4")), Conc = gl(3, 3, length = 90, la = c(0.10, 0.25, 0.50))) Rep <- apply(data.frame(Solo = rep(1:5, each = 18), Sol = rep(rep(c(0.5, 1), each = 9), 5), Conc = rep(rep(c(0.4, 0.9, 1.5), each = 3), 10)), 1, sum) fat_data$Rep <- Rep + rnorm(length(Rep),0,0.5) adi_data <- data.frame(Solo = gl(5, 3, la = LETTERS[1:5]), Sol = rep("Água", 15)) Rep <- rep(1:5, each = 3) adi_data$Rep <- Rep + rnorm(length(Rep),0,0.5) Muito obrigado. -- Alisson Lucrecio da Costa

1 0

Re: [R-br] Emparelhamento de arquivos
by Éder Comunello 09 Mai '15

09 Mai '15

Desculpe, Mauro! Troque xts2 por xts3: merge.xts(xts1, xts2, join="inner") # Data: # numeric(0) # # Index: # integer(0) merge.xts(xts1, xts3, join="inner") # Open High Low Close Open.1 High.1 Low.1 Close.1 # 2007-01-04 50.42096 50.42096 50.26414 50.33236 50.42096 50.42096 50.26414 50.33236 # 2007-01-05 50.37347 50.37347 50.22103 50.33459 50.37347 50.37347 50.22103 50.33459 # 2007-01-06 50.24433 50.24433 50.11121 50.18112 50.24433 50.24433 50.11121 50.18112 Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 8 de maio de 2015 23:42, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > Por que eu não estou conseguindo rodar este seu exemplo? > > > > >

1 0

Emparelhamento de arquivos
by Edilson Santos 08 Mai '15

08 Mai '15

Caros usuários, Preciso emparelhar dados de dois arquivos "xts", por data, sendo que em cada um dos arquivos somente algumas datas são coincidentes. São essas que eu gostaria de emparelhadas. Alguém idéia como poderia fazer. Utilizando o código abaixo não tive sucesso: OS dois arquivos: Shanghai.xts e vale.xts ShVale=merge.xts(Shanghai.xts,vale.xts,all = FALSE,suffixes=NULL, join="outer",retside=TRUE,retclass="xts") ShVale = na.exclude(object=ShVale) o arquivo resultante (ShVale) tem não tem dados, embora existam datas coincidentes entre os dois arquivos originais. Edilson Floriano dos Santos

2 6

Re: [R-br] Obtendo efeitos INDIVIDUAIS de múltiplas variáveis qualitativas.
by Rafael Costa 08 Mai '15

08 Mai '15

Prezados colegas, Esclareço que antes de clicar em DOWNLOAD, por favor desmarquem a opção"Use our download manager and get recommended downloads" para baixar o arquivo diretamente sem problemas. Qualquer contribuição será muito bem vinda! Em 4 de maio de 2015 22:51, Rafael Costa <rafael_costa(a)caen.ufc.br> escreveu: > Mauro, obrigado pela atenção. Segue em anexo o arquivo. > > Em 4 de maio de 2015 21:18, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > >> Poderia me enviar este arquivo anexado, pois meu antivírus não deixa >> baixar! >> >> >> >> Prezados colegas, >> >> >> >> Possuo dados de um questionário e quero estimar o efeito individual de >> cada variável explicativa (todas são qualitativas) sobre a variável >> dependente (contínua). Contudo o padrão é considerar os coeficientes >> estimados como a diferença entre o o grupo de referência (valor estimado do >> intercepto) e o grupo do coeficiente. Cada variável qualitativa diz >> respeito a uma característica de uma determinada atividade e a variável >> contínua é o tempo tomado para executar tal atividade. Destaco que o nível >> de referência de cada fator diz respeito ao caso em que nenhuma das opções >> para aquele fator foi marcado. Os dados estão em " >> http://www.datafilehost.com/d/c7f0d342". Não os coloquei no CMR, pois >> ainda não consegui (peço minhas sinceras desculpas), mas espero que isto >> não seja empecilho.Também não coloquei apenas uma amostra dos dados, pois >> estava havendo problemas de matriz singular. >> >> >> >> 1º (e principal) Problema - Para obter o efeito individual dos níveis de >> cada fator, considerei que o grupo de referência possui efeito zero e fiz >> os seguintes passos: >> >> >> >> #Como o arquivo não está sendo carregado no CMR, supõe-se aqui que ele >> foi baixado da internet e já se encontra carregado no R. >> >> #Utilizarei uma regressão quantílica, por isso segue o pacote abaixo. >> >> install.packages(quantreg) >> >> library(quantreg) >> >> #Transformando as variáveis categóricas em objetos únicos >> >> p_1 = table(1:length(tabela1.1$p1),as.ordered(tabela1.1$p1)) >> >> p_21 = table(1:length(tabela1.1$p21),as.factor(tabela1.1$p21)) >> >> p_22 = table(1:length(tabela1.1$p22),as.factor(tabela1.1$p22)) >> >> p_23 = table(1:length(tabela1.1$p23),as.factor(tabela1.1$p23)) >> >> p_24 = table(1:length(tabela1.1$p24),as.factor(tabela1.1$p24)) >> >> p_25 = table(1:length(tabela1.1$p25),as.factor(tabela1.1$p25)) >> >> p_34 = table(1:length(tabela1.1$p34),as.ordered(tabela1.1$p34)) >> >> p_5 = table(1:length(tabela1.1$p5),as.ordered(tabela1.1$p5)) >> >> p_6 = table(1:length(tabela1.1$p6),as.ordered(tabela1.1$p6)) >> >> p_7 = table(1:length(tabela1.1$p7),as.ordered(tabela1.1$p7)) >> >> p_8 = table(1:length(tabela1.1$p8),as.ordered(tabela1.1$p8)) >> >> p_9 = table(1:length(tabela1.1$p9),as.ordered(tabela1.1$p9)) >> >> #Regressão do modelo sem intercepto, mas também considerando que o grupo >> de referência tem valor 0, haja vista que o grupo de referência significa >> que nenhum dos fatores foi marcado (se nenhum item foi marcado, considero >> que o tempo tomado para executar a atividade é praticamente zero). >> >> qrModel=rq(data=tabela1.1, pontoefetivo ~ 0 + p_1[,-1] + p_21[,-1] + >> p_22[,-1] + p_23[,-1] + p_24[,-1] + p_25[,-1] + p_34[,-1] + p_5[,-1] + >> p_6[,-1] + p_7[,-1] + p_8[,-1] + p_9[,-1], tau=0.5) >> >> summary(qrModel) >> >> >> >> A minha ideia foi que, dado que o grupo de referência é zero, o >> coeficiente estimado é justamente o efeito individual do nível escolhido da >> variável. A minha ideia está certa? Se não, o que faço para obter esses >> efeitos individuais? >> >> >> >> 2º Problema: Supondo que fiz tudo certinho acima, os fatores ordenados >> não possuem efeitos crescentes [Ver summary(qrModel)]. Mas não deveriam >> possuir? Se deveriam, o que faço para garantir tal efeito? >> >> >> >> 3º Problema: Mas uma vez supondo que tudo está correto, esperava que >> todos os efeitos sobre o tempo de execução das atividades fossem não >> negativos. Ou estou errado a respeito disso? Se estou certo, o que faço >> para garantir que todos os valores sejam não negativos. >> >> >> >> Desde já agradeço a todos e aguardo qualquer ajuda. >> <http://www.datafilehost.com/get.php?file=c7f0d342>-- >> >> >> >> >> *Rafael Carneiro da Costa Doutorando em Economia Centro de Pós-Graduação >> em Economia - CAEN Universidade Federal do Ceará - UFC* >> >> >> ------------------------------ >> [image: Avast logo] <http://www.avast.com/> >> >> Este email foi escaneado pelo Avast antivírus. >> www.avast.com >> >> > > > -- > > > > *Rafael Carneiro da CostaDoutorando em EconomiaCentro de Pós-Graduação em > Economia - CAENUniversidade Federal do Ceará - UFC* > -- *Rafael Carneiro da CostaDoutorando em EconomiaCentro de Pós-Graduação em Economia - CAENUniversidade Federal do Ceará - UFC*

1 0

Criar mapa com Image
by Wagner Wolff 08 Mai '15

08 Mai '15

Olá pessoal, quero gerar um mapa em que a variável em questão é obtida por meio de uma função, na qual as variáveis de entrada são dados espacializados em um grid. Entretanto estou me deparando com este erro: Error in image.default(x, y, z) : increasing 'x' and 'y' values expected Vocês teriam uma solução? Segue o código exemplo Qp <- function(mu,sigma,area,p){ Qp<-(exp(mu+(sigma*qnorm(1-p))))*area } grid <- pred_grid(range(Longitude),range(Latitude), by=0.5) gr <- locations.inside(grid, limite) x <- gr$Var1 y <- gr$Var2 z <- Qp(kc.mu.ANO$predict,kc.sigma.ANO$predict,1,0.9) image(x,y,z) contour(z, add=T, nlev=5)

2 3

R erro Rstudio
by Andre Oliveira 08 Mai '15

08 Mai '15

Pessoal,não uso Windows, mas hoje precisei de usar e encontrei este problema ao iniciar o Rstudio. Alguém poderia ajudar resolver este erro? É sobre permissão mas não sei resolver! R version 3.2.0 (2015-04-16) -- "Full of Ingredients" Copyright (C) 2015 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. Error installing package: Erro: ERROR: no permission to install to directory 'C:/Program Files/R/R-3.2.0/library' Error installing package: Erro: ERROR: no permission to install to directory 'C:/Program Files/R/R-3.2.0/library' 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

5 7

Banco de dados
by Andre Oliveira 07 Mai '15

07 Mai '15

Pessoal,gostaria de deixa aqui meus agradecimentos nas sugestões do grupo sobre banco de dados públicos e no problema do Rstudio não citarei ninguém para não cometer injustiça de esquecimento. No Rstudio de fato era permissão que o aluno da escola adicionou permissão daquela pasta do Rstudio. Só para constar voltarei para linux semana que vem hj foi algo pontual. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espírito Santo. IFES

1 0

Sequência e repetição de série temporal (minutos)
by Aline Lipsky 07 Mai '15

07 Mai '15

Boa tarde prezados colegas, Eu preciso criar um fator dena série temporal entre 0 a 10 minutos, com intervalo de 15 segundos e repetir essa sequência 21 vezes em ordem crescente. Para isso usei o seguinte script: library ( chron ) sort (rep (seq ( chron ( times = c("00:00:00") ), chron ( times = c("00:10:00") ), chron ( times = c("00:00:15"))), 21 )) Mas quando faço isso os valores que retornam não são o intervalo de tempo descrito.Consigo apenas criar a sequência sem repetição com o seguinte script: seq ( chron ( times = c ( "00:00:00" ) ), chron ( times = c ( "00:10:00" ) ), chron ( times = c ( "00:00:15" ))) Alguém saberia o que posso fazer? Devo usar outro pacote? Grata pela atenção. Atenciosamente, Aline Lipsky

2 3