- R-br - DINF UFPR

[OFF] Linguagem Julia.
by Pedro Rafael 15 Abr '12

15 Abr '12

Linguagem Julia, que já ouviu falar e quem já a usa. Ouvi rumores de pessoas da estatística falando que ela pretende ser um "R" mais otimizado, isso é verdade? Vocês acreditam que ela venha a ser muito popular no meio estatístico? Link: http://julialang.org/ -- Saudações, Pedro Rafael Diniz Marinho.

2 3

Webinar: R for SAS and SPSS Users
by Mauro Sznelwar 15 Abr '12

15 Abr '12

Repassaram para mim, vejam se tem alguma utilidade. Webinar: R for SAS and SPSS Users This free webinar will provide a general introduction to R and also give insight into the differences between R, SPSS and SAS. The session will then offer a brief tutorial on how to create simple Bar Plots in R. The webinar will be sponsored by The Orange County R User Group and Lion Data Systems, LLC. If this webinar's start time does not work for you, please email info(a)liondatasystems.com to let us know a start time that would be better for your country. We can choose different start times for our future R webinars based on audience demand. You can also tweet about the event using a message similar to the following: #Rstats for #SAS and #SPSS users webinar with @LionDataSystems liondatasystems.com/events #in You and your colleagues can register for this webinar at: liondatasystems.com/events Webinar Date: May 9, 2012 Webinar Time: 9 am Pacific Time Webinar Duration: About 1 hour Please email info(a)liondatasystems.com with questions.

1 0

Problema pra criar Função - duvidas com rowSums
by Augusto Ribas 14 Abr '12

14 Abr '12

Ola pessoal, estou com problemas pra deixar uma função generica o suficiente, acho que provavelmente a solução deve ser bem simples, mas não estou conseguindo solucionar alguem poderia me dar uma ajudinha? Minha intenção a partir de uma matriz de network entre hospedeiros e parasitos conseguir calcular quantos parasitos morrem ao um ou mais hospedeiros serem extintos. em baixo eu dou um exemplo de matriz e a função que num funciona direito para quando so 1 hospedeiro morre. #vamos supor uma matrix onde 1 indica que aquele parasita(linha) vive #dentro do hospedeiro tal(coluna) #são 6 parasitos e 6 hospedeiros dados<-matrix(c(1,1,1,0,0,0, 0,0,1,0,0,0, 0,0,1,0,0,0, 0,0,0,1,0,0, 0,0,0,1,1,0, 0,0,0,0,0,1),byrow=T,ncol=6,nrow=6, dimnames=list(paste("P",1:6,sep=""),paste("H",1:6,sep=""))) dados #Separando os hospedeiros hospedeiros<-dimnames(dados)[[2]] hospedeiros #Sorteando uma ordem de extinção para os hospedeiros extinto<-sample(hospedeiros,ncol(dados),replace=F) extinto #função pra ver qts parasitas extinguem dados a morte dos hospedeiros extintos extcas<-function(dados,extintos) { interações<-dimnames(dados)[[1]][which(rowSums(dados[,extintos]>0)>0)] spextintas<-sum(rowSums(dados[interações,]>0)<=length(extintos)) return(spextintas) } #ela funciona mais ou menos #se todos os hospedeiros morrem, todos os parasitos morrem extcas(dados,extinto) #se morre 3 hospedeiros, os 3 primeiros morre nesse caso... extcas(dados,extinto[1:3]) #o mesmo acima pra 2 extcas(dados,extinto[1:2]) #Agora qd tento pra ver so pra 1 hospedeiro não funciona extcas(dados,extinto[1:1]) #o problema ta no rowSums #eu to usando esse comando pra ver qts vezes cada parasito interage com #os hospedeiros rowSums(dados[,extinto]>0) #a questão é até 2 ele funciona rowSums(dados[,extinto[1:2]]>0) #mas no 1 ele nao funciona rowSums(dados[,extinto[1:1]]>0) #pq ai nao fica uma matriz fica um vetor #então tem como eu faze o R entender isso dados[,extinto[1:1]]>0 #como uma matriz de 1 coluna e 6 linhas? Ai acho que o rowSums funcionaria #ou não sei se existe alguma solução mais interessante pra chegar no #mesmo resultado #mas tem que ser uma matrix que nao perca a informação do nome da coluna nem da linha pelo que entendo #minha meta final é fazer um grafico onde no eixo x é o numero de #especies extintas acumulado de hospedeiros e no Y qts parasitas #extinguiram juntos #então se essa função funcionar eu faço um loop pra cada número de #hospedeiros mortos, qts parasitas morrem juntos. #vetor pra receber o resultado resultado<-rep(NA,6) resultado #loop for(i in 2:6) { resultado[i]<-extcas(dados,extinto[1:i]) } resultado #grafico plot(resultado~c(1:6),type="b",xlab="Número de hospedeiros extintos", ylab="Número de parasitos extintos") #eu quero chegar nesse grafico, que fica faltando o resultado pro primeiro #hospedeiro morto. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Dados como fator para numérico
by Sérgio Henrique almeida da silva ju 14 Abr '12

