- R-br - DINF UFPR

Re: [R-br] xtabs ou equivalente
by Jose Claudio Faria 22 Jul '11

22 Jul '11

Test! -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

PLOT SERIE TEMPORAL
by Enéias Façanha 22 Jul '11

22 Jul '11

PessoALL, Utilizando os dados abaixo, como eu faço para plotar essas duas séries, "MILLE_META" e "MILE_RECEITA", no mesmo gráfico, inserindo os meses no eixo x? MILLE_META MILLE_RECEITA JAN 339.155,42 275.649,07 FEV 320.427,98 267.601,66 MAR 333.703,99 279.127,36 ABR 318.293,90 184.393,12 MAI 342.795,95 262.414,92 JUN 328.074,05 268.703,78 Obrigado, -- Enéias Façanha

2 1

Minimização Condicionada
by Luciano Ramos Gonçalves 22 Jul '11

22 Jul '11

Gostaria de saber qual é a melhor função para realizar uma minimização com restrições: nlm(), optim(), optimize() ou constrOptim(), ou outra? Abaixo vai a função objeto, e as restrições descritas. Qualquer dica já ajuda. O meu problema maior foi aplicar as restrições aos métodos de otimização. Obs.: Salvo engano, tudo abaixo é plenamente reproduzível. # Tenho a função abaixo. LogLik <- function(parametros,dados){ alpha <- parametros[1] delta <- parametros[2] epsilon.c <- parametros[3] epsilon.v <- parametros[4] mu <- parametros[5] B <- dados[,2] S <- dados[,3] logl <- sum(-epsilon.c-epsilon.v+((max(c(B,S))+min(c(B,S)))/2)* (log((epsilon.v/(mu+epsilon.v)))+ log((epsilon.c/(mu+epsilon.c))))+ B*log(mu+epsilon.c)+S*log(mu+epsilon.v))+ sum( log( (alpha)*(1-delta)*exp(-mu)* ((((epsilon.v/(mu+epsilon.v))))^ (S-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.c/(mu+epsilon.c))^ (-((max(c(B,S))+min(c(B,S)))/2)))+ (alpha)*(delta)*exp(-mu)* ((((epsilon.c/(mu+epsilon.c))))^ (B-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.v/(mu+epsilon.v))^ (-((max(c(B,S))+min(c(B,S)))/2)))+ (1-alpha)* ((epsilon.c/(mu+epsilon.c))^ (B-((max(c(B,S))+min(c(B,S)))/2)))* ((epsilon.v/(mu+epsilon.v))^ (S-((max(c(B,S))+min(c(B,S)))/2))) )) return(logl) } # Supondo os valores BS <- cbind(1,500,501) theta <- c(0.4,0.5,40,40,50) # Quero minimizar a função: LogLik(theta,BS) # Com as seguintes restrições, que eu não consegui implementar (as restrições são # aos valores de theta, ou "parametros" na função, em ordem): # 0 <= alpha <= 1 ("alpha" é um valor entre 0 e 1) # 0 <= delta <= 1 ("delta" é um valor entre 0 e 1) # epsilon.c > 20 ("epsilon.c" é maior que 20) # epsilon.v > 20 ("epsilon.v" é maior que 20) # mu > 30 ("mu" é maior que 30) Grato, Luciano.

1 0

[Dúvida] Salvando arquivo CSV sem aspas nas strings
by Pedro Rafael 21 Jul '11

21 Jul '11

Pessoal estou salvando um data.frame com o objeto write.table(objeto, file = "nome_do_arquivo.csv", sep = ";"). Contudo quero que as strings não contenham aspas. Como fazer para que o csv não contenham nas strings? Saudações, Pedro

3 2

xtabs ou equivalente
by Vinicius Brito Rocha 21 Jul '11

21 Jul '11

Senhores, estou usando o xtabs, no intuito de fazer alto parecido com o que se faz em uma tabela dinamica tenho um data.frame com dois fatores e uma variavel resposta continua estou fazendo algo do tipo: xtabs(y~A+B,data=base) tenho como resposta uma tabela com os níveis de A e B e uma estatística dentro de cada célula. Minha pergunta é: Como faço para mudar a estatística do xtabs? média, quantidade, mediana,.. não tem um parâmetro FUN dentro do xtabs. tem algum pacote com uma função equivalente que faça isso? Abs. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional 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.

