- R-br - DINF UFPR

[Dúvida] Programação paralella e distribuida em R.
by Pedro Rafael 31 Ago '12

31 Ago '12

Pessoal estou começando a estudar programação paralela e distribuida no R. Gostaria de saber se estas tarefas podem ser feitas com o R e se os resultados obtidos são interessantes. Estou lendo um pouco sobre OpenMP e OpenMPI, alguem trabalha com essa API utilizando R? [ ], Pedro Rafael Diniz Marinho.

5 5

R na plataforma Debian
by Kate Dresh 31 Ago '12

31 Ago '12

Olá pessoal, Alguém poderia me ajudar no uso do R na plataforma Debian. Tenho instalado na minha máquina, e estava rodando normalmente há alguns dias, porém quando fui trabalhar com hoje não estou conseguindo abrir, pois dar uma mensagem de falha de segmentação. Atenciosamente, Kátia

2 1

Re: [R-br] Mais códigos HTML e função readlines
by Paulo Nogueira 30 Ago '12

30 Ago '12

Opa, só agora reparei que estávamos trocando email fora da lista... Segue para deixar por lá. abs Em 28 de agosto de 2012 17:58, Paulo Nogueira <paulons(a)gmail.com> escreveu: > Olá Augusto. > As páginas simples, que contém o html puro são simples de importar e > trabalhar no R, mas tem algumas outras bem mais complicadas. > Já fiz muita pesquisa por "Web Scrapping R" e suas variantes e ainda > não consegui resolver problemas como esse proposto no tópico. > > Nesse caso, não acho que a página acesse o conteúdo da tabela através > de outros links, mas acredito que esse conteúdo seja "dinâmico" e > gerado por funções javascript desenvolvidas para o site. > Se você utiliza o Chrome ou Firefox, clique na tabela com o botáo > direito do mouse e selecione "Inspecionar elemento". Note que o html > agora está lá, mas porque provavelmente o browser chama essas funções > em tempo real, ou algo assim. > > Preciso fazer o R funcionar assim, mas não sei se é possível. Talvez > tenha de mudar de ferramenta, o que me colocaria na estaca quase zero > novamente. > > Abraços, > Paulo Nogueira Starzynski > > Em 28 de agosto de 2012 15:19, Augusto Ribas <ribas.aca(a)gmail.com> escreveu: >> Rapaz, eu não entendo de web design. Mas eu tava tentando mexer nessa >> pagina que vc queria. >> Eu fiz uns scrapper, mas a paginas que procurava era muito mais simples. >> Essa pagina procura em outros sites, e monta a tabela dela, por isso a >> pagina que vc faz download so tem links. >> >> Mas se vc fizer uma busca por Screen Scrap vai ver algum material interessante. >> >> Olhe esse exemplo: >> http://franklincenterhq.org/2068/simple-r-screen-scraping-example/ >> >> O cara baixa a pagina, dai baixa os links da tabela e ai retira a >> informação que deseja. >> >> Eu não consegui replicar o exemplo dele no site que vc precisa, mas se >> tiver mais sorte que eu poste na lista o seu resultado :) >> Mas acho que vai ser algo bem similar a esse exemplo que vc vai precisar. >> >> Em 28 de agosto de 2012 00:00, Paulo Nogueira <paulons(a)gmail.com> escreveu: >>> Pessoal, não sei se a lista se presta a esclarecer dúvidas como essa >>> minha, mas pelo menos estou usando R para resolver os problemas. rs >>> >>> Ultimamente tenho trabalhado bastante com determinadas buscas na web >>> para coletar dados e estruturá-los. Porém, como meus conceitos de web >>> se restringem ao conhecimento de usuário, e não de desenvolvedor, >>> acabo tendo alguns problemas. >>> Vou exemplificar uma das minhas dificuldades com um problema prático: >>> >>> Utilizar readlines() e os pacores RCurl, XML, etc, para tratar o >>> código fonte é algo que já faço até que bem. Mas problemas surgem >>> quando o conteúdo que eu preciso extrair de informação não vem no >>> código fonte? >>> Vejam esse exemplo: >>> http://www.skyscanner.com.br/passagens-aereas/saoa/slz/121002/130409/tarifa… >>> >>> Após carregar toda a busca, o site SkyScanner exibe os preços >>> encontrados para o trecho São Paulo - São Luis. >>> [1] No navegador Firefox ou Chrome, ao clicar em "Exibir Código >>> Fonte", não aparecem as informações do vôos, apenas o código que >>> estrutura a página. >>> [2] Ao clicar em "Inspecionar Elemento" eu consigo analisar o código >>> fonte do meu interesse. >>> [3] Ao utilizar readlines passando o a url acima o R importa o código >>> fonte [1] e não o [2], com as informações de interesse. >>> >>> Imagino que isso tenha a ver com as funções de javascript e coisas desse tipo... >>> Alguém sabe contornar essa situação com o R? Se não, o que devo >>> estudar para agilizar a busca dessa solução? Estou meio perdido nessa >>> parte e agradeço ajudas de qualquer tipo. rs >>> >>> Peço desculpas se por acaso isso for muito off-topic. >>> Abraços, >>> Paulo Nogueira >>> _______________________________________________ >>> 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. >> >> >> >> -- >> Grato >> Augusto C. A. Ribas >> >> Site Pessoal: http://augustoribas.heliohost.org >> Lattes: http://lattes.cnpq.br/7355685961127056

3 3

Compilar arquivo .rnw
by Andre Oliveira 30 Ago '12

30 Ago '12

Pessoas, boa noite. Estou com um problema que ao tentar compilar meu arquivo .rnw tenho o erro abaixo. R CMD Sweave apostila.rnw . . . . . . Loading required package: snow26 : echo keep.source term verbatim Loading required package: fdth Erro em driver$finish(drobj) : the output file 'apostila_R.tex' has disappeared Calls: <Anonymous> -> <Anonymous> Execução interrompida Alguma sugestão? Já reinstalei snow e fdth e não resolveu. André Oliveira Souza

2 1

ggplot2
by Vinicius Brito Rocha 30 Ago '12

30 Ago '12

Senhores, preciso montar um gráfico sobre o resultado de uma clusterização, gostaria de usar o R+ggplot2. Alguém tem um código semi pronto para plotar as duas componentes principais de forma elegante com ggplot2? Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

[Off-topic] Problema codificação arquivo .tex
by Andre Oliveira 30 Ago '12

30 Ago '12

Boa tarde, meus arquivos .tex não estão compilando acentos corretamente por exemplo André compila AndrÃ@. Como posso resolver isto aqui no meu Ubuntu 12.04. obrigado André Oliveira Souza

4 6

Ajuda em ilustração animada!
by Ivan Bezerra Allaman 28 Ago '12

28 Ago '12

