- R-br - DINF UFPR

encontrar resposta no tempo anterior
by Maurício Lordêlo 27 Fev '13

27 Fev '13

Saudações a todos desta lista. Estou precisando de uma ajuda no sentido de definir uma nova variável em função de outra já existente. No exemplo ilustrativo, tenho uma variável resposta binária (denominada de "resp") observada em quatro indivíduos (id) e medida em quatro momentos distintos (time). Preciso desta resposta no tempo anterior (inicialmente de primeira e segunda ordem). A resposta prévia de primeira ordem (resp_ant1) indica o valor da resposta daquele indivíduo no tempo imediatamente anterior. De forma análoga, tem-se a resposta prévia de segunda ordem (resp_ant2). Neste exemplo com poucas observações fica "fácil" visualizar estas respostas prévias e criar estas novas variáveis. Alguém teria uma sugestão de uma forma mais simples para encontrar as variáveis "resp_ant1" e "resp_ant2"? Desde já agradeço, Maurício id=rep(1:4, each = 4) time=rep(1:4, 4) resp=c(1,1,0,1,0,0,0,1,0,1,1,1,0,1,0,1) resp_ant1=c(NA,1,1,0,NA,0,0,0,NA,0,1,1,NA,0,1,0) resp_ant2=c(NA,NA,1,1,NA,NA,0,0,NA,NA,0,1,NA,NA,0,1)

2 3

Vaga para Professor Visitante UFMG
by Carolina Pena 27 Fev '13

27 Fev '13

Bom dia, Gostaria de divulgar vaga de Professor Visitante na UFMG. As inscrições já estão abertas. Perfil do candidato: Graduação em Turismo, Estatística,Economia e Administração e Doutorado em Turismo, Estatística, Economia ou áreas afins (Área 1) https://www2.ufmg.br/ concursos/Concursos/Professor- Substituto-Visitante/ Editais352 Abraços, Carolina Pena

1 0

Dados de solos do Brasil
by Alisson Lucrecio 26 Fev '13

26 Fev '13

Caros Colegas R-br, Boa tarde. Eu estou precisando de dados de solos do Brasil. Alguém sabe onde posso encontrar um lista ou planilha com esses dados? Obrigado. Alisson Lucrécio da Costa

1 0

Pesquisar elementos dentro de uma variável
by Victor Eduardo 26 Fev '13

26 Fev '13

Pessoal, tem como eu pesquisar dentro de uma variável no R apenas algumas siglas ou números? Tenho um banco de dados com mais de 5000 observações e estou interessado apenas na região nordeste, ou seja, preciso fazer uma busca no meu banco de dados para que pegue somente essa região. A variável está dessa maneira: [1] "EBA2927400001" "EDF5300100001" "EES3205300001" "EMS5002700001" "EMS5002700002" [6] "EPI2211000001" "EPR4106900002" "ERJ3304551111" "ERS" "ERS4314900001" [11] "ESC4205400001" "ESP3550309999" "ETO1721000001" "MAC1200200001" "MAC1200200002" [16] "MAC1200300001" "MAC1200340001" "MAC1200400001" "MAC1200430001" "MAL2700100001" No caso, pegaria BA, AC, AL e os demais estados do nordeste. Como posso fazer uma busca dentro dessa variável para me retornar apenas as siglas que quero? Abraços e ótima semana!

4 4

formato do output
by alanarocha＠sapo.pt 26 Fev '13

26 Fev '13

bom dia, o meu quadro das estatísticas do output esta de uma forma estranha: > # Estatística por ano > porano<-apply(dados[,-1],2,function(x)c(summary(x),0)[1:7]) > desviopadrão<-sd(dados[,-1], na.rm= TRUE) Warning message: sd(<data.frame>) is deprecated. Use sapply(*, sd) instead. > variância<-var(dados[,-1], na.rm= TRUE) > amplitude<-apply(dados[,-1],2,range, na.rm= TRUE) > rownames(porano) <- c('Minimo','Primeiro > quartil','Mediana','Média','Terceiro Quartil','Máximo','NAs') > porano<-rbind(porano,'Desviopadrão','Variância','Amplitude') > colnames(porano) = 2008:2012 > desviopadrão X2008 X2009 X2010 X2011 X2012 5.492936 4.433807 4.226329 3.919534 4.968671 > variância X2008 X2009 X2010 X2011 X2012 X2008 30.172348 15.315341 12.590909 13.140152 8.883523 X2009 15.315341 19.079545 13.514205 9.153409 5.238636 X2010 12.590909 13.514205 18.329545 11.909091 5.261364 X2011 13.140152 9.153409 11.909091 15.047348 8.272727 X2012 8.883523 5.238636 5.261364 8.272727 13.528409 > amplitude X2008 X2009 X2010 X2011 X2012 [1,] 61 64 68 65 57 [2,] 82 83 84 85 85 > porano 2008 2009 2010 2011 Minimo "61" "64" "68" "65" Primeiro quartil "70" "71" "72.25" "75" Mediana "74" "75" "75.5" "77" Média "72.88" "74.09" "75.68" "76.97" Terceiro Quartil "76" "76" "79" "78.75" Máximo "82" "83" "84" "85" NAs "4" "3" "3" "3" "Desviopadrão" "Desviopadrão" "Desviopadrão" "Desviopadrão" "Variância" "Variância" "Variância" "Variância" "Amplitude" "Amplitude" "Amplitude" "Amplitude" 2012 Minimo "57" Primeiro quartil "75" Mediana "78" Média "77.08" Terceiro Quartil "81" Máximo "85" NAs "0" "Desviopadrão" "Variância" "Amplitude" > > eu quero da forma: 2008 2009 2010 2011 2012 name.est valor de cada ano sem aspas exemplo: mínimo valor2008 valor2009 valor2010 valor2011 valor2012 desviopadrão valor2008 valor2009 valor2010 valor2011 valor2012 o que poderá estar a acontecer? Cumprimentos Ana Rocha

