Fevereiro 2012 - R-br

intervalo de confiança em modelo de regressão
by Marcelo Cardoso mello 16 Fev '12

16 Fev '12

Gostaria de estimar um ic de um modelo lm. fazendo: predict(model, interval="confidence") Com esta função eu tenho os intervalos em cada valor predito. Mas ... existe uma função que me dá o intervalo para qualquer valor predito? por exemplo: x<-c(12,16,19,18,22) y<-c(40,36,45,49,52) model<-lm(y~x) predict(model, interval="confidence") f<-function(x){coef(model)[1]+coef(model)[2]*x} f(12) Tenho o inetrvalo para este ponto. Como faço para obter o intervalo de f(14) por exemplo? Tem como obter um intervalo médio para todo o modelo, ou seja, uma quantidade para mais e para menor média do modelo?

2 1

[NEWS] Mudanças nas versões do R e melhora de performance
by Fernando Mayer 16 Fev '12

16 Fev '12

Não sei se todos já estão sabendo, mas li isso hoje e achei importante repassar para a lista: Sobre as mudanças de releases do R: R 2.14 to be released on October 31; R 2.13 patch on September 13 http://blog.revolutionanalytics.com/2011/09/r-214-to-be-released-on-october… Brian Ripley on The R Development Process http://blog.revolutionanalytics.com/2011/08/brian-ripley-on-the-r-developme… Sobre a melhora de performance: With byte compiler, R 2.14 will be even faster http://blog.revolutionanalytics.com/2011/08/with-byte-compiler-r-214-will-b… []s, --- Fernando Mayer Dept. Ecologia e Zoologia - ECZ/CCB Universidade Federal de Santa Catarina - UFSC URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com

4 4

Erro: 'Sigma' is not positive definite
by Junior Beleti 16 Fev '12

16 Fev '12

Boa noite a todos. Estou trabalhando com a função mvrnorm, passando como parâmetros mvrnorm(n=1, matriz1[1:5202],matriz2[5202:5202]), sendo as matrizes (matriz1 e matriz2) matrizes hipotéticas, onde estão mostrando suas dimensões. Em experimentos realizados, estou dividindo as matrizes realizando quatro chamas ao invés de uma como a seguir: 1) mvrnorm(n=1, matriz1[1:2601], matriz2[1:2601, 1:2601]) 2) mvrnorm(n=1, matriz1[1:2601], matriz2[2602:5202, 1:2601]) 3) mvrnorm(n=1, matriz1[2602:5202], matriz2[1:2601, 2602:5202]) 4) mvrnorm(n=1, matriz1[2602:5202], matriz2[2602:5202, 2602:5202]) O problema está nas chamadas 2 e 3, nas quais estou recebendo seguinte mensagem : " 'Sigma' is not positive definite " Alguém pode me auxiliar quanto a essa mensagem? Obrigado, Carlos Beleti.

4 4

falha na instalação pacote maps no ubuntu
by Adriano S. Melo 15 Fev '12

15 Fev '12

Prezados colegas, Não estou conseguindo instalar o pacote 'maps'. Uso ubuntu 10.04 e R 2.14.1. Consigo instalar outros pacotes normalmente, mas não o 'maps'. Já tentei baixar os arquivos manualmente (inclusive de versões anteriores), mas sempre tenho erro, aparentemente na compilação. A mensagem de erro segue abaixo. Fico grato por qualquer ajuda. Abraços, Adriano S. Melo Dep. Ecologia - UFG install.packages('maps') Installing package(s) into ‘/home/aa/R/i486-pc-linux-gnu-library/2.14’ (as ‘lib’ is unspecified) --- Please select a CRAN mirror for use in this session --- Carregando a interface Tcl/Tk ... pronto tentando a URL ' http://www.vps.fmvz.usp.br/CRAN/src/contrib/maps_2.2-5.tar.gz' Content type 'application/x-gzip' length 1428846 bytes (1.4 Mb) URL aberta ================================================== downloaded 1.4 Mb * installing *source* package ‘maps’ ... ** package ‘maps’ successfully unpacked and MD5 sums checked ** libs ** arch - gcc -std=gnu99 -g -O2 Gmake.c -o Gmake gcc -std=gnu99 -g -O2 Lmake.c -o Lmake Converting world to world2 f convert.awk < world.line > world2.line /bin/bash: f: comando não encontrado make: [world2.line] Erro 127 (ignorado) make county.L state.L usa.L nz.L world.L world2.L italy.L france.L state.vbm.L state.carto.L make[1]: Entrando no diretório `/tmp/RtmpOHSZmi/R.INSTALL79b85cf8d6ef/maps/src' ./Lmake 0 s b county.line county.linestats ../inst/mapdata/county.L ./Lmake 0 s b state.line state.linestats ../inst/mapdata/state.L ./Lmake 0 s b usa.line usa.linestats ../inst/mapdata/usa.L ./Lmake 0 s b nz.line nz.linestats ../inst/mapdata/nz.L ./Lmake 0 s b world.line world.linestats ../inst/mapdata/world.L ./Lmake 0 s b world2.line world2.linestats ../inst/mapdata/world2.L Cannot read left and right at line 1 make[1]: ** [world2.L] Erro 1 make[1]: Saindo do diretório `/tmp/RtmpOHSZmi/R.INSTALL79b85cf8d6ef/maps/src' make: ** [ldata] Erro 2 ERROR: compilation failed for package ‘maps’ * removing ‘/home/aa/R/i486-pc-linux-gnu-library/2.14/maps’ The downloaded packages are in ‘/tmp/RtmpSDbTyV/downloaded_packages’ Mensagens de aviso perdidas: In install.packages("maps") : installation of package ‘maps’ had non-zero exit status