Bom dia senhores! Estou fazendo uma ilustração animada para demonstrar que a média é influenciada por valores extremos com o seguinte código: png(file='exemplos/med0%1d.png',width=500,height=250) x1 <- list() medias <- list() for(i in 1:10){ x <- c(7,10,13.4,15.3,17.6,19.5,20.6,21.1) y <- rep(0,8) x1[[i]] <- c(x,i+21.1) medias[[i]] <- mean(x1[[i]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x1[[i]][9],0,col='red',lwd=7,cex=0.1) points(medias[[i]],-0.2,col='blue',lwd=10,pch=2) } dev.off() png(file='exemplos/med1%1d.png',width=500,height=250) x2 <- list() medias2 <- list() for(j in seq(10,1,-1)){ x2[[j]] <- c(12.3-j,x) medias2[[j]] <- mean(x2[[j]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x1[[j]][9],0,col='red',lwd=7,cex=0.1) points(medias[[j]],-0.2,col='blue',lwd=10,pch=2) } dev.off() png(file='exemplos/med2%1d.png',width=500,height=250) for(j in 1:10){ x2[[j]] <- c(12.3-j,x) medias2[[j]] <- mean(x2[[j]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x2[[j]][1],0,col='red',lwd=7,cex=0.1) points(medias2[[j]],-0.2,col='blue',lwd=10,pch=2) } dev.off() É possível juntar isso em um único código? Desde já grato por qualquer ajuda. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 2

Converter list em data.frame
by Thiago Veloso 28 Ago '12

28 Ago '12

Olá pessoal, Uma dúvida básica: possuo uma lista de objetos que contam a frequência de valores (count) em classes (value). Nessa lista, cada índice ([[1]], [[2]] etc) corresponde a um passo de tempo. Vejam abaixo uma amostra: > head(resd) [[1]] ! essa é a primeira unidade de tempo. value count [1,] 1 267 ! na classe '1' foram encontradas 267 ocorrências, e assim por diante [2,] 6 86 [3,] 7 21 [4,] 8 7 [5,] 9 3683 [6,] 10 537 [7,] 11 135 [8,] 12 17 [9,] 13 74 [10,] NA 2228 [[2]] ! essa é a segunda unidade de tempo. value count [1,] 1 2257 ! na classe '1' foram encontradas 2257 ocorrências, e assim por diante [2,] 2 57 [3,] 3 289 [4,] 4 2 [5,] 6 88 [6,] 7 13 [7,] 8 32 [8,] 9 1137 [9,] 10 756 [10,] 11 99 [11,] 12 22 [12,] 13 67 [13,] 14 1 [14,] 15 7 [15,] NA 2228 [[3]] value count [1,] 1 2170 [2,] 2 61 [3,] 3 316 [4,] 4 3 [5,] 6 74 [6,] 7 1 [7,] 8 12 [8,] 9 781 [9,] 10 1100 [10,] 11 66 [11,] 12 160 [12,] 13 46 [13,] 14 9 [14,] 15 28 [15,] NA 2228 [[4]] value count [1,] 1 2174 [2,] 2 56 [3,] 3 300 [4,] 4 3 [5,] 6 62 [6,] 8 4 [7,] 9 564 [8,] 10 1317 [9,] 11 75 [10,] 12 113 [11,] 13 42 [12,] 14 85 [13,] 15 32 [14,] NA 2228 [[5]] value count [1,] 1 2174 [2,] 2 58 [3,] 3 301 [4,] 4 5 [5,] 6 54 [6,] 8 4 [7,] 9 469 [8,] 10 1393 [9,] 11 83 [10,] 12 76 [11,] 13 41 [12,] 14 136 [13,] 15 33 [14,] NA 2228 [[6]] value count [1,] 1 2185 [2,] 2 52 [3,] 3 300 [4,] 4 4 [5,] 6 52 [6,] 8 3 [7,] 9 397 [8,] 10 1457 [9,] 11 83 [10,] 12 74 [11,] 13 38 [12,] 14 146 [13,] 15 36 [14,] NA 2228 Como converter essa lista em um data.frame, mantendo o índice, e salvar em um arquivo ascii? Agradece pelas dicas, Thiago.

4 5

[Off-topic] - Códigos R - Finanças
by Gabriel Bruno de Lemos 27 Ago '12

27 Ago '12

Boa noite! Para quem trabalha com finanças no R (e/ou Matlab) segue um link que acredito ser útil http://sfb649.wiwi.hu-berlin.de/quantnet/index.php?p=start Este link é do site da Editora Springer. Eles reuniram em um único lugar vários códigos do R e Matlab dos livros que lançaram, tais como: Applied Quantitative Finance, Statistical Tools for Finance and Insurance, Applied Multivariate Statistical Analysis e outros. Abs Gabriel Gabriel Bruno de Lemos Mestrando em Estatística e Experimentação Agronômica Esalq / Usp skype: gb_lemos msn: gb_lemos(a)hotmail.com (19) 8212-6999

1 0

Matriz de proximidade no R
by Gilbert Queiroz 27 Ago '12

27 Ago '12

É possível determinar a Matriz de Proximidade no R? Se sim, como? Abs. Gilbert.

3 2

Re: [R-br] Converter list em data.frame
by Diego Bilski 27 Ago '12

27 Ago '12

você pode usar a função ldply() do pacote plyr > dframe <- ldply(resd, data.frame) mas antes é preciso atribuir nomes aos índices da sua lista. Diego > Message: 2 > Date: Mon, 27 Aug 2012 05:37:12 -0700 (PDT) > From: Thiago Veloso <thi_veloso(a)yahoo.com.br> > To: Lista R <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Converter list em data.frame > Message-ID: > <1346071032.86187.YahooMailNeo(a)web121903.mail.ne1.yahoo.com> > Content-Type: text/plain; charset=iso-8859-1 > > Olá pessoal, > > Uma dúvida básica: possuo uma lista de objetos que contam a frequência de > valores (count) em classes (value). Nessa lista, cada índice ([[1]], [[2]] > etc) corresponde a um passo de tempo. Vejam abaixo uma amostra: > > > head(resd) > [[1]] ! essa é a primeira unidade de tempo. > value count > [1,] 1 267 ! na classe '1' foram encontradas 267 ocorrências, e > assim por diante > [2,] 6 86 > [3,] 7 21 > [4,] 8 7 > [5,] 9 3683 > [6,] 10 537 > [7,] 11 135 > [8,] 12 17 > [9,] 13 74 > [10,] NA 2228 > > [[2]] ! essa é a segunda unidade de tempo. > value count > [1,] 1 2257 ! na classe '1' foram encontradas 2257 ocorrências, e > assim por diante > [2,] 2 57 > [3,] 3 289 > [4,] 4 2 > [5,] 6 88 > [6,] 7 13 > [7,] 8 32 > [8,] 9 1137 > [9,] 10 756 > [10,] 11 99 > [11,] 12 22 > [12,] 13 67 > [13,] 14 1 > [14,] 15 7 > [15,] NA 2228 > > [[3]] > value count > [1,] 1 2170 > [2,] 2 61 > [3,] 3 316 > [4,] 4 3 > [5,] 6 74 > [6,] 7 1 > [7,] 8 12 > [8,] 9 781 > [9,] 10 1100 > [10,] 11 66 > [11,] 12 160 > [12,] 13 46 > [13,] 14 9 > [14,] 15 28 > [15,] NA 2228 > > [[4]] > value count > [1,] 1 2174 > [2,] 2 56 > [3,] 3 300 > [4,] 4 3 > [5,] 6 62 > [6,] 8 4 > [7,] 9 564 > [8,] 10 1317 > [9,] 11 75 > [10,] 12 113 > [11,] 13 42 > [12,] 14 85 > [13,] 15 32 > [14,] NA 2228 > > [[5]] > value count > [1,] 1 2174 > [2,] 2 58 > [3,] 3 301 > [4,] 4 5 > [5,] 6 54 > [6,] 8 4 > [7,] 9 469 > [8,] 10 1393 > [9,] 11 83 > [10,] 12 76 > [11,] 13 41 > [12,] 14 136 > [13,] 15 33 > [14,] NA 2228 > > [[6]] > value count > [1,] 1 2185 > [2,] 2 52 > [3,] 3 300 > [4,] 4 4 > [5,] 6 52 > [6,] 8 3 > [7,] 9 397 > [8,] 10 1457 > [9,] 11 83 > [10,] 12 74 > [11,] 13 38 > [12,] 14 146 > [13,] 15 36 > [14,] NA 2228 > > Como converter essa lista em um data.frame, mantendo o índice, e salvar > em um arquivo ascii? > > Agradece pelas dicas, > Thiago. > >

1 0

dúvida sobre o pacote sendmailR
by Luiz Henrique da Conceicao Leal 27 Ago '12

27 Ago '12

Boa tarde pessoal Estou tendando enviar um e-mail através do R (conforme script abaixo) e está ocorrendo o seguinte problema: library(sendmailR) from = sprintf("<lhleal(a)inmetro.gov.br>") to = "<aprodrigues(a)inmetro.gov.br>" subject = "Hello from R" msg = "It works!" sendmail(from, to, subject, msg,control=list(smtpServer=" webmail.inmetro.gov.br")) *Erro em waitFor(code) : SMTP Error: 5.7.1 Client was not authenticated * Algúem sabe como resolver? Desde já agradeço Luiz

3 3

OFF - Emacs
by Sérgio Henrique almeida da silva ju 25 Ago '12

25 Ago '12

Olá Sou iniciante no emacs e segui o tutorial do Walmes ( http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/) porém antes de modificar o meu .emacs as minhas janelas de Buffers ficavam lado a lado e agora fica em baixo. Como faço para elas ficarem lado a lado, pois prefiro dessa forma? Além disso, como faço para mudar a cor da seleção de um texto, pois com o fundo escuro a cor atual fica quase imperceptível. Obrigado! -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 3

variáveis diferentes resoluções espaciais
by Heloíse Pavanato 25 Ago '12

25 Ago '12

Boa tarde a todos, Tenho algumas variáveis ambientais de diferentes resoluções espaciais (1°,1/60°,1/20°) e gostaria de utilizar a variável de escala mais grosseira para padronizar as outras. Alguém sabe que pacote/função posso usar para fazer isso? Obrigada, Heloise

1 0

Re: [R-br] comparação multipla para dados desbalanceados
by Ivan Bezerra Allaman 24 Ago '12

24 Ago '12

Quanto ao teste de Tukey, se você pegar o teste original, não há nenhuma restrição comparar médias com repetições diferentes. Agora, quanto a pegar a soma de quadrados do erro correta para tal teste, existe uma longa discussão aqui no grupo quanto a isso e eu não irei entrar neste mérito. O pacote TukeyC, agricolae entre outros desempenha perfeitamente a comparação para dados desbalanceados. No entanto, atenção!!!! As funções do pacote TukeyC utilizam internamente a função "aov", ou seja, estará sendo utilizado a soma de quadrados do tipo II por default. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

comparação multipla para dados desbalanceados
by Fernando Antonio de souza 24 Ago '12

24 Ago '12

Caros amigos, Estou necessitando realizar um teste de comparação de médias para dados oriundos de um experimento em delineamento inteiramente casualizado com arranjo fatorial, desbalanceado. Li que para arranjos desbalanceados pode-se utilizar o teste tukey-Kramer. O pacote "agricolae" possui a função HSD.test no entanto não diz se esta função faz o tukey-kramer ou somente tukey (há diferenças entre elas?). Então encontrei o pacote tuckeyC que faz o teste de tukey para dados desbalanceados entretanto não conheço este pacote e gostaria de alguma sugestão sobre o mesmo. Minha outra dúvida é sobre o uso do test SNK. Ele só é permitido para dados balanceados? Pois pacotes como o expeDes e o agricolae só funcionam para ele. Ou isto é uma limitação do pacote? Alguém poderia me sugerir alguma literatura para que eu possa me inteirar sobre o tema? att,

1 0

Join my network on LinkedIn
by Lucas Petri Damiani via LinkedIn 24 Ago '12

24 Ago '12

LinkedIn ------------ Lucas Petri Damiani requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Lucas Petri Damiani http://www.linkedin.com/e/-33g73l-h69bpfvf-2x/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Lucas Petri Damiani http://www.linkedin.com/e/-33g73l-h69bpfvf-2x/rso/200604584/hmbf/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-h69bpfvf-2x/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

[Dúvida] Comentários multiplos.
by Pedro Rafael 23 Ago '12

23 Ago '12

Pessoal vocês sabem se existe alguma forma de fazer comentários múltiplos no R? [ ], Pedro Rafael Diniz Marinho.

6 11

como identificar pontos no mapa
by Natalia Martins 23 Ago '12

23 Ago '12

Boa noite pessoal, estou tentando identificar uns pontos em um mapa. Os comandos que utilizei são os seguintes: geo = data.frame(dados$idestacao,dados$latitude,dados$longitude) head(geo) dados.idestacao dados.latitude dados.longitude 1 1 -24.3209 -54.0111 2 1 -24.3209 -54.0111 3 1 -24.3209 -54.0111 4 1 -24.3209 -54.0111 5 1 -24.3209 -54.0111 6 1 -24.3209 -54.0111 dadosgeo <- as.geodata(geo,coords.col=c(3,2),data.col=1) points(dadosgeo) meso = readShapePoly("E:/natalia/41ME2500G.shp") estado = readShapePoly("E:/natalia/41UF2500G.shp") polygon(estado@polygons[[1]]@Polygons[[1]]@coords) polygon(meso@polygons[[4]]@Polygons[[1]]@coords) identify(dadosgeo) mas quando peço a identificação aparece o seguinte erro: 'x' is a list, but does not have components 'x' and 'y' Não entendi onde estou errando, no entanto quando pedi pra ver meus dadosgeo vi que eles tem apenas as idestacao, as quais deveriam ser identificada no mapa. Muito obrigada. Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

3 6

Mensagem de erro no teste exato de Fisher
by Emerson Cotta Bodevan 23 Ago '12

23 Ago '12

Prezados, considere os dados *>teste<-matrix(c(3,8,15,75,1,12,0,2,18,49,27,56),nrow=4)* ao executar o comando *> fisher.test(teste)* aparece a mensagem abaixo. >Erro em fisher.test(teste) : FEXACT error 6. >LDKEY is too small for this problem. >Try increasing the size of the workspace. Alguém sabe o que significa? Agradeço qualquer ajuda. Abs, Emerson. -- "A educação não muda o mundo. A educação muda as pessoas. As pessoas mudam o mundo." * Paulo Freire.*

4 4

Re: [R-br] [STAT-MATH] NÚMEROS PRIMOS NO R
by Leonard de Assis 22 Ago '12

22 Ago '12

Olá, não sei um algoritmo rápido, mas o caminho se passa por: 1. Testar até o último inteiro menor ou igual à raíz quadrada do número em questão (Acredito ser impossível existir um número divisível por algum número maior que este limite) 2. Ir armazenando os números primos encontrados e testar somente dentro desta lista (Lembrei da fatoração, geralmente utilizamos números primos quando fazemos isto) Acredito que devam existir mais maneiras de otimizar, teria que debruçar sobre o problema com calma []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com Em 22/08/2012 11:35, Francisco Eduardo de Castro escreveu: > > Bom dia, > > Alguém sabe como fazer um programa no R que encontre os 10 mil > primeiros números primos, cuja velocidade de execução da tarefa seja > de até 15 segundos? > > att, > > Francisco > __._,_.___ > <mailto:dumat84@yahoo.com.br?subject=Res%3A%20N%DAMEROS%20PRIMOS%20NO%20R> > | através de email > <mailto:STAT-MATH@yahoogrupos.com.br?subject=Res%3A%20N%DAMEROS%20PRIMOS%20NO%20R> > | Responder através da web > <http://br.groups.yahoo.com/group/STAT-MATH/post;_ylc=X3oDMTJxaTNpbmhvBF9TAz…> > | Adicionar um novo tópico > <http://br.groups.yahoo.com/group/STAT-MATH/post;_ylc=X3oDMTJlMmFqbGU5BF9TAz…> > > Mensagens neste tópico > <http://br.groups.yahoo.com/group/STAT-MATH/message/25756;_ylc=X3oDMTM2YTRuZ…> > (1) > Atividade nos últimos dias: > > * Novos usuários > <http://br.groups.yahoo.com/group/STAT-MATH/members;_ylc=X3oDMTJmbGdkY2poBF9…> > 17 > > Visite seu Grupo > <http://br.groups.yahoo.com/group/STAT-MATH;_ylc=X3oDMTJlZmIwa2JrBF9TAzk3NDk…> > > Quem procura, acha. Encontre a pessoa ideal online > <http://global.ard.yahoo.com/SIG=15m70od66/M=758712.15012700.14784528.129601…> > > ------------------------------------------------------------------------ > Ainda dá tempo de encontrar companhia para passar o Dia dos Namorados! > <http://global.ard.yahoo.com/SIG=15mi8vfse/M=830568.15128937.14797398.129601…> > > ------------------------------------------------------------------------ > Ofertas HP > <http://global.ard.yahoo.com/SIG=15mjg3tka/M=833022.15256589.14900301.127961…> > > > Yahoo! Grupos > <http://br.groups.yahoo.com/;_ylc=X3oDMTJkdGozbmUwBF9TAzk3NDkwNDM1BGdycElkAz…> > > Trocar para: SÃ³ Texto > <mailto:STAT-MATH-traditional@yahoogrupos.com.br?subject=Mudar%20Formato%20de%20Envio:%20Tradicional>, > Resenha DiÃ¡ria > <mailto:STAT-MATH-digest@yahoogrupos.com.br?subject=Envio%20de%20email:%20Resenha> > . Sair do grupo > <mailto:STAT-MATH-unsubscribe@yahoogrupos.com.br?subject=Sair%20do%20grupo> > . Termos de uso <http://br.yahoo.com/info/utos.html> > . > > __,_._,___

3 3

Re: [R-br] Microdados PNAD2009
by Roney Fraga 22 Ago '12

22 Ago '12

Dificuldade em baixar os dados do site do IBGE ou de importar o arquivo .txt para o R? Att Roney

1 0

Microdados PNAD2009
by Maria Cristina Nascimento 22 Ago '12

22 Ago '12

Estou com dificuldades em baixar os microdados da PNAD 2009. Mesmo seguindo o manual disponível no site do IBGE e aplicando as dicas dos usuários, o "R" apresenta erros. Será que alguém pode me ajudar ? Obrigada, Cristina

2 1

[Dúvida] Pacote com rotinas de cálculo numérico.
by Pedro Rafael 22 Ago '12

22 Ago '12

Pessoal quais bons pacotes para realização de cálculos numéricos em R, em especial rotinas para diferenciação e integração unidimensional e multidimensional? Me parece que no package:stats já existe torina de integração numérica, contudo quais outros pacotes vocês conhecem para esta tarefa? O R2CUBA é um bom pacote, ou seja, está bem implementado para tais tarefas? [ ], Pedro Rafael Diniz Marinho.

1 0

ajuda looop
by Heloíse Pavanato 22 Ago '12

22 Ago '12

Luciano e Tito, Obrigada pela ajuda, mas infelizmente não deu certo pois eu tenho duas matrizes diferentes, e a sintaxe que vocês me passaram serviria para vetores. Eu usei um pouco diferente e deu certo: coord <- coord1 <- coord2 <- matrix(NA,480,216) for(i in 1:480){ for(j in 1:216){ if (lat[i,j] < latmax & lat[i,j] > latmin) { if(lon[i,j] < lonmax & lon[i,j] > lonmin) coord[i,j] <- 1 coord1[i,j] <- lat[i,j] coord2[i,j] <- lon[i,j] } } } Agora a minha resposta são matrizes de NA e 1 quando a lat/lon está contida na área do polígono. Obrigada pelas dicas, Abc, Heloise --- Oc. Heloise Pavanato Pós-graduação em Oceanogra Biológica Instituto de Oceanografia Universidade Federal do Rio Grande

1 0

Re: [R-br] Reparar matrizes não positivas definidas
by Ivan Bezerra Allaman 22 Ago '12

22 Ago '12

Talvez em http://statsadventure.blogspot.com.br/2011/11/dealing-with-non-positive-def… \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Ajuda com entrada de argumento
by Clécio da Silva Ferreira 22 Ago '12

22 Ago '12

Caros, estou terminando de montar uma library em modelos assimétricos para o R. Gostaria de que a entrada fosse no mesmo formato de modelos lineares (lm) e glm, ou seja, y~x1+x2+...+xp (com os dados carregados, obviamente). Alguém saberia como fazer isso? Att, xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Prof. Clécio da Silva Ferreira x x Depto Estatística - UFJF x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

4 4

Ajuda com função
by clara luz 22 Ago '12

22 Ago '12

Olá grupo! Estou trabalhando em uma função, mas sou uma iniciante no R, se vcs puderem me dar dicas sobre como torna-lá mais rápida ou simples, seria ótimo!! A função visa alocar diferentes valores ao redor de uma célula focal pertencente a uma matriz, ou seja, a partir de uma célula focal, preencher os vizinhos primários (os 8 vizinhos mais próximos), secundários (os próximos 16 vizinhos ao redor), e assim por diante, com diferentes valores de uma função. aloj<-read.table("2012_aloj.txt") #matriz de '0' e '1', cada '1' será minha célula focal (CF) ep<-expand.grid(1:nrow(aloj),1:ncol(aloj)) #função do R que expande uma matriz, gera nova matriz de 2 colunas, mostrando pra cada célula sua posição x na col 1 e y na col 2 gram<-which(aloj==1,arr.ind = TRUE) #encontra '1's, seleciona as células focais (CFs), suas posições em uma nova matriz 'gram' pbtp<-array(0,dim=c(nrow(aloj),ncol(aloj),nrow(gram))) # array de 3 dimensões onde serão guardadas os valores dos vizinhos de cada célula focal (VZ), cada matriz é referente a uma CF, uma célula com '1', no exemplo são cerca de 250 CFs e 250 matrizes for (k in 1:nrow(gram)) { # analisa cada CF, cada 'k' é referente a uma CF, no caso cerca de 250 'k's x<-gram[k,1] y<-gram[k,2] # x,y é a posição da CF na matriz 'aloj' ps<-numeric(d) #criar um vetor temporário para alocar os valores gerados no 'for p', valores que serão alocados nos vizinhos, um valor pra VZ primário (VZ1), outro pra secundário (VZ2). #d é um argumento variável, pode ser d<-5 for(p in 1:d){ # dentro do for k, para gerar os valores dos VZ ps[p]<-2.718^(-(0.1*p)) #calcula o valor pra cada cél viz, cada 'p' gera um 'ps', esse 'ps' será o valor dos VZ, sendo que ps1=>VZ1, ps2=>VZ2,... #meu problema maior está aqui, a idéia é colocar cada valor ps[p] no seu local adequado (ps1=>VZ1, ps2=>VZ2,...) s<-which(ep[,1]>(x-(p+0.1)) & ep[,1]<(x+(p+0.1)) & ep[,2]>(y-(p+0.1)) & ep[,2]<(y+(p+0.1))) # s é um vetor que seleciona as células ao redor da minha CF, para p=1 encontra os VZ1, porém para p=2 encontra os VZ2 e também os VZ1, ele seleciona todo o intervalo... #cada 'p' gera um 's' diferente, p=1, 's' com 9 elementos (8 VZ1 e 1 CF), p=2, 's' com 25 elementos (16 VZ2, 8 VZ1 e 1 CF).... for (n in 1:length(s)){ #dentro do for p, para colocar os valores ps[d] nos VZ certos v<-ep[s[n],] t<-pbtp[v[1,1],v[1,2],k] if(t==0){ #só analisa as células com '0' pbtp[v[1,1],v[1,2],k]<-ps[p] } # fecha if } #fecha for n } #fecha for p } #fecha for k Até funciona, mas está muito lento, para o exemplo acima, com 250 CF, são cerca de 30s, mas vou usar isso para até 20.000 CF... se alguém souber como tornar a função mais eficiente, seria ótimo, ou mesmo indicar alguma função do R mesmo que já faça isso. Aloj em anexo. Valeu galera!! Clara Luz B. Sant'Anna Bióloga, Bacharel e Licenciada Mestranda em Ecologia - Unicamp Laboratório de Biogeografia da Conservação e Mudanças Climáticas- UFG (19) 9233.1062

2 1

Reparar matrizes não positivas definidas
by Henrique Leme Felizatti 21 Ago '12

21 Ago '12

Caros, estou trabalhando com alguns conjuntos de dados que me retornam matrizes de correlação que são não positivas definidas (por questões numéricas, não teóricas). Devido a isso, não é possível inverter essas matrizes o que inviabiliza a utilização de cópulas para os estudos de simulação. Notei que Rebonato (1999), Grubisic and Pietersz (2003), Kercheval (2009) e and Rapisardo (2006) apresentaram maneiras diferentes de lidar com o assunto. Vocês saberiam me dizer se existe algum pacote ou função no R que trata deste problema ou implementa alguma das metodologias citadas nas referências anteriores ? Grato desde já Henrique

2 1

ajuda loop
by Heloíse Pavanato 21 Ago '12

21 Ago '12

Boa tarde a todos, Estou com um problema em um loop. Tenho duas matrizes de 480 por 216 com latitudes e longitudes e quero selecionar dentre estas somente as que estão dentro de determinado limite (latLim e lonLim), que são pares. Estou usando: lat <- read.table(file="lat.txt") lon <- read.table(file="lon.txt") library(R.matlab) # para ler .mat sx <- readMat('SX.mat',maxLength=NULL, fixNames=TRUE,verbose=F,sparseMatrixClass='SparseM') # limite de longitude lonSX <- sx$SX[3,1,1] lonlim <- as.numeric(lonSX$Lon[1:130]) # limite de latitude latSX <- sx$SX[4,1,1] latlim <- as.numeric(latSX$Lat[1:130]) #### latA <- lonA <- matrix(NA,480,216) for(i in 1:(dim(lat)[1])){ for(j in 1:(dim(lon)[2])){ if (c(lat[i,j],lon[i,j]) > min(latLim,lonLim) # seleciona a lat/lon só para a área do polígono & c(lat[i,j],lon[i,j]) < max(latLim,lonLim)) latA[i,j] <- lat[i,j] lonA[i,j] <- lon[i,j] } } Mas a sintaxe não está correta pois eu não estou sabendo dizer que cada par latLim,lonLim é um delimitador de lat e lon. Se alguém puder me dar uma mãozinha, por favor. Os arquivos estão em anexo. Obrigada, Heloise --- Oc. Heloise Pavanato Pós-graduação em Oceanogra Biológica Instituto de Oceanografia Universidade Federal do Rio Grande

3 2

Re: [R-br] desempenho de um modelo linear a partir de valores previstos e outros dados
by Gilenio Borges Fernandes 20 Ago '12

20 Ago '12

Prezado Pedro Emmanuel A. A. do Brasil Talvez a função cv.lm do pacote DAAG possa lhe ajudar Vejas os códigos abaixo #cv.lm {DAAG} R Documentation #Cross-Validation for Linear Regression #Description #This function gives internal and cross-validation measures of #predictive accuracy for ordinary linear regression. The data #are randomly assigned to a number of ‘folds’. Each fold is removed, #in turn, while the remaining data is used to re-fit the regression #model and to predict at the deleted observations. # x <- c(30,20,60,80,40,50,60,30,70,60) y <- c(73,50,128,170,87,108,135,69,148,132) a=data.frame(x,y);attach(a) mod<-lm(y ~ x) mod summary(mod) anova(mod) model.frame(mod) predict(mod, se.fit = TRUE) library(DAAG) cv.lm(df = a, form.lm = formula(y ~ x), m=3, dots = FALSE, seed=29, plotit=TRUE, printit=TRUE) Boa sorte Gilenio Fernandes Em 9 de agosto de 2012 10:36, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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: Como estimar R2 de modelo linear a partir de previsões em > novos dados? (Ivan Bezerra Allaman) > 2. Re: Erro nls: matriz gradiente singular com estimativas de > parâmetros iniciais (Ivan Bezerra Allaman) > 3. Correlação de pearson (Marcelo Claro de Souza) > 4. Re: Correlação de pearson (FHRB Toledo) > 5. Re: Correlação de pearson (Marcelo Claro de Souza) > 6. bibliog para estudo, preferencialmente em R. (Cleber N.Borges) > 7. Re: Correlação de pearson (Paulo Justiniano) > 8. Re: Erro nls: matriz gradiente singular com estimativas de > parâmetros iniciais (Gustavo Dias Azevedo) > 9. Re: Como estimar R2 de modelo linear a partir de previsões em > novos dados? (Pedro Emmanuel Alvarenga Americano do Brasil) > 10. Re: Como estimar R2 de modelo linear a partir de previsões em > novos dados? (Fernando Colugnati) > 11. Re: Como estimar R2 de modelo linear a partir de previsões em > novos dados? (Pedro Emmanuel Alvarenga Americano do Brasil) > 12. Regressão (kaue veras) > 13. Re: Erro nls: matriz gradiente singular com estimativas de > parâmetros iniciais (Gustavo Dias Azevedo) > 14. Re: Regressão (Walmes Zeviani) > 15. Re: Erro nls: matriz gradiente singular com estimativas de > parâmetros iniciais (Benilton Carvalho) > 16. Re: Regressão (kaue veras) > 17. Re: Como estimar R2 de modelo linear a partir de previsões em > novos dados? (Fernando Colugnati) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 8 Aug 2012 08:13:25 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Como estimar R2 de modelo linear a partir de > previsões em novos dados? > Message-ID: > <1344438805.53229.YahooMailNeo(a)web161803.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia Pedro! > > Tu concordas comigo que a partir do momento que ajustastes um modelo por > meio de uma amostra, o modelo irá fazer uma estimativa da variável resposta > independentemente de qual amostra você utilize (desde que os pontos estejam > dentro do intervalo no qual o modelo foi ajustado é claro) com aquela > precisão no qual foi construído o modelo. Se vc utilizar um modelo ajustado > com a amostra w1 e depois usar o modelo para fazer estimativas com a > amostra w2, estas estimativas foram estimadas com a precisão dada no > primeiro ajuste. Se você realmente, quer avaliar o ajuste do modelo feito > com a amostra w1 e depois com a amostra w2 é simples, basta ajustar um > modelo feito com a amostra w1 e depois ajustar outro modelo feito com a > amostra w2 e comparar os R2, embora não vejo muito sentido nisso, pois é > claro, que os valores serão diferentes sempre, pois é aquela velha > história, se retirarmos 'n' amostras de uma população e retirarmos de cada > amostra a > média, estas médias serão diferentes obviamente pelo simples processo de > amostragem. > > Abraço! > > (S,f,P) > Allaman > > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120808/5d962bb8/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 8 Aug 2012 08:18:33 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Erro nls: matriz gradiente singular com > estimativas de parâmetros iniciais > Message-ID: > <1344439113.80614.YahooMailNeo(a)web161805.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Se você indicou exatamente os valores não deveria dar erro!! Perceba que > primeiro vc nos informa o seguinte modelo: > > y = K*x(^a)*w(^b)/(z^c) > > e depois na função você informa outro modelo: > > modelo <- > nls(y~K*(x^a)*(w^b)*(z^c),data=dados,start=list(K=1,a=0.33,b=0.66,c=-1.33)) > > Provavelmente o erro deve ser este!! > > (S,f,P) > Allaman > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120808/ff73f34f/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 8 Aug 2012 08:25:20 -0700 (PDT) > From: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > To: Rbr <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Correlação de pearson > Message-ID: > <1344439520.35595.YahooMailNeo(a)web160106.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Estou rodando minhas correlações de Pearson segundo modelo abaixo: > > a=c(1,2,3,4,5,6) > b=c(10,13,7,6,8,9) > c=c(12,15,16,22,30,15) > d=c(5,6,7,8,9,10) > x=cbind(a,b,c,d) > x=data.frame(x) > cor(x,method='pearson') > É possível determinar o número de casas decimais que será exibido na > matriz de correlações? > Muito obrigado. > > > Marcelo Claro de Souza > > Biologist, PhD student in Plant Biology > Institute of Bioscience - UNESP, Brazil > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120808/7b4d9d17/atta… > > > > ------------------------------ > > Message: 4 > Date: Wed, 8 Aug 2012 12:35:04 -0300 > From: FHRB Toledo <fernandohtoledo(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Marcelo Claro de Souza > <marcelo_claro(a)yahoo.com.br> > Subject: Re: [R-br] Correlação de pearson > Message-ID: > <CAN55XP7rHnW0QrOKuGa8FfYQfkn-Lri3uunTd= > oPaPj9wHzPvw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Acho que não é exatamente o que você quer, mas... > > round(cor(x, method = 'pearson'), dig = 4) > > 2012/8/8 Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > > > Estou rodando minhas correlações de Pearson segundo modelo abaixo: > > > > a=c(1,2,3,4,5,6) > > b=c(10,13,7,6,8,9) > > c=c(12,15,16,22,30,15) > > d=c(5,6,7,8,9,10) > > x=cbind(a,b,c,d) > > x=data.frame(x) > > cor(x,method='pearson') > > > > É possível determinar o número de casas decimais que será exibido na > > matriz de correlações? > > Muito obrigado. > > > > Marcelo Claro de Souza > > Biologist, PhD student in Plant Biology > > Institute of Bioscience - UNESP, Brazil > > > > > > _______________________________________________ > > 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/20120808/2053a3cb/atta… > > > > ------------------------------ > > Message: 5 > Date: Wed, 8 Aug 2012 08:37:11 -0700 (PDT) > From: Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > To: FHRB Toledo <fernandohtoledo(a)gmail.com> > Cc: Rbr <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Correlação de pearson > Message-ID: > <1344440231.27198.YahooMailNeo(a)web160106.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Olá Fernando, > Acho que não fui bem claro na pergunta, mas era isso mesmo. > Muito obrigado. > > > Marcelo Claro de Souza > > Biologist, PhD student in Plant Biology > Institute of Bioscience - UNESP, Brazil > > > > > ________________________________ > De: FHRB Toledo <fernandohtoledo(a)gmail.com> > Para: r-br(a)listas.c3sl.ufpr.br; Marcelo Claro de Souza < > marcelo_claro(a)yahoo.com.br> > Enviadas: Quarta-feira, 8 de Agosto de 2012 12:35 > Assunto: Re: [R-br] Correlação de pearson > > > Acho que não é exatamente o que você quer, mas... > round(cor(x, method = 'pearson'), dig = 4) > > > 2012/8/8 Marcelo Claro de Souza <marcelo_claro(a)yahoo.com.br> > > Estou rodando minhas correlações de Pearson segundo modelo abaixo: > > > > > >a=c(1,2,3,4,5,6) > >b=c(10,13,7,6,8,9) > >c=c(12,15,16,22,30,15) > >d=c(5,6,7,8,9,10) > >x=cbind(a,b,c,d) > >x=data.frame(x) > >cor(x,method='pearson') > > > >É possível determinar o número de casas decimais que será exibido na > matriz de correlações? > >Muito obrigado. > > > > > >Marcelo Claro de Souza > > > >Biologist, PhD student in Plant Biology > >Institute of Bioscience - UNESP, Brazil > > > > > >_______________________________________________ > >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/20120808/63478da2/atta… > > > > ------------------------------ > > Message: 6 > Date: Wed, 08 Aug 2012 12:46:19 -0300 > From: "Cleber N.Borges" <klebyn(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] bibliog para estudo, preferencialmente em R. > Message-ID: <502289CB.2000600(a)yahoo.com.br> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Prezados, > Boa Tarde, > > Peço por gentileza indicações de Referencias Bibliograficas > para 02 tópicos de estatística, que são: > > 1) Intervalo de confiança de uma prevalência; > 2) Análise não paramétrica para respostas subjetivas/discretas. > > Acredito que o primeiro seja referente à Epidemologia, porém nunca > estudei nada a respeito. > Caso alguém puder confirmar que o termo 'Prevalência' seja usual, > ficaria grato. > Eu sou químico de formação, porém tenho uns 30 livros de estatística > e nunca me deparei com esses 02 tópicos. > > Agradeço antecipadamente por qualquer ajuda. > Obrigado. > > Cleber N.Borges > > > > > ------------------------------ > > Message: 7 > Date: Wed, 8 Aug 2012 13:54:15 -0300 (BRT) > From: Paulo Justiniano <paulojus(a)leg.ufpr.br> > To: Rbr <r-br(a)listas.c3sl.ufpr.br>, Marcelo Claro de Souza > <marcelo_claro(a)yahoo.com.br> > Subject: Re: [R-br] Correlação de pearson > Message-ID: <alpine.DEB.2.00.1208081353380.7502(a)pataxo.est.ufpr.br> > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > options()$digits > format() > prettynum() > > e funcoes relacionadas nas documentacoes destas podem fazer o que deseja > > > > > On Wed, 8 Aug 2012, Marcelo Claro de Souza wrote: > > > Estou rodando minhas correlações de Pearson segundo modelo abaixo: > > > > a=c(1,2,3,4,5,6) > > b=c(10,13,7,6,8,9) > > c=c(12,15,16,22,30,15) > > d=c(5,6,7,8,9,10) > > x=cbind(a,b,c,d) > > x=data.frame(x) > > cor(x,method='pearson') > > > > É possível determinar o número de casas decimais que será exibido na > matriz de correlações? > > Muito obrigado. > > > > Marcelo Claro de Souza > > Biologist, PhD student in Plant Biology > > Institute of Bioscience - UNESP, Brazil > > > > > > > > ------------------------------ > > Message: 8 > Date: Wed, 8 Aug 2012 14:25:42 -0300 > From: Gustavo Dias Azevedo <gustavoazevedo(a)id.uff.br> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Erro nls: matriz gradiente singular com > estimativas de parâmetros iniciais > Message-ID: > <CAG2s6y_7ERWh6pfK5D70_ZP4JRwRXm05SdE4r= > EwMWAYNNiORw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Foi apenas um erro meu de digitação. O modelo correto encontra-se como > escrito na função: (y~K*(x^a)*(w^b)*(z^c). > > Em 8 de agosto de 2012 12:18, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br>escreveu: > > > Se você indicou exatamente os valores não deveria dar erro!! Perceba que > > primeiro vc nos informa o seguinte modelo: > > > > y = K*x(^a)*w(^b)/(z^c) > > > > e depois na função você informa outro modelo: > > > > modelo <- nls(y~K*(x^a)*(w^b)* > > (z^c),data=dados,start=list(K=1,a=0.33,b=0.66,c=-1.33)) > > > > Provavelmente o erro deve ser este!! > > > > (S,f,P) > > Allaman > > * > > * > > \begin{signature} > > <<>>= > > Prof. Dr. Ivan Bezerra Allaman > > Universidade Estadual de Santa Cruz > > Departamento de Ciências Exatas e Tecnológicas > > Ilhéus/BA - Brasil > > Fone: +55 73 3680-5596 > > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > > @ > > \end{signature} > > > > _______________________________________________ > > 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/20120808/249cec6a/atta… > > > > ------------------------------ > > Message: 9 > Date: Wed, 8 Aug 2012 18:52:20 -0300 > From: Pedro Emmanuel Alvarenga Americano do Brasil > <emmanuel.brasil(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Como estimar R2 de modelo linear a partir de > previsões em novos dados? > Message-ID: > <CAFfGvy+cVsO=j= > g6j5p3tVJB5-DB9WW-ajB6SJhqU5ZFtqwT6w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Ivan, > > Entendi o seu comentario mas isso nao me serve. Eu nao quero outro modelo. > O que eu quero e saber se a qualidade das previsoes no futuro sem mantem > como no ajuste inicial. > Caso sejam piores, o modelo inicial necessita de calibracao. > No pacote rms ha as fucoes val.prob e val.surv que o fazem para modelos > logisticos e para os modelos de sobrevivencia,mas mao encontrei para os > mpdelos lineares. > > Se eu vou utilizar esses modelos para prever eventos em pacientes que serao > avaliados no futuro eu gostaria de saber o quanto esse modelo e bom para > esse fim. Por isso me interessa as previsoes no w2 muito mais o que as > previsoes no w1 pelo mesmo modelo. > > Pedro Brasil > via Android (:)= > Em 08/08/2012 12:13, "Ivan Bezerra Allaman" <ivanalaman(a)yahoo.com.br> > escreveu: > > > Bom dia Pedro! > > > > Tu concordas comigo que a partir do momento que ajustastes um modelo por > > meio de uma amostra, o modelo irá fazer uma estimativa da variável > resposta > > independentemente de qual amostra você utilize (desde que os pontos > estejam > > dentro do intervalo no qual o modelo foi ajustado é claro) com aquela > > precisão no qual foi construído o modelo. Se vc utilizar um modelo > ajustado > > com a amostra w1 e depois usar o modelo para fazer estimativas com a > > amostra w2, estas estimativas foram estimadas com a precisão dada no > > primeiro ajuste. Se você realmente, quer avaliar o ajuste do modelo feito > > com a amostra w1 e depois com a amostra w2 é simples, basta ajustar um > > modelo feito com a amostra w1 e depois ajustar outro modelo feito com a > > amostra w2 e comparar os R2, embora não vejo muito sentido nisso, pois é > > claro, que os valores serão diferentes sempre, pois é aquela velha > > história, se retirarmos 'n' amostras de uma população e retirarmos de > cada > > amostra a média, estas médias serão diferentes obviamente pelo simples > > processo de amostragem. > > > > Abraço! > > > > (S,f,P) > > Allaman > > > > * > > * > > \begin{signature} > > <<>>= > > Prof. Dr. Ivan Bezerra Allaman > > Universidade Estadual de Santa Cruz > > Departamento de Ciências Exatas e Tecnológicas > > Ilhéus/BA - Brasil > > Fone: +55 73 3680-5596 > > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > > @ > > \end{signature} > > > > _______________________________________________ > > 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/20120808/53fe561d/atta… > > > > ------------------------------ > > Message: 10 > Date: Wed, 8 Aug 2012 20:37:15 -0300 > From: Fernando Colugnati <fcolugnati(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Cc: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Como estimar R2 de modelo linear a partir de > previsões em novos dados? > Message-ID: > <CAOLzYt0-4qFRTNRz-XPYuQAagvG11mwt_t= > uhDxMAeG5WkRruA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Pedro, me desculpe, mas conceitualmente isso não faz sentido do ponto de > vista frequentista que vc está empregando. No máximo vc conseguirá bandas > de confiança para extrapolações e interpolações a partir dos seus dados, e > verificar se as "previsões com novos dados" (que não são previsões no > sentido que vc está querendo) caem dentro destas bandas. o que isso > significa, não sei ao certo. Na verdade técnicas como Análise Discriminante > e modelos de clasificação utilizam este tipo de abordagem como validação do > modelo, a chamada Crossvalidation, mas mesmo lá, são feitas apenas medidas > de "acerto" nas classificações, dado que se sabe o estado real das > observações desta nova amostra (eg: Doente e Não Doente). > > "quero e saber se a qualidade das previsoes no futuro sem mantem como no > ajuste inicial" . Isso não vai acontecer, principalmente se seus novos > dados estiverem em uma amplitude diferente de observação (algo que o Ivan > já apontou no email dele). > > Este seu raciocínio me parece muito mais algo Bayesiano....aliás, modelos > de regressão para prognóstico de pacientes é algo muito pouco preconizado, > vide literatura (Bland, Altman, Greenland, Rothman, etc...). > > Abs > > > Em 8 de agosto de 2012 18:52, Pedro Emmanuel Alvarenga Americano do Brasil > < > emmanuel.brasil(a)gmail.com> escreveu: > > > Ivan, > > > > Entendi o seu comentario mas isso nao me serve. Eu nao quero outro > modelo. > > O que eu quero e saber se a qualidade das previsoes no futuro sem mantem > > como no ajuste inicial. > > Caso sejam piores, o modelo inicial necessita de calibracao. > > No pacote rms ha as fucoes val.prob e val.surv que o fazem para modelos > > logisticos e para os modelos de sobrevivencia,mas mao encontrei para os > > mpdelos lineares. > > > > Se eu vou utilizar esses modelos para prever eventos em pacientes que > > serao avaliados no futuro eu gostaria de saber o quanto esse modelo e > bom > > para esse fim. Por isso me interessa as previsoes no w2 muito mais o que > as > > previsoes no w1 pelo mesmo modelo. > > > > Pedro Brasil > > via Android (:)= > > Em 08/08/2012 12:13, "Ivan Bezerra Allaman" <ivanalaman(a)yahoo.com.br> > > escreveu: > > > >> Bom dia Pedro! > >> > >> Tu concordas comigo que a partir do momento que ajustastes um modelo por > >> meio de uma amostra, o modelo irá fazer uma estimativa da variável > resposta > >> independentemente de qual amostra você utilize (desde que os pontos > estejam > >> dentro do intervalo no qual o modelo foi ajustado é claro) com aquela > >> precisão no qual foi construído o modelo. Se vc utilizar um modelo > ajustado > >> com a amostra w1 e depois usar o modelo para fazer estimativas com a > >> amostra w2, estas estimativas foram estimadas com a precisão dada no > >> primeiro ajuste. Se você realmente, quer avaliar o ajuste do modelo > feito > >> com a amostra w1 e depois com a amostra w2 é simples, basta ajustar um > >> modelo feito com a amostra w1 e depois ajustar outro modelo feito com a > >> amostra w2 e comparar os R2, embora não vejo muito sentido nisso, pois é > >> claro, que os valores serão diferentes sempre, pois é aquela velha > >> história, se retirarmos 'n' amostras de uma população e retirarmos de > cada > >> amostra a média, estas médias serão diferentes obviamente pelo simples > >> processo de amostragem. > >> > >> Abraço! > >> > >> (S,f,P) > >> Allaman > >> > >> * > >> * > >> \begin{signature} > >> <<>>= > >> Prof. Dr. Ivan Bezerra Allaman > >> Universidade Estadual de Santa Cruz > >> Departamento de Ciências Exatas e Tecnológicas > >> Ilhéus/BA - Brasil > >> Fone: +55 73 3680-5596 > >> E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > >> @ > >> \end{signature} > >> > >> _______________________________________________ > >> 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. > > > > > > -- > Fernando A.B. Colugnati > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120808/363fbe30/atta… > > > > ------------------------------ > > Message: 11 > Date: Thu, 9 Aug 2012 09:27:36 -0300 > From: Pedro Emmanuel Alvarenga Americano do Brasil > <emmanuel.brasil(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Como estimar R2 de modelo linear a partir de > previsões em novos dados? > Message-ID: > < > CAFfGvyLB1A6mJkMEDJzGFpW_zw+RxBZ9NnaWiDAZjNR0qAfaeQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Amigos de R, > > Conceitualmente falando, a unica coisa que eu poderia comentar é essa ideia > não é minha. Ha uma literatura extensa a respeito de modelos para predição. > > Vide > > http://books.google.com.br/books?id=kHGK58cLsMIC&printsec=frontcover&dq=cli… > > ou > > > http://books.google.com.br/books?id=kMyXEJEtFmkC&printsec=frontcover&dq=dia… > > Dentro de modelos de predição clinica ou ferramentas de suporte a decisão, > ha conceitos como desenvolvimento ou ajuste de modelos, calibração ou > penalização de modelos, validação interna, e validação externa. Essa ultima > pode ser validação cruzada, temporal, ou simultanea propriamente dita. > > Imaginem que há um modelo em que estima a probabilidade de um sujeito ser > portador de uma doença de transmissão aérea. Com esse modelo o médico pode > decidir se o paciente ficará em isolamento respiratorio ou não na chegada > do paciente no hospital. Mas percebam que o modelo foi desenvolvido em > outros pacientes. O que interessa nesse momento é se para este paciente em > particular o modelo funcionará razoavelmente. O que representa essa ideia é > a validação externa, e isso vai alem representação da população pela > amostra ou de características da amostra. Geralmente o desempenho dos > modelos em validação externa é pior que o desempenho no ajuste inicial. > Isso é o superajuste, que parece ser muito frequente em modelos para > predição, e é por isso que se recomenda a calibração de modelos para uso em > predição. Essas ideias possuem analogia com o aprendizado, validação e > simulação de redes neurais. > > Mais uma vez a pergunta inicial. Eu encontrei funções que fazem a > estimativa de desempenho de modelo em valores previstos pelos modelos em > outros dados para logisticos e sobrevivencia, mas não encontrei para > modelos lineares. Alguem poderia me indicar alguma forma de estimar o > desempenho de um modelo linear a partir de valores previstos e outros > dados? > > Abraço forte, > > Dr. Pedro Emmanuel A. A. do Brasil > Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 > Instituto de Pesquisa Clínica Evandro Chagas > Fundação Oswaldo Cruz > Rio de Janeiro - Brasil > Av. Brasil 4365, > CEP 21040-360, > Tel 55 21 3865-9648 > email: pedro.brasil(a)ipec.fiocruz.br > email: emmanuel.brasil(a)gmail.com > > ---Apoio aos softwares livres > www.zotero.org - gerenciamento de referências bibliográficas. > www.broffice.org ou www.libreoffice.org - textos, planilhas ou > apresentações. > www.epidata.dk - entrada de dados. > www.r-project.org - análise de dados. > www.ubuntu.com - sistema operacional > > > > Em 8 de agosto de 2012 20:37, Fernando Colugnati <fcolugnati(a)gmail.com > >escreveu: > > > Pedro, me desculpe, mas conceitualmente isso não faz sentido do ponto de > > vista frequentista que vc está empregando. No máximo vc conseguirá bandas > > de confiança para extrapolações e interpolações a partir dos seus dados, > e > > verificar se as "previsões com novos dados" (que não são previsões no > > sentido que vc está querendo) caem dentro destas bandas. o que isso > > significa, não sei ao certo. Na verdade técnicas como Análise > Discriminante > > e modelos de clasificação utilizam este tipo de abordagem como validação > do > > modelo, a chamada Crossvalidation, mas mesmo lá, são feitas apenas > medidas > > de "acerto" nas classificações, dado que se sabe o estado real das > > observações desta nova amostra (eg: Doente e Não Doente). > > > > "quero e saber se a qualidade das previsoes no futuro sem mantem como no > > ajuste inicial" . Isso não vai acontecer, principalmente se seus novos > > dados estiverem em uma amplitude diferente de observação (algo que o Ivan > > já apontou no email dele). > > > > Este seu raciocínio me parece muito mais algo Bayesiano....aliás, modelos > > de regressão para prognóstico de pacientes é algo muito pouco > preconizado, > > vide literatura (Bland, Altman, Greenland, Rothman, etc...). > > > > Abs > > > > > > Em 8 de agosto de 2012 18:52, Pedro Emmanuel Alvarenga Americano do > Brasil > > <emmanuel.brasil(a)gmail.com> escreveu: > > > > Ivan, > >> > >> Entendi o seu comentario mas isso nao me serve. Eu nao quero outro > >> modelo. O que eu quero e saber se a qualidade das previsoes no futuro > sem > >> mantem como no ajuste inicial. > >> Caso sejam piores, o modelo inicial necessita de calibracao. > >> No pacote rms ha as fucoes val.prob e val.surv que o fazem para modelos > >> logisticos e para os modelos de sobrevivencia,mas mao encontrei para os > >> mpdelos lineares. > >> > >> Se eu vou utilizar esses modelos para prever eventos em pacientes que > >> serao avaliados no futuro eu gostaria de saber o quanto esse modelo e > bom > >> para esse fim. Por isso me interessa as previsoes no w2 muito mais o > que as > >> previsoes no w1 pelo mesmo modelo. > >> > >> Pedro Brasil > >> via Android (:)= > >> Em 08/08/2012 12:13, "Ivan Bezerra Allaman" <ivanalaman(a)yahoo.com.br> > >> escreveu: > >> > >>> Bom dia Pedro! > >>> > >>> Tu concordas comigo que a partir do momento que ajustastes um modelo > por > >>> meio de uma amostra, o modelo irá fazer uma estimativa da variável > resposta > >>> independentemente de qual amostra você utilize (desde que os pontos > estejam > >>> dentro do intervalo no qual o modelo foi ajustado é claro) com aquela > >>> precisão no qual foi construído o modelo. Se vc utilizar um modelo > ajustado > >>> com a amostra w1 e depois usar o modelo para fazer estimativas com a > >>> amostra w2, estas estimativas foram estimadas com a precisão dada no > >>> primeiro ajuste. Se você realmente, quer avaliar o ajuste do modelo > feito > >>> com a amostra w1 e depois com a amostra w2 é simples, basta ajustar um > >>> modelo feito com a amostra w1 e depois ajustar outro modelo feito com a > >>> amostra w2 e comparar os R2, embora não vejo muito sentido nisso, pois > é > >>> claro, que os valores serão diferentes sempre, pois é aquela velha > >>> história, se retirarmos 'n' amostras de uma população e retirarmos de > cada > >>> amostra a média, estas médias serão diferentes obviamente pelo simples > >>> processo de amostragem. > >>> > >>> Abraço! > >>> > >>> (S,f,P) > >>> Allaman > >>> > >>> * > >>> * > >>> \begin{signature} > >>> <<>>= > >>> Prof. Dr. Ivan Bezerra Allaman > >>> Universidade Estadual de Santa Cruz > >>> Departamento de Ciências Exatas e Tecnológicas > >>> Ilhéus/BA - Brasil > >>> Fone: +55 73 3680-5596 > >>> E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > >>> @ > >>> \end{signature} > >>> > >>> _______________________________________________ > >>> 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. > >> > > > > > > > > -- > > Fernando A.B. Colugnati > > > > > > > > _______________________________________________ > > 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/20120809/65b7b375/atta… > > > > ------------------------------ > > Message: 12 > Date: Thu, 9 Aug 2012 09:57:01 -0300 > From: kaue veras <kaueveras(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Regressão > Message-ID: > < > CAENvtmjB1QJK+_k3kyUNSNpU1N2JgQnXWzWeF4CzLYMuv2cD9g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia caros colegas, > > Poderiam me dar uma ajuda? Eu gostaria de realizar dois tipos de > regressões, linear e múltipla, podem me orientar como posso realiza-las? > Através de algum pacote ou função? Acredito que não seja muito complicado, > até já realizei uma vez, mas não estou me lembrando como fiz. > > Desde já agradeço. > > Atenciosamente, > > Kauê P. Veras da Cunha > Estatística / 8º Período - UERJ > MSN: kaueveras(a)hotmail.com > E-mail: kaueveras(a)gmail.com > Cel: 8254-9601 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120809/8faeea48/atta… > > > > ------------------------------ > > Message: 13 > Date: Thu, 9 Aug 2012 10:01:29 -0300 > From: Gustavo Dias Azevedo <gustavoazevedo(a)id.uff.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Erro nls: matriz gradiente singular com > estimativas de parâmetros iniciais > Message-ID: > < > CAG2s6y8BLQZYO5GNVM3n84pGd0KmTLCvv0TW+rYoMCQjK1UqRw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Alguém tem ideia de como resolver este problema? Será que o problema está > no modelo? > > Em 8 de agosto de 2012 14:25, Gustavo Dias Azevedo < > gustavoazevedo(a)id.uff.br > > escreveu: > > > Foi apenas um erro meu de digitação. O modelo correto encontra-se como > > escrito na função: (y~K*(x^a)*(w^b)*(z^c). > > > > Em 8 de agosto de 2012 12:18, Ivan Bezerra Allaman < > > ivanalaman(a)yahoo.com.br> escreveu: > > > >> Se você indicou exatamente os valores não deveria dar erro!! Perceba que > >> primeiro vc nos informa o seguinte modelo: > >> > >> y = K*x(^a)*w(^b)/(z^c) > >> > >> e depois na função você informa outro modelo: > >> > >> modelo <- nls(y~K*(x^a)*(w^b)* > >> (z^c),data=dados,start=list(K=1,a=0.33,b=0.66,c=-1.33)) > >> > >> Provavelmente o erro deve ser este!! > >> > >> (S,f,P) > >> Allaman > >> * > >> * > >> \begin{signature} > >> <<>>= > >> Prof. Dr. Ivan Bezerra Allaman > >> Universidade Estadual de Santa Cruz > >> Departamento de Ciências Exatas e Tecnológicas > >> Ilhéus/BA - Brasil > >> Fone: +55 73 3680-5596 > >> E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > >> @ > >> \end{signature} > >> > >> _______________________________________________ > >> 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/20120809/014cc132/atta… > > > > ------------------------------ > > Message: 14 > Date: Thu, 9 Aug 2012 10:03:04 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Regressão > Message-ID: > <CAFU= > EkaDC8FVLvGvgTHBxbD5sffug_by394BNNOAo0pm40C6ww(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Kaue, > > http://zoonek2.free.fr/UNIX/48_R/09.html > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120809/f1143ff7/atta… > > > > ------------------------------ > > Message: 15 > Date: Thu, 9 Aug 2012 14:13:08 +0100 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Erro nls: matriz gradiente singular com > estimativas de parâmetros iniciais > Message-ID: > < > CAO-arWMi_+55rRzm2gJLXdTbmPY9umDYoSayoyZ84QMvFRxv_Q(a)mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > > CMR.... > > > ------------------------------ > > Message: 16 > Date: Thu, 9 Aug 2012 10:21:02 -0300 > From: kaue veras <kaueveras(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Regressão > Message-ID: > < > CAENvtmg43HxT8j-HVZJmMF5b242-iT87zJEvCftoQmYk9WT7Mg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Muito obrigado Walmes, muito útil. > > Att, > > Kauê P. Veras da Cunha > Estatística / 8º Período - UERJ > MSN: kaueveras(a)hotmail.com > E-mail: kaueveras(a)gmail.com > Cel: 8254-9601 > > 2012/8/9 Walmes Zeviani <walmeszeviani(a)gmail.com> > > > Kaue, > > > > http://zoonek2.free.fr/UNIX/48_R/09.html > > > > À disposição. > > Walmes. > > > > > ========================================================================== > > Walmes Marques Zeviani > > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 > W) > > Departamento de Estatística - Universidade Federal do Paraná > > fone: (+55) 41 3361 3573 > > VoIP: (3361 3600) 1053 1173 > > e-mail: walmes(a)ufpr.br > > twitter: @walmeszeviani > > homepage: http://www.leg.ufpr.br/~walmes > > linux user number: 531218 > > > ========================================================================== > > > > _______________________________________________ > > 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/20120809/c3c6d5cb/atta… > > > > ------------------------------ > > Message: 17 > Date: Thu, 9 Aug 2012 10:36:50 -0300 > From: Fernando Colugnati <fcolugnati(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Como estimar R2 de modelo linear a partir de > previsões em novos dados? > Message-ID: > <CAOLzYt2qZhUq13hGsd5qXpGz8zwJ7tPF4cckP78= > QJo2NUJ4NQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bem, a sugestão que me parece mais pertinente é ver o que é utilizado para > estes modelos logísticos e de sobrevivência, e adapte para os lineares, > pensando que os outros dois modelos também são lineares a partir da > definição apropriada da família de distribuição e função link que vc > utiliza. A extensão para mim pareceria natural, dentro do framework de > GLMs. > > O R deve fornecer algo parecido em pacotes para CART, Redes Neurais, etc... > > Outra sugestão, veja como algumas fórmulas foram propostas para por exemplo > Filtração Glomerular para classificação do estágio doenças renais crônicas > (a partir da Creatinina), previsão de % de gordura por soma de pregas, > Framinghan, etc...são todos modelos deste tipo que vc parece buscar, e que > profissionais de saúde usam no dia e dia, e por mais críticas que possamos > ter como estatísticos a estes modelos, eles funcionam nos serviços. > > Mas certamente, a ideia do R2 sugerida, para mim não faz sentido. > > > Abs > > Em 9 de agosto de 2012 09:27, Pedro Emmanuel Alvarenga Americano do Brasil > < > emmanuel.brasil(a)gmail.com> escreveu: > > > Amigos de R, > > > > Conceitualmente falando, a unica coisa que eu poderia comentar é essa > > ideia não é minha. Ha uma literatura extensa a respeito de modelos para > > predição. > > > > Vide > > > http://books.google.com.br/books?id=kHGK58cLsMIC&printsec=frontcover&dq=cli… > > > > ou > > > > > > > http://books.google.com.br/books?id=kMyXEJEtFmkC&printsec=frontcover&dq=dia… > > > > Dentro de modelos de predição clinica ou ferramentas de suporte a > decisão, > > ha conceitos como desenvolvimento ou ajuste de modelos, calibração ou > > penalização de modelos, validação interna, e validação externa. Essa > ultima > > pode ser validação cruzada, temporal, ou simultanea propriamente dita. > > > > Imaginem que há um modelo em que estima a probabilidade de um sujeito ser > > portador de uma doença de transmissão aérea. Com esse modelo o médico > pode > > decidir se o paciente ficará em isolamento respiratorio ou não na chegada > > do paciente no hospital. Mas percebam que o modelo foi desenvolvido em > > outros pacientes. O que interessa nesse momento é se para este paciente > em > > particular o modelo funcionará razoavelmente. O que representa essa > ideia é > > a validação externa, e isso vai alem representação da população pela > > amostra ou de características da amostra. Geralmente o desempenho dos > > modelos em validação externa é pior que o desempenho no ajuste inicial. > > Isso é o superajuste, que parece ser muito frequente em modelos para > > predição, e é por isso que se recomenda a calibração de modelos para uso > em > > predição. Essas ideias possuem analogia com o aprendizado, validação e > > simulação de redes neurais. > > > > Mais uma vez a pergunta inicial. Eu encontrei funções que fazem a > > estimativa de desempenho de modelo em valores previstos pelos modelos em > > outros dados para logisticos e sobrevivencia, mas não encontrei para > > modelos lineares. Alguem poderia me indicar alguma forma de estimar o > > desempenho de um modelo linear a partir de valores previstos e outros > dados? > > > > Abraço forte, > > > > Dr. Pedro Emmanuel A. A. do Brasil > > Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 > > Instituto de Pesquisa Clínica Evandro Chagas > > Fundação Oswaldo Cruz > > Rio de Janeiro - Brasil > > Av. Brasil 4365, > > CEP 21040-360, > > Tel 55 21 3865-9648 > > email: pedro.brasil(a)ipec.fiocruz.br > > email: emmanuel.brasil(a)gmail.com > > > > ---Apoio aos softwares livres > > www.zotero.org - gerenciamento de referências bibliográficas. > > www.broffice.org ou www.libreoffice.org - textos, planilhas ou > > apresentações. > > www.epidata.dk - entrada de dados. > > www.r-project.org - análise de dados. > > www.ubuntu.com - sistema operacional > > > > > > > > Em 8 de agosto de 2012 20:37, Fernando Colugnati <fcolugnati(a)gmail.com > >escreveu: > > > > Pedro, me desculpe, mas conceitualmente isso não faz sentido do ponto de > >> vista frequentista que vc está empregando. No máximo vc conseguirá > bandas > >> de confiança para extrapolações e interpolações a partir dos seus > dados, e > >> verificar se as "previsões com novos dados" (que não são previsões no > >> sentido que vc está querendo) caem dentro destas bandas. o que isso > >> significa, não sei ao certo. Na verdade técnicas como Análise > Discriminante > >> e modelos de clasificação utilizam este tipo de abordagem como > validação do > >> modelo, a chamada Crossvalidation, mas mesmo lá, são feitas apenas > medidas > >> de "acerto" nas classificações, dado que se sabe o estado real das > >> observações desta nova amostra (eg: Doente e Não Doente). > >> > >> "quero e saber se a qualidade das previsoes no futuro sem mantem como no > >> ajuste inicial" . Isso não vai acontecer, principalmente se seus novos > >> dados estiverem em uma amplitude diferente de observação (algo que o > Ivan > >> já apontou no email dele). > >> > >> Este seu raciocínio me parece muito mais algo Bayesiano....aliás, > modelos > >> de regressão para prognóstico de pacientes é algo muito pouco > preconizado, > >> vide literatura (Bland, Altman, Greenland, Rothman, etc...). > >> > >> Abs > >> > >> > >> Em 8 de agosto de 2012 18:52, Pedro Emmanuel Alvarenga Americano do > >> Brasil <emmanuel.brasil(a)gmail.com> escreveu: > >> > >> Ivan, > >>> > >>> Entendi o seu comentario mas isso nao me serve. Eu nao quero outro > >>> modelo. O que eu quero e saber se a qualidade das previsoes no futuro > sem > >>> mantem como no ajuste inicial. > >>> Caso sejam piores, o modelo inicial necessita de calibracao. > >>> No pacote rms ha as fucoes val.prob e val.surv que o fazem para modelos > >>> logisticos e para os modelos de sobrevivencia,mas mao encontrei para os > >>> mpdelos lineares. > >>> > >>> Se eu vou utilizar esses modelos para prever eventos em pacientes que > >>> serao avaliados no futuro eu gostaria de saber o quanto esse modelo e > bom > >>> para esse fim. Por isso me interessa as previsoes no w2 muito mais o > que as > >>> previsoes no w1 pelo mesmo modelo. > >>> > >>> Pedro Brasil > >>> via Android (:)= > >>> Em 08/08/2012 12:13, "Ivan Bezerra Allaman" <ivanalaman(a)yahoo.com.br> > >>> escreveu: > >>> > >>>> Bom dia Pedro! > >>>> > >>>> Tu concordas comigo que a partir do momento que ajustastes um modelo > >>>> por meio de uma amostra, o modelo irá fazer uma estimativa da variável > >>>> resposta independentemente de qual amostra você utilize (desde que os > >>>> pontos estejam dentro do intervalo no qual o modelo foi ajustado é > claro) > >>>> com aquela precisão no qual foi construído o modelo. Se vc utilizar um > >>>> modelo ajustado com a amostra w1 e depois usar o modelo para fazer > >>>> estimativas com a amostra w2, estas estimativas foram estimadas com a > >>>> precisão dada no primeiro ajuste. Se você realmente, quer avaliar o > ajuste > >>>> do modelo feito com a amostra w1 e depois com a amostra w2 é simples, > basta > >>>> ajustar um modelo feito com a amostra w1 e depois ajustar outro modelo > >>>> feito com a amostra w2 e comparar os R2, embora não vejo muito sentido > >>>> nisso, pois é claro, que os valores serão diferentes sempre, pois é > aquela > >>>> velha história, se retirarmos 'n' amostras de uma população e > retirarmos de > >>>> cada amostra a média, estas médias serão diferentes obviamente pelo > simples > >>>> processo de amostragem. > >>>> > >>>> Abraço! > >>>> > >>>> (S,f,P) > >>>> Allaman > >>>> > >>>> * > >>>> * > >>>> \begin{signature} > >>>> <<>>= > >>>> Prof. Dr. Ivan Bezerra Allaman > >>>> Universidade Estadual de Santa Cruz > >>>> Departamento de Ciências Exatas e Tecnológicas > >>>> Ilhéus/BA - Brasil > >>>> Fone: +55 73 3680-5596 > >>>> E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > >>>> @ > >>>> \end{signature} > >>>> > >>>> _______________________________________________ > >>>> 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. > >>> > >> > >> > >> > >> -- > >> Fernando A.B. Colugnati > >> > >> > >> > >> _______________________________________________ > >> 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. > > > > > > -- > Fernando A.B. Colugnati > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120809/8b10f55f/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 18, assunto 8 > **************************************** > -- Gilenio Borges Fernandes Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

2 1

Empilhar data-frames com colunas diferentes
by Paulo Nogueira 20 Ago '12

20 Ago '12

Boa tarde grupo! Estou com um problema que acredito ser bastante simples de resolver, porém, pesquisei e ainda não encontrei alternativa prática. A grosso modo, quero empilhar bases para formar uma única. O problema é que elas não possuem o mesmo layout. Todas possuem algumas colunas em comum e outras não. Quando uma base não possuir determinada coluna gostaria de deixar NA nos registros oriundos dessa base. Para quem conhece SAS, procuro o resultado de fazer simplesmente: DATA base_final; SET Base1 Base2; RUN; No R ainda não consigo utilizar a função rbind com o mesmo resultado. Segue exemplo para reprodução. # Exemplo d1<-data.frame(index=1,a=1,b=91)d2<-data.frame(index=1,b=92,a=2,c=200)rbind.data.frame(d1,d2)Erro em rbind.data.frame(d1, d2) : números de colunas dos argumentos não correspondem O resultado que desejo é index a b c 1 1 91 NA 1 2 92 200 O rbind troca a posição das colunas baseada nos nomes. Isso ajuda. Mas não sei como fazer funcionar quando algumas colunas estão apenas em algumas bases, o que no exemplo acontece com a coluna c da base d2. Como vocês fariam? Abraços, Paulo Nogueira Starzynski

3 4

Re: [R-br] Regressão Gama
by Ivan Bezerra Allaman 20 Ago '12

20 Ago '12

Neste site existem "n" gráficos com intervalos de confiança. Provavelmente seja possível adaptar para o seu caso. http://ridiculas.wordpress.com/ \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 2

Re: [R-br] TABELAS POR GRUPO
by Ivan Bezerra Allaman 20 Ago '12

20 Ago '12

?fdth ?descr \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

TABELAS POR GRUPO
by Edson Lira 20 Ago '12

20 Ago '12

Bom dia,estou usando esta rotina para analisar o banco postado abaixo. (Veja no link) g<-read.csv2("C:\\Consulta\\HEMOAM\\Elene_fisio\\gest.csv") g[g==""]<-NA names(g) tabelas<-function(x){ qtde<-table(x) perc<-100*prop.table(table(x)) cbind(qtde,perc) } options(digits=3) #imc g$IMC_f<-rep(NA,length(g$IMC)) summary(g$IMC) g$IMC_f[g$IMC<18.5]<-'Baixo peso' g$IMC_f[18.5<=g$IMC &g$IMC<25]<-'Normal' g$IMC_f[25<=g$IMC &g$IMC<30]<-'Pré-obesidade' g$IMC_f[30<=g$IMC &g$IMC<35]<-'Obesidade - I' g$IMC_f[35<=g$IMC &g$IMC<40]<-'Obesidade - II' g$IMC_f[g$IMC>=40]<-'Obesidade - III' #Idade g$ida<-rep(NA,length(g$Idade)) summary(g$Idade) g$ida[18<=g$Idade &g$Idade<25]<-'1. 18 a 25' g$ida[25<=g$Idade &g$Idade<32]<-'2. 25 a 32' g$ida[32<=g$Idade &g$Idade<39]<-'3. 32 a 39' apply(g[,c(-1,-2,-3,-7,-11,-14,-15,-16)],2,tabelas) http://www.datafilehost.com/download-64b7d14d.html Gostaria de obter um resulado semelhante a rotina tabelas, só que estratificado por grupo(Controle, Tratado). Alguém tem uma sugestão? Edson Lira Estatístico Manaus-Amazonas

1 0

Re: [R-br] Somar vendas por semana e por mês
by Hugo Sena 18 Ago '12

18 Ago '12

Boas. Obrigado pela dica, penso que a função aggregate() resolve mas não consigo utilizar. Tenho os seguintes dados: *Data* *Qtd* 28-10-2003 1 28-10-2003 1 31-10-2003 1 05-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 11-11-2003 1 E o que pretendo é somar as quantidades em cada dia de forma a ficar assim: Data Qtd 28-10-2003 2 31-10-2003 1 05-11-2003 1 11-11-2003 6 Podem dar uma ajuda? E depois é suposto exportar esta agregação para um ficheiro? Como? Já procurei no google mas ainda não consegui fazer no meu caso em particular. Obrigado. Cumps.

1 0

Re: [R-br] Anova em faixas
by Juliano Farias 18 Ago '12

18 Ago '12

Srs, Grato pela ajuda. Éder o link que você repassou tem exemplos com dois fatores. O que eu tenho é com três fatores. Em Inglês o delineamento é chamado de strip split plot. Segue o quadro da análise da variância para o meu tipo de experimento. Não consegui montar esse quadro da ANOVA no R. CV GL QM F(SobHo) B b-1 V1 V1/V5 A a-1 V2 V2/V5 C c-1 V3 V3/V5 A*C (a-1)(c-1) V4 V4/V5 Erro(ac) (ac-1)(b-1) V5 V5/V11 D d-1 V6 V6/V7 Erro(d) (b-1)(d-1) V7 V7/V11 A*D (a-1)(d-1) V8 V8/V11 C*D (c-1)(d-1) V9 V9/V11 A*C*D (a-1)(c-1)(d-1) V10 V10/V11 Erro (ac-1)(d-1)(b-1) V11 ----- Atenciosamente, Juliano Em 16 de agosto de 2012 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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: Regressão Gama (Ivan Bezerra Allaman) > 2. Re: Regressão Gama (Walmes Zeviani) > 3. Anova em faixas (Juliano Farias) > 4. Re: Anova em faixas (Eder David Borges da Silva) > 5. Re: Anova em faixas (Listeiro 037) > 6. Re: Anova em faixas (Eder David Borges da Silva) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Aug 2012 08:20:07 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <1345044007.77714.YahooMailNeo(a)web161804.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Neste site existem "n" gráficos com intervalos de confiança. Provavelmente > seja possível adaptar para o seu caso. > > http://ridiculas.wordpress.com/ > > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/b1f6625c/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 15 Aug 2012 13:31:28 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <CAFU=EkYGjAuNnRnBD7fV_5+n4Yd0QS2NM3E+pZ6CyKX= > Ez+U1A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Além da sugestão do Ivan, como trata-se de um GLM, recomendo a leitura > dessa matéria > > > http://freakonometrics.blog.free.fr/index.php?post/2011/11/01/Fisher,-delta… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/9188fb40/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 15 Aug 2012 21:32:40 -0300 > From: Juliano Farias <julianofarias(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Anova em faixas > Message-ID: > < > CAEer2Pik5nco26+F88tkCr02nAx0-kH3geR7GfF01p33+Zwz7g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Srs, > > Preciso de uma ajuda. > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > problema está no delineamento que foi em faixas. > Tenho o quadro da análise da variança para esse tipo de delineamento, porém > não consegui montar a programação e não achei um pacote para esse tipo de > delineamento. > O experimento foi realizado em blocos (B) ao acaso com quatro repetições em > esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com os > dois níveis do fator C foram casualizados nas parcelas principais, e os > dois níveis do fator D foram casualizados em faixas perpendiculares às > parcelas principais. > > Grato pela atenção. > > > -- > Atenciosamente, > > Juliano Ricardo Farias > Eng. Agrônomo > Doutorando em Entomologia > Lab. de Resistência de Artrópodes a Pesticidas > ESALQ/USP > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/1097b154/atta… > > > > ------------------------------ > > Message: 4 > Date: Thu, 16 Aug 2012 07:35:30 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > <CALKkXXru-3qZThcWiKNAeHAN= > TMYnzF2uyDMLE4Mj1RAyZk47w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Juliano, > Acho que o script neste link pode te ajudar, > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > Att > Éder > > Em 15 de agosto de 2012 21:32, Juliano Farias <julianofarias(a)gmail.com > >escreveu: > > > Srs, > > > > Preciso de uma ajuda. > > > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > > problema está no delineamento que foi em faixas. > > Tenho o quadro da análise da variança para esse tipo de delineamento, > > porém não consegui montar a programação e não achei um pacote para esse > > tipo de delineamento. > > O experimento foi realizado em blocos (B) ao acaso com quatro repetições > > em esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com > > os dois níveis do fator C foram casualizados nas parcelas principais, e > > os dois níveis do fator D foram casualizados em faixas perpendiculares às > > parcelas principais. > > > > Grato pela atenção. > > > > > > -- > > Atenciosamente, > > > > Juliano Ricardo Farias > > Eng. Agrônomo > > Doutorando em Entomologia > > Lab. de Resistência de Artrópodes a Pesticidas > > ESALQ/USP > > > > _______________________________________________ > > 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/20120816/3100052d/atta… > > > > ------------------------------ > > Message: 5 > Date: Thu, 16 Aug 2012 09:27:05 -0300 > From: Listeiro 037 <listeiro_037(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: <375172.55957.bm(a)smtp146.mail.mud.yahoo.com> > Content-Type: text/plain; charset=UTF-8 > > Olá. > > Há como colocar shebang nesse script de R abaixo? > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > Juliano, > > Acho que o script neste link pode te ajudar, > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > Att > > Éder > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > <julianofarias(a)gmail.com>escreveu: > > > > -- > "There are three kinds of lies: lies, damn lies, and statistics." > Benjamin Disraeli > > > ------------------------------ > > Message: 6 > Date: Thu, 16 Aug 2012 10:44:36 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > < > CALKkXXqS7WV+sdjDoWbdMyZgHZsOf_qapeZbeXL+1qMmEnpirQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > ?Shebang > > Em 16 de agosto de 2012 09:27, Listeiro 037 <listeiro_037(a)yahoo.com.br > >escreveu: > > > Olá. > > > > Há como colocar shebang nesse script de R abaixo? > > > > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > > > Juliano, > > > Acho que o script neste link pode te ajudar, > > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > > Att > > > Éder > > > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > > <julianofarias(a)gmail.com>escreveu: > > > > > > > -- > > "There are three kinds of lies: lies, damn lies, and statistics." > > Benjamin Disraeli > > _______________________________________________ > > 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/20120816/c74a1216/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 18, assunto 21 > ***************************************** > -- Atenciosamente, Juliano Ricardo Farias Eng. Agrônomo Doutorando em Entomologia Lab. de Resistência de Artrópodes a Pesticidas ESALQ/USP

2 1

Re: [R-br] Digest R-br, volume 18, assunto 21
by Juliano Farias 17 Ago '12

17 Ago '12

Listeiro 037, Pelo que entendi do shebang, sim pode colocar no script do R. Obrigado, Juliano. Em 16 de agosto de 2012 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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: Regressão Gama (Ivan Bezerra Allaman) > 2. Re: Regressão Gama (Walmes Zeviani) > 3. Anova em faixas (Juliano Farias) > 4. Re: Anova em faixas (Eder David Borges da Silva) > 5. Re: Anova em faixas (Listeiro 037) > 6. Re: Anova em faixas (Eder David Borges da Silva) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Aug 2012 08:20:07 -0700 (PDT) > From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> > To: R Brasil <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <1345044007.77714.YahooMailNeo(a)web161804.mail.bf1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Neste site existem "n" gráficos com intervalos de confiança. Provavelmente > seja possível adaptar para o seu caso. > > http://ridiculas.wordpress.com/ > > > > \begin{signature} > <<>>= > Prof. Dr. Ivan Bezerra Allaman > Universidade Estadual de Santa Cruz > Departamento de Ciências Exatas e Tecnológicas > Ilhéus/BA - Brasil > Fone: +55 73 3680-5596 > E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com > @ > \end{signature} > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/b1f6625c/atta… > > > > ------------------------------ > > Message: 2 > Date: Wed, 15 Aug 2012 13:31:28 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman > <ivanalaman(a)yahoo.com.br> > Subject: Re: [R-br] Regressão Gama > Message-ID: > <CAFU=EkYGjAuNnRnBD7fV_5+n4Yd0QS2NM3E+pZ6CyKX= > Ez+U1A(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Além da sugestão do Ivan, como trata-se de um GLM, recomendo a leitura > dessa matéria > > > http://freakonometrics.blog.free.fr/index.php?post/2011/11/01/Fisher,-delta… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/9188fb40/atta… > > > > ------------------------------ > > Message: 3 > Date: Wed, 15 Aug 2012 21:32:40 -0300 > From: Juliano Farias <julianofarias(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Anova em faixas > Message-ID: > < > CAEer2Pik5nco26+F88tkCr02nAx0-kH3geR7GfF01p33+Zwz7g(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Srs, > > Preciso de uma ajuda. > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > problema está no delineamento que foi em faixas. > Tenho o quadro da análise da variança para esse tipo de delineamento, porém > não consegui montar a programação e não achei um pacote para esse tipo de > delineamento. > O experimento foi realizado em blocos (B) ao acaso com quatro repetições em > esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com os > dois níveis do fator C foram casualizados nas parcelas principais, e os > dois níveis do fator D foram casualizados em faixas perpendiculares às > parcelas principais. > > Grato pela atenção. > > > -- > Atenciosamente, > > Juliano Ricardo Farias > Eng. Agrônomo > Doutorando em Entomologia > Lab. de Resistência de Artrópodes a Pesticidas > ESALQ/USP > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120815/1097b154/atta… > > > > ------------------------------ > > Message: 4 > Date: Thu, 16 Aug 2012 07:35:30 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > <CALKkXXru-3qZThcWiKNAeHAN= > TMYnzF2uyDMLE4Mj1RAyZk47w(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Juliano, > Acho que o script neste link pode te ajudar, > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > Att > Éder > > Em 15 de agosto de 2012 21:32, Juliano Farias <julianofarias(a)gmail.com > >escreveu: > > > Srs, > > > > Preciso de uma ajuda. > > > > Preciso fazer uma ANOVA de um experimento com três fatores, porém o > > problema está no delineamento que foi em faixas. > > Tenho o quadro da análise da variança para esse tipo de delineamento, > > porém não consegui montar a programação e não achei um pacote para esse > > tipo de delineamento. > > O experimento foi realizado em blocos (B) ao acaso com quatro repetições > > em esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com > > os dois níveis do fator C foram casualizados nas parcelas principais, e > > os dois níveis do fator D foram casualizados em faixas perpendiculares às > > parcelas principais. > > > > Grato pela atenção. > > > > > > -- > > Atenciosamente, > > > > Juliano Ricardo Farias > > Eng. Agrônomo > > Doutorando em Entomologia > > Lab. de Resistência de Artrópodes a Pesticidas > > ESALQ/USP > > > > _______________________________________________ > > 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/20120816/3100052d/atta… > > > > ------------------------------ > > Message: 5 > Date: Thu, 16 Aug 2012 09:27:05 -0300 > From: Listeiro 037 <listeiro_037(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: <375172.55957.bm(a)smtp146.mail.mud.yahoo.com> > Content-Type: text/plain; charset=UTF-8 > > Olá. > > Há como colocar shebang nesse script de R abaixo? > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > Juliano, > > Acho que o script neste link pode te ajudar, > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > Att > > Éder > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > <julianofarias(a)gmail.com>escreveu: > > > > -- > "There are three kinds of lies: lies, damn lies, and statistics." > Benjamin Disraeli > > > ------------------------------ > > Message: 6 > Date: Thu, 16 Aug 2012 10:44:36 -0300 > From: Eder David Borges da Silva <eder(a)leg.ufpr.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Anova em faixas > Message-ID: > < > CALKkXXqS7WV+sdjDoWbdMyZgHZsOf_qapeZbeXL+1qMmEnpirQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > ?Shebang > > Em 16 de agosto de 2012 09:27, Listeiro 037 <listeiro_037(a)yahoo.com.br > >escreveu: > > > Olá. > > > > Há como colocar shebang nesse script de R abaixo? > > > > > > Em Thu, 16 Aug 2012 07:35:30 -0300 > > Eder David Borges da Silva <eder(a)leg.ufpr.br> escreveu: > > > > > Juliano, > > > Acho que o script neste link pode te ajudar, > > > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > > > Att > > > Éder > > > > > > Em 15 de agosto de 2012 21:32, Juliano Farias > > > <julianofarias(a)gmail.com>escreveu: > > > > > > > -- > > "There are three kinds of lies: lies, damn lies, and statistics." > > Benjamin Disraeli > > _______________________________________________ > > 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/20120816/c74a1216/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 18, assunto 21 > ***************************************** > -- Atenciosamente, Juliano Ricardo Farias Eng. Agrônomo Doutorando em Entomologia Lab. de Resistência de Artrópodes a Pesticidas ESALQ/USP

2 1

Re: [R-br] dúvida - modelo [OFF-TOPIC]
by Fernando Colugnati 17 Ago '12

17 Ago '12

Só um pitaco na modelagem do ponto de vista fisiológico...o fato do cara não fazer atividade física não implica ele em ter Gasto Energético 0 !!!!!! Se vc utilizar esta abordagem vai ter uma subestimação não real das médias ajustadas e uma superestimação de efeitos. Neste caso vcs deveriam ter no mínimo o gasto energético basal! Se vcs estiverem medindo o gasto energético pelo IPAQ (um questionário, que na verdade fornece METs, uma medida proxy de gasto energético) sugiro ver outros trabalhos que adequam este problema. Estes 0s devem ser pessoas que não responderam ao questionário por não fazerem atividades, mas isso não quer dizer que não gastam energia....neste caso, estes caras não deveriam entrar na modelagem. Outro ponto, mesmo medido como valores inteiros, também não faz sentido os dados seguirem uma distribuição Geométrica...na verdade vcs devem ter uma discretização deste dado, mas veja a definição da Geométrica e verá que não faz sentido. Nâo assumo uma distribuição pq os dados parecem com aquela distribuição, mas sim pelo processo que pode ter gerado aqueles dados, mas isso é pessoal...há escolas que fazem desta forma. Desculpem se dei pitacos desnecessários, mas pela breve descrição do problema isso me chamou a atenção! Abs Em 16 de agosto de 2012 22:47, Sérgio Henrique <sergio.edfisica(a)gmail.com>escreveu: > Oi > > Obrigado por sempre ajudar Benilton. > > Deixa eu tentar explicar o objetivo do trabalho. > > Na verdade que vai trabalhar com esses dados é um amigo no doutorado. > > A variável resposta é Atividade Física e ele quer modelar, na verdade, > seria gasto energético por AF, os 0 são indivíduos que não praticam AF, ou > seja, tem gasto 0, acima desse valor, são gastos energéticos de atividades > físicas quaisquer... > Os dados não são normais e quando tento normalizar com boxcox não roda por > conta dos 0. > Quando fiz um histograma dessa variável, me pareceu que os mesmos > apresentam distribuição geométrica, parei por aí... > > Obrigado mais uma vez. > > Enviado via iPod > > Em 16/08/2012, às 22:11, Benilton Carvalho <beniltoncarvalho(a)gmail.com> > escreveu: > > > modelagem envolve muito mais que os dados com os quais vc esta' > > trabalhando... na verdade, o processo gerador dos dados e' q e' > > relevante... vc precisa entender a "historia" dos dados... > > > > sobre os seus dados, pergunte-se: 1) eles sao sempre inteiros? 2) se > > sim, eles representam o numero de tentativas (de alguma coisa) ate' o > > primeiro sucesso (ou falha)? 3) se sim, cada uma dessas tentativas tem > > a mesma chance de sucesso (ou falha)? > > > > se sua resposta for "nao" para qq uma das perguntas acima, e' pouco > > provavel que uma geometrica seja a resposta q vc procura... > > > > b > > > > 2012/8/17 Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com > >: > >> Olá > >> > >> Estou modelando um conjunto de dados e a variável resposta parece ter > >> distribuição geométrica, qual é a melhor forma de modelar esses dados? > >> > >> Abraços > >> > >> > >> dput(dados) > >> structure(list(AFtotal = structure(c(103L, 76L, 58L, 105L, 125L, > >> 163L, 106L, 146L, 227L, 187L, 5L, 117L, 147L, 39L, 232L, 2L, > >> 67L, 150L, 46L, 118L, 18L, 171L, 157L, 226L, 80L, 96L, 142L, > >> 162L, 128L, 152L, 109L, 10L, 108L, 102L, 47L, 208L, 122L, 78L, > >> 176L, 3L, 196L, 40L, 209L, 119L, 179L, 29L, 1L, 81L, 190L, 24L, > >> 204L, 28L, 92L, 136L, 104L, 85L, 91L, 15L, 141L, 13L, 11L, 121L, > >> 112L, 210L, 31L, 6L, 70L, 156L, 114L, 68L, 107L, 112L, 211L, > >> 167L, 1L, 215L, 88L, 213L, 194L, 136L, 220L, 199L, 136L, 9L, > >> 136L, 118L, 38L, 192L, 42L, 62L, 155L, 19L, 79L, 48L, 138L, 172L, > >> 134L, 231L, 231L, 216L, 136L, 195L, 130L, 101L, 217L, 205L, 174L, > >> 36L, 111L, 140L, 106L, 99L, 22L, 218L, 170L, 7L, 8L, 136L, 19L, > >> 62L, 25L, 123L, 1L, 137L, 220L, 97L, 52L, 169L, 190L, 216L, 1L, > >> 53L, 172L, 100L, 79L, 86L, 201L, 61L, 198L, 17L, 218L, 67L, 106L, > >> 193L, 1L, 138L, 164L, 200L, 190L, 133L, 129L, 74L, 186L, 72L, > >> 41L, 177L, 48L, 113L, 50L, 229L, 178L, 106L, 124L, 75L, 14L, > >> 148L, 218L, 191L, 19L, 219L, 30L, 1L, 71L, 112L, 172L, 21L, 221L, > >> 16L, 203L, 27L, 69L, 143L, 82L, 207L, 223L, 32L, 54L, 60L, 59L, > >> 1L, 97L, 149L, 63L, 132L, 212L, 214L, 1L, 1L, 165L, 77L, 20L, > >> 73L, 26L, 168L, 160L, 233L, 139L, 1L, 188L, 95L, 175L, 206L, > >> 93L, 144L, 135L, 56L, 6L, 189L, 197L, 45L, 30L, 35L, 1L, 4L, > >> 136L, 154L, 89L, 65L, 184L, 225L, 230L, 151L, 1L, 173L, 43L, > >> 44L, 33L, 182L, 34L, 153L, 1L, 1L, 90L, 106L, 110L, 98L, 1L, > >> 219L, 1L, 1L, 116L, 49L, 64L, 57L, 131L, 185L, 145L, 17L, 94L, > >> 1L, 180L, 1L, 120L, 159L, 51L, 172L, 181L, 183L, 7L, 87L, 106L, > >> 228L, 55L, 126L, 157L, 218L, 158L, 115L, 133L, 222L, 84L, 166L, > >> 1L, 37L, 1L, 1L, 202L, 66L, 12L, 1L, 83L, 23L, 110L, 233L, 231L, > >> 161L, 106L, 224L, 127L), .Label = c("0,00", "1.031,00", "1.035,00", > >> "1.042,50", "1.064,00", "1.080,00", "1.120,00", "1.155,00", "1.171,50", > >> "1.177,50", "1.179,00", "1.200,00", "1.204,50", "1.210,00", "1.218,00", > >> "1.220,00", "1.230,00", "1.233,00", "1.260,00", "1.275,00", "1.287,00", > >> "1.290,00", "1.314,00", "1.320,00", "1.356,00", "1.359,50", "1.400,00", > >> "1.410,00", "1.418,00", "1.440,00", "1.466,00", "1.480,00", "1.500,00", > >> "1.566,00", "1.572,00", "1.590,00", "1.611,00", "1.652,00", "1.680,00", > >> "1.722,00", "1.734,00", "1.740,00", "1.752,00", "1.800,00", "1.818,00", > >> "1.821,00", "1.833,00", "1.836,00", "1.854,00", "1.872,00", "1.900,00", > >> "1.920,00", "1.938,00", "10.080,00", "10.230,00", "10.557,00", > >> "10.680,00", "11.238,00", "12.000,00", "12.270,00", "12.420,00", > >> "120,00", "13.644,00", "13.860,00", "14.590,00", "150,00", "165,00", > >> "180,00", "2.002,50", "2.060,00", "2.078,00", "2.100,00", "2.120,00", > >> "2.157,00", "2.320,00", "2.332,50", "2.376,00", "2.382,00", "2.430,00", > >> "2.457,00", "2.475,00", "2.504,00", "2.520,00", "2.526,00", "2.532,00", > >> "2.556,00", "2.572,00", "2.580,00", "2.616,00", "2.640,00", "2.645,00", > >> "2.655,00", "2.706,00", "2.739,00", "2.745,00", "2.767,50", "2.784,00", > >> "2.798,00", "2.820,00", "2.841,00", "2.895,00", "2.937,00", "2.940,00", > >> "2.965,00", "21.515,00", "240,00", "247,50", "259,00", "26.862,00", > >> "264,00", "273,00", "297,00", "3.019,00", "3.035,00", "3.036,00", > >> "3.282,00", "3.306,00", "3.360,00", "3.497,00", "3.504,00", "3.510,00", > >> "3.540,00", "3.600,00", "3.615,00", "3.660,00", "3.710,00", "3.840,00", > >> "3.892,50", "3.900,00", "3.960,00", "3.977,00", "3.996,00", "300,00", > >> "318,00", "330,00", "360,00", "372,00", "396,00", "4.020,00", > >> "4.086,00", "4.092,00", "4.100,00", "4.140,00", "4.158,00", "4.200,00", > >> "4.219,50", "4.331,00", "4.359,00", "4.554,00", "4.618,00", "4.680,00", > >> "4.770,00", "407,50", "410,00", "420,00", "424,00", "480,00", > >> "495,00", "5.052,00", "5.160,00", "5.340,00", "5.622,00", "5.688,00", > >> "5.715,00", "5.745,00", "5.832,00", "5.940,00", "5.958,00", "510,00", > >> "516,00", "537,00", "540,00", "565,50", "570,00", "575,00", "576,00", > >> "585,00", "6.097,50", "6.285,00", "6.495,00", "6.556,00", "6.576,00", > >> "6.720,00", "6.746,00", "6.828,00", "639,00", "645,00", "655,00", > >> "66,00", "660,00", "678,00", "690,00", "693,00", "696,00", "7.056,00", > >> "7.074,00", "7.200,00", "7.236,00", "7.245,00", "7.500,00", "7.600,00", > >> "7.680,00", "7.860,00", "7.920,00", "720,00", "738,00", "754,00", > >> "754,50", "780,00", "795,00", "8.268,00", "8.467,50", "8.592,00", > >> "8.755,00", "8.820,00", "80,00", "801,00", "810,00", "812,00", > >> "840,00", "876,00", "882,00", "895,50", "9.120,00", "9.150,00", > >> "9.310,00", "9.508,00", "924,00", "930,00", "933,00", "960,00", > >> "99,00", "990,00"), class = "factor")), .Names = "AFtotal", class = > >> "data.frame", row.names = c(NA, > >> -299L)) > >> > >> > >> -- > >> Sérgio Henrique Almeida da Silva Junior > >> Doutorando em Epidemiologia em Saúde Pública > >> Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > >> http://lattes.cnpq.br/1611345552843383 > >> Tel: (21) 94429486/78101651 id: 123*20942 > >> > >> > >> > >> _______________________________________________ > >> 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 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. > -- Fernando A.B. Colugnati

3 3

Re: [R-br] Utilizando DEA com R
by Roney Fraga 17 Ago '12

17 Ago '12

Cleyton Fábio, Já utilizei o pacote Benchmarking, segue link com alguns exemplos: http://r-br.2285057.n4.nabble.com/R-br-Analise-Envoltoria-de-Dados-DEA-td37… utilizei os dados do livro: EVALUACION DE LA EFICIENCIA MEDIANTE EL ANALISIS ENVOLVENTE DE DATOS dos autores Vicente Coll Serrano e Olga M Blasco Blasco de 2006. Esse livro esta disponível no google livros. Todos os cálculos tem exemplos no excel/calc, isso ajuda a entender o método e ainda permite comparar os resultados apresentados pelo livro com os realizados por você no excel/calc e no R. Gosto de utilizar o pacote Benchmarking porque tenho o livro Bogetoft e Otto (2011) Benchmarking with DEA, SFA, and R. Como você pode notar, esse livro é dos mesmos autores do pacote, acho o texto didático. Outro livro que gosto é Coelli, Rao, O'Donnell, Battese (2005) Introduction to Efficiency and Productivity Analysis 2ed. Atenciosamente Roney

3 2

Matriz de proximidade
by Gilbert Queiroz 17 Ago '12

17 Ago '12

Pessoal, Alguém sabe se é possível obter a matriz de proximidade pelo R e em qual pacote se faz isto? Abs. Gilbert.

1 0

dúvida - modelo
by Sérgio Henrique almeida da silva ju 17 Ago '12

17 Ago '12

Olá Estou modelando um conjunto de dados e a variável resposta parece ter distribuição geométrica, qual é a melhor forma de modelar esses dados? Abraços dput(dados) structure(list(AFtotal = structure(c(103L, 76L, 58L, 105L, 125L, 163L, 106L, 146L, 227L, 187L, 5L, 117L, 147L, 39L, 232L, 2L, 67L, 150L, 46L, 118L, 18L, 171L, 157L, 226L, 80L, 96L, 142L, 162L, 128L, 152L, 109L, 10L, 108L, 102L, 47L, 208L, 122L, 78L, 176L, 3L, 196L, 40L, 209L, 119L, 179L, 29L, 1L, 81L, 190L, 24L, 204L, 28L, 92L, 136L, 104L, 85L, 91L, 15L, 141L, 13L, 11L, 121L, 112L, 210L, 31L, 6L, 70L, 156L, 114L, 68L, 107L, 112L, 211L, 167L, 1L, 215L, 88L, 213L, 194L, 136L, 220L, 199L, 136L, 9L, 136L, 118L, 38L, 192L, 42L, 62L, 155L, 19L, 79L, 48L, 138L, 172L, 134L, 231L, 231L, 216L, 136L, 195L, 130L, 101L, 217L, 205L, 174L, 36L, 111L, 140L, 106L, 99L, 22L, 218L, 170L, 7L, 8L, 136L, 19L, 62L, 25L, 123L, 1L, 137L, 220L, 97L, 52L, 169L, 190L, 216L, 1L, 53L, 172L, 100L, 79L, 86L, 201L, 61L, 198L, 17L, 218L, 67L, 106L, 193L, 1L, 138L, 164L, 200L, 190L, 133L, 129L, 74L, 186L, 72L, 41L, 177L, 48L, 113L, 50L, 229L, 178L, 106L, 124L, 75L, 14L, 148L, 218L, 191L, 19L, 219L, 30L, 1L, 71L, 112L, 172L, 21L, 221L, 16L, 203L, 27L, 69L, 143L, 82L, 207L, 223L, 32L, 54L, 60L, 59L, 1L, 97L, 149L, 63L, 132L, 212L, 214L, 1L, 1L, 165L, 77L, 20L, 73L, 26L, 168L, 160L, 233L, 139L, 1L, 188L, 95L, 175L, 206L, 93L, 144L, 135L, 56L, 6L, 189L, 197L, 45L, 30L, 35L, 1L, 4L, 136L, 154L, 89L, 65L, 184L, 225L, 230L, 151L, 1L, 173L, 43L, 44L, 33L, 182L, 34L, 153L, 1L, 1L, 90L, 106L, 110L, 98L, 1L, 219L, 1L, 1L, 116L, 49L, 64L, 57L, 131L, 185L, 145L, 17L, 94L, 1L, 180L, 1L, 120L, 159L, 51L, 172L, 181L, 183L, 7L, 87L, 106L, 228L, 55L, 126L, 157L, 218L, 158L, 115L, 133L, 222L, 84L, 166L, 1L, 37L, 1L, 1L, 202L, 66L, 12L, 1L, 83L, 23L, 110L, 233L, 231L, 161L, 106L, 224L, 127L), .Label = c("0,00", "1.031,00", "1.035,00", "1.042,50", "1.064,00", "1.080,00", "1.120,00", "1.155,00", "1.171,50", "1.177,50", "1.179,00", "1.200,00", "1.204,50", "1.210,00", "1.218,00", "1.220,00", "1.230,00", "1.233,00", "1.260,00", "1.275,00", "1.287,00", "1.290,00", "1.314,00", "1.320,00", "1.356,00", "1.359,50", "1.400,00", "1.410,00", "1.418,00", "1.440,00", "1.466,00", "1.480,00", "1.500,00", "1.566,00", "1.572,00", "1.590,00", "1.611,00", "1.652,00", "1.680,00", "1.722,00", "1.734,00", "1.740,00", "1.752,00", "1.800,00", "1.818,00", "1.821,00", "1.833,00", "1.836,00", "1.854,00", "1.872,00", "1.900,00", "1.920,00", "1.938,00", "10.080,00", "10.230,00", "10.557,00", "10.680,00", "11.238,00", "12.000,00", "12.270,00", "12.420,00", "120,00", "13.644,00", "13.860,00", "14.590,00", "150,00", "165,00", "180,00", "2.002,50", "2.060,00", "2.078,00", "2.100,00", "2.120,00", "2.157,00", "2.320,00", "2.332,50", "2.376,00", "2.382,00", "2.430,00", "2.457,00", "2.475,00", "2.504,00", "2.520,00", "2.526,00", "2.532,00", "2.556,00", "2.572,00", "2.580,00", "2.616,00", "2.640,00", "2.645,00", "2.655,00", "2.706,00", "2.739,00", "2.745,00", "2.767,50", "2.784,00", "2.798,00", "2.820,00", "2.841,00", "2.895,00", "2.937,00", "2.940,00", "2.965,00", "21.515,00", "240,00", "247,50", "259,00", "26.862,00", "264,00", "273,00", "297,00", "3.019,00", "3.035,00", "3.036,00", "3.282,00", "3.306,00", "3.360,00", "3.497,00", "3.504,00", "3.510,00", "3.540,00", "3.600,00", "3.615,00", "3.660,00", "3.710,00", "3.840,00", "3.892,50", "3.900,00", "3.960,00", "3.977,00", "3.996,00", "300,00", "318,00", "330,00", "360,00", "372,00", "396,00", "4.020,00", "4.086,00", "4.092,00", "4.100,00", "4.140,00", "4.158,00", "4.200,00", "4.219,50", "4.331,00", "4.359,00", "4.554,00", "4.618,00", "4.680,00", "4.770,00", "407,50", "410,00", "420,00", "424,00", "480,00", "495,00", "5.052,00", "5.160,00", "5.340,00", "5.622,00", "5.688,00", "5.715,00", "5.745,00", "5.832,00", "5.940,00", "5.958,00", "510,00", "516,00", "537,00", "540,00", "565,50", "570,00", "575,00", "576,00", "585,00", "6.097,50", "6.285,00", "6.495,00", "6.556,00", "6.576,00", "6.720,00", "6.746,00", "6.828,00", "639,00", "645,00", "655,00", "66,00", "660,00", "678,00", "690,00", "693,00", "696,00", "7.056,00", "7.074,00", "7.200,00", "7.236,00", "7.245,00", "7.500,00", "7.600,00", "7.680,00", "7.860,00", "7.920,00", "720,00", "738,00", "754,00", "754,50", "780,00", "795,00", "8.268,00", "8.467,50", "8.592,00", "8.755,00", "8.820,00", "80,00", "801,00", "810,00", "812,00", "840,00", "876,00", "882,00", "895,50", "9.120,00", "9.150,00", "9.310,00", "9.508,00", "924,00", "930,00", "933,00", "960,00", "99,00", "990,00"), class = "factor")), .Names = "AFtotal", class = "data.frame", row.names = c(NA, -299L)) -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

4 3

Anova em faixas
by Juliano Farias 17 Ago '12

17 Ago '12

Srs, Preciso de uma ajuda. Preciso fazer uma ANOVA de um experimento com três fatores, porém o problema está no delineamento que foi em faixas. Tenho o quadro da análise da variança para esse tipo de delineamento, porém não consegui montar a programação e não achei um pacote para esse tipo de delineamento. O experimento foi realizado em blocos (B) ao acaso com quatro repetições em esquema fatorial em faixas 2 x 2 x 2. Os dois níveis do fator A, com os dois níveis do fator C foram casualizados nas parcelas principais, e os dois níveis do fator D foram casualizados em faixas perpendiculares às parcelas principais. Grato pela atenção. -- Atenciosamente, Juliano Ricardo Farias Eng. Agrônomo Doutorando em Entomologia Lab. de Resistência de Artrópodes a Pesticidas ESALQ/USP

3 4

Somar vendas por semana e por mês
by Hugo Sena 17 Ago '12

17 Ago '12

Boas! Vou fazer análise de séries temporais e preciso de colocar as vendas de alguns produtos tendo: em cada linha as vendas totais por semana e também tendo em cada linha, as vendas totais por mês. Tenho os dados em Excel, alguma ajuda para isto ser rápido e automático?. Cumps.

2 1

Regressão Gama
by Luiz Henrique da Conceicao Leal 15 Ago '12

15 Ago '12

Bom dia pessoal Construi um regressão gama e gostaria de construir um intervalo de confiança (no gráfico) para reta estimada. Alguém tem alguma sugestão? modelo=glm(y~x,family=Gamma(link = "identity")) plot(x,y,pch=20) abline(modelo,col="blue",lwd=2) Desde já agradeço Luiz Henrique Leal

1 0

Duvida no TAPPLY
by Edson Lira 15 Ago '12

15 Ago '12

Estou usando esta rotina e obtenho tapply(acj$Quantidade,list(Laboratório=acj$lab,Mês=acj$Mês),sum,na.rm=T)#Laboratórios E obtenho esta saída abaixo, que me dá o quantitativo de exames mensal, para obter a soma nas linhas e a soma nas colunas,o que devo acrescdentar a rotina? Mês Laboratório Jan Fev Mar Abr Mai Jun Jul Ago Set Out Bioquímica 24927 25529 24589 23736 29551 30175 28760 29220 24291 22189 Estudo de Hemoglobinopatias 300 226 232 201 293 360 209 243 233 208 H. Pardini 546 479 436 360 433 403 405 455 357 310 HEMATOLOGIA 8764 8877 8398 7969 9389 9360 9264 9585 8145 7548 HORMÔNIOS 1137 1276 1307 1271 1579 1522 1479 1519 1257 1286 IMUNOLOGIA ERITROCITÁRIA 1961 1942 1890 1820 2132 2206 1996 1992 1736 1564 PARASITOLOGIA/URINÀLISE 4365 4265 3477 3526 4183 4581 4581 4582 3898 3667 HEMOSTASIA 1174 987 1025 928 1266 1167 1257 1134 792 732 IMUNOLOGIA 1839 1570 1428 2220 1715 1522 1266 1490 1349 1253 SOROLOGIA 3327 3116 2863 3431 4437 4426 4294 4134 2834 1993 [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 4

Utilizando DEA com R
by Cleyton Fábio 15 Ago '12

15 Ago '12

Ola gente, tudo bem? Estou trabalhando em um projeto onde quero analisar a eficiência de municípios brasileiros utilizando Analise Envoltória de Dados. A ideia é implementar uma solução(100% baseada em software livre) que facilite este trabalho, conheci o R e descobri que existe alguns pacotes como o DEA e o Benchmarking, porem ainda estou tendo um pouco de dificuldades em utilizar a aplicação. Gostaria de saber se alguem tem alguma documentação sobre a aplicação de DEA ou Benchmarking, poderia ser artigos, monografias, tutoriais e a possibilidade de disponibilizar. O trabalho é a minha monografia! grato, Cleyton Fábio

5 4

Duvidas com spDists
by Julianna Trindade 15 Ago '12

15 Ago '12

Pessoal, tenho alguns pontos num shapefile e preciso saber as distâncias entre todos os pontos. Consegui fazer isso utilizando spDists mas a matriz que resulta não identifica os pontos pela coluna de atributo que eu preciso, apresenta apenas uma matriz quadrada utilizando 1,1; 1,2...2,1;2,2...e etc pra definir as linhas e as colunas e eu preciso que ao invés de 1,1 apareçam os atributos de uma coluna específica. Isso seria possível?? Além disso, qual unidade é utilizada nas distancias?Km?metros? O codigo eh o seguinte: rm(list=ls(all=TRUE)) library(base) library(maps) library(sp) library(maptools) library(spdep) library(classInt) library(RColorBrewer) library(shape) library(SDMTools) library(latticeExtra) pontos = readShapePoints("c:/users/Julianna/Documents/SecSaude/hospitais_produtividade.shp") proj4string(pontos) <- CRS("+init=epsg:4291") distancias = spDists(pontos) Resultando: pontosSpatialPointsDataFrame[1] distancias23x23 double matrix -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

2 2

[OFF] Biblioteca C para matrizes.
by Pedro Rafael 14 Ago '12

14 Ago '12

Pessoal gostaria de saber qual uma boa biblioteca livre de C para trabalhar com operações matriciais como inversas e produtos de matrizes. Qualquer dica será interessante. [ ] Pedro Rafael Diniz Marinho.

3 28

Vaga para Estágio - Desenvolvimento em R - Porto Alegre
by Fernando Henrique Ferraz Pereira da Rosa 14 Ago '12

14 Ago '12

PROGRAMA DE ESTÁGIO O Programa de estágio da Souza Cruz desenvolve jovens talentos, agregando valor à sua formação acadêmico-profissional. Ao final do Programa, os estagiários têm a oportunidade de participar do Programa de Trainee ou serem indicados para outros desafios na Empresa. Vaga para os cursos de Estatística / Economia / Engenharia de Produção / Processamento de Dados / Ciência da Computação e correlatos Bolsa auxílio no valor de R$1.053 (6h) e R$702 (4h) , auxílio transporte, restaurante, plano de saúde (opcional) e seguro de vida. Área de lotação: Product Centre – Suporte Estatístico. A área provê suporte estatístico para as diversas áreas da empresa (laboratórios, fábricas, desenvolvimento de produtos, entre outros). Resumo de atividades: atender às necessidades de manutenção e desenvolvimento de softwares e ferramentas estatísticas da área e dos seus clientes. A demanda por uso e aplicação de conhecimentos estatísticos será balanceada de acordo com a área de formação do candidato. A posição poderá ser assumida por alguém sem conhecimento prévio em estatística, caso no qual será dado mais enfoque ao desenvolvimento de software e programação, ou caso seja preenchida por alguém com formação em estatística, haverá maior enfoque nos métodos estatísticos. Áreas de atuação: Planos amostrais Planejamento de Experimentos Programas intra e interlaboratoriais Relatórios Gerenciais (MIS) Otimização de Processos Controle Estatístico de Processo / 6-sigma PRÉ-REQUISITOS Universitários cursando o penúltimo ou último ano de graduação. Bons conhecimentos em informática (preferencialmente deve conhecer ao menos uma linguagem de programação). DESEJÁVEL Conhecimento em R e/ou SAS, Excel VBA, lógica de programação. LOCAL DE TRABALHO Unidade Fabril de Cachoeirinha (Região Metropolitana de Porto Alegre) INSCRIÇÕES Enviar CV para fernando.henrique.rosa(a)souzacruz.com.br -- "Though this be randomness, yet there is structure in't." Fernando H Rosa - Statistician http://www.fernandohrosa.com.br / http://www.feferraz.net - Estatística, Matemática e Computação AprendaAlemao.net - Seu ponto de partida para melhorar seu Alemão! @fhrosa

1 0