14 Abr '12

Olá Estou trabalhando com alguns dados numéricos e que o R esta considerando como fator. > dput(dados) structure(list(NU_LESOES = c(12, 2, 1, 2, 2, 4, 2, 12, 12, 9, 1, 5, 6, 7, 1, 8, 10, 3, 4, 11)), .Names = "NU_LESOES", row.names = c(NA, -20L), data_types = "C", class = "data.frame") >str(dados2) 'data.frame': 20 obs. of 1 variable: $ NU_LESOES: Factor w/ 12 levels "1","10","15",..: 12 2 1 2 2 4 2 12 12 9 ... - attr(*, "data_types")= chr "C" >table(dados2$NU_LESOES) 1 10 15 2 20 3 30 5 7 8 80 NA 3 4 1 2 1 1 1 1 1 1 1 3 Porém quando uso o comando: dados2$NU_LESOES=as.numeric(dados2$NU_LESOES) os dados mudam > dados2$NU_LESOES=as.numeric(dados2$NU_LESOES) > table(dados2$NU_LESOES) 1 2 3 4 5 6 7 8 9 10 11 12 3 4 1 2 1 1 1 1 1 1 1 3 O que pode está acontecendo? -- 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

3 10

Cálculo de vértices à partir do centróide
by ASANTOS 14 Abr '12

14 Abr '12