1 0

Resolvendo uma distribuição Uniforme SIMPLES.
by Jackeline Bonetti Campos 25 Fev '13

25 Fev '13

Boa noite! Preciso resolver uma distribuição Uniforme informando sua densidade e sua E(x) porém sem utilizar os comandos prontos do R. como por ex dunif. No caso, vou ter que dar valores para a e b (limite inferior e superior) a=2 b=3 , por ex. E terei que colocar um comando para que o usuário insira o valor de x. sendo que a fdp ou f(x)=1(sqrt(b-a)) para a<=x<=b e f(x)=0 , caso contrário. Ou seja, pra x<=a || x>=b e a E(x)= b+a(sqrt(2)) Portanto tenho que fazer um comando para descobrir a f(x) e a E(x) utilizando somente comandos como while, if, else,for, etc. Alguém consegue me ajudar? Estou desesperada já, rs. Att, Jackeline.

1 0

Re: [R-br] Janela de gráficos do R não mostro o gráfico (Ubuntu 12.04)
by Andre Oliveira 25 Fev '13

25 Fev '13

O meu é. R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) Aqui o problema continua! André Oliveira Souza ________________________________ De: David Valentim Dias <dvdscripter(a)gmail.com> Para: Andre Oliveira <andreolsouza(a)yahoo.com.br> Enviadas: Segunda-feira, 25 de Fevereiro de 2013 15:53 Assunto: Re: [R-br] Janela de gráficos do R não mostro o gráfico (Ubuntu 12.04) Andre, Aqui o problema já desapareceu. R version 2.15.2 (2012-10-26) Platform: x86_64-unknown-linux-gnu (64-bit) Abraços Em 25 de fevereiro de 2013 11:07, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: Não sei o motivo, mas isto apareceu de uma hora para outra no Ubuntu 12.04. > >André Oliveira Souza > > >________________________________ > De: David Valentim Dias <dvdscripter(a)gmail.com> >Para: r-br(a)listas.c3sl.ufpr.br; Andre Oliveira <andreolsouza(a)yahoo.com.br> >Enviadas: Domingo, 24 de Fevereiro de 2013 15:05 >Assunto: Re: [R-br] Janela de gráficos do R não mostro o gráfico (Ubuntu 12.04) > > > >Notei o mesmo problema no Gnome shell, parece que as decorações usando composition de ambos os sistemas deixam o R sem saber que a device esta pronta para plotar os dados. Isso não ocorre no Rstudio que é baseado no Qt. > > > > >Em 24 de fevereiro de 2013 14:35, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: > >Boa tarde, >>estou com um pequeno problema que é o seguinte. >> >> >> >> >>x=rnorm(9999) >>hist(x) >> >> >> >> >>Este gráfico não aparece na janela do R ele fica em branco só consigo ver o gráfico se maximizar a janela. Alguém tem alguma sugestão de resolver este problema? >> >> >>obrigado >> >>André Oliveira Souza >>_______________________________________________ >>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. >> > > >-- > >Biostatistician Consultant and Cogito Programmer's http://www.facebook.com/CogitoConsultoriaEstatistica >Curriculum: http://lattes.cnpq.br/7541377569511492 > > -- Biostatistician Consultant and Cogito Programmer's http://www.facebook.com/CogitoConsultoriaEstatistica Curriculum: http://lattes.cnpq.br/7541377569511492

1 0

ordenar linhas, separar elementos e contar elementos
by Manoel Nascimento 25 Fev '13

25 Fev '13

Olá! Estou precisando fazer uma operação que deve ser supersimples, mas que como bom novato em R não sei fazer sem me embananar todo. Tenho uma tabela com as seguintes colunas: Municípios População..2010. Ranking.populacional PIB.2006..R..1.000. PIB.2007..R..1.000. PIB.2008..R..1.000. PIB.2009..R..1.000. PIB.2010..R..1.000. Ranking.PIB PIB.per.capita..R.. Ranking.PIB.per.capita Ela mostra os municípios brasileiros, sua população em 2010, o ranking populacional de cada município, seu PIB entre 2006 e 2010, seu ranking de PIB, o PIB per capita e seu ranking de acordo com o PIB per capita. As operações que preciso fazer são as seguintes: a) Ordenar as cidades em ordem decrescente, do maior para o menor PIB; b) Separar as mil cidades com maior PIB, e as mil cidades com menor PIB; c) Em cada um destes dois subconjuntos (cidades com maior PIB e cidades com menor PIB), contar quantas cidades há com mais e com menos de 100 mil habitantes. Como faço isto? []'s Manolo