4 8

R is Hot!
by Alejandro C. Frery 21 Jul '11

21 Jul '11

Boa noite, Recomendo a leitura da matéria "R is Hot!" http://t.co/MFs3Bz9 Alejandro -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

Incluir formula em um função
by Leo Bastos 20 Jul '11

20 Jul '11

Olá pessoal, Tenho algumas funções que eu gostaria de incluir uma formula como parametro de entrada. Da mesma forma que as funcoes lm e glm. Alguem poderia me indicar um exemplo, ou indicar um artigo/tutorial/link/whatever para que eu possa incluir uma formula com entrada de uma funcao? Por exemplo, suponha que eu tenho um data frame com um vetor de resposta y e algumas covariáveis. E suponha q eu quero uma função q retorne uma estimativa dos coeficientes de uma regressão linear. (sim, eu sei q o lm faz isso, mas eu gostaria de um exemplo de função onde eu possa usar uma formula como parametro de entrada). Uma função com um vetor y e uma matrix X é trivial, mas como fazer usando uma formula como parametro de entrada? Agradeço antecipadamente qq ajuda. Abraço, Leo Bastos

2 1

Fwd: [Dúvida] Ambiente para programação em R
by Pedro Rafael 19 Jul '11

19 Jul '11

Eu já conhecia o RKward e é bem interessante. Recentemente conheci o RStudio e acredito que ele está sendo mantido com mais regularidade ou estou enganado? Pedro ---------- Mensagem encaminhada ---------- De: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> Data: 18 de julho de 2011 08:51 Assunto: Re: [R-br] [Dúvida] Ambiente para programação em R Para: "jcfaria [via R-br]" <ml-node+3675065-198716963-223914(a)n4.nabble.com> Eu já conhecia o RKward e é bem interessante. Recentemente conheci o RStudio e acredito que ele está sendo mantido com mais regularidade ou estou enganado? Pedro Em 18 de julho de 2011 08:44, jcfaria [via R-br] < ml-node+3675065-198716963-223914(a)n4.nabble.com> escreveu: > Eu > > uso o vim. > > > Eu geralmente > > programo com uma janela para o vim e duas rodando o R: uma para as > > páginas de ajuda e outra para o código. > > Gustavo, > > Um opção muito interessante (eu prefiro assim) é habilitar o > vim/vim-r-plugin para abrir os arquivos de help do R dentro do próprio > vim (<leader>rh) em uma nova tab (para alternar entre tabs: gt (next > Tab) gT (prior Tab)) . > > Vantagens: > 1. Fica bem mais fácil para ler o help (sintaxe); > 2. Fica bem mais fácil de navegar no help (marcas, atalhos do vim, etc); > 3. Dá para interagir com o R (com todos os recursos) a partir do vim, > por exemplo, executar linha por linha os exemplos para entender a > função que está no help; > 4. As tabs vão se abrindo sequencialmente (de uma vc pode abrir outra, etc) > > > Teste colocando a linha abaixo no seu .vimrc: > > let vimrplugin_vimpager = "tabnew" " Open help always in new tab > > Abs, > -- > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > Jose Claudio Faria > Estatistica - prof. Pleno > UESC/DCET/Brasil > [hidden email] <http://user/SendEmail.jtp?type=node&node=3675065&i=0> > ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ > _______________________________________________ > R-br mailing list > [hidden email] <http://user/SendEmail.jtp?type=node&node=3675065&i=1> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Duvida-Ambiente-para-programacao-em-… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

4 3

[Dúvida] Ambiente para programação em R
by Pedro Rafael 18 Jul '11

18 Jul '11

Pessoal qual um bom ambiente para gerenciar projetos bem como para programação em R no Linux? Pedro

4 4

Produzindo um Codigo Minimo Reproduzivel (CMR)
by Benilton Carvalho 18 Jul '11

18 Jul '11

Na esperanca de que melhores CMRs sejam empregados na lista, adaptei um material que descreve como apropriadamente criar um CMR. Para detalhes, visitem: http://www.leg.ufpr.br/doku.php/software:rbr e sigam ate' o final da pagina. Ou, para uma versao "resumida": http://gist.github.com/1088208 b

2 1