Fevereiro 2013 - R-br

Re: [R-br] grafico de função de autocorrelação
by Ivan Bezerra Allaman 22 Fev '13

22 Fev '13

library(tseries) y = log(AirPassengers) acf(diff(y)) library(plotrix) draw.ellipse(0.5,0.2,0.5,0.2) (S,f,P) Allaman

4 3

grafico de função de autocorrelação
by Natalia Martins 22 Fev '13

22 Fev '13

Boa tarde prezados, Fiz alguns gráficos (FAC e FACP) no R e eu preciso fazer uns círculos nas autocorrelações que estao fora do intervalo. Eu salvei como jpg com 100% de qualidade e fiz isso no paint e depois as salvei como pdf, mas a qualidade está horrivel. Alguem saberia me dizer se eu posso fazer esses circulo no R para manter a qualidade? Grata -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

3 3

cluster 5MM
by Vinicius Brito Rocha 22 Fev '13

22 Fev '13

Pessoal, tenho uma base com 5 milhoes de registros e 20 variáveis. preciso clusteriza-la. alguém pode me ajudar. Sera que o clara do pacote cluster resolve? dentre os pacotes bigmemory e biganalytics tem funçoes de cluster? Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* 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 4

Re: [R-br] aggregate
by Ivan Bezerra Allaman 22 Fev '13

22 Fev '13

Não complique o fácil!! da=data.frame(A=rep(LETTERS[1:3],3), B=rep(letters[1:3],3), y=rnorm(9), x=rnorm(9), z=rnorm(9)) aggregate(da[,-c(1:2)], list(A=da[['A']], B=da[['B']]), function(x) c('Média'=mean(x), 'Desvio'=sd(x), 'CV'=sd(x)/mean(x)*100)) (S,f,P) Allaman

2 1

aggregate
by Edson Lira 21 Fev '13

21 Fev '13

Boa tarde a todos, tenho um bd com estas variáveis > names(ccq) [1] "Grupo" "EA" "EB" "EC" "ED" "EE" "EF" "EG" "EH" "EI" "EJ" "PA" "PB" "PC" "PD" "PE" "PF" "PG" "PH" "PI" "PJ" "AA" [23] "AB" "AC" "AD" "AE" "AF" "AG" "AH" "AI" "AJ" "Nome" "Turno" "cargo" "Sexo" "te" "tc" "t_ccq" "esc" "Tipo" Criei uma funçao: eda<-function(x){ média<-mean(x,na.rm=T) desvio<-sd(x,na.rm=T) cv<-desvio/média*100 cbind(média,desvio,cv) } Como faço para obter estas medidas de uma única vez para todas variáveis numérias por exemplo por turno que tem 3 níveis: 1,2,Com. Fiz no aggregate e me dar esta saída(pequena amostra) > ag<-aggregate(. ~ Turno, data = ccq, eda) > ag Turno Grupo.1 Grupo.2 Grupo.3 EA.1 EA.2 EA.3 EB.1 EB.2 EB.3 EC.1 EC.2 EC.3 ED.1 ED.2 ED.3 EE.1 EE.2 EE.3 EF.1 EF.2 EF.3 EG.1 EG.2 EG.3 EH.1 EH.2 EH.3 EI.1 EI.2 1 1 1 0 0 8.0 1.4 17.8 7.83 2.31 29.51 7.8 1.7 21.3 8.48 1.20 14.17 8.3 1.7 20.0 9.13 0.97 10.60 9.00 1.21 13.40 8.43 1.62 19.19 7.6 1.5 2 2 1 0 0 8.2 1.1 13.3 9.44 0.88 9.34 8.9 1.4 15.3 8.78 0.83 9.49 8.8 1.3 14.8 8.78 1.72 19.55 8.67 1.66 19.13 9.56 0.88 9.23 8.6 1.5 3 Com 1 0 0 7.8 1.1 14.2 8.11 1.50 18.52 8.0 1.2 15.4 8.33 1.21 14.51 8.2 1.7 20.3 9.15 0.95 10.37 9.22 0.93 10.12 8.81 1.18 13.36 8.0 1.3 EI.3 EJ.1 EJ.2 EJ.3 PA.1 PA.2 PA.3 PB.1 PB.2 PB.3 PC.1 PC.2 PC.3 PD.1 PD.2 PD Vejam que por exemplo, EA.1 EA.2 EA.3 são respectivamente a média, desvio e cv. Como melhorar esta saída, dando nomes a EA.1, EA.2,EA.3 Talvez um loop? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

amostragem 3D dados espaciais
by Samuel Carvalho 21 Fev '13

