Outubro 2011 - R-br

[OFF] Morre John McCarthy, pioneiro da Inteligência Artificial e pai do Lisp.
by Pedro Rafael 25 Out '11

25 Out '11

Morre John McCarthy, pioneiro da Inteligência Artificial e pai do Lisp. O matemático e pesquisador da área de ciências da computação John McCarthy, responsável pela criação do termo “inteligência artificial” e inventor da linguaguem de programação Lisp, morreu no último domingo (23) em Palo Alto, na região da Califórnia conhecida como Vale do Silício. http://g1.globo.com/tecnologia/noticia/2011/10/morre-john-mccarthy-pioneiro… -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Leitura arquivos .DB
by Eder David Borges da Silva 25 Out '11

25 Out '11

Bom dia pessoal, Gostaria de pedir ajudar a vocês na leitura de arquivos .DB (database), tentei funções do pacote RODBC, require(RODBC) con0 <- odbcConnectDbase('http://www.leg.ufpr.br/~eder/E2D00007.DB') con1 <- odbcConnectDbase('http://www.leg.ufpr.br/~eder/M00777600m000.DB') porem estou em windows 64 Bit, preciso ler neste sistema operacional, e da msn de erro pelo driver ser 32 bit e outras coisa s mais, desta forma alguém tem uma sugestão para leitura destes arquivos? Alternativa a essa? OBS: os arquivos estão no caminho indicado acima. Grato pela atenção

6 27

order
by Leonard de Assis 25 Out '11

25 Out '11

Boa noite a todos, Estava escrevendo um email gigantesco, detalhando um problema que eu estava enfrentando em uma simulação aqui, mas no meio do mesmo eu descobri um jeito de fazer o que pretendo, então sobrou pra mim o passo 2 da mesma: Conseguir um meio menos "marreta" de fazer. em meu caso, Menos "marreta" implica em algo mais rápido que o que já tenho, pois apesar de ser em um dataframe pequeno, tenho que fazer este passo algumas milhares de vezes em meu algorítmo. então vamos lá: eu preciso adicionar em um dataframe, uma coluna informando a ordem de cada um dos eventos que ocorreram. Eu atualmente faço da seguinte forma: G[order(-G$A,-G$B),]$Pos <- 1:length(G$V) do jeito que está, funciona, minha intenção é algo mais eficiente, se possível. -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

1 0

Regressão polinomial. Uma discussão teórica!
by Ivan Bezerra Allaman 25 Out '11

25 Out '11

Bom tarde senhores!! Me deparei com uma situação que em outrora não havia questionado tal situação, seguindo então o que a literatura indicava. A situação ao qual me refiro, é que ao desdobrar a soma de quadrados de regressão em termos lineares, quadráticos e cúbicos, os três termos foram significativos. FATORES GL SQ QM F value p value Silag:ordered(Tempo):FosForm11.L 1 16.26 16.26 33.97 0.0000 Silag:ordered(Tempo):FosForm11.Q 1 0.22 0.22 0.46 0.5000 Silag:ordered(Tempo):FosForm11.C 1 10.51 10.51 21.96 0.0000 Residuals 108 51.68 0.48 A literatura nos recomenda escolher o termo significativo de maior grau. Sinceramente, ao pararmos para pensar, porque não escolher o termo com a maior SQ, no caso de significância é claro dos dois termos? Afinal, o SQ nos informa o quanto da variação se deve aquele fator não? Esse pensamento, podemos ter no caso de um cálculo a mão certo! Mais hoje em dia, em que temos o valor do "p value", não seria mais prudente escolhermos o termo com o maior "p value"? que consequentemente, é claro, terá o maior SQ? Gostaria de saber a opinião dos colegas a respeito disso. Abraço a todos! (S,f,P) Allaman \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}

6 7

Instalando versão mais nova do R no Ubuntu
by Pedro Rafael 24 Out '11

24 Out '11

