- R-br - DINF UFPR

Intervalos de credibilidade e HPD
by Edmundo Caetano 27 Dez '12

27 Dez '12

Prezados, Estou com uma dúvida pretendo estabelecer intervalos de credibilidade ou HPD para alguns parâmetros. Alguém pode me ajudar. Abrços Edmundo Caetano "Ninguém é tão sábio que nada tenha para aprender, nem tão tolo que nada tenha para ensinar." Blaise Pascal ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quarta-feira, 26 de Dezembro de 2012 12:00 Assunto: Digest R-br, volume 22, assunto 35 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. Re: Script S100 revisão 1 (Rodrigo Sant'Ana) ---------------------------------------------------------------------- Message: 1 Date: Wed, 26 Dec 2012 00:15:09 -0200 From: "Rodrigo Sant'Ana" <rodrigo.gringo(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Script S100 revisão 1 Message-ID: <CACpsV-wF1SV6EKoBb4A6RgAeTG8MQG3XOBP8LQyA_Vb-4+Pr+w(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Prezado Mauro, se fizeres uma pesquisa nos últimos posts da lista, verás que o arquivo "recorte.txt" está anexado a um dos posts anteriores do Sr. Hélio. Ok? Feliz Natal a você e a todos da lista, Boas festas. cordialmente, _______________ Rodrigo Sant'Ana Oceanógrafo Em 24 de dezembro de 2012 04:59, Mauro Sznelwar <sznelwar(a)uol.com.br>escreveu: > ** > > Alguém teria o arquivo recorte.txt que faz parte do Script? > > Pessoal, > > seguindo a proposta do Sr. Benilton, criei um gist publico para > otimização desta discussão. > Assim podemos maximizar a troca de informações e experiências sobre estas > análises, e quaisquer > outras análises vindouras, sem criar uma rotina de emails cheios de anexos > para a lista da R-BR. > > Acho que a ideia do grupo de estudos super válida, e ela se tornará ainda > melhor se utilizarmos algo > como o gist para otimizar nosso trabalho. Fica a dica, e o link para > darmos prossegmento nesta > discussão. > > link: https://gist.github.com/4262669 > > Como primeiro passo, a titulo de exemplo, fiz uma pequena proposta de > alteração no script original > proposto pelo Sr. Hélio, simplesmente alterei a forma de carregar a base > de dados s100 prevista no > pacote geoR. > > Quem nunca utilizou o gist, é simples, para inserir novas alterações no > script, sejam novas análises ou mesmo > otimizações no código, basta clicar no botão "Revisions", após clicar em > "edit", realizar as inserções, > alterações, sugestões, etc, e clicar em "Update Gist". > Todas as revisões ficaram salvas e para visualiza-las bastará clicar no > botão "Revisions" novamente. > > Peço desculpas caso tenha sido muito invasivo, mas acredito que este seja > um caminho interessante para > adotarmos, assim como inúmeros outros exemplos sugeridos anteriormente. > > Boa noite a todos, > > _______________ > Rodrigo Sant'Ana > Oceanógrafo > > > > . > > > _______________________________________________ > 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/20121226/a3909095/atta…> ------------------------------ _______________________________________________ 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 22, assunto 35 *****************************************

1 0

Experimento de mistura
by Ari Clecius 24 Dez '12

24 Dez '12

Boa tarde, estou tentando realizar um planejamento de mistura de três componentes onde um dos componentes não pode passar de 0.05 ou 5%, estou usando o pacote qualityTools, algoém poderia me dizer como configurar os limites? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)87687786 (085)33669042

1 1

R para Linux
by Gilbert Queiroz 23 Dez '12

23 Dez '12

Pessoal, R para Linux, roda legal??? Abs.

4 3

Re: [R-br] R para Linux
by Sérgio Henrique 23 Dez '12

23 Dez '12

Eu gosto bastante do emacs! Enviado via iPad Em Dec 23, 2012, às 8:29 AM, Rodolfo Timoteo da Silva <zhushazang(a)gmail.com> escreveu: Além do RStudio (rstudio.org) tem o rkward (rkward.sourceforge.net) que podem te auxiliar no uso. Att Em 22-12-2012 20:14, FHRB Toledo escreveu: Roda! 2012/12/22 Gilbert Queiroz <gilbert_queiroz(a)yahoo.com.br> > Pessoal, > R para Linux, roda legal??? > Abs. > > _______________________________________________ > 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. _______________________________________________ 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