2 7

Janela de gráficos do R não mostro o gráfico (Ubuntu 12.04)
by Andre Oliveira 25 Fev '13

25 Fev '13

Boa tarde, estou com um pequeno problema que é o seguinte. x=rnorm(9999) hist(x) Este gráfico não aparece na janela do R ele fica em branco só consigo ver o gráfico se maximizar a janela. Alguém tem alguma sugestão de resolver este problema? obrigado André Oliveira Souza

2 2

números em percentagem no R
by alanarocha＠sapo.pt 25 Fev '13

25 Fev '13

Bom dia, alguem me pode dizer como indicar ao R que os valores do banco de dados são com pontos como casa decimal e percentagem e os cabeçalhos tem de ser 2008 2009 2010 etc porque o que tem de aparecer no grafico final é essa formatação. vou mostrar o codigo que fiz: clinica 2008 2009 2010 2011 2012 1 PBarca 70% 72% 72% 76% 77% 2 Braga 74% 75% 78% 76% 73% 3 Fafe 74% 75% 74% 72% 76% 4 Maia 70% 75% 71% 76% 80% 5 VNGaia 80% 83% 84% 84% 81% 6 SMFeira 73% 75% 74% 78% 75% setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep="%") head(dados) # Estatística por ano perc<-100*prop.table(table(dados)) porano<-apply(dados[,-1],2,function(x) c(summary(x),0)dados) rownames(porano) <- c(perc,'Min.','1st Qu.','Median','Mean','3rd Qu.','Max.','NAs') porano perc boxplot(dados[,-1]) preciso que apareça no no summary() e gráfico a formatação que esta nos dados acima Como poderei fazer por causa dos cabeçalhos serem numeros? Cumprimentos Ana

1 0

Cálculo de probabilidade em variáveis contínuas
by Jackeline Bonetti Campos 25 Fev '13

25 Fev '13

Olá! Sou novata no R. e ficarei muito grata se alguém puder ajudar. Quero implementar um programa que realize cálculos de probabilidade e momentos em variáveis aleatorias discretas (P(X <x); P(X = x); V AR[X]). Seria dado opções de funções ao usuário e este daria os valores de x. E o programa usaria de uma interface (um menu) com as opções... Porém, a única coisa que já consegui fazer foi o menu. Os outros comandos estão complicados. Fico no aguardo. Att, Jackeline.

2 2

Redes de cocitação
by Roney Fraga 24 Fev '13

24 Fev '13

Caros, Pretendo fazer uma rede de citação e co-citação, a mesma terá mais de 150000 nós ou (vértices, indivíduos, etc.). Estou buscando as informações no Web of Science. Para exemplificar fiz uma busca de 50 artigos, e exportei para o seguinte endereço: https://docs.google.com/document/d/14w8lQgAwZsmwZ-xGbV5vy4Z3lzvIkFa6Bj44GHi… minha dificuldade é deixar esses dados da seguinte forma, https://docs.google.com/spreadsheet/ccc?key=0AnxuqxFeBpZodG9waVNPYmdNV0V5ZW… uma matriz quadrada relacionando todos os atores. Atualmente utilizo o software VantagePoint para importar esses dados e criar uma matriz de adjacência. Considerando as várias dificuldades de uso do VantagePoint, preço, velocidade de processamento, acesso, etc, busco uma alternativa para importar os dados do web of science e criar uma matriz de adjacência pelo R, ou outro software acessível. primeira pergunta: existe alguma função do R para importar esses dados (exportados do web of science) e criar uma matriz de adjacência? segunda pergunta: é possível criar uma rede entre 150 e 250 mil nós usando o R? terceira pergunta: qual pacote (igraph, sna, statnet) lida melhor com redes grandes? Att -- Roney

3 3

Tópicos
by Jackeline Bonetti Campos 24 Fev '13

24 Fev '13

1 0

GRANDES RELATÓRIOS NO R
by Tiago Souza Marçal 23 Fev '13

23 Fev '13

Bom dia pessoal, eu fiz uma analise no r que gerou um relatório muito grande e eu não estou conseguindo todas as analises. Como faço para ver todas as análises? Att. Tiago.

2 2

média de tempo de grupos
by Fátima Lima Paula 22 Fev '13

22 Fev '13

Prezados, tenho um banco com causas de readmissões e tempos entre as internações. Queria criar uma coluna com o tempo médio de cada grupo de causas. Alguém pode me ajudar? Obrigada. nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano","Paula","Clea") tempo=c(68,3,14,89,7,50,30,5) causa=c("circ","resp","osteo","resp","circ","circ","osteo","osteo") banco1=cbind(nome,tempo,causa) banco1=data.frame(banco1) banco1