3 3

R homologação uso comercial
by Vinicius Brito Rocha 15 Fev '12

15 Fev '12

Senhores, estou dando entrada no processo de homologação do R na empresa em que trabalho. Estou sob alguns questionamentos e gostaria da ajuda de vcs. uma delas é: Pré-requisitos de Hardware e Software? alguém por acaso sabe? -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 3

maximin latin hypercube
by Edmar Alino 15 Fev '12

15 Fev '12

Prezados, Tenho um sistema tipo black-box com 10 dimensoes, isto e, y=f(xi), onde i=1,2,..10. E preciso treinar uma rede neural. Para tal gostaria de gerar os pontos internos de maneira organizada e ate otimizada. Parece que a melhor solucao seja o maximin latin hypercube (LHD), no entanto nao sei como gerar os pontos. Pelo site www.spacefillingdesings.nl consegue-se gerar no maximo 100 pontos, considerando 10 variaveis. Existe alguma solucao no R? Alguem pode me ajudar? Saliento que nao obrigatoriamente a geracao dos pontos tem que ser do tipo maximin LHD, embora seja interessante que os pontos sejam bem distribuidos e nao colapsando. Grato desde ja! -- *Edmar Alino (12) 8111-4699*

1 0

Duvida - Calculo do percentual de erro de predicao (percentage prediction error)
by Humberto Hazin 15 Fev '12

15 Fev '12

Prezados, Estou com uma dúvida que pode parecer besteira, porém gostarai da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo glm? Meus dados tem 500 linhas e 4 colunas (anexo). O modelo é: m1<- glm(lgcpue=fy+fFF+fq+fA, family=gaussian, data=teste) Meu objetivo é calcular o percentual de desvio (bias) de um modelo por ano via box plot para saber o quanto ele esta sobre ou sub estimando minhas predições com a variável fFF e sem a variável fFF, por exemplo! todos as variaveis são fatores Agradeço desde já Humberto Hazin

1 0

[Dúvida] Problema na instalação: rgdal
by Pedro Rafael 14 Fev '12

14 Fev '12

Pessoal alguém sabe me informar como corrigir este problema na instalação do pacote rgdal? > install.packages("rgdal") Installing package(s) into ‘/usr/local/lib/R/site-library’ (as ‘lib’ is unspecified) tentando a URL 'http://www.vps.fmvz.usp.br/CRAN/src/contrib/rgdal_0.7-8.tar.gz' Content type 'application/x-gzip' length 1427851 bytes (1.4 Mb) URL aberta ================================================== downloaded 1.4 Mb * installing *source* package ‘rgdal’ ... ** package ‘rgdal’ successfully unpacked and MD5 sums checked gdal-config: gdal-config ./configure: line 2013: gdal-config: command not found Error: gdal-config not found The gdal-config script distributed with GDAL could not be found. If you have not installed the GDAL libraries, you can download the source from http://www.gdal.org/ If you have installed the GDAL libraries, then make sure that gdal-config is in your path. Try typing gdal-config at a shell prompt and see if it runs. If not, use: --configure-args='--with-gdal-config=/usr/local/bin/gdal-config' with appropriate values for your installation. ERROR: configuration failed for package ‘rgdal’ * removing ‘/usr/local/lib/R/site-library/rgdal’ The downloaded packages are in ‘/tmp/RtmpcOSraQ/downloaded_packages’ Mensagens de aviso perdidas: In install.packages("rgdal") : installation of package ‘rgdal’ had non-zero exit status -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

6 12

Legenda do eixo x em um gráfico de controle.
by Juliana Freitas de Mello e Silva 14 Fev '12

14 Fev '12

Boa tarde! Gostaria de saber se é possível mudar os valores do eixo 'x' em um gráfico de controle de qualidade do tipo "xbar.one". Minha intenção seria, ao invés de uma seqüência de 1 a 30, por exemplo, colocar as datas referentes ao valores utilizados. Espero ter sido clara. Obrigada!

2 2

ide no Ubuntu
by rodrigo koblitz 13 Fev '12

13 Fev '12

Buenos, estou começando no R e, acostumado com outros programas, estou tendo alguns problemas com o shell do Linux. Basicamente a seta de voltar casas não funciona, o pior é a seta pra cima e para baixo para locomover entre comandos já digitados. Procurei algo na net, mas não vi nada que resolvesse. Alguém tem alguma sugestão, abs, Koblitz

5 4