(sem assunto)
by Raul Carlos de Mello 22 Dez '12

22 Dez '12

Boa tarde a todos a lista. Estou fazendo leitura de dados do site da bovespa e melhores e maiores da exame e tenho o seguinte problema que gostaria de ajuda: - Faço a teste da existência do site com o comando; - Ok, tudo funciona sem problema - Após a verificação existência da página busco pelos dados contidos nas tabelas. Entretanto nem todas a tabelas existem dados. Como posso fazer, ao utilizar "readHTMLTable", para evitar que o programa seja interrompido em decorrência de um erro pela tabela não existir. -- Att, Prof. Adm. Raul Carlos de Mello http://lattes.cnpq.br/3784085036882680

1 0

Re: [R-br] Cálculo do percentage prediction error em ZINB
by Humberto Hazin 21 Dez '12

21 Dez '12

Opa Leonard, Muito obrigado vou dar uma olhada! Sim estou usando o pacote pscl para rodar os modelos ZINB! Se eu conseguir resolver o problema coloco na lista. Cara se não me engano tive esse problema também a um tempo atrás, veja se código abaixo para ver se resolve o teu problema. Só uma observação, com esse código não precisas rodar o pacote pscl. https://www.dropbox.com/s/8l480on83dvm3m7/zeroinflCorreto.R Um abraço, Humberto Enviado do Email do Windows De: Leonard de Assis Enviado: ‎19‎ de ‎dezembro‎ de ‎2012 ‎22‎:‎04 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Cálculo do percentage prediction error em ZINB Humberto vc está usando o pscl? Eu estou no momento com um problema envolvendo ZINB, só que o meu é um pouco mais chato, pois meu modelo é multivariado (E pra complicar, está dando que uma das dimensões é singular, hehehe) Sugiro ler os tutoriais, vc acha um monte via rseek. Esse aqui é um deles (De antemão, não tem exatamente o que vc precisa, mas tem muita coisa) http://www.ats.ucla.edu/stat/r/dae/zinbreg.htm []s Leonard de Assis http://about.me/ldeassis Em 19/12/2012 17:44, Humberto Hazin escreveu: Prezados, Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano Eu iniciei da seguinte forma: dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/v7vruztus11nnbi/try.csv original<-read.csv(”original.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv # Criar Fatores para o banco original original$area<-rep(1,nrow(original)) original$area[original$lat<=-15&original$lat>=-45]<-2 original$A=as.factor(original$area) summary(original) original$S<-factor(original$id.group.fleet) original$Q<-factor(original$Quarter) original$Y<-factor(original$Year) original$logEf<-log(original$effort1) summary(original) # Criar Fatores para o banco try try$area<-rep(1,nrow(try)) try$area[try$lat<=-15&try$lat>=-45]<-2 try$A=as.factor(try$area) summary(try) try$S<-factor(try$id.group.fleet) try$Q<-factor(try$Quarter) try$Y<-factor(try$Year) try$logEf<-try$effort1*0 summary(try) zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negbin",data=original) pred<-predict(zeroinfl.fit_original,newdata=try,type="response") Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error Alguem poderia me ajudar! Desde já Agradeço Humberto Enviado do Email do Windows _______________________________________________ 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.

4 4

Junção de dados no mclapply
by Junior Beleti 21 Dez '12

21 Dez '12

