- R-br - DINF UFPR

Pedindo ajuda para simulação em R
by Mauro Sznelwar 21 Ago '11

21 Ago '11

1) Este algoritmo da distribuição geométrica(p) Conhecendo a cdf: F(x)=1-(1-p)^x+1 F(x-1)=1-(1-p)^x Usando o método da transformada inversa F(x-1)<R<=F(x) 1-(1-p)^x<R<=1-(1-p)^x+1 (1-p)^x+1<R-1<=1-(1-p)^x ln(1-R)/ln(1-p)-1<=X<ln(1-R)/ln(1-p) X=[(ln(1-R)/ln(1-p))-1] 2)Gerar a distribuição de Poisson através do processo de Poisson que pode ser simulado gerando variáveis aleatórias exponenciais até sua soma exceder 1. Lembrar o método da inversa para distribuição exponencial X1= -1/A(ln(Ri)) Então A1+....+Ax<=1<A1+A2+....+Ax+Ax+1 torna-se -somatória de i=1 até x(1/A(ln(Ri)<=1< -somatória(i=1 até x+1(1/lambda(ln(Ri)) somatória de i=1 até x((ln(Ri)>=lambda>somatória(i=1 até x+1(ln(Ri)) ln(produtória(i=1 até x(Ri)>= -lambda>ln(produtória(i=1 até x(Ri) (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Isto leva ao algoritmo da rejeição Usando (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Passo 1: N=0 ; P=1 Passo 2: Gere número aleatório Rn-1 e substitua P por P*Rn-1 Passo 3: Se P<e^-lambda, então aceita N=n, senão rejeita o n, soma n por um e retorne ao passo 2. 3) Gerar a distribuição de Pareto que é definida por f(x|alfa)=alfa*x^-alfa-1 sobre (1,infinito). Mostre que pode ser gerada como potencia -1/alfa de uma variável uniforme. Eu tentei fazer a inversa mas não sei se é isto: a funçao de distribuição deve ser F(x)=1-x^-alfa e ao fazer a tranformada inversa: log(1-x)=log(u)/-alfa 4)Gerar a distribuição Beta onde U e V são i.i.d e a distribuição é (U^1/alfa)/((U^1/alfa)+V^(1/Beta)), condicional em U^1/alfa+V^1/B<=1.Tem que ser gerado com a potencia -1/alfa de uma uniforme. Desde já agradeço [edit] -------------------------------------------------------------------------------- _______________________________________________ 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.

2 2

[Dúvida] Salvar gráfico com boa qualidade
by Pedro Rafael 21 Ago '11

21 Ago '11

Pessoal como salvar um gráfico no R com uma ótima qualidade para que o mesmo possa ser ampliado sem prejudicar a qualidade da imagem. Existe formas de salvar a imagem e definir o tamanho da figura? Saudações, Pedro Rafael

4 6

OFF - Topic - procurando um artigo
by Pedro Emmanuel Alvarenga Americano do Brasil 20 Ago '11

20 Ago '11

Amigos de R, Será que alguns de voces poderia ter acesso ou tem o pdf desse artigo aqui? http://www.tandfonline.com/doi/abs/10.1080/03610910701212181#preview , Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 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

3 3

Função tapply considera comprimentos deiferentes com NA
by Alexandre Santos 20 Ago '11

20 Ago '11

Bom Dia Pessoal, Estava tentando tirar à média de 3 fatores para densidade de insetos, sendo 4, 16 e 64 insetos e cada fator com 10 repetições, porém em 4 insetos tenho 3 NA' s e quando tento usar tapply() a função dá erro obviamente, pois são comprimentos diferentes, sendo: > medias<-tapply(na.omit(n.provas),densidade,mean) Erro em tapply(na.omit(n.provas), densidade, mean) : argumentos devem ter o mesmo comprimento > medias 16prey 4prey 64prey 4.0 NA 12.4 Alguém saberia sugerir outra função, sendo que já tentei aggregate() também sem sucesso, Obrigado, Alexandre

3 3

R simulação: Entendendo as distribuições de probabilidade
by Jose Claudio Faria 20 Ago '11

20 Ago '11

Uma das coisas que considero mais interessantes ao adotar o R em um cursos introdutórios de Probabilidade e Estatística, é possibilitar que os estudante entendam (efetivamente) o que são as distribuições de probabilidades (inf. paramétrica) que eles estão começando a usar para tomar decisões em situação de incerteza. (Em outra palavras,: jogar a mágica na lata de lixo!) Uma das que mais gosto, que pode vir a ser útil para outras pessoas está abaixo: (Estou socializado: sujeito a críticas e sugestões...) # Simulação da origem da distribuição F (Snedecor) # Autor: J.C.Faria # Tamanho das amostras nn <- 05 # Tamanho da amostra do numerador nd <- 10 # Tamanho da amostra do denominador # Características da população mu <- 10 sigma <- 2 # Normal padrão # mu <- 0 # sigma <- 1 # Número de pares de amostras n <- 1e4 # Simulação s2n <- apply(matrix(rnorm(n*nn, mu, sigma), ncol=nn), 1, var) s2d <- apply(matrix(rnorm(n*nd, mu, sigma), ncol=nd), 1, var) require(fdth) plot(fdt(s2n/s2d, start=0, end=10, h=.05), type='d', xlim=c(0, 10), ylim=c(0, 0.8)) curve(df(x, nn-1, nd-1), col='darkgreen', add=TRUE, lwd=2) A conclusão óbvia a que se chega: A distribuição F descreve como se distribui a razão entre duas estimativas da variância de uma população normal! Portanto: Pode ser usada para decidir (sob as devidas condições) se duas estimativas da variância podem, ou não, serem consideradas como provenientes de uma mesma população normal. (Que é o uso mais básico da distribuição F, ou seja, teste de precisão). Gostaria de ver a galera desentocar esse tipo de script e postar aqui na lista! (Tenho mais uns poucos nessa linha) Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

5 10

funcao para extrair o ano de uma variável data
by Paula Mendes Luz 19 Ago '11

19 Ago '11

oi pessoal acredito que deva existir uma funcao para isso... se alguem souber me avisa por favor! tenha uma variavel que é uma data (as.Date) no R e gostaria de extrair o ano. como posso fazer isso? obrigada, paula.

2 2

BATCH
by FHRB Toledo 19 Ago '11

19 Ago '11

Senhores, Gostaria de saber se existe uma forma elegante de rodar um script em BATCH passando argumentos como opção na linha de comando? Explico: tenho um script que faz uma tarefa "X", para tal no script tenho alguns argumentos especificados do tipo arg1 <- valorA. Entretanto queria rodar o mesmo escript variando o argumento parao valorB por exemplo. Na minha ignorância estou editando o script no valorB do arg e rodando denovo. Pois bem, existe uma forma melhor e menos tosca de se fazer isso? desde já grato, FH

2 2

Pedindo ajuda para simulação
by Mauro Sznelwar 19 Ago '11

19 Ago '11

1) Como para simular este algoritmo da distribuição geométrica(p)? Conhecendo a cdf: F(x)=1-(1-p)^x+1 F(x-1)=1-(1-p)^x Usando o método da transformada inversa F(x-1)<R<=F(x) 1-(1-p)^x<R<=1-(1-p)^x+1 (1-p)^x+1<R-1<=1-(1-p)^x ln(1-R)/ln(1-p)-1<=X<ln(1-R)/ln(1-p) X=[(ln(1-R)/ln(1-p))-1] 2)Gerar a distribuição de Poisson através do processo de Poisson que pode ser simulado gerando variáveis aleatórias exponenciais até sua soma exceder 1. Lembrar o método da inversa para distribuição exponencial X1= -1/A(ln(Ri)) Então A1+....+Ax<=1<A1+A2+....+Ax+Ax+1 torna-se -somatória de i=1 até x(1/A(ln(Ri)<=1< -somatória(i=1 até x+1(1/lambda(ln(Ri)) somatória de i=1 até x((ln(Ri)>=lambda>somatória(i=1 até x+1(ln(Ri)) ln(produtória(i=1 até x(Ri)>= -lambda>ln(produtória(i=1 até x(Ri) (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Isto leva ao algoritmo da rejeição Usando (produtória(i=1 até x(Ri)>=e^ -lambda>(produtória(i=1 até x(Ri) Passo 1: N=0 ; P=1 Passo 2: Gere número aleatório Rn-1 e substitua P por P*Rn-1 Passo 3: Se P<e^-lambda, então aceita N=n, senão rejeita o n, soma n por um e retorne ao passo 2. [edit]

1 0

Imprimir mais de um data.frame em um único arquivo
by Carolina Pena 19 Ago '11

19 Ago '11

Olá, Estou com o seguinte problema: gostaria de imprimir três data.frames em um mesmo arquivo csv. Eu gostaria de uma saída no seguinte formato: Titulo data.frame 1: Dados de Belo Horizontes Codigo Risco 1 2742 0.0004 2 2267 0.0006 3 2554 0.0007 4 2963 0.0008 Titulo data.frame 2: Dados de Minas Gerais Codigo Risco 1 4743 0.0004 2 4200 0.0003 3 4359 0.0009 4 4555 0.0001 Titulo data.frame 2: Dados do Brasil Codigo Risco 1 5743 0.0006 2 5200 0.0008 3 5359 0.0002 4 5555 0.0004 Eu pensei em colocar tudo em uma lista, mas até agora não consegui encontrar nenhum comando que cria arquivos a partir de uma lista. Eu costumo utilizar o write.table, mas nesse caso teria que criar um arquivo para cada data.frame. Obrigada! Abraços, Carolina

2 1

[OFF-Topic] Agenda de webnários sobre Mendeley
by Walmes Zeviani 18 Ago '11

18 Ago '11

Pessoal, Acredito que nem todos os interessados assistiram o webnário de terça-feira sobre Mendeley anunciado pelo Fernando Toledo. Bem, eu assisti e achei muito conveniente divulgar a agenda disponível no site para os próximos webnários. Recomendo que tomem conhecimento desde aplicativo. Para se cadastrar para os próximos webnários visite o link e escolha uma data. https://mendeley.webex.com/mw0306ld/mywebex/default.do?siteurl=mendeley 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 ==========================================================================

1 0