- R-br - DINF UFPR

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