21 Fev '13

Boa noite caros membros da r-br Por gentileza, alguém saberia me dizer como amostrar um grid espacial em 3D? Tentei usar a função spsample do pacote sp, porém a mesma me parece que amostra apenas grid's em duas dimensoes. Obg, Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

3 2

função summary das clinicas
by alanarocha＠sapo.pt 21 Fev '13

21 Fev '13

Bom dia, O exemplo do Rodrigo é optimo mas o problema é que ao todo são 37 clinicas... outra duvida eu preciso que as estatisticas do meu banco de dados sejam feitas não clinica a clinica mas sim de ano a ano para todas. o que tenho de modificar na função summary() para aparecer por exemplo mindoismiloito mindoismilnove etc por baixo os valores minimos dos minimos de cada ano. ou então summary normal e depois imprimir na primeira linha os cabeçalhos (nomes) e por coluna o nome da estatistica e a frente os resultados das estatisticas de cada ano. exemplo: doismiloito doismilnove doismildez doismilonze doismildoze min val val val val val max val val val val val eu não me tinha esplicado bem. Cumprimentos Ana Rocha

2 1

PerformanceAnalytics
by Marcos Vagner Macedo 21 Fev '13

21 Fev '13

Boa noite, Alguém com experiencia no pacote PerformanceAnalytics que pode me ajudar com algumas dúvidas? Um abraço, Marcos Macedo

2 1

função summary quase resolvido
by alanarocha＠sapo.pt 20 Fev '13

20 Fev '13