Pessoal estou tentando instalar a versão mais recente do R (2.13.2) no Ubuntu 11.10 (64bist) e adicionei o repositório deb http://cran.stat.ucla.edu/bin/linux/ubuntu oneiric/ e em seguida dei um update. O seguinte erro apareceu na hora de atualizar o repositório. *W: Erro GPG: http://cran.stat.ucla.edu oneiric/ Release: As assinaturas a seguir não puderam ser verificadas devido à chave pública não estar disponível: NO_PUBKEY 51716619E084DAB9* Após esse erro tentei fazer como o site do R manda, fui atrás de validar a chave e fiz: *gpg --keyserver keyserver.ubuntu.com --recv-key E084DAB9* Após rodar o comando acima a falha continua. O que pode ser? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 17

ByVal, ByRef
by Leonard de Assis 23 Out '11

23 Out '11

Já existe no R alguma forma de se passar parâmetros uma função explicitando se estes são por valor ou por referência? -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

2 1

Erro ao instalar o pacote qtbase
by Luís Gustavo 23 Out '11

23 Out '11

Olá pessoal, estou tentando instalar o pacote "qtbase", seguindo os comandos no site http://qtinterfaces.r-forge.r-project.org/. Porém, quando chego na etapa de "R CMD INSTALL qtbase", vai instalando tudo ok, e já no final ocorre o erro logo abaixo. Lembrando que estou usando linux em uma maquina virtual. Como posso resolver este problema? Obrigado. Install the project... -- Install configuration: "" -- Up-to-date: /home/luis/qtbase/src/./qtbase.so make[1]: Saindo do diretório `/home/luis/qtbase/src-build' installing to /home/luis/R/i686-pc-linux-gnu-library/2.11/qtbase/libs ** R ** demo ** inst ** preparing package for lazy loading ** help *** installing help indices ** building package indices ... ** testing if installed package can be loaded Error : .onLoad failed in loadNamespace() for 'qtbase', details: call: FUN(c("operator!=", "operator&", "operator*", "operator+", "operator-", error: não foi possível encontrar a função "list2env" ERROR: loading failed * removing ‘/home/luis/R/i686-pc-linux-gnu-library/2.11/qtbase’ -- Luís Gustavo Silva e Silva

2 2

plotar erro padrão usando "arrows"
by Rodrigo Plei 21 Out '11

21 Out '11

Prezados, Como plotar o erro padrão em um gráfico de médias usando a função "arrows"? Agradeço desde logo a atenção. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

3 6

Dúvida em Cluster!
by Elton Gean Araújo Araújo 21 Out '11

21 Out '11

Boa tarde prezados colegas da lista. É a primeira vez que trabalho com analise de cluster e me deparei com um problema. Eu quero agrupar os países segundo os indices de satisfação (Sat), mais quando rodo o gráfico, o eixo x não sai com os nomes dos países, mais sim a identificação só (1,2,3...). Alguém poderia me dizer o que está faltando para aparecer os nomes dos países? Segue o script que estou usando e os dados em anexo. hc <- hclust(dist(model.matrix(~-1 + Sat, dados)) , method= "single") plot(hc,hang=-0.1). Desde já, agradeço qualquer ajuda. Abraços, Elton.

3 3

substitui loop por funcoes do R
by Samuel Carvalho 21 Out '11

21 Out '11

Caros membros do grupo, Alguém poderia por gentileza me auxiliar com o código a seguir? No caso a idéia é substituir o loop que calcula a diferença entre as observações dentro de uma mesma parcela sendo que a primeira observação de cada parcela deve ter valor zerado. Para aqueles mais familiarizados estou calculando a taxa de incremento para sucessivas medições. Segue o código ###RCMR dados<-data.frame(parcela=rep(1:5,each=3), y=runif(15)) dados dados$taxa<-0 for(i in 1:(nrow(dados)-1)) { if(dados$parcela[i] == dados$parcela[i+1]) { dados$taxa[i+1] <- dados$y[i+1]-dados$y[i] } } dados #A seguir tento reproduzir o código acima, porém feito incorretamente dados$taxa<-0 dados$taxa[which(diff(dados$parcela)==0)] <- diff(dados$y) #feito de maneira incorreta dados Abraços Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 2