Fevereiro 2013 - R-br

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

14 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