Boa tarde pessoal, Bom descobri onde estou errando, quando faço: for(i in 1:(length(xcc))){ for(j in 1:(length(ycc))){ O loop vai pegar a primeira coordenada de xcc e fazer todas as combinações com ycc. Na verdade, eu queria uma ajuda para o loop e fazer com que ele pegue apenas o primeiro valor de xcc e comparar apenas com o primeiro valor de ycc, seguido de xcc[2] com ycc[2] ate xcc[11] com ycc[11], meu script completo é: xcc<-c(672036,672332,672674,672079,671386,671364,670753,671854,669986,670860,671366) ycc<-c(6783650,6784105,6783470,6782354,6782416,6782418,6782199,6783073,6781853,6781791,6781911) res=NULL for(i in 1:(length(xcc))){ for(j in 1:(length(ycc))){ xa<-xcc[i]-20 ya<-ycc[j]+20 v1<-cbind(xa,ya) xb<-xcc[i]+20 yb<-ycc[j]+20 v2<-cbind(xb,yb) xc<-xcc[i]+20 yc<-ycc[j]-20 v3<-cbind(xc,yc) xd<-xcc[i]-20 yd<-ycc[j]-20 v4<-cbind(xd,yd) lim.p<-rbind(v1,v2,v3,v4,v1) write.table(lim.p, file ="parc.mod", col.names=T, row.names=F,append=T) result=read.table("parc.mod",header=F) }} Obrigado -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 2

Para quem eh fa~ de "bleeding-edge tech" (em termos de R)
by Benilton Carvalho 14 Abr '12

14 Abr '12

As promessas sobre as quais comentei em algum ponto no passado (de que o R-core iria incorporar suporte a objetos mais longos que 2^31-1) comecam a ser cumpridas.... Eu imagino que no R-2.16.x eh possivel que jah tenhamos alguma dessa funcionalidade (lembrem-se que o release schedule do R mudou e agora eh apenas 1x por ano)... Enfim, o texto abaixo informa as mudancas que ja' foram feitas no codigo-fonte. b \subsection{LONG VECTORS}{ \itemize{ \item There are the beginnings of support for vectors longer than \eqn{2^{31}-1}{2^31 - 1} elements on 64-bit platforms. This applies to raw, logical, integer, double, complex and character vectors, as well as lists. (Elements of character vectors remain limited to \eqn{2^{31}-1}{2^31 - 1} bytes.) \item What can be done with such vectors is currently very limited, and most operations will return the error \sQuote{long vectors not supported yet}. They can be serialized and unserialized, \code{identical()} and \code{object.size()} work and means can be computed. \item \code{dist()} can produce dissimilarity objects for more than 65536 rows (but for example \code{hclust()} cannot process such objects). \item \code{serialize()} to a raw vector is no longer limited in size on 64-bit platforms. \item The C-level function \code{R_alloc} can now allocate \eqn{2^{34}}{2^34} or more bytes on 64-bit platforms. } }

2 1

Tinn-R com Windows 7 64 bit, R 2.15.0
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Usuários do Tinn-R, Testemunho que praticamente não usei Windows nos últimos 2/3 anos. Nesse tempo estava ajudando o Jakson Aquino no projeto do Vim-R-plugin (que considero hoje uma das melhores alternativas de GUI para usuários intermediários e avançados de computador e R). Pois bem, esses dias decidi portar alguns recursos gráficos de um pacote que roda apenas sob Windows e com o R em 32 bit para todas as plataformas suportadas pelo R. Precisava então estudar o código fonte do pacote no Windows com o R rodando 32 bit. Meu velho Windows Vista 32 bit estava bem pesado, ele veio num notebook 64 bit que havia comprado a alguns anos. Instalei o Windows 7 64 bit, instalei o R (32/64) e o Tinn-R (confesso que estava com saudades...). Estou trabalhando a 2 semanas no projeto do pacote e gostaria de informar que nunca tinha usado esses programas nessa plataforma. Testemunho que o Tinn-R me surpeendeu... muito... ele vem atravessando bravamente a evolução do Windows e do R. Para quem usa Windows, acredito que ele ainda tem muito contribuir, e por muitos anos, pois sua plataforma de desenvolvimento (Object Pascal) foi/é ainda muito consistente. Bom, sugestões de melhorias são sempre bem vindas: Em 2 ou 3 meses deverá sair uma nova versão corrigindo umas coisinhas pequenas que não atrapalham seu uso efetivo. Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Ajuda para evitar um loop com elegância (agora correto)
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Pessoal, Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível: # R help m <- matrix(runif(24), ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(m) # média das colunas mr <- rowMeans (m) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 7

SciELO Books
by Benilton Carvalho 14 Abr '12

14 Abr '12

Uma sugestao de um modo legal de distribuicao de livros (agradecam tbm `a licenca Creative Commons): http://books.scielo.org/ e algo um pouco mais "especifico": http://books.scielo.org/id/p2qh6 benilton

4 3

Ajuda para evitar loop com elegância
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Pessoal, Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível: m <- matrix(1:24, ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(matrixdata) # média das colunas mr <- rowMeans (matrixdata) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 1