3 3

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 21 Fev '13

21 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 20 Fev '13

20 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

Step GLM - remover valores ausentes?
by Elisa Seyboth 20 Fev '13

20 Fev '13

Olá! Estou tentando fazer uma seleção de GLMs através da função step(). Porém, ocorre o seguinte erro: Erro em step(M1) : Número de linhas em uso mudou: remover valores ausentes? Acho que ele poderia estar relacionado aos NAs presentes na tabela, mas eles já haviam sido eliminados da tabela a ser utilizada no primeiro modelo (M1) com a função na.omit(). Teriam alguma sugestão? Desde já agradeço, Elisa. Oc. Elisa Seyboth Mestranda em Oceanografia Biológica Laboratório de Tartarugas e Mamíferos Marinhos Universidade Federal do Rio Grande - FURG Tel.: (53) 32336537 / 81093713

2 2

função sumary do R
by alanarocha＠sapo.pt 20 Fev '13

20 Fev '13

Boa tarde, eu estou a trabalhar com as estatisticas descritivas no R. a leitura do ficheiro não deu erros mas a função sumary não faz nada. código: setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") head(dados) Clinica X2008 X2009 X2010 X2011 X2012 1 PBarca 70% 72% 72% 76% 77% 2 Braga 74% 75% 78% 76% 73% 3 Fafe 74% 75% 74% 72% 76% 4 Maia 70% 75% 71% 76% 80% 5 VNGaia 80% 83% 84% 84% 81% 6 SMFeira 73% 75% 74% 78% 75% summary(dados) Clinica X2008 X2009 X2010 X2011 X2012 Abrantes: 1 74% : 5 75% : 8 74% : 4 78% : 6 81% : 7 Almada : 1 : 4 76% : 5 79% : 4 76% : 5 76% : 4 Alverca : 1 69% : 3 71% : 4 : 3 75% : 4 80% : 4 Amadora : 1 70% : 3 : 3 72% : 3 : 3 75% : 3 APDP : 1 71% : 3 68% : 2 73% : 3 77% : 3 77% : 3 Barreiro: 1 76% : 3 70% : 2 77% : 3 72% : 2 79% : 3 (Other) :31 (Other):16 (Other):13 (Other):17 (Other):14 (Other):13 porque aparece o X nas no nome das variaveis? e não percebo porque o sumary apresenta isto. Cumprimentos Ana Rocha

4 3

Testando frequências de distribuição normal usando o pacote MASS
by ASANTOS 19 Fev '13

19 Fev '13

Bom tarde Pessoal, Estou tentando testar frequências de distribuição normal para um vetor com 100,000 valores com o pacote MASS, o problema é que quando faço o teste de Qui quadrado entre a distribuição empírica o teste não funciona: Erro em chisq.test(x = table(emp), y = table(est)) : 'x' e 'y' devem ter ao menos 2 níveis Quando faço para binomial negativa funciona, mas para distribuição normal não. Alguém teria uma outra solução, segue CRM: require(MASS) res<-rnorm(100000)### Com 100,000 valores não posso usar shapiro.test ( k <- fitdistr(res,"normal")) par <- k$estimate (md <- par[1])#Média (ep <- par[2])#sd (N <- length(res)) est <-N*dnorm(res,mean=md,sd=ep) ## Estimativas ##Teste Chi - Quadrado para valiar a significância da distribuição fecdf <- ecdf(res) ###ecdf- Função de distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) chisq.test(x=table(emp),y=table(est))## Teste de Chi quadrado Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 3

Conexão em postgres
by Edson Lira 19 Fev '13

19 Fev '13

Boa Tarde, Gostaria de saber como faço para criar um conexão com o postgres. segue o exemplo que fiz require(RPostgreSQL) drv <- dbDriver("PostgreSQL") con <- dbConnect(drv,user="consulta",password="idoctor@)!)",dbname="iDoctor_DW") mas essa conexão é local, gostaria de acessar o proprio servidor de uma outra máquina, como faço? Edson Lira Estatístico Manaus-Amazonas

1 0

tabela muito grande
by Edson Lira 18 Fev '13

18 Fev '13

Bom dia a todos. Estou analisando um banco que necessito fazer a seguinte rotina: GEN<-table(m$LOCALDENAS,m$MUNICIPIOR,m$GENOTIPAGE,m$REG) ftable(addmargins(GEN, c(1,2,3,4), FUN=list(Total = sum))) Que dá essa saída: D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 CAREIRO A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 CASTANHO A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 COARI A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 CODAJAS A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 EIRUNEPE A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 IPIXUNA A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 ITAITUBA A 0 0 0 0 0 0 0 0 0 0 0 0 A sáida é muito grande e não consigo visualizar na tela de saída do R. Como eliminar os zeros. Edson Lira Estatístico Manaus-Amazonas

3 4