Bom dia a todos. Estou realizando experimentos para execução paralela. Para isso criei uma função que realiza a multiplicação de matrizes na sua forma digamos "primitiva", ou seja, como na matemática mesmo, realizando a multiplicação linha por coluna. Sei que existe o operador %*% que já realiza a multiplicação matricial, porém, como já disse antes, estou apenas realizando experimentos. Segue código: m1 <- matrix(0,3,3) #criando matriz m1[1,1]<-1 m1[1,2]<-2 m1[1,3]<-3 m1[2,1]<-4 m1[2,2]<-1 m1[2,3]<-4 m1[3,1]<-2 m1[3,2]<-3 m1[3,3]<-2 m2 <- matrix(0,3,3) #criando matriz m2[1,1]<-2 m2[1,2]<-1 m2[1,3]<-2 m2[2,1]<-3 m2[2,2]<-4 m2[2,3]<-1 m2[3,1]<-2 m2[3,2]<-3 m2[3,3]<-3 m3 <- matrix(0,3,3) #criando matriz numlinha<-3 numcoluna<-3 N<-3 #funcao que realiza a multiplicação myf = function(N,mat1,mat2,num,x,y){ num + (mat1[x,N] * mat2[N,y]) } #forma sequencial (já testada e que está correta) aux <-0 for(i in 1:numlinha){ for(j in 1:numcoluna){ for(k in 1:N){ aux <-myf(k,m1,m2,aux,i,j) } m3[i,j] <- aux aux <-0 } } ## resultado da execucao: ## matriz m1 multiplicada pela matriz m2 que gera a matriz m3 > m1 [,1] [,2] [,3] [1,] 1 2 3 [2,] 4 1 4 [3,] 2 3 2 > m2 [,1] [,2] [,3] [1,] 2 1 2 [2,] 3 4 1 [3,] 2 3 3 > m3 [,1] [,2] [,3] [1,] 14 18 13 [2,] 19 20 21 [3,] 17 20 13 Se eu fosse utilizar a função mclapply, ao invés das tres repetições (for), bastaria colocar: require(parallel) aux <- 0 for(i in 1:numlinha){ mclapply(1:numcoluna,function(h){ for(k in 1:N){ aux <-myf(k,m1,m2,aux,i,h) } m3[i,h] <- aux aux <-0 },mc.cores=2) } O problema é que não consigo armazenar os valores em cada posição da matriz (m3[i,h] <- aux), ou seja, ele não salva na matriz de destino m3, visto que isso ocorre dentro da chamado ao mclapply. Tentei utilizar o "unlist" mas não obtive sucesso. Alguem pode me auxiliar com isso? Att, Junior Beleti.

2 1

ler arquivo excel 2007 com o R (64bit)
by Alexandro (Yahoo) 20 Dez '12

20 Dez '12

#Saudações a todos da lista R-Br, utilizei o comando abaixo no R (64bit) para abrir um arquivo excel 2007 require(RODBC);cnct = odbcConnectExcel2007(file.choose()) #escolhi o arquivo de interesse e apareceu a mensagem: #Carregando pacotes exigidos: RODBC #Mensagens de aviso perdidas: #1: In odbcDriverConnect(con, tabQuote = c("[", "]"), ...) : # [RODBC] ERROR: state IM002, code 0, message [Microsoft][ODBC Driver Manager] Nome da fonte de dados não encontrado e nenhum driver padrão especificado #2: In odbcDriverConnect(con, tabQuote = c("[", "]"), ...) : # ODBC connection failed #Uso processador i5 4GB de memória e windows 7 (64bit). #Quando executei o mesmo comando no R (32bit), funcionou perfeitamente, isto é o que me deixa intrigado! #o que pode ter ocorrido e o que fazer para ter este problema solucionado? #será que este comando no R (64bit) pode não estar implementado? #Muito Obrigado ... #Alexandro

9 10

Cálculo do percentage prediction error em ZINB
by Humberto Hazin 19 Dez '12

19 Dez '12

Prezados, Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano Eu iniciei da seguinte forma: dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv try<-read.csv(”try.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/v7vruztus11nnbi/try.csv original<-read.csv(”original.csv”, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv # Criar Fatores para o banco original original$area<-rep(1,nrow(original)) original$area[original$lat<=-15&original$lat>=-45]<-2 original$A=as.factor(original$area) summary(original) original$S<-factor(original$id.group.fleet) original$Q<-factor(original$Quarter) original$Y<-factor(original$Year) original$logEf<-log(original$effort1) summary(original) # Criar Fatores para o banco try try$area<-rep(1,nrow(try)) try$area[try$lat<=-15&try$lat>=-45]<-2 try$A=as.factor(try$area) summary(try) try$S<-factor(try$id.group.fleet) try$Q<-factor(try$Quarter) try$Y<-factor(try$Year) try$logEf<-try$effort1*0 summary(try) zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negbin",data=original) pred<-predict(zeroinfl.fit_original,newdata=try,type="response") Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error Alguem poderia me ajudar! Desde já Agradeço Humberto Enviado do Email do Windows

2 1

problemas com R-br listas
by Fernando A. Souza 19 Dez '12

19 Dez '12

Caros colegas estou com um problema Gostaria de mudar a forma como eu recebo as postagem na lista. Geralmente recebo vários e-mails de todos os que postam. Recentemente um amigo fez o cadastro na lista e ele geralmente recebe um único e-mail com a listas das mensagens de um dado período. Gostaria de receber desta forma tambem. Pois evitaria a sobrecarca de minha caixa de e-mail. Algúem sabe como se faz isto? att -- Fernando Antônio de Souza Zootecnista, Msc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

3 3