Setembro 2015 - R-br

precisão da operação
by Cleber N.Borges 22 Set '15

22 Set '15

> 12345678912345678+1 [1] 12345678912345680 > 1.2345678912345678+0 [1] 1.2345678912345679 Prezados, pq exatamente o fato acima apresentado acontece? se a capacidade de representação dos números estão na faixa de 1.xxe-308 a 1.xxe+308 agradeço antecipadamente pela atenção cleber --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

4 3

Automatizar a inclusão de váriaveis
by Nicolay Cunha 22 Set '15

22 Set '15

Prezados, Gostaria de ajuda para automatizar a inclusão de variáveis dependentes em um modelo de RDA. Usarei o valor de F como critério para seleção do número de variáveis Y a serem usadas posteriormente, entretanto, preciso de uma forma automatizada para incluir as 59 colunas que tenho em minha matriz. Basicamente o que preciso fazer é calcular RDA s com o número crescente de Ys. Abaixo um CMR. # library(vegan) y <- matrix(sample(rnorm(3540)+10), 60) x <- rep(c("A", "B", "C"), each = 20) # Aqui gero o modelo com somente as colunas 1 e 2 da matriz y m1 <- capscale((y[,1:2]) ~ x, dist="bray") anova(m1) # agora incluindo a terceira coluna m2 <- capscale((y[,1:3]) ~ x, dist="bray") anova(m2) # adicionando a quarta coluna m3 <- capscale((y[,1:4]) ~ x, dist="bray") anova(m3) # e assim até incluir as 59 colunas na RDA m59 <- capscale((y[,1:59]) ~ x, dist="bray") anova(m59) Feito isso, preciso acessar os valores de F gerados para cada modelo para fazer a comparação. Alguém saberia me ajudar com isso? Muito Obrigado, Nicolay

1 1

Re: [R-br] Digest R-br, volume 57, assunto 13
by Carlos Matos 21 Set '15

21 Set '15

Prezados: Talvez os rankings de linguagens de programação abaixo possam ajudar: http://www.tiobe.com/index.php/content/paperinfo/tpci/index.html http://redmonk.com/sogrady/2015/07/01/language-rankings-6-15/ At.te, Em 20/9/2015 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. Estísticas da evolução do uso do R (Andre Oliveira) > 2. Re: Estísticas da evolução do uso do R (Walmes Zeviani) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 20 Sep 2015 12:08:15 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Estísticas da evolução do uso do R > Message-ID: > <964026555.429902.1442750895740.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Pessoas bom dia,me preparo para apresentação sobre o R. Gostaria de mostrar um gráfico oficial acerca do uso do R em reação a outros software ao longo do tempo. Gostaria de mostrar a evolução do número de bibliotecas ao longo da tempo e por último quantas bibliotecas oficiais existem disponíveis hoje nos repositórios. Alguém poderia indicar ou sugerir como posso conseguir tais informações? Julgo que estas informações seriam interessantes para a divulgação na apresentação. > André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150920/d859cfb7/atta…> > > ------------------------------ > > Message: 2 > Date: Sun, 20 Sep 2015 10:40:56 -0300 > From: Walmes Zeviani <walmeszeviani(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, Andre > Oliveira <andreolsouza(a)yahoo.com.br> > Subject: Re: [R-br] Estísticas da evolução do uso do R > Message-ID: > <CAFU=EkZoynY+0623ZuUsz50yALxg=WjpPxWWu0O-Rgqw=6akvA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Nessa exposição de números eu incluiria as listas de discussão. O markmail > reúne as oficiais. O gráfico mostra o número de mensagens. > > http://r-project.markmail.org/search/ > > Você pode explorar as 24 listas disponíveis e os mecanismos de buscas. > Lembrar que a queda na troca de mensagens depois de 2010 não representa > redução de popularidade mas sim que o histórico de mensagens disponível > retroalimenta a lista, ou seja, os novos usuários precisam postar menos > pois o histórico é rico, existem menos problemas novos, além disso muitos > novos usuários usam o stackoverflow. > > A nossa lista pode ser considerada também. Pegue os números em > http://listas.inf.ufpr.br/pipermail/r-br/. O tamanho do zip reflete o > tráfego. > > Google com "R popularity softwareA softwareB" a adicione na busca o nome de > outros softwares para restringir as busca. > > À disposição. > Walmes. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150920/7ba77fa4/atta…> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > 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 57, assunto 13 > ***************************************** -- ------------------------------------------------------------------------ *Carlos Alberto Oliveira de Matos* Professor Assistente Doutor Unesp - Campus de Itapeva Telefone: (15) 3524-9100 Ramal 9129 *Opção do Windows 8.1: Atualize seu computador sem afetar os arquivos.* *Resultado: perda de emails desde 26/5/14 até 19/3/2015* ****

2 1

Estísticas da evolução do uso do R
by Andre Oliveira 21 Set '15

21 Set '15

Pessoas bom dia,me preparo para apresentação sobre o R. Gostaria de mostrar um gráfico oficial acerca do uso do R em reação a outros software ao longo do tempo. Gostaria de mostrar a evolução do número de bibliotecas ao longo da tempo e por último quantas bibliotecas oficiais existem disponíveis hoje nos repositórios. Alguém poderia indicar ou sugerir como posso conseguir tais informações? Julgo que estas informações seriam interessantes para a divulgação na apresentação. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

5 4

somar 1+1 mas com R e C
by Cleber N.Borges 21 Set '15

21 Set '15