[Dúvida] Paralelizar código R.
by Pedro Rafael 17 Fev '13

17 Fev '13

Pessoal não entendo muito de computação paralela. Quais os pacotes mais utilizados em R com esse fim? Só devo paralelizar o código quando tenho um cluster de computadores conectados ou faz sentido também distribuir parte do processo com processos escravos em meu notebook por exemplo? Alguém poderia me dar uma ideia sobre esse assunto? Quem tiver algum material falando sobre computação paralela em R poderia enviar para meu email? Obrigado desde já. [ ], Pedro Rafael Diniz Marinho.

2 3

Letras gregas: como utilizar menos linhas de comandos....
by andrebvs＠bol.com.br 16 Fev '13

16 Fev '13

1 0

Concurso Substituto - Mat. Aplicada - FEIS/UNESP
by Alan R Panosso 16 Fev '13

16 Fev '13

Caros amigos, Encaminho para a divulgação dessa oportunidade. Abraços e a todos e muito obrigado. Alan De: Berenice Damasceno [mailto:berenice@mat.feis.unesp.br] Enviada em: sexta-feira, 15 de fevereiro de 2013 16:22 Para: Berenice Damasceno Assunto: Confirmando período de inscrição: Concurso Substituto - Mat. Aplicada Boa tarde! Peço sua colaboração na divulgação de concurso para professores emergenciais, na Área de Matemática Aplicada, Computação e Estatística, no Departamento de Matemática da Faculdade de Engenharia de Ilha Solteira (FEIS) - UNESP. Inscrições de 18 a 22/02/2013 na Seção de Comunicações da FEIS 02 (duas) vagas para o conjunto de disciplinas "Bioestatística"; "Técnicas Experimentais com Animais"; "Técnicas e Análises Experimentais"; "Estatística e Probabilidade"; "Probabilidade e Estatística". Informações sobre o concurso poderão ser obtidas acessando: http://www.feis.unesp.br/#!/concursos/departamento-de-matematica/ (a partir de 18/02/13) (o edital de abertura de inscrições será publicado no DOE de SP amanhã, dia 16/02/13) Grata pela colaboração, cordialmente, Berenice. ************************************************************************* Profª. Berenice C. Damasceno Departamento de Matemática - Faculdade de Engenharia de Ilha Solteira - FEIS - UNESP Alameda Rio de Janeiro, 266 15 385 - 000 - Ilha Solteira - SP.-Brasil E_mail: berenice(a)mat.feis.unesp.br ou damascenoberenice(a)yahoo.com.br http://www.mat.feis.unesp.br/berenice_damasceno.php Tel.: +55 (18) 37 43 10 91 Fax : +55 (18) 37 43 11 55 Cel.: +55 (15) 97 15 67 88 *************************************************************************

1 0

Contrastes multivariados
by Tiago Souza Marçal 15 Fev '13

15 Fev '13

Boa noite pessoal, estou utilizando a MANOVA para analisar um experimento fatorial com duas características: f1<-c(rep("A",20),rep("B",20)) f2<-rep(c(gl(4,5)),2) res1<-c(0,0, 20, 0, 20, 20, 0, 0, 20, 40, 30, 40, 40, 20, 40, 30, 20, 30, 30, 50, 40, 40, 40, 70, 60, 40, 40, 40, 60, 20, 30, 30, 50, 40, 40, 50, 30, 40, 30, 30)res2<-c(0,0,1,0,1,1,0,0,1,1.75,1.333333,1,1.5,1,1.75,1.333333,1.5,1.666667,1.333333,1.666667,1,1,1.25,1,1.166667,1.25,1.25,1.25,1.333333,1,2.333333,1,1.2,1,1,2.2,1.333333,1.5,1.333333,1)f1<-factor(f1)f2<-factor(f2) y<-cbind(res1,res2) require(mvoutlier) aq.plot(y, alpha=0.05) # O que devo fazer com os outlines multivariados? mav<-manova(y~f1*f2) summary(mav) # Como testar os contrastes multivariados de interesse no R? Ex: Nível A de f1 vs Nível B de f1 dentro do nível um de f2 # Com fazer a análise Stepdown no R? O R possui outros testes multivariados implementados? Desde já eu agradeço pela ajuda. Att. Tiago.

1 0

parcela subdividida
by Ze Henrique 15 Fev '13

15 Fev '13

Boa noite a todos. Gostaria de saber se alguem saberia como faço para calcular os teste de bartlett, Lilliefors para parcela subdividida. Desde já obrigado a todos. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

2 2

Dúvida
by ana paula coelho madeira 15 Fev '13

15 Fev '13

Boa tarde, É possível fazer análise de dados no R de um delineamento rotacional? Como faço?AgradecidaAna Paula

2 1

Histograma de frequência com intervalos de classe
by Adriano Montaguti 15 Fev '13

15 Fev '13

