Abril 2013 - R-br - DINF UFPR

Alocação proporcional no R
by Matheus Monteiro 29 Abr '13

29 Abr '13

Ainda trabalhando com os 50 municipios mais populosos do Brasil, divido a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo os 40 restantes. Preciso obter uma amostra com reposição de 10 municípios usando o alocação proporcional e calcular a estimativa da média e variancia. Alguém pode ajudar? pop = c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569)

2 6

R 3.0 Ubuntu 13.04
by Frederico Alvarez 29 Abr '13

29 Abr '13

Boa tarde a todos. Estou tentando instalar o R 3.0.0 no ubuntu 13.10 via apt. Mas acontece que no repositório da UFPR não há o pacote r-base-core para 64 bit. Porém, há dois pacotes para 32 bits. Alguém sabe como posso resolver esse problema? Saudações. -- Frederico Steinmetz Alvarez Engenheiro de Produção Elétrica Mestrando em Energia Eólica Laboratório Didático Pedagógico de Fluidos - LDPFLU Departamento de Engenharia Mecânica - DEMEC Universidade Federal de Pernambuco - UFPE (81) 8885-1105 (Oi) (48) 8448-5994 (Oi) (48) 9143-5141 (Vivo)

9 11

Alocação Proporcional
by Matheus Monteiro 29 Abr '13

29 Abr '13

Ainda trabalhando com os 50 municipios mais populosos do Brasil, divido a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo os 40 restantes. Preciso obter uma amostra com reposição de 10 municípios usando o alocação proporcional e calcular a estimativa da média e variancia. Alguém pode ajudar? Estou com dúvida aqui, pois preciso encontrar uma maneira de alocar proporcionalmente a população inteira e não apenas o estrato 2 pop = c(11376685,6390290,2710968,2648532,2500532,2395785, 1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610, 1016128,953393,867067,830231,817590,805397,801746,774886, 742428,680496,668877,654786,643603,619746,619536,613815, 600692,587701,568099,561329,526338,525225,515707,491807, 483821,474596,474596,472300,460062,446911,442701,433158, 425169,424948,422569) # order pop <- sort(pop, T) # Alocação sub1 <- pop[(1:10)] sub2 <- pop[(11:50)] # Alocação proporcional DÚVIDA prop <- sub2/sum(sub2) set.seed(51) sub3 <- sample(sub2, 10, replace=TRUE, prob=prop)

2 1

Reduzir objeto auxiliares...
by andrebvs＠bol.com.br 29 Abr '13

29 Abr '13

2 1

(sem assunto)
by Antonio Silva 29 Abr '13

29 Abr '13

Olá Eu estou tentado utilizar o teste de Qui-quadrado para verificar a dependência entre a distribuição de duas espécies, como indicado no livro do Jean Valetin (Ecologia Numérica 2a ed, pg 49). Eu tenho um dataframe com dados de presença/ausência de duas espécies em um determinado número de unidades amostrais A equação indicada no livro é: X2= p*(|ad-bc|-p/2)^2 / ((a+b)*(c+d)*(a+c)*(b+d)) Onde a é o número de dupla presença (1-1), b o número de 1-0, c o número de 0-1 e d o número de 0,0 (dupla ausência) p é a+b+c+d Há alguma função no R para este cálculo? Eu não consegui utilizar a função chisq.test para este cálculo. Obrigado! Antonio Olinto

4 4

Estratificação no R
by Matheus Monteiro 29 Abr '13

29 Abr '13

Estou trabalhando com os 50 municipios mais populosos do Brasil, foi pedido para dividir a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo o restante. Preciso obter uma amostra de 5 municípios de cada estrato e calcular a estimativa da media amostral e variancia. Sei que preciso instalar o pacote gtools mas não estou conseguindo. Não estou conseguindo, alguém pode ajudar? Por enquanto está assim: pop = c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569) Npop=length(pop) Npop EstI=c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014) EstII=c(1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569) nEstI=length(EstI) nEstII=length(EstII)

5 8

Sequência no R...
by andrebvs＠bol.com.br 29 Abr '13

29 Abr '13

2 2

Manipulação de dados diários
by Marcos Vagner Macedo 28 Abr '13

28 Abr '13

Boa tarde, Estou com dificuldade na utilização de dados no formato zoo. Estou usando o comando abaixo. z <- read.zoo("dados.csv", header = TRUE, sep = ",", format = "%m-%d-%Y") Está aparecendo o seguinte erro. Erro em read.zoo("dados.csv", header = TRUE, sep = ",", format = "%m-%d-%Y") : index has 4016 bad entries at data rows: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 ... > Vocês poderiam me ajudar? Além disso, gostaria de saber como pegar sub amostras desse banco de dados. A primeira amostra deve ser um período reduzido do banco de dados (OBS:. dados diários - dias uteis) A segunda amostra deve ser algumas séries no período reduzido. Como o banco é grande o ideal é que os filtros seja feitos de uma forma que seja possível automatizar o processo. Exemplo: 5 séries . Periodo 10/01/2000 a 28/12/2012 Analise: 3 séries Periodo 08/01/2012 a 06/08/2012

2 1

Join my network on LinkedIn
by Sérgio Henrique Almeida da Silva Junior via LinkedIn 28 Abr '13

28 Abr '13

LinkedIn ------------ Sérgio Henrique Almeida da Silva Junior requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Sérgio Henrique Almeida da Silva Junior http://www.linkedin.com/e/-33g73l-hg1l8cl7-5x/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Sérgio Henrique Almeida da Silva Junior http://www.linkedin.com/e/-33g73l-hg1l8cl7-5x/rso/250437145/r-Yw/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-hg1l8cl7-5x/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Ifelse para duas condições
by Fátima Lima Paula 27 Abr '13

27 Abr '13

Preciso fazer um ifelse para duas condições. Exemplo nome=c("Ricardo","Maria","Tereza","3edro","José", "Germano") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima") idade=c(68,67,89,78,65,67) idade2=c(68,67,90,77,65,70) ano=c(1944,1945,1923,1954,1947,1945) ano2=c(1953,1945,1923,1951,1947,1956) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209) banco1=cbind(nome,sobrenome,idade,idade2,ano,ano2,dataint) banco1=data.frame(banco1) banco1 banco1$ano=as.character(banco1$ano) banco1$ano2=as.character(banco1$ano2) banco1$idade=as.character(banco1$idade) banco1$idade2=as.character(banco1$idade2) Quero criar uma coluna que me diga sim se as idades e os anos estiverem iguais e não caso contrário. Fiz assim, banco$teste=if(banco1$idade==banco1$idade2 & banco1$ano==banco1$ano2)sim else nao Naturalmente, deu erro. Alguém pode me ajudar, por favor? -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 2