Prezados, Gostaria de uma ajuda em como fazer a soma de 1+1 mas usando a interface entre R e C. a chamada da função .Call faz travar e fechar o R. Agradeço por qualquer ajuda nisso Cleber # um CMR (do conceito) abaixo: ################################################ > texto_code <-" + #include <stdlib.h> + #include <R.h> + #include <Rdefines.h> + + void testef( double *k, double *res ) + { + res = k + 1; + } + " > > sink('teste.c') > cat( texto_code ) > sink() > > system("R CMD SHLIB teste.c") cygwin warning: MS-DOS style path detected: C:/R/etc/x64/Makeconf Preferred POSIX equivalent is: /cygdrive/c/R/etc/x64/Makeconf CYGWIN environment variable option "nodosfilewarning" turns off this warning. Consult the user's guide for more details about POSIX paths: http://cygwin.com/cygwin-ug-net/using.html#using-pathnames gcc -m64 -I"C:/R/include" -DNDEBUG -I"d:/RCompile/r-compiling/local/local320/include" -O2 -Wall -std=gnu99 -mtune=core2 -c teste.c -o teste.o gcc -m64 -shared -s -static-libgcc -o teste.dll tmp.def teste.o -Ld:/RCompile/r-compiling/local/local320/lib/x64 -Ld:/RCompile/r-compiling/local/local320/lib -LC:/R/bin/x64 -lR > dyn.load( "teste.dll" ) > is.loaded( 'testef' ) [1] TRUE > > k=1 > ############################ ########## # ( .Call('testef', k=as.double(k), res=as.double(0) ) ) # # aqui o R trava e fecha !!!!!!!!!!!!!!!!!!!!!! --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 1

Ajuda usuários linux
by Fernando Antonio de souza 21 Set '15

21 Set '15

Caros amigos, Sei que a lista é restrita a questões do R, mas como aqui há usuários linux pode ser q alguém me ajude. ja procurei na NET em várias comunidades linux e não consegui nenhuma solução. Recentemente formato meu computador para o sistema Debian jessie com KDE. Após a formatação tudo estava correndo bem. Porém após uma atualização automática (acredito q foi isso) minha tecla CTR e FN não funciona. O comando xev no terminal não retorna nenhum evento. Já reformatei e desta vez instalou gnome e continua a mesma coisa. Possuo um Dell vostro 3550. Alguém sabe como resolver.

4 4

Elipse circundados pontos e cálculo do centro de massa
by ASANTOS 21 Set '15

21 Set '15

Boa tarde Listeiros, Tenho dois conjuntos de pontos A e B (dados$Class) e suas coordenadas (dados$x e dados$y), gostaria de saber se é possível desenhar uma elipse que circunde todos os pontos A e outra elipse que circunde todos os pontos B, posteriormente gostaria de calcular o centro de massa para cada conjunto de pontos A e B, isso é possível? Segue CRM: #Criação do data frame dados<-NULL #Coordenadas x dados$x<-c(0.350,0.487,0.637,0.775,0.825,0.087,0.237,0.400,0.575,0.737,0.062,0.212) #Coordenadas y dados$y<-c(0.025,0.087,0.050,0.025,0.125,0.187,0.150,0.162,0.212,0.237,0.362,0.337) #Classificação dados$Class<-c("A","A","A","A","A","A","B","B","B","B","B","B") #Visualização plot(dados$x,dados$y) text(dados$x,dados$y+0.01,labels=dados$Class,col="red") # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

3 3

OPLS
by marcelo claro de souza 18 Set '15

18 Set '15

Olá pessoal, Preciso rodar umas análises de OPLS (Orthogonal Partial Least Squares) entretanto parece que não existe nenhum pacote no R para rodar esse tipo de análise. Alguém teria um script ou então poderia me indicar algum pacote para o R? Muito obrigado Abraço Marcelo -- *Dr. Marcelo Claro de Souza* Postdoctoral researcher FCFRP - USP Department of Pharmacy Laboratory of Pharmacognosy Ribeirão Preto, São Paulo, Brazil Av. do Café s/n CEP 14040-230

2 1

Link de transmissão ao vivo pelo YouTube do Meetup em R (IME-USP) 18 de setembro
by Luis Benites Sánchez 18 Set '15

18 Set '15

Caros, Link de transmissão ao vivo pelo Youtube do Meetup de R no IME-USP. http://youtu.be/EDXuYbXTl5c Um abraço Luis Benites Sánchez Organizador do Meetup em R -- Luis B. Sánchez Ph.D Candidate Dept. of Statistics-IME, University at São Paulo (USP) São Paulo, BRAZIL lsanchez(a)ime.usp.br Skype: lbenitesanchez Twitter: @lbenitesanchez www.luis-benites.com/

1 0

Escala de um eixo
by Carise Schmidt 17 Set '15

17 Set '15

Boa tarde, Por gentileza, se alguém puder me auxiliar.. Desejo plotar um gráfico de incidência para três localidades ao longo de um período de tempo (anos). Contudo, não estou conseguindo alterar a escala de forma que todos os anos apareçam no gráfico. Da forma como estou plotando os anos aparecem apenas de 2 em 2. Estou usando os dados abaixo: AnoBRSCCC20014,41857818,73841140,917820024,52434718,5791112,655620036,326457 19,4572784,2161320047,37766619,5902581,1693320058,17768923,0116690,396152006 7,84705922,9260271,5683320075,38251819,0436833,8531720085,9479118,53753 82,0774320096,79231719,5792460,8548420106,29338719,5248957,7562320117,539273 22,146469,7660720127,01532323,9061876,69847 E os comandos são os seguintes: >taxa <- read.table("clipboard",header=TRUE,dec=",") >matplot(taxa$Ano,cbind(taxa$BR,taxa$SC,taxa$CC),lty=2,col="black",lwd=1,type="b",pch=c(1,2,4),xlab="Ano", ylab="Incidência (X 100.000 habitantes)",cex.axis=0.5) Carise Schmidt

4 4