Boa tarde a todos. Estou tendo aulas de probabilidade e estatística nesse semestre e preciso de fazer alguns histogramas no R. Basicamente eu tenho essa pequena amostra: X<-c( 4.5 ,4.7 ,3.2 ,6.3 ,6.4 ,7.1 ,7.2 ,5.1 ,3.7 ,4.2 ,16.2 ,15.9 ,17.3 ,5.1 , 4.6, 8.4, 9.0, 3.6, 3.2, 6.3, 18.4, 3.4, 9.1, 5, 3.6, 2.6, 10.1, 5.6, 6.3) Eu estou tentando construir um histograma de freqüência com intervalos de classe para uma variável contínua. Usei essa função aqui: hist(X,breaks=c(2.6,5.9,9.2,12.5,15.8,19.1),main="Chefes de família do Bairro Alto segundo renda mensal, em salários mínimos - Curitiba - 2002", ylab="Nº de Chefes",xlab="Renda mensal, em salário mínimos", col="red",xlim=c(2.6,19.1),sub="Fonte: Dados fictícios",cex.sub=0.9) Problemas: A ultima classe deveria estar com frequência igual a 5, mas pelo visto conta somente 4. No eixo x eu não quero que ele fique mostrando os valores em escala "5 - 10 - 15" e sim o Ponto médio de cada classe (4.25, 7.55, 10.85, 14.15, 17.45) e que fique centralizado com a largura da coluna. Dados: Aqui estão os intervalos de classe e as frequências: *Tabela 7 - Chefes de família do Bairro Alto segundo renda mensal, em salários mínimos - Curitiba - 2002* * Renda mensal em salários mínimos* *Freqüência* 2,6 |¾ 5,9 15 5,9 |¾ 9,2 9 9,2 |¾ 12,5 1 12,5 |¾ 15,8 0 15,8 |¾| 19,1 5 Total 30 *Fonte:* Dados Fictícios Amplitude de classe = 3.3 Número de classes = 5 Agradeço desde já. Adriano.

2 5

Formatação de hora
by Edson Lira 14 Fev '13

14 Fev '13

Boa tarde a todos, estou trabalhando com um banco que duas variáveis são: a hora de entrada e saída, e o banco que foi criado não diferencia am e pm. Quando a hora é 13:20, foi digitado 01:20, como diferenciar esse horário? Se facilitar não tenho atendimento das 19:00 até as 06:00 hs. Preciso calcular o tempo de atendimento, isso já resolvi usando o difftime. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 3

Big Data
by Vinicius Brito Rocha 14 Fev '13

14 Fev '13

Pessoal, vcs tem recomendações de material sobre Big Data e o R? Andei procurando no google, mas confesso que estou tendo dificuldades. 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.

5 4

plotar nome de município e quantidade num mapa
by Alexandro (Yahoo) 14 Fev '13

14 Fev '13

#Olá,tudo bem! #Gostaria de saber como plotar o nome de município e uma quantidade num mapa: #Tenho uma planilha com os 645 municípios do Estado de São Paulo, das quais como CMR, coloco 5 deles num data.frame: cidade=c("São Paulo","Diadema","São Bernardo do Campo","Santo André","São Caetano") valor=c(6000,450,2500,2000,1000) saopaulo=data.frame(cidade,valor) saopaulo #Tenho os arquivos dbf, sbn, sbx, shp e shx e uso os seguintes comandos: ### 1-Carregando arquivo shape library(maptools); sp = readShapePoly("C:/r/shapefile-sp-ibge/35mu500gc.shp") ### 2-Carregando planilha com o nome do município e valor require(RODBC);cnct = odbcConnectExcel2007(file.choose()) saopaulo = sqlQuery(cnct, "select * from \"Plan1$\" ");odbcClose(cnct);attach(saopaulo);saopaulo sp.dados <- attr(sp, "data") sp.dados$indice <- 1:dim(sp.dados)[1] sp.dados <- merge(sp.dados, saopaulo) sp.dados <- sp.dados[order(sp.dados$indice), ] ### 3-Atribuindo cores para cada intervalo intervalo = c(0,1,500,1000,5000,9999) cat3 <- findInterval(saopaulo$valor,intervalo ) options(scipen=10) table(cat3) cores3 <- c("white","yellow","blue","green","red","black") ### 4-Plotando o mapa #par(mar=c(0,0,0,0)) windows(width=19, height=9, xpos=0, ypos=0,pointsize = 13,rescale="fixed") plot(sp, col=cores3[cat3]);legend("bottomleft",leglabs( c(0,1,500,1000,5000,9999), "<", ">"),fill=cores3, bty="n",text.width=10) ### 5-Colocando números sp.polgns<-attr(sp,"polygons") sp.cntrd<-lapply(sp.polgns,slot,"labpt") text(as.character(saopaulo$valor),cex=0.5) sp.x<-sapply(sp.cntrd,function(x) x[1]) sp.y<-sapply(sp.cntrd,function(x) x[2]) text(sp.x,sp.y,as.character(saopaulo$valor),cex=0.3) #Dúvidas: #1- Como colocar os nomes dos municípios (saopaulo$cidade) em cima dos números no item 5 (para que não se sobreponham)? #2- Quais são os comandos para se plotar uma determinada região. Ex: Municípios do ABC ou Litoral Paulista (somente municípios de interesse no mapa)? Muito obrigado! Felicidades... Alexandro Vieira Lopes

