- R-br - DINF UFPR

Re: [R-br] Representação usando expression()
by Alexandre Santos 25 Ago '11

25 Ago '11

Obrigado Benilton, Resolveu o meu problema e ficaram muito mais elegante os argumentos também, Abraço, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122 Tomei a liberdade de simplificar um pouco sua expressao (e' equivalente ao modelo que vc descreve): txt = expression(frac(1, 1+e^{0.9588-0.0841*x+0.0010*x^2})) plot(1, type='n') legend(1, 1, txt, bty='n') b

1 0

Representação usando expression()
by Alexandre Santos 25 Ago '11

25 Ago '11

Bom Bom dia pessoal, Gostaria de saber s é possível fazer com que em: text(33,0.9,expression(italic("y")==italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x"))/(1+italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x")))),cex=0.8) A equação aparecesse no gráfico como: ("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x")) sobre (1+italic("exp")(-"0.9588"+"0.0841"*italic("*x")-"0.0010"*italic("*x") e não separado por / , isso é possível? Obrigado, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

2 1

Otimizar tempo de processamento
by Tiago Vieira 25 Ago '11

25 Ago '11

Estou tentando rodar o programa de simulação de percolação (rede de bethe), mas esta demandando um tempo de processamento muito alto, pois há preciso de muitas iterações. O computador que estou utilizando tem processador i3 e 4gb de memória. Alguém pode me dar uma dica de como posso diminuir o tempo de processamento? Abaixo segue o código: programa <- function(f,n,m) { # f: número máximo de filhos # n: número de gerações # m: número de iterações S <- matrix(0,101,2) for (p in seq(0,1,by=0.01)) { pp <- 0 for (g in 1:m) { y <- 1 for (k in 1:n) { y <- sum(rbinom(y,f,p)) if (y == 0) break } if (y > 0) pp <- pp + 1 } S[p/0.01+1,1] <- p S[p/0.01+1,2] <- pp/m } return(S) } prog <- programa(2,100,1000) plot(prog, xlab = "p", ylab = "Probabilidade de Percolação") Desde já agradeço pela atenção! Tiago Vieira.

6 12

Problema de memória: MATRIZES GIGANTES
by Márcio Luís Moreira de Souza 24 Ago '11

24 Ago '11

Gente... Tô com outro problema. Agora é em relação à memória do R. Não tô conseguindo gerar matrizes grandes (no meu caso são gigantes, na casa dos milhões de linhas por dezenas de colunas). Alguma dica?? Dá pra aumentar a memória utilizada no R (tenho um processador I5)?? E em computadores mais "fraquinhos", como faço?? Abraços!! -- Prof. Márcio Souza, M.Sc. (DE/UFJF) www.profmsouza.blogspot.com http://lattes.cnpq.br/4578008002785684

5 10

Invitation to connect on LinkedIn
by Marcelo Costa via LinkedIn 24 Ago '11

24 Ago '11

LinkedIn ------------ Marcelo Costa requested to add you as a connection on LinkedIn: ------------------------------------------ Fábio, I'd like to add you to my professional network on LinkedIn. - Marcelo Accept invitation from Marcelo Costa http://www.linkedin.com/e/-33g73l-grqxst4r-3t/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View invitation from Marcelo Costa http://www.linkedin.com/e/-33g73l-grqxst4r-3t/vhj1xG42iBdUFK5d_HMEGL4iFMC51… ------------------------------------------ DID YOU KNOW LinkedIn can help you find the right service providers using recommendations from your trusted network? Using LinkedIn Services, you can take the risky guesswork out of selecting service providers by reading the recommendations of credible, trustworthy members of your network. http://www.linkedin.com/e/-33g73l-grqxst4r-3t/svp/inv-25/ -- (c) 2011, LinkedIn Corporation

1 0

Re: [R-br] [OFF] Web Framework usando R
by Pedro Rafael 24 Ago '11

24 Ago '11

Aproveitando o mesmo tópico, qual seria uma boa linguagem para a pessoa aprender fora o R para trabalhar com estatística? Em 24 de agosto de 2011 17:51, Pedro Rafael <pedro.rafael.marinho(a)gmail.com>escreveu: > Aproveitando o mesmo tópico, qual seria uma boa linguagem para a pessoa > aprender fora o R para trabalhar com estatística? > > > > 2011/8/24 beniltoncarvalho [via R-br] < > ml-node+3766456-521461530-223914(a)n4.nabble.com> > > algo alem do que ha' descrito na documentacao? >> >> http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-Web-Interfaces >> >> E, aqui, algo que possa ser relevante (eh um bom comeco, IMHO): >> >> >> http://nsaunders.wordpress.com/2010/04/19/getting-your-web-application-and-… >> >> b >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=3766456&i=0> >> 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-OFF-Web-Framework-usando-R-tp3766338… >> 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

1 0

Factors no R
by Leandro Marino 24 Ago '11

24 Ago '11

Caros, Em algum ponto do tempo houve uma mudança na função factor. ### no passado: R.version # _ #platform i386-pc-mingw32 #arch i386 #os mingw32 #system i386, mingw32 #status #major 2 #minor 9.0 #year 2009 #month 04 #day 17 #svn rev 48333 #language R #version.string R version 2.9.0 (2009-04-17) set.seed(1) x <- round(runif(150,0,10),0) table(x) #x # 0 1 2 3 4 5 6 7 8 9 10 # 4 14 15 13 16 19 16 18 15 14 6 table(as.character(factor(x,levels=c(0:10),labels=c(rep('Tipo 1',3), rep('Tipo 2',4),rep('Tipo 3',2),'Tipo 4','Tipo 5')))) #Tipo 1 Tipo 2 Tipo 3 Tipo 4 Tipo 5 # 33 64 33 14 6 R.version # _ #platform i386-pc-mingw32 #arch i386 #os mingw32 #system i386, mingw32 #status #major 2 #minor 13.1 #year 2011 #month 07 #day 08 #svn rev 56322 #language R #version.string R version 2.13.1 (2011-07-08) set.seed(1) x <- round(runif(150,0,10),0) table(x) #x # 0 1 2 3 4 5 6 7 8 9 10 # 4 14 15 13 16 19 16 18 15 14 6 table(as.character(factor(x,levels=c(0:10),labels=c(rep('Tipo 1',3), rep('Tipo 2',4),rep('Tipo 3',2),'Tipo 4','Tipo 5')))) # #Tipo 1 Tipo 2 Tipo 3 Tipo 4 Tipo 5 # 33 64 33 14 6 #Mensagens de aviso perdidas: #In `levels<-`(`*tmp*`, value = c("Tipo 1", "Tipo 1", "Tipo 1", "Tipo 2", : # não se permite mais níveis duplicados em fatores Ele até faz mas com mensagem de erro. Como poderia fazer isto sem ter o erro!? Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

3 6

[OFF] Web Framework usando R
by Pedro Rafael 24 Ago '11

24 Ago '11

Pessoal não sei se concordam mas a informática "atualmente" foca muito em Web. Muitas aplicações vem sendo feita para Web. Muitas empresas monitoram suas estatísticas de forma online via shmartphones, tablets via Web. Existe algum Framework para desenvolvimento de aplicações Web que utilize o R como linguagem principal? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

3 2

(Resolvido) aov() e anova() no agricolae
by Alexandre Santos 24 Ago '11

24 Ago '11

Obrigado Ivan, Aumentou meu grau de confiança nas decisões que estou tomando, Abraço, Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122 Alexandre, Até onde eu sei, o HSD.test é utilizado para erros normais. Mais você pode utilizar sem problemas o seu modelo glm com a função glht para comparar as médias. Não vejo problemas também, você utilizar tais análises para uma variável e tais análises para outra variável, desde que esteja justificada o uso de tais metodologias. Allaman (S,f,P) \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-5076 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. Alexandre dos Santos Engenheiro Florestal, MSc. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 3829-5122

1 0

ajuda para um scrip
by Mauro Sznelwar 24 Ago '11

24 Ago '11

Tenho este algoritmo p=1, N=0, c=e^-lambda Repita N=N+1 gere Ui atualize p<c.e^-lambda Faça X=N-1 Estou fazendo esta rotina abaixo, mas o resultado não está consistente! Qual o melhor jeito? p=1 n=0 lambda=2 c=exp(-lambda) for(n in 1:100){ u=runif(n) p=p*u p<c } x=n-1

4 10