Olá, agora já me aparecem as estatísticas. Mas as clinicas não estão na mesma ordem das que estão no ficheiro. há alguma forma de as coloccar por ordem e em cada linha o nome da clinica e os valores dos 5 anos? setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") head(dados) Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") head(dados) summary(dados) > setwd("G:/SIN/Users/Ana Rua/indicadores_ano") > > dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") Warning message: In read.table("indicadores_08to12.csv", header = TRUE, sep = ";") : incomplete final line found by readTableHeader on 'indicadores_08to12.csv' > head(dados) [1] Clinica..2008..2009..2010..2011..2012..PBarca.0.70.0.72.0.72.0.76.0.77.Braga.0.74.0.75.0.78.0.76.0.73.Fafe.0.74.0.75.0.74.0.72.0.76.Maia.0.70.0.75.0.71.0.76.0.80.VNGaia.0.80.0.83.0.84.0.84.0.81.SMFeira.0.73.0.75.0.74.0.78.0.75.Viseu.0.71.0.68.0.78.0.76.0.76.Guarda.0.69.0.64.0.69.0.74.0.81.Covilhã.0.70.0.71.0.77.0.78.0.79.Coimbra.0.80.0.77.0.77.0.77.0.77.Abrantes.0.81.0.83.0.82.0.80.0.81.Entroncamento.0.68.0.74.0.73.0.75.0.75.Santarém.0.76.0.75.0.73.0.82.0.83.TVedras.0.77.0.69.0.75.0.78.0.81.VFXira.0.74.0.71.0.76.0.79.0.79.Alverca.0.79.0.76.0.80.0.84.0.80.Amadora.0.69.0.70.0.73.0.75.0.75.Lumiar.0.69.0.71.0.72.0.75.0.76.DAD.....0.80.Restelo.0.72.0.71.0.70.0.75.0.74.CUF.0.64.0.66.0.68.0.73.0.71.SAMS.0.82.0.80.0.83.0.85.0.81.APDP.....0.57.Montijo.0.74.0.76.0.71.0.74.0.85.Barreiro.0.77.0.78.0.79.0.77.0.78.Almada.0.75.0.75.0.72.0.76.0.79.Seixal.0.76.0.75.0.77.0.81.0.81.Setúbal.0.71.0.76.0.79.0.79.0.77.Portalegre.0.74.0.78.0.81.0.81.0.76.Évora.0.82.0.81.0.81.0.78.0.81.Grãndola.0.76.0.76.0.7 8.0.78.0.80.Beja.0.71.0.76.0.74.0.73.0.69.Tavira.0.63.0.70.0.70.0.65.0.71.Faro.0.61.0.74.0.79.0.77.0.78.Portimão.0.63.0.75.0.79.0.78.0.82.Funchal..0.68.0.74.0.72.0.74.Machico.....0.73..Clinica..2008....2009 [2] X2010.2011 [3] X2012.PBarca.0.70.0.72.0.72.0.76.0.77.Braga.0.74.0.75.0.78.0.76.0.73.Fafe.0.74.0.75.0.74.0.72.0.76.Maia.0.70.0.75.0.71.0.76.0.80.VNGaia.0.80.0.83.0.84.0.84.0.81.SMFeira.0.73.0.75.0.74.0.78.0.75.Viseu.0.71.0.68.0.78.0.76.0.76.Guarda.0.69.0.64.0.69.0.74.0.81.Covilhã.0.70.0.71.0.77.0.78.0.79.Coimbra.0.80.0.77.0.77.0.77.0.77.Abrantes.0.81.0.83.0.82.0.80.0.81.Entroncamento.0.68.0.74.0.73.0.75.0.75.Santarém.0.76.0.75.0.73.0.82.0.83.TVedras.0.77.0.69.0.75.0.78.0.81.VFXira.0.74.0.71.0.76.0.79.0.79.Alverca.0.79.0.76.0.80.0.84.0.80.Amadora.0.69.0.70.0.73.0.75.0.75.Lumiar.0.69.0.71.0.72.0.75.0.76.DAD.....0.80.Restelo.0.72.0.71.0.70.0.75.0.74.CUF.0.64.0.66.0.68.0.73.0.71.SAMS.0.82.0.80.0.83.0.85.0.81.APDP.....0.57.Montijo.0.74.0.76.0.71.0.74.0.85.Barreiro.0.77.0.78.0.79.0.77.0.78.Almada.0.75.0.75.0.72.0.76.0.79.Seixal.0.76.0.75.0.77.0.81.0.81.Setúbal.0.71.0.76.0.79.0.79.0.77.Portalegre.0.74.0.78.0.81.0.81.0.76.Évora.0.82.0.81.0.81.0.78.0.81.Grãndola.0.76.0.76.0.78.0.78.0.80.Beja.0.71.0.76.0.74.0 .73.0.69.Tavira.0.63.0.70.0.70.0.65.0.71.Faro.0.61.0.74.0.79.0.77.0.78.Portimão.0.63.0.75.0.79.0.78.0.82.Funchal..0.68.0.74.0.72.0.74.Machico.....0.73.. <0 rows> (or 0-length row.names) > > Type 'contributors()' for more information and Error: unexpected string constant in "Type 'contributors()'" > 'citation()' on how to cite R or R packages in publications. Error: unexpected symbol in "'citation()' on" > > Type 'demo()' for some demos, 'help()' for on-line help, or Error: unexpected string constant in "Type 'demo()'" > 'help.start()' for an HTML browser interface to help. Error: unexpected 'for' in "'help.start()' for" > Type 'q()' to quit R. Error: unexpected string constant in "Type 'q()'" > setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") head(dados) summary(dados) > setwd("G:/SIN/Users/Ana Rua/indicadores_ano") > > dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") > head(dados) clinica doismiloito doismilnove doismildez doismilonze doismildoze 1 PBarca 0.70 0.72 0.72 0.76 0.77 2 Braga 0.74 0.75 0.78 0.76 0.73 3 Fafe 0.74 0.75 0.74 0.72 0.76 4 Maia 0.70 0.75 0.71 0.76 0.80 5 VNGaia 0.80 0.83 0.84 0.84 0.81 6 SMFeira 0.73 0.75 0.74 0.78 0.75 > > Type 'contributors()' for more information and Error: unexpected string constant in "Type 'contributors()'" > 'citation()' on how to cite R or R packages in publications. Error: unexpected symbol in "'citation()' on" > > Type 'demo()' for some demos, 'help()' for on-line help, or Error: unexpected string constant in "Type 'demo()'" > 'help.start()' for an HTML browser interface to help. Error: unexpected 'for' in "'help.start()' for" > Type 'q()' to quit R. Error: unexpected string constant in "Type 'q()'" > > setwd("G:/SIN/Users/Ana Rua/indicadores_ano") > > dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") > head(dados) clinica doismiloito doismilnove doismildez doismilonze doismildoze 1 PBarca 0.70 0.72 0.72 0.76 0.77 2 Braga 0.74 0.75 0.78 0.76 0.73 3 Fafe 0.74 0.75 0.74 0.72 0.76 4 Maia 0.70 0.75 0.71 0.76 0.80 5 VNGaia 0.80 0.83 0.84 0.84 0.81 6 SMFeira 0.73 0.75 0.74 0.78 0.75 > summary(dados) clinica doismiloito doismilnove doismildez Abrantes: 1 Min. :0.6100 Min. :0.6400 Min. :0.6800 Almada : 1 1st Qu.:0.7000 1st Qu.:0.7100 1st Qu.:0.7225 Alverca : 1 Median :0.7400 Median :0.7500 Median :0.7550 Amadora : 1 Mean :0.7288 Mean :0.7409 Mean :0.7568 APDP : 1 3rd Qu.:0.7600 3rd Qu.:0.7600 3rd Qu.:0.7900 Barreiro: 1 Max. :0.8200 Max. :0.8300 Max. :0.8400 (Other) :31 NA's :4 NA's :3 NA's :3 doismilonze doismildoze Min. :0.6500 Min. :0.5700 1st Qu.:0.7500 1st Qu.:0.7500 Median :0.7700 Median :0.7800 Mean :0.7697 Mean :0.7708 3rd Qu.:0.7875 3rd Qu.:0.8100 Max. :0.8500 Max. :0.8500 NA's :3 > > setwd("G:/SIN/Users/Ana Rua/indicadores_ano") Error: unexpected '>' in ">" > > Error: unexpected '>' in ">" > > dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") Error: unexpected '>' in ">" > Warning message: Error: unexpected symbol in "Warning message" > In read.table("indicadores_08to12.csv", header = TRUE, sep = ";") : Error: unexpected symbol in "In read.table" > incomplete final line found by readTableHeader on 'indicadores_08to12.csv' Error: unexpected symbol in " incomplete final" > > head(dados) Error: unexpected '>' in ">" > [1] > Clinica..2008..2009..2010..2011..2012..PBarca.0.70.0.72.0.72.0.76.0.77.Braga.0.74.0.75.0.78.0.76.0.73.Fafe.0.74.0.75.0.74.0.72.0.76.Maia.0.70.0.75.0.71.0.76.0.80.VNGaia.0.80.0.83.0.84.0.84.0.81.SMFeira.0.73.0.75.0.74.0.78.0.75.Viseu.0.71.0.68.0.78.0.76.0.76.Guarda.0.69.0.64.0.69.0.74.0.81.Covilhã.0.70.0.71.0.77.0.78.0.79.Coimbra.0.80.0.77.0.77.0.77.0.77.Abrantes.0.81.0.83.0.82.0.80.0.81.Entroncamento.0.68.0.74.0.73.0.75.0.75.Santarém.0.76.0.75.0.73.0.82.0.83.TVedras.0.77.0.69.0.75.0.78.0.81.VFXira.0.74.0.71.0.76.0.79.0.79.Alverca.0.79.0.76.0.80.0.84.0.80.Amadora.0.69.0.70.0.73.0.75.0.75.Lumiar.0.69.0.71.0.72.0.75.0.76.DAD.....0.80.Restelo.0.72.0.71.0.70.0.75.0.74.CUF.0.64.0.66.0.68.0.73.0.71.SAMS.0.82.0.80.0.83.0.85.0.81.APDP.....0.57.Montijo.0.74.0.76.0.71.0.74.0.85.Barreiro.0.77.0.78.0.79.0.77.0.78.Almada.0.75.0.75.0.72.0.76.0.79.Seixal.0.76.0.75.0.77.0.81.0.81.Setúbal.0.71.0.76.0.79.0.79.0.77.Portalegre.0.74.0.78.0.81.0.81.0.76.Évora.0.82.0.81.0.81.0.78.0.81.Grãndola.0.76.0.76.0 > .78.0.78.0.80.Beja.0.71.0.76.0.74.0.73.0.69.Tavira.0.63.0.70.0.70.0.65.0.71.Faro.0.61.0.74.0.79.0.77.0.78.Portimão.0.63.0.75.0.79.0.78.0.82.Funchal..0.68.0.74.0.72.0.74.Machico.....0.73..Clinica..2008....2009 Error: unexpected '[' in "[" > [2] X2010.2011 Error: unexpected '[' in "[" > [3] > X2012.PBarca.0.70.0.72.0.72.0.76.0.77.Braga.0.74.0.75.0.78.0.76.0.73.Fafe.0.74.0.75.0.74.0.72.0.76.Maia.0.70.0.75.0.71.0.76.0.80.VNGaia.0.80.0.83.0.84.0.84.0.81.SMFeira.0.73.0.75.0.74.0.78.0.75.Viseu.0.71.0.68.0.78.0.76.0.76.Guarda.0.69.0.64.0.69.0.74.0.81.Covilhã.0.70.0.71.0.77.0.78.0.79.Coimbra.0.80.0.77.0.77.0.77.0.77.Abrantes.0.81.0.83.0.82.0.80.0.81.Entroncamento.0.68.0.74.0.73.0.75.0.75.Santarém.0.76.0.75.0.73.0.82.0.83.TVedras.0.77.0.69.0.75.0.78.0.81.VFXira.0.74.0.71.0.76.0.79.0.79.Alverca.0.79.0.76.0.80.0.84.0.80.Amadora.0.69.0.70.0.73.0.75.0.75.Lumiar.0.69.0.71.0.72.0.75.0.76.DAD.....0.80.Restelo.0.72.0.71.0.70.0.75.0.74.CUF.0.64.0.66.0.68.0.73.0.71.SAMS.0.82.0.80.0.83.0.85.0.81.APDP.....0.57.Montijo.0.74.0.76.0.71.0.74.0.85.Barreiro.0.77.0.78.0.79.0.77.0.78.Almada.0.75.0.75.0.72.0.76.0.79.Seixal.0.76.0.75.0.77.0.81.0.81.Setúbal.0.71.0.76.0.79.0.79.0.77.Portalegre.0.74.0.78.0.81.0.81.0.76.Évora.0.82.0.81.0.81.0.78.0.81.Grãndola.0.76.0.76.0.78.0.78.0.80.Beja.0.71.0.76.0.74 > .0.73.0.69.Tavira.0.63.0.70.0.70.0.65.0.71.Faro.0.61.0.74.0.79.0.77.0.78.Portimão.0.63.0.75.0.79.0.78.0.82.Funchal..0.68.0.74.0.72.0.74.Machico.....0.73.. Error: unexpected '[' in "[" > <0 rows> (or 0-length row.names) Error: unexpected '<' in "<" > > Error: unexpected '>' in ">" > > Type 'contributors()' for more information and Error: unexpected '>' in ">" > Error: unexpected string constant in "Type 'contributors()'" Error: unexpected symbol in "Error: unexpected string" > > 'citation()' on how to cite R or R packages in publications. Error: unexpected '>' in ">" > Error: unexpected symbol in "'citation()' on" Error: unexpected symbol in "Error: unexpected symbol" > > Error: unexpected '>' in ">" > > Type 'demo()' for some demos, 'help()' for on-line help, or Error: unexpected '>' in ">" > Error: unexpected string constant in "Type 'demo()'" Error: unexpected symbol in "Error: unexpected string" > > 'help.start()' for an HTML browser interface to help. Error: unexpected '>' in ">" > Error: unexpected 'for' in "'help.start()' for" Error: unexpected string constant in "Error: unexpected 'for'" > > Type 'q()' to quit R. Error: unexpected '>' in ">" > Error: unexpected string constant in "Type 'q()'" Error: unexpected symbol in "Error: unexpected string" > > Error: unexpected '>' in ">" setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") head(dados) clinica doismiloito doismilnove doismildez doismilonze doismildoze 1 PBarca 0.70 0.72 0.72 0.76 0.77 2 Braga 0.74 0.75 0.78 0.76 0.73 3 Fafe 0.74 0.75 0.74 0.72 0.76 4 Maia 0.70 0.75 0.71 0.76 0.80 5 VNGaia 0.80 0.83 0.84 0.84 0.81 6 SMFeira 0.73 0.75 0.74 0.78 0.75 summary(dados) clinica doismiloito doismilnove doismildez Abrantes: 1 Min. :0.6100 Min. :0.6400 Min. :0.6800 Almada : 1 1st Qu.:0.7000 1st Qu.:0.7100 1st Qu.:0.7225 Alverca : 1 Median :0.7400 Median :0.7500 Median :0.7550 Amadora : 1 Mean :0.7288 Mean :0.7409 Mean :0.7568 APDP : 1 3rd Qu.:0.7600 3rd Qu.:0.7600 3rd Qu.:0.7900 Barreiro: 1 Max. :0.8200 Max. :0.8300 Max. :0.8400 (Other) :31 NA's :4 NA's :3 NA's :3 doismilonze doismildoze Min. :0.6500 Min. :0.5700 1st Qu.:0.7500 1st Qu.:0.7500 Median :0.7700 Median :0.7800 Mean :0.7697 Mean :0.7708 3rd Qu.:0.7875 3rd Qu.:0.8100 Max. :0.8500 Max. :0.8500 NA's :3 cumprimentos Ana

3 2

Soundex - Marcelino
by Fátima Lima Paula 20 Fev '13

20 Fev '13

Olá Daniel Marcelino. Estou escrevendo meu artigo e preciso saber se existe algum endereço em que se possa obter o Soundex-BR que você criou e me passou há algum tempo. Eu gostaria de citá-lo. Obrigada Fátima

1 0