2 1

Re: [R-br] Resultado de prcomp()
by Ivan Bezerra Allaman 13 Fev '13

13 Fev '13

Jovem, Pela sua dúvida aconselho que leia um material teórico sobre componentes principais. Provavelmente verás que seu problema já está resolvido. (S,f,P) Allaman

2 1

Resultado de prcomp()
by Carlos Andrade 13 Fev '13

13 Fev '13

Amigos da lista, Estou praticando análise de componentes principais para reduzir o número de variáveis para uma análise de cluster. Encontrei um exercício e algumas apostilas do R para análise multivariada. Seguem os dados do exercício, que pede para usar 2 componentes principais: > Lista1MB213 city bread hamburger butter apples tomatoes 1 Anchorage 70.9 135.6 155.0 63.9 100.1 2 Atlanta 36.4 111.5 144.3 53.9 95.9 3 Baltimore 28.9 108.8 151.0 47.5 104.5 4 Boston 43.2 119.3 142.0 41.1 96.5 5 Buffalo 34.5 109.9 124.8 35.6 75.9 6 Chicago 37.1 107.5 145.4 65.1 94.2 7 Cincinnati 37.1 118.1 149.6 45.6 90.8 8 Cleveland 38.5 107.7 142.7 50.3 93.2 9 Dallas 35.5 116.8 142.5 62.4 90.7 10 Detroit 40.8 108.8 140.1 39.7 96.1 11 Honolulu 50.9 131.7 154.4 65.0 93.9 12 Houston 35.1 102.3 150.3 59.3 84.5 13 KansasCity 35.1 99.8 162.3 42.6 87.9 14 LosAngeles 36.9 96.2 140.4 54.7 79.3 15 Milwaukee 33.3 109.1 123.2 57.7 87.7 16 Minneapolis 32.5 116.7 135.1 48.0 89.1 17 NewYork 42.7 130.8 148.7 47.6 92.1 18 Philadelphia 42.9 126.9 153.8 51.9 101.5 19 Pittsburgh 36.9 115.4 138.9 43.8 91.9 20 StLouis 36.9 109.8 140.0 46.7 79.0 21 SanDiego 32.5 84.5 145.9 48.5 82.3 22 SanFrancisco 40.0 104.6 139.1 59.2 81.9 23 Seattle 32.2 105.4 136.8 54.0 88.6 24 Washington 31.8 116.7 154.8 57.6 86.6 dados <- Lista1MB213[,2:6] standardise <- as.data.frame(scale(dados)) lanche.pca <- prcomp(standardise) > summary(lanche.pca) Importance of components: PC1 PC2 PC3 PC4 PC5 Standard deviation 1.5591 0.9654 0.9134 0.7252 0.52613 Proportion of Variance 0.4862 0.1864 0.1669 0.1052 0.05536 Cumulative Proportion 0.4862 0.6726 0.8394 0.9446 1.00000 A minha dúvida é como distribuir as cinco variáveis no dois componentes PC1 e PC2 que juntos explicam 67,26% da variância. Como seria a equação de combinação linear? Há muito tempo li em um tutorial do SPSS a solução, mas não lembro onde encontrá-lo. Agradeço antecipadamente qualquer ajuda. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

3 3

LEG Ufpr
by Gilbert Queiroz 13 Fev '13

13 Fev '13

Pessoal, Alguem sabe dizer o que houve com a página do LEG Ufpr??? Está fora desde domingo, 10 Fev...Abs. Gilbert.

2 1

Converter dias em anos
by Fátima Lima Paula 13 Fev '13

13 Fev '13

Olá pessoal, como passar a diferença entre datas (idade) em dias para anos? Já havia feito isso em algum lugar, mas não encontro. datanasc=c(19250405,19300809,19400402) dataint=c(20080302,20090506,20080104) datanasc=as.character(datanasc,"%Y%m%d") dataint=as.character(dataint,"%Y%m%d") datanasc=as.Date(datanasc,"%Y%m%d") dataint=as.Date(dataint,"%Y%m%d") idade=dataint-datanasc idade Time differences in days [1] 30282 28760 24748 Agora preciso ter a idade de cada um em anos, com números inteiros. Alguém pode me ajudar? Obrigada Fátima

3 7

Selecionar vários processadores
by Alexandre Santos 13 Fev '13

13 Fev '13

Bom dia pessoal, Estou querendo reduzir o tempo de execução de alguns cálculos utilizando o uso de vários processadores em paralelo, uma vez que o R utiliza somente um núcleo. Pesquisei na internet e encontrei vários pacotes como: multicore, parallel, foreach, doMC e outros, no entanto, eles tem apenas algumas funções em paralelo, como por exemplo o pacote foreach e o doMC, onde você substitui o for{} dos loops por foreach{} e consegue obter especificando o numero de núcleos que deseja utilizar(registerDoMC(cores=4)), cálculos mais rápidos. No entanto, gostaria de saber se existe uma maneira de especificar o número de núcleos que deseja utilizar e o R trabalhar em paralelo para todas as funções, mesmo a for{}? Meu pc é um core i5 com 8MB de RAM com SO Ubuntu 11.04. Obrigado, Alexandre

