- R-br - DINF UFPR

dúvidas em outputs
by alanarocha＠sapo.pt 22 Nov '12

22 Nov '12

Olá, eu tenho um output da função summary de 4 variaveis ou seja tenho 4 valores de medias diferentes medianas diferentes desvios padroes diferentes. eu ando a procura de algo que me esplique por exemplo qual dos 4 vaores de medias é imais sgnificativo, e o mesmo relativamente, a mediana quantis desvio padrão variancia. obrigada desde já. vai aqui meu exemplo para ser mais facil. clínicas X1s09 X1s10 X1s11 X1s12 1 Ponte da Barca 70.5 73.9 75.7 81.0 2 Braga 72.5 79.9 77.3 73.4 3 Fafe 74.7 74.7 71.1 74.4 4 Maia 76.7 72.7 77.3 78.1 5 VNGaia 84.6 84.3 83.4 81.8 6 SMFeira 73.0 73.8 76.7 79.4 7 Viseu 67.4 75.7 75.6 74.5 8 Guarda 70.6 65.4 68.7 78.9 9 CovilhÆ 76.7 76.0 76.4 81.3 10 Coimbra 77.5 77.9 77.0 77.7 11 Abrantes 84.0 81.2 82.9 79.1 12 Entroncamento 74.4 73.5 74.1 72.0 13 Santar‚m 75.2 69.0 80.6 83.3 14 TVedras 70.6 74.1 81.1 78.5 15 VFXira 70.7 76.7 78.3 79.7 16 Alverca 74.3 80.9 84.4 80.2 17 Amadora 69.5 70.5 74.7 75.5 18 Lumiar 71.4 72.3 73.2 76.9 19 Restelo 71.2 69.7 72.5 74.0 20 CUF 64.2 65.5 71.6 66.5 21 SAMS 78.8 82.9 84.0 82.3 22 APDP 0.0 0.0 0.0 57.6 23 Montijo 77.3 71.4 69.8 84.3 24 Barreiro 78.6 78.6 77.8 77.6 25 Almada 76.5 73.1 74.8 78.2 26 Seixal 78.6 75.4 80.9 80.0 27 Set£bal 76.3 79.8 78.2 74.8 28 Portalegre 77.0 80.0 80.1 77.7 29 vora 81.2 79.6 80.9 81.8 30 Grƒndola 75.9 78.4 80.2 78.8 31 Beja 73.9 76.4 75.0 67.8 32 Tavira 73.0 66.5 66.9 75.1 33 Faro 70.6 78.5 79.1 77.6 34 PortimÆo 71.9 80.2 79.0 80.6 dim(dados) [1] 34 5 #estatísticas descritivas sapply(dados[,2:4],mean) média: X1s09 X1s10 X1s11 X1s121s09 72.332 73.19 1 74.685 74.685 sapply(dados[,2:5],median) medianas: X1s09 X1s10 X1s11 X1s12 74.35 75.55 77.15 78.15 sapply(dados[,2:5],min) mínimo: X1s09 X1s10 X1s11 X1s12 0.0 0.0 0.0 57.6 sapply(dados[,2:5],max) máximo: X1s09 X1s10 X1s11 X1s12 84.6 84.3 84.4 84.3 sapply(dados[,2:5],sd) desvio padrão: X1s09 X1s10 X1s11 X1s12 13.501958 13.796999 13.890047 5.241185 sapply(dados[,2:5],var) variância: X1s09 X1s10 X1s11 X1s12 182.30286 190.35719 192.93341 27.47002 sapply(dados[,2:5],quantile) quantis: X1s09 X1s10 X1s11 X1s12 0% 0.000 0.00 0.000 57.600 25% 70.825 72.40 74.250 74.875 50% 74.350 75.55 77.150 78.150 75% 76.925 79.35 80.175 80.150 100% 84.600 84.30 84.400 84.300 diff(range(dados$X1s09, na.rm=T)) # amplitude amplitude total X1s09 X1s10 X1s11 X1s12 84.6 diff(range(dados$X1s10, na.rm=T)) # amplitude [1] 84.3 diff(range(dados$X1s11, na.rm=T)) # amplitude [1] 84.4 diff(range(dados$X1s12, na.rm=T)) # amplitude [1] 26.7 Ana Rocha

1 0

Criando Rotina
by Giovane Pretto Azambuja 22 Nov '12

22 Nov '12

Olá, pessoal Sou um pouco novato em R, mas vamos tentar. Tenho de fazer a análise estatística, pela função summary(), de uma pancada de objetos. Gostaria de saber se existe uma maneira em que o R faça essas análises, e as armazene em um objeto diferente para cada análise, de maneira automática. Tipo uma rotina. Creio que consigo fazer isso com o comando for, porém como são vários objetos não me veio em mente como organiza-lo. OBS: Não precisa formular a rotina inteira, apenas dê uma luz. Agradeço quem puder ajudar Giovane P. Azambuja Enviado do Email do Windows

3 2

Enc: Problemas com rgl!
by Ivan Bezerra Allaman 22 Nov '12

22 Nov '12

Boa tarde Daniel! Tenho a seguinte mensagem como camando: Lendo listas de pacotes... Pronto Construindo árvore de dependências Lendo informação de estado... Pronto E: Você deve colocar algumas URIs 'source' em seu sources.list O que isso quer dizer eu não sei, só sei que tenho a seguinte conteúdo no meu sources.list: deb http://mint-packages.c3sl.ufpr.br/ maya main upstream import deb http://archive.ubuntu.com/ubuntu/ precise main restricted universe multiverse deb http://archive.ubuntu.com/ubuntu/ precise-updates main restricted universe multiverse deb http://security.ubuntu.com/ubuntu/ precise-security main restricted universe multiverse deb http://archive.canonical.com/ubuntu/ precise partner deb http://packages.medibuntu.org/ precise free non-free # deb http://archive.getdeb.net/ubuntu precise-getdeb apps # deb http://archive.getdeb.net/ubuntu precise-getdeb games #Repositório do R deb http://cran-r.c3sl.ufpr.br/bin/linux/ubuntu precise/ deb http://cran-r.c3sl.ufpr.br/ precise-backports main restricted universe Inté, \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-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature} ________________________________ De: Daniel C Bezerra <danielcbezerra(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Enviadas: Sábado, 8 de Setembro de 2012 15:57 Assunto: Re: [R-br] Problemas com rgl! Tente rodar no terminal: sudo apt-get build-dep r-cran-rgl On Sat, Sep 8, 2012 at 3:30 PM, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> wrote: Boa tarde senhores! > > >Ao rodar um gráfico em 3D com a rgl percebi que estava muito lento, travando. Então descobri que não estava instalado o driver de vídeo da minha máquina. Ao instalar e tentar usar a rgl obtive erros. Então removi o pacote que estava dentro do usr blablab e tentei reinstalar o pacote sem sucesso com o seguinte erro: > > >* installing *source* package ‘rgl’ ... >** package ‘rgl’ successfully unpacked and MD5 sums checked >checking for gcc... gcc -std=gnu99 >checking whether the C compiler works... yes >checking for C compiler default output file name... a.out >checking for suffix of executables... >checking whether we are cross compiling... no >checking for suffix of object files... o >checking whether we are using the GNU C compiler... yes >checking whether gcc -std=gnu99 accepts -g... yes >checking for gcc -std=gnu99 option to accept ISO C89... none needed >checking how to run the C preprocessor... gcc -std=gnu99 -E >checking for gcc... (cached) gcc -std=gnu99 >checking whether we are using the GNU C compiler... (cached) yes >checking whether gcc -std=gnu99 accepts -g... (cached) yes >checking for gcc -std=gnu99 option to accept ISO C89... (cached) none needed >checking for libpng-config... yes >configure: using libpng-config >configure: using libpng dynamic linkage >checking for X... libraries , headers >checking GL/gl.h usability... yes >checking GL/gl.h presence... yes >checking for GL/gl.h... yes >checking GL/glu.h usability... yes >checking GL/glu.h presence... yes >checking for GL/glu.h... yes >checking for glEnd in -lGL... yes >checking for gluProject in -lGLU... yes >checking for freetype-config... yes >configure: using Freetype and FTGL >configure: creating ./config.status >config.status: creating src/Makevars > >. >. >. > >Error : .onLoad failed in loadNamespace() for 'rgl', details: > call: dyn.load(file, DLLpath = DLLpath, ...) > error: impossível carregar objeto compartilhado '/usr/local/lib/R/site-library/rgl/libs/rgl.so': > libGL.so.1: cannot open shared object file: No such file or directory >Erro: loading failed >Execução interrompida >ERROR: loading failed >* removing ‘/usr/local/lib/R/site-library/rgl’ > >The downloaded source packages are in > ‘/tmp/RtmpS6wZHV/downloaded_packages’ >Mensagens de aviso perdidas: >In install.packages("rgl") : > installation of package ‘rgl’ had non-zero exit status > > >Já pesquisei em vários lugares, instalar várias coisas e o erro persiste. Grato se alguém puder ajudar. Utilizo o mint 13 maya ubuntu. > > >Abraços! > > > >\begin{signatu ><<>>= >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-5596 >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. >

3 3

Regressão linear
by Narede Golaco 22 Nov '12

22 Nov '12

Amigos gostaria de saber como se analisa uma regressão linear apartir do test T ? Agradeço a ajuda e qualquer exemplo, assim como a indicação de material, pois ja coloquei no google " teste T e regressão linear '' e não aparece nada de relevante que avalie a regressão com o teste T. Agradeço desde ja a ajuda.

4 6

Comparar graficamente densidades estimadas
by Matheus Monteiro 21 Nov '12

21 Nov '12

Obtive as densidades aproximadas da 't' e da 'normal', gostaria de compara-las através de graficos dessas densidades estimadas. Qual seria o melhor comando para isso? Uso o que está entre parênteses ou o outro? Obrigado. *t * > fitdistr(x[,1],"t") m s df -0.001874414 0.009122738 3.522064232 ( 0.001726835) ( 0.001862692) ( 2.234067399)* **normal* fitdistr(x[,1],"normal") mean sd 0.0001040816 0.0128201843 (0.0018314549) (0.0012950342) x são dados simulados.

2 4

duvida
by Hélio Gallo Rocha 21 Nov '12

21 Nov '12

Boa noite, Espero não estar enviando este e-mail para todos, mas sim para algum admin da lista. Não estou conseguindo configurar para receber todas as msg separadamente e não resumos. Grato pela atenção Hélio

2 1

Re: [R-br] INLA Geoestatística
by Elias T. Krainski 21 Nov '12

21 Nov '12

Oi Taynãna, DIC negativo pode ocorrer com distribuições contínuas sem problemas. O mesmo vale para log-Verossimilhança. Veja sobre DIC em http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/dicpage.shtml ### exemplo de log-verossimilhança negativa y1 <- rnorm(100, 0, 100) ### log-verossimilhança em 0, 100 sum(dnorm(y1, 0, 100, log=TRUE)) ### exemplo de log-verossimilhança positiva y2 <- rnorm(100, 0, 0.01) ### log-verossimilhança em 0, 0.01 sum(dnorm(y2, 0, 0.01, log=TRUE)) ### para distribuições discretas, o ### log-verossimilhança será sempre negativa y3 <- rpois(100, 3) ### log lik sum(dpois(y3, 3, log=TRUE)) att. Elias T. Krainski > >Ei, Elias. >Agora sim, consegui rodar, mas o DIC está dando negativo, com qualquer das estratégias utilizadas (gaussiana, laplace ou laplace simplificada)... >Percebeu isso? > >Abraço, Taynãna > > >Taynãna, fiz > inla.upgrade(testing=TRUE) >atualizou para a versao de 5 de novembro >e consegui executar todo o script SwissRainfall.R sem problemas... > > > > Att. >Elias T. Krainski > >>Oi, Elias. >> >> >>O Havard Rue respondeu dizendo p tentar o comando inla.upgrade(testing=TRUE), mas não deu certo. >>A versão que estou usando é de 5 de novembro. >> >> >>Abraço, Taynãna >> >>Oi Taynãna, >> >>qual a versão do seu INLA. >> inla.version() >>Eu consegui executar sem problemas na versão de 18 de agosto >> >> >> >>Elias T. Krainski >> >>> >>> >>>Pessoal, >>> >>>Estou tentando reproduzir um exemplo de modelagem espacial no INLA com dados de geoestatística, do artigo >>> >>>de Blangiardo et al (2012). O código em R está disponível e o exemplo 3 (SwissRainfall.R). >>> >>>Ao rodar o modelo aparece a seguinte mensagem: >>> >>>> mod = inla(formula, >>>+ data=inla.stack.data(stack, spde=spde),family="gaussian", >>>+ control.predictor=list(A=inla.stack.A(stack), compute=TRUE), >>>+ >>control.compute=list(cpo=TRUE, dic=TRUE), >>>+ keep=FALSE, verbose=TRUE) >>>Erro em .local(x, i, j, ..., value) : >>> not-yet-implemented 'Matrix[<-' method >>>Calls: inla -> inla.predictor.section -> [<- -> [<- -> .local >>> >>>Alguém tem idéia do quê possa estar acontecendo? >>> >>>Att, >>>Taynãna >>>Pesquisadora em Saúde Pública - ENSP/FIOCRUZ >

2 1

Geoestatística: pacotes para o R [iniciantes]
by Ligia Souza 21 Nov '12

21 Nov '12

Boa tarde, Estou procurando, novamente, indicações. Desta vez de pacotes do R para usar geoestatística, de preferência para iniciantes hehehe Abraços, -- *LIGIA PEREIRA DE SOUZA* *Instituto Brasileiro para Educação e Negócios - IBEN * * * *Graduanda em Ecologia* * * *UNESP - Rio Claro/SP* *Lattes - http://lattes.cnpq.br/3635101224851131<https://wwws.cnpq.br/cvlattesweb/PKG_MENU.menu?f_cod=729600D125BF53DC3FC4A0…> * * *

3 2

problemas da correlação resolvidos
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

muito obrigado Daniel Ana

2 1

matriz de correlação
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

Olá Daniel, auentão quer dizer q vou ter de fazer varios cor.test pra relacionar X1s09 primeiro com X1s10 depois com Xs11 e depois com X1s12 e assim sucessivamente? Ana

2 1