5 6

Comparação entre duas populações com distribuição de Poisson
by Alexandre Santos 12 Fev '13

12 Fev '13

Bom dia pessoal, Tenho duas populações com distribuição de Poisson e gostaria de saber qual a abordagem mais correta para comprar essa duas populações, se empregando análise de variância sobre um GLM ajustado ou teste de Chi quadrado sendo que os dois testes fornecem o mesmo resultado, sendo o CRM: # Dados y<- c(mapply(rpois, lambda=c(5,5), MoreArgs=list(n=15)))##Criação da variável resposta trat <- gl(2,15) ##Criação da variável resposta #Usando GLM m.1 <- glm(y~trat,family="poisson") anova.m.1 <- anova(m.1, test="Chi") anova.m.1 # #Usando teste Chi^2 chisq.test(y,correct=TRUE) # Obrigado, Alexandre

2 2

BLUP
by Ze Henrique 10 Fev '13

10 Fev '13

Bom dia a todos. Gostaria de saber se alguem conhece como trabalhar com BLUP no R, pacote ou algo que faça este tipo de estatística. Desde já agradeço a ajuda. Obrigado. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

2 2

Método de Newton Rhapson...
by andrebvs＠bol.com.br 09 Fev '13

09 Fev '13

2 9

Join my network on LinkedIn
by Juliana Passos via LinkedIn 09 Fev '13

09 Fev '13

LinkedIn ------------ Juliana Passos requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Juliana Passos http://www.linkedin.com/e/-33g73l-hcytrjre-4i/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Juliana Passos http://www.linkedin.com/e/-33g73l-hcytrjre-4i/rso/230788258/wFJc/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hcytrjre-4i/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Programando com C e R
by ctrucios 08 Fev '13

08 Fev '13

Boa tarde pessoal, Estou tentando utilizar C e R juntos, fiz 3 funções em C e chamo elas do R. Mas estive testando e testando e deu este erro. Na verdade as vezes da este erro e as vezes trava de vez. Não sei se é pela forma em que estou programando em C (na verdade estou bem enferrujado e estou utilizando ponteiros) deixo aqui o código quem puder de uma olhadinha. [image: Imagem inline 2] Na verdade é um mesmo código que faço três vezes (assim que descobri o erro). O problema aparece quando rodo mais de dois vezes (por isso que aqui copiei o codigo 3 vezes) source("volatnGarch.r") source("sbootGarch.r") source("forebootGarch.r") yb = rnorm(10000) d = 1 p=1 q=1 garchOrder = c(p,q) Coeff_b = c(0.01,0.1,0.85) aux1 = Coeff_b[1]/(1-Coeff_b[2]-Coeff_b[3]) aux3 = cumsum(c(1,1)) C = sbootGarch(aux1,yb,Coeff_b,garchOrder) A = volatnGarch(yb,Coeff_b,garchOrder) e_b = rnorm(5) n.ahead = length(e_b) B = forebootGarch(yb,e_b,A[[1]],Coeff_b,garchOrder) rm(list=ls(all=TRUE)) source("volatnGarch.r") source("sbootGarch.r") source("forebootGarch.r") yb = rnorm(10000) d = 1 p=1 q=1 garchOrder = c(p,q) Coeff_b = c(0.01,0.1,0.85) aux1 = Coeff_b[1]/(1-Coeff_b[2]-Coeff_b[3]) aux3 = cumsum(c(1,1)) C = sbootGarch(aux1,yb,Coeff_b,garchOrder) A = volatnGarch(yb,Coeff_b,garchOrder) e_b = rnorm(5) n.ahead = length(e_b) B = forebootGarch(yb,e_b,A[[1]],Coeff_b,garchOrder) rm(list=ls(all=TRUE)) source("volatnGarch.r") source("sbootGarch.r") source("forebootGarch.r") yb = rnorm(10000) d = 1 p=1 q=1 garchOrder = c(p,q) Coeff_b = c(0.01,0.1,0.85) aux1 = Coeff_b[1]/(1-Coeff_b[2]-Coeff_b[3]) aux3 = cumsum(c(1,1)) C = sbootGarch(aux1,yb,Coeff_b,garchOrder) A = volatnGarch(yb,Coeff_b,garchOrder) e_b = rnorm(5) n.ahead = length(e_b) B = forebootGarch(yb,e_b,A[[1]],Coeff_b,garchOrder) rm(list=ls(all=TRUE)) as funções em C seguem no seguinte link. https://www.dropbox.com/sh/5f8bpx6pt2cgam6/fcuvUtFili/Codigos *Carlos Trucíos Maza ** : ctrucios(a